| 查看: 735 | 回复: 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 ] |
回复此楼» 猜你喜欢
2026年循环经济功能材料国际会议（ICFMCE 2026）
已经有0人回复
-
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有246人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有5人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
-
光学工程学硕调剂信息
已经有26人回复
-
欢迎加入课题组
已经有0人回复
-
散金币，求好运，祝面上顺利！
已经有34人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助用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/99
-
专业技术开发及第三方检测
+1/87
-
物理学 调剂
+1/84
-
华南师范大学（211）博士招生- 电子、自动化、机械、生物学、物理相关专业
+2/42
-
1
+1/37
-
2026年重庆工商大学管理科学与工程专业学术硕士招生调剂，风景总在柳暗花明处
+1/36
-
关于碳量子点抗菌的研究
+1/35
-
纺织科学与工程、材料化工方向招收研究生
+1/28
-
欢迎交流咨询
+1/28
-
新加坡国立大学药学系化学生物学课题组招PhD
+1/28
-
医学327分专硕调剂
+1/16
-
北京某研究院生物化学与分子生物学相关专业申请考核制博士招生（春季）
+1/15
-
2026年西南科技大学功能涂层课题组简介
+1/12
-
课题组招收环境及相关专业调剂硕士研究生（欢迎优秀学生加入）
+1/6
-
上海第二工业大学-朱大海课题组招生（过线就能调剂）
+1/5
-
重庆大学诚招2026年生物材料方向博士生
+1/4
-
静电纺丝：工艺原理、形貌调控与典型应用
+1/4
-
江苏南京林业大学-高分子功能材料和器件-招2名硕士研究生
+1/3
-
课题组招收环境及相关专业调剂硕士研究生（欢迎优秀学生加入）
+1/2
-
085404 270求调剂，b区求导师收留接受跨专业
+1/1
ahbiao(金币+10): 谢谢~~ 2011-03-30 08:10:16
|
本帖内容被屏蔽 |
2楼2011-03-29 10:10:20
3楼2011-03-30 08:12:07
ahbiao(金币+40): 2011-03-31 08:34:11
|
本帖内容被屏蔽 |
4楼2011-03-30 15:03:23
恩，谢谢！ 5楼2011-03-31 08:34:35
