| 查看: 1112 | 回复: 15 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】这段程序是什么语言写的? 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 ] |
回复此楼» 猜你喜欢
40+中年人失业后，还能找到好工作吗？
已经有11人回复
-
用NBS溴化时溶剂除了CCl4外还有什么
已经有6人回复
-
申请科研项目纯靠实力能成功吗？
已经有11人回复
-
极限拉扯，如何抉择
已经有15人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
4楼2009-10-29 20:41:37
3楼2009-10-29 20:18:38
mickk
铁杆木虫 (职业作家)
- 应助: 11 (小学生)
- 金币: 9204.8
- 散金: 1002
- 红花: 3
- 帖子: 4609
- 在线: 643.4小时
- 虫号: 179848
- 注册: 2006-02-09
- 性别: GG
- 专业: 光子与光电子器件
5楼2009-10-30 07:28:21
6楼2009-10-30 08:30:07
200