| 查看: 749 | 回复: 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 ] |
回复此楼» 猜你喜欢
070200求调剂，一志愿某211，288分
已经有19人回复
-
法国博士后职位
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有174人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有0人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有8人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有9人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助用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人回复
» 抢金币啦！回帖就可以得到:
-
长江大学化学与环境工程学院2026年硕士研究生调剂公告（理学、工学均可）
+5/770
-
江西科师大招收与湖南大学联合培养学硕1名
+1/488
-
武汉纺织大学材料学院士兵计划调剂信息
+1/96
-
长江大学化学与环境工程学院2026年硕士研究生调剂公告（理学、工学均可）
+1/90
-
招收化学、海洋科学、材料、环境专业调剂
+3/40
-
中国石油大学（北京）理学院2026年硕士研究生接收调剂生：物理、化学、材料等
+1/37
-
国家级人才课题组招研究生 （催化方向）调剂
+1/36
-
中国林科院木材所硕士调剂
+1/6
-
齐齐哈尔大学李莉课题组诚招2026级考研调剂生（学硕和专硕）
+1/6
-
王金凤教授团队
+1/5
-
渤海大学材料与环境工程学院发光及透明陶瓷材料、地聚物吸波混凝土方向接收研究生调剂
+1/5
-
烟台大学化学化工学院班庆福课题组招收调剂生
+1/5
-
欢迎调剂攀枝花学院化学工程专业，资源与环境专业，线上复试
+1/5
-
齐齐哈尔大学李莉课题组诚招2026级考研调剂生（学硕和专硕）
+1/5
-
海南医科大学智能医学课题组招收生物信息学方向专业硕士研究生！
+1/3
-
吉林工程技术师范学院 交叉学科研究院（08/07） 招收调剂
+1/3
-
找好工作来我这-广东唯一石化院校——资源与环境专硕招生
+1/3
-
上海交通大学化学化工学院张智涛课题组诚聘博士后
+1/2
-
福建理工大学材料学院招收专业代码08开头学术型硕士研究生
+1/2
-
招生
+1/2
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
