24小时热门版块排行榜    

CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知
查看: 703  |  回复: 15
当前主题已经存档。

cypemail

铁虫 (小有名气)

[交流] 【求助】这段程序是什么语言写的? C么?【已完结】

r=20;
alkoi=1.162-0.342/r;
for(pp=1;pp<1000;pp+=1){
p=pp/1000;
for(aa0=0;abs(alkoi-aa0)>0.001; ){
aa0=alkoi;
a0=0.5*(r+aa0);
a11=0.414*aa0*p-0.707*aa0-r*p;
a2=p^2*(0.5*r+0.0858*aa0)-p*0.293*aa0;
a1=-a11/(2*a0);
z11=a1-(a1^2-a2/a0)^0.5;
pe=0.293*p+0.707*z11;
if(r>1){
if(pe<(0.342-0.342/r)/ln(r)){alkoi=0.82+pe*ln(r)};
};
};
if(z11<0){z11=0};
px=p;
py1=z11;
pyx=p;
data1_a[pp]=px;
data1_b[pp]=pyx;
data1_c[pp]=py1;
};
alkoi=1.162-0.342/r;
for(pp=1;pp<1000;pp+=1){
p=pp/1000;
for(aa0=0;abs(alkoi-aa0)>0.001; ){
aa0=alkoi;
a0=0.5*(r+aa0);
a11=0.414*aa0*p-0.707*aa0-r*p;
a2=p^2*(0.5*r+0.0858*aa0)-p*0.293*aa0;
a1=-a11/(2*a0);
z22=a1+(a1^2-a2/a0)^0.5;
pe=0.293*p+0.707*z22;
if(r>1){
if(pe>1-0.342*(1-1/r)/ln(r)){alkoi=0.82+(1-pe)*ln(r)};
};
};
if(z22>0.999999){z11=0.999999};
py2=z22;
data1_d[pp]=py2;
};

[ Last edited by 余泽成 on 2009-11-3 at 09:57 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aitao

至尊木虫 (正式写手)

★ ★
小木虫(金币+0.5):恭喜抢沙发,给个红包
nono2009(金币+1,VIP+0):鼓励新虫!欢迎常来程序软件版交流。 10-29 20:31
nono2009(金币+0,VIP+0):请参阅4楼的comments. 10-30 07:28
据我所知是的
2楼2009-10-29 20:03:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cypemail

铁虫 (小有名气)

引用回帖:
Originally posted by aitao at 2009-10-29 20:03:
据我所知是的

请问您了解 求解非线性方程的Ridder's method么?
3楼2009-10-29 20:18:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

magic7004

金虫 (职业作家)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
余泽成(金币+1,VIP+0):谢谢参与交流! 10-29 22:19
看语法貌似不是C
引用回帖:
if(pe<(0.342-0.342/r)/ln(r)){alkoi=0.82+pe*ln(r)};
};
};
if(z11<0){z11=0};
px=p;

C里面貌似}后面没有;的
流氓不可怕,可怕的是流氓有文化,有文化又BH的流氓无敌~~!
4楼2009-10-29 20:41:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mickk

铁杆木虫 (职业作家)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
余泽成(金币+1,VIP+0):排除法,可行。欢迎常来程序软件版! 10-30 07:44
排除fortran,matlab,mathematics
5楼2009-10-30 07:28:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

magic7004

金虫 (职业作家)


小木虫(金币+0.5):给个红包,谢谢回帖交流
搂住想要做什么呢?

这段程序看起来很清楚的,搂住需要什么语言就按什么语言的语法改写一下不就行了?
流氓不可怕,可怕的是流氓有文化,有文化又BH的流氓无敌~~!
6楼2009-10-30 08:30:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aitao

至尊木虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by cypemail at 2009-10-29 20:18:

请问您了解 求解非线性方程的Ridder's method么?

这个不太了解...
7楼2009-10-30 11:44:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhuhongaaa

金虫 (正式写手)

扔到C,C++,java编译应该都可以通过的。。。
分子模拟的主页http://varmilion.tk/
8楼2009-10-30 11:52:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaoxiaoqi

木虫 (著名写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
余泽成(金币+1,VIP+0):感谢交流,欢迎常来程序软件版! 10-30 17:49
排除语法错误
猜测作者的意思来看,应该不是java,x^y在java中写作math.pow(x,y);^在java里是逻辑运算符
                             也不会是c,
原因有二:除非作者定义函数lnx=logx,c语言中求自然对数用logx
              }; 在c里面一般是算作语法错误的。
个人观点,仅供参考。
9楼2009-10-30 15:32:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cypemail

铁虫 (小有名气)


余泽成(金币+1,VIP+0):谢谢交流,欢迎常来程序软件版! 10-30 17:49
引用回帖:
Originally posted by zhaoxiaoqi at 2009-10-30 15:32:
排除语法错误
猜测作者的意思来看,应该不是java,x^y在java中写作math.pow(x,y);^在java里是逻辑运算符
                             也不会是c,
原因有二:除非作者定义函数lnx=logx,c语言中求自然对数用l ...

这个是别人附的求解非线性方程的代码, 我看不懂怎么求解的。如果是c的话,我放在turbo c的程序里运行就能出结果了么?
10楼2009-10-30 16:19:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cypemail 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见