| 查看: 853 | 回复: 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论文润色/翻译怎么收费?
已经有101人回复
-
南科大活性流体和软物质课题组诚招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人回复
» 抢金币啦！回帖就可以得到:
-
北京理工大学集成电路学院-柔性电子器件与智造研究所-2027级博士招生
+1/480
-
发榜前的“心理过山车”：你经历了吗？
+5/305
-
别研究filecode了，祝都有一个好结果：）
+2/94
-
1,3,2-Dioxaphospholane, 2-(2,2,2-trifluoroethoxy)-
+1/87
-
3-羟基丙腈（CAS：109-78-4）
+1/87
-
3,4-二(三甲基硅氧基)-1-丁烯（CAS：89543-84-0）
+1/87
-
1,2-双(三甲基硅氧基)苯（CAS：5075-52-5）
+1/85
-
招聘机械领域专利代理人
+1/80
-
请问谁有2025最新影响因子Excel表？金币求助！
+1/33
-
浙江大学医学院附属第一医院招聘科研助理
+1/28
-
哈工大（深圳）何自开团队诚招化学工程与技术专业博士生
+1/27
-
哈工大（深圳）何自开团队诚招2027级化学工程与技术专业博士生
+1/27
-
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/25
-
中山大学计算数学方向招2027届考核制博士生
+1/16
-
双一流高校招收博士（现代技术管理/环境技术经济管理）及硕士研究生（长期有效）
+1/13
-
中国石油大学（华东）材料科学与工程学院/新颖能源材料智能设计和开发方向 特聘副教授
+2/13
-
瑞典林雪平大学 高峰课题组 诚聘 氮化物基无机材料方向 博士后1名
+1/7
-
福建师范大学海峡柔性电子学院-2027年9月入学博士-考核制-电催、能源
+1/7
-
天津大学张小岩教授团队招收2027级推免生（光电子能源、热管理材料与器件、微纳传感）
+1/5
-
夏尽秋来，蓄力沉淀，静待硕果成文
+1/4
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
