24小时热门版块排行榜    

查看: 1115  |  回复: 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的回帖

nono2009

超级版主 (文学泰斗)

No gains, no pains.

优秀区长优秀区长优秀区长优秀区长优秀版主

★ ★
余泽成(金币+1,VIP+0):辛苦了! 10-31 11:19
jjdg(金币+1,VIP+0):感谢解释! 11-4 10:10
伪代码(pseudocode),又称为虚拟代码,是高层次描述算法的一种方法。它不是一种现实存在的编程语言(已经出现了类似伪代码的语言,参见Nuva);它可能综合使用多种编程语言中语法、保留字,甚至会用到自然语言。

它以编程语言的书写形式指明算法的职能。相比于程序语言(例如Java, C++,C, Delphi 等等)它更类似自然语言。它是半形式化、不标准的语言。我们可以将整个算法运行过程的结构用接近自然语言的形式(这里,你可以使用任何一种你熟悉的文字,中文,英文 等等,关键是你把你程序的意思表达出来)描述出来. 使用伪代码, 可以帮助我们更好的表述算法, 不用拘泥于具体的实现.

人们在用不同的编程语言实现同一个算法时意识到,他们的实现(注意:这里是实现,不是功能)很不同。尤其是对于那些熟练于不同编程语言的程序员要理解一个(用其他编程语言编写的程序的)功能时可能很难,因为程序语言的形式限制了程序员对程序关键部分的理解。这样伪代码就应运而生了。

当考虑算法功能(而不是其语言实现)时,伪代码常常得到应用。计算机科学在教学中通常使用虚拟码,以使得所有的程序员都能理解。

(http://zh.wikipedia.org/wiki/%E4%BC%AA%E4%BB%A3%E7%A0%81)
13楼2009-10-31 10:54:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 16 个回答

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的回帖
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[职场人生] 40+中年人失业后,还能找到好工作吗? +11 xqy2009 2026-10-13 13/650 2026-10-15 21:32 by bear2007
[基金申请] 申请科研项目纯靠实力能成功吗? 200+7 David1982 2026-10-13 11/550 2026-10-15 10:30 by zyqchem
[硕博家园] 极限拉扯,如何抉择 +11 Isak0320 2026-10-12 15/750 2026-10-14 21:55 by 猫鱼不是鱼
[论文投稿] 二审4个多月了还是Awaiting Reviewer Scores 10+3 纳斯达克ML 2026-10-11 4/200 2026-10-14 17:30 by babero
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +8 wangyl_123 2026-10-10 13/650 2026-10-14 15:12 by zuesli
[教师之家] 课题组会帮你通过考核么? +4 好运来~9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
信息提示
请填处理意见