| 查看: 756 | 回复: 4 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
forcal 代码运行后出错
|
|||
|
运行了别人的Forcal代码，出现了“Forcal表达式运行错误！ 实数表达式名称：main 所在模块：-1 实数函数名：fcopt::SimOpt 错误代码：3”为什么？能不能帮忙调试下，指出问题所在！！谢谢~~ 代码：!using["sys","IMSL","fcopt"]; f(t,HCp,dHCp : R : k1,a,E1, Acid,T)= { R=8.3144, dHCp=-k1*Acid^a*exp[-E1/(R*T)]*HCp }; 优化(_k1, _a, _E1 : i, j, s, t1, t2, HCp, IDO : k1,a,E1, Acid,T, hf, tArray, max, TempArray,AcidArray,maxHCp)= { k1=_k1, a=_a, E1=_E1, s=0, j=0, (j IDO=1, t1=tArray.A(0,0), HCp=tArray.A(0,j+1), i=1, (i IVPRK(&IDO, hf, &t1, t2, 1e-6, 0, &HCp), s=s+[HCp-tArray.A(i,j+1)]^2, i++ }, IVPRK(3, hf, &t1, t2+0.1, 1e-6, 0, &HCp), j++ }, s }; main(:i, k1, a, E1, min : hf, tArray, max, TempArray,AcidArray,maxHCp)= { ClearImslErr(), //清空IMSL的错误输出 ERSET(0, 0, 0), //关闭IMSL错误输出 hf=HFor("f"  ,max=10, maxHCp=8,//实验数据组数 tArray=new{rtoi(real_s),rtoi(max),rtoi(maxHCp+1)}.free().io::arrayns{ //存放实验数据 " 0 14.277 12.892 13.880 12.859 13.526 15.058 14.705 14.042 5 14.057 12.721 12.986 12.576 13.138 11.872 14.530 12.709 10 13.418 11.721 11.350 11.616 12.309 9.015 13.831 12.133 20 11.691 9.383 10.122 9.069 10.231 7.119 12.748 9.872 30 9.220 7.622 7.795 7.321 7.870 5.223 10.619 5.949 40 7.875 6.287 6.781 5.880 6.358 4.451 8.678 4.741 60 6.554 4.687 6.101 4.076 5.310 3.162 8.043 3.527 80 5.776 4.946 4.416 3.784 2.915 2.284 5.270 3.196 100 4.646 3.586 3.434 2.716 2.502 2.070 4.072 2.489 120 2.998 3.101 1.985 2.458 2.370 1.933 3.866 2.989 " }, TempArray=new{rtoi(real_s),rtoi(maxHCp)}.free().io::arrayns{"413.15 413.15 423.15 423.15 433.15 433.15 443.15 443.15"}, AcidArray=new{rtoi(real_s),rtoi(maxHCp)}.free().io::arrayns{"0.7 0.9 0.5 0.7 0.3 0.5 0.1 0.3"}, min=1.1, k1=0.3, a=1.2, E1=1e5, //提供一组合适的初值 i=SimOpt[HFor("优化" , optstep,0, optexpand,1.5, optcontract,0.5, optmax,5000 : &k1, &a, &E1, &min], //求n维极值的单形调优法printff{"\r\n实际迭代次数={1,i}, k1={2,r}, a={3,r}, E1={4,r}, 目标函数终值={5,r}\r\n",i, k1, a, E1, min}, ERSET(0, 2, 2), ERSET(0, 1, 0) //打开IMSL错误输出 }; [ Last edited by ahbiao on 2011-3-27 at 21:41 ] |
回复此楼» 猜你喜欢
B区学生调剂-兰州交通大学材料科学与工程学院
已经有12人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有13人回复
-
物理学I论文润色/翻译怎么收费?
已经有168人回复
-
桂林理工大学物理学专业招收调剂
已经有18人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有6人回复
-
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有19人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助用MS计算声子谱出错的问题
已经有15人回复
- fortran转换CHGCAR出错
已经有4人回复
- VASP并行计算出错
已经有9人回复
- Gaussian 计算拉曼光谱最后时刻出错
已经有14人回复
- 急！关于nsfc申请代码的问题
已经有6人回复
- 得不到纯产物，哪出错了，亲们帮忙看看啊
已经有3人回复
- 高斯09中的L502出错的问题
已经有5人回复
- 关于CASTEP 中计算时出错
已经有6人回复
- 出错 warning#2 #180
已经有7人回复
- 选错申请代码了，影响大吗
已经有29人回复
- MS一打开就出错，任何任务都不能运行
已经有12人回复
- 博士后基金申请在打印转化成word时出错
已经有13人回复
» 抢金币啦！回帖就可以得到:
-
上海交通大学张万斌院士团队招聘科研助理和博士候选人
+1/81
-
网爆某大型玻璃厂家
+1/68
-
万人领军
+3/47
-
遵义医科大学老师招收2026年微生物与生化药学学术型硕士调剂学生（3名）
+1/39
-
坐标郑州，征女友
+1/36
-
江苏盐城工学院材料与化工专硕还有调剂名额，招生08开头专硕或学硕
+1/31
-
广西科技大学化工专业现仍有较多调剂指标，欢迎符合条件的同学申请调剂
+1/29
-
瑞典 Linköping University 诚聘 氮化物基无机材料方向 博士后1名
+1/23
-
招聘学术型博士--噬菌体方向
+1/21
-
黑龙江科技大学 化工招收调剂 07理学 08工学 09农学 均可。
+3/18
-
坐标广州，征女友
+2/14
-
南京林业大学-国家级青年人才团队 招2026级申请考核制博士 （有机合成化学方向）
+1/12
-
北京理工大学（珠海）招收智能科学与技术方向26级博士
+1/7
-
非全博士求助
+1/6
-
青岛科技大学高分子学院硕士招生
+1/5
-
从上海旅行归来，想分享一下我的真实感受，每一句话都是真情实感
+1/3
-
科研助理
+1/3
-
研究助理招聘 | 中山大学光电材料与技术全国重点实验室--集成信息光电子研究团队
+1/2
-
小身材大能量，易度质量流量计在雾化器测量中的精准控制
+1/2
-
中国药科大学招聘药物化学方向博士后（长期有效）
+1/1
ahbiao(金币+40): 2011-03-31 08:34:11
|
本帖内容被屏蔽 |
4楼2011-03-30 15:03:23
ahbiao(金币+10): 谢谢~~ 2011-03-30 08:10:16
|
本帖内容被屏蔽 |
2楼2011-03-29 10:10:20
3楼2011-03-30 08:12:07
恩，谢谢！ 5楼2011-03-31 08:34:35
