| 查看: 717 | 回复: 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 ] |
回复此楼» 猜你喜欢
所感
已经有3人回复
-
要不要辞职读博？
已经有7人回复
-
不自信的我
已经有11人回复
-
北核录用
已经有3人回复
-
实验室接单子
已经有3人回复
-
磺酰氟产物，毕不了业了！
已经有8人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有10人回复
-
26申博（荧光探针方向，有机合成）
已经有4人回复
-
论文终于录用啦！满足毕业条件了
已经有26人回复
-
2026年机械制造与材料应用国际会议 （ICMMMA 2026）
已经有4人回复
★ ★
小木虫(金币+0.5):恭喜抢沙发，给个红包
nono2009(金币+1,VIP+0):鼓励新虫！欢迎常来程序软件版交流。 10-29 20:31
nono2009(金币+0,VIP+0):请参阅4楼的comments. 10-30 07:28
小木虫(金币+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
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
6楼2009-10-30 08:30:07
7楼2009-10-30 11:44:25
zhuhongaaa
金虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 396.9
- 红花: 2
- 帖子: 607
- 在线: 23小时
- 虫号: 748787
- 注册: 2009-04-14
- 性别: GG
- 专业: 高分子物理与高分子物理化
8楼2009-10-30 11:52:29
9楼2009-10-30 15:32:52
10楼2009-10-30 16:19:26