| 查看: 698 | 回复: 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 ] |
回复此楼» 猜你喜欢
三甲基碘化亚砜的氧化反应
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
孩子确诊有中度注意力缺陷
已经有12人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
-
论文投稿，期刊推荐
已经有4人回复
-
硕士和导师闹得不愉快
已经有13人回复
4楼2009-10-29 20:41:37
3楼2009-10-29 20:18:38
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