| 查看: 855 | 回复: 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 ] |
回复此楼» 猜你喜欢
如何从铁电相到顺电相。
已经有1人回复
-
各位大佬，求一份最新Wien2k版本作为学习用途
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有218人回复
-
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
已经有42人回复
-
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有24人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助用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/453
-
1,3,2-Dioxaphospholane, 2-(2,2,2-trifluoroethoxy)-
+1/87
-
3-羟基丙腈（CAS：109-78-4）
+1/86
-
3,4-二(三甲基硅氧基)-1-丁烯（CAS：89543-84-0）
+1/86
-
三(三甲基硅基)硼酸酯（CAS：4325-85-3）
+1/84
-
四(三甲基硅基)赤藓醇醚（CAS：18547-29-0 ）
+1/81
-
香港城市大学郑星课题组诚聘博士后 (生态学、环境科学、植物学、大气科学背景)
+1/78
-
别研究filecode了，祝都有一个好结果：）
+2/72
-
【坐标北京，诚征男友，结婚对象】
+5/42
-
中国石油大学（华东）计算机科学与技术学院/智能化新能源系统设计和开发方向 特聘副教
+2/28
-
中国石油大学（华东）材料科学与工程学院/新颖能源材料智能设计和开发方向 特聘副教授
+2/23
-
求相关合成条件
+1/16
-
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/16
-
海南大学化学院—功能分子器件团队2027博士招生+博士后招聘
+1/14
-
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
+1/9
-
中山大学李国强课题组招收博士生（柔性电子方向）
+1/5
-
海南大学化学院-科研助理招聘 （可读博，可与华中科技大学，天津大学联培）+博后招聘
+1/5
-
夏尽秋来，蓄力沉淀，静待硕果成文
+1/4
-
中国科学技术大学环境学院招聘招生
+1/4
-
悉尼科技大学（QS 第87）数据科学/AI招收2027年秋季入学全奖及CSC博士生
+1/1
3楼2011-03-30 08:12:07
ahbiao(金币+10): 谢谢~~ 2011-03-30 08:10:16
|
本帖内容被屏蔽 |
2楼2011-03-29 10:10:20
ahbiao(金币+40): 2011-03-31 08:34:11
|
本帖内容被屏蔽 |
4楼2011-03-30 15:03:23
恩，谢谢！ 5楼2011-03-31 08:34:35
