| 查看: 622 | 回复: 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 |
回复此楼» 猜你喜欢
导师想让我从独立一作变成了共一第一
已经有8人回复
-
博士读完未来一定会好吗
已经有23人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有11人回复
-
读博
已经有4人回复
-
JMPT 期刊投稿流程
已经有4人回复
-
心脉受损
已经有5人回复
-
Springer期刊投稿求助
已经有4人回复
-
小论文投稿
已经有3人回复
-
申请2026年博士
已经有6人回复
5楼2016-10-18 19:46:25
2楼2016-10-18 16:55:07
3楼2016-10-18 18:44:51
4楼2016-10-18 19:24:11