| 查看: 869 | 回复: 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 ] |
» 猜你喜欢
0856材料化工调剂 总分330
已经有13人回复
330分求调剂
已经有5人回复
291求调剂
已经有3人回复
086000生物与医药调剂
已经有8人回复
352分-085602-一志愿985
已经有3人回复
291求调剂
已经有4人回复
309求调剂
已经有4人回复
0856求调剂
已经有13人回复
0703化学
已经有11人回复
329求调剂
已经有7人回复
7楼2009-10-30 11:44:25
3楼2009-10-29 20:18:38

4楼2009-10-29 20:41:37
mickk
铁杆木虫 (职业作家)
- 应助: 11 (小学生)
- 金币: 9196.3
- 散金: 1002
- 红花: 3
- 帖子: 4608
- 在线: 643.3小时
- 虫号: 179848
- 注册: 2006-02-09
- 性别: GG
- 专业: 光子与光电子器件
5楼2009-10-30 07:28:21














回复此楼