| 查看: 841 | 回复: 6 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
求助，matlab定义的循环函数错在哪里 已有1人参与
|
||
大神帮我看看，这个函数的定义错在哪里了。 function concen= LibrConcen(P,T ) x1=40; x2=70; A0=-2.00755; B0=124.937; A1=0.16976; B1=-7.7165; A2=-0.003133362; B2=0.152286; A3=0.0000197668; B3=-0.0079509; C=7.05; D=-1603.54; E=-104095.5; delta=D^2-4*E*(C-log10(P)); delta=sqrt(delta); temref=-2*E/(D+delta)-273.15; while (1) x3=0.5*(x1+x2); f1=B0+temref*A0-T+(B1+temref*A1)*x1+(B2+temref*A2)*x1^2+(B3+temref*A3)*x1^3; f2=B0+temref*A0-T+(B1+temref*A1)*x2+(B2+temref*A2)*x2^2+(B3+temref*A3)*x2^3; f3=B0+temref*A0-T+(B1+temref*A1)*x3+(B2+temref*A2)*x3^2+(B3+temref*A3)*x3^3; if (f1*f3>0) x1=x3; else x2=x3; end if(abs(f3)<=10^(-6)) break else continue end end concen=x3/100; end |
回复此楼» 猜你喜欢
40+中年人失业后，还能找到好工作吗？
已经有14人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有5人回复
-
申请科研项目纯靠实力能成功吗？
已经有12人回复
-
用NBS溴化时溶剂除了CCl4外还有什么
已经有6人回复
-
极限拉扯，如何抉择
已经有15人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
7楼2016-10-19 08:55:15
2楼2016-10-18 16:55:07
3楼2016-10-18 18:44:51
4楼2016-10-18 19:24:11
10