| 查看: 868 | 回复: 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 ] |
回复此楼» 猜你喜欢
309求调剂
已经有3人回复
-
070300化学354求调剂
已经有4人回复
-
330分求调剂
已经有3人回复
-
299求调剂
已经有8人回复
-
一志愿北京理工大学本科211材料工程294求调剂
已经有6人回复
-
300求调剂，材料科学英一数二
已经有8人回复
-
招收生物学/细胞生物学调剂
已经有5人回复
-
070305高分子化学与物理 304分求调剂
已经有7人回复
-
289求调剂
已经有13人回复
-
一志愿哈尔滨工业大学材料与化工方向336分
已经有9人回复
15楼2009-11-01 20:09:51
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
20