马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
学习1stopt不久,自己也算了一些东西,很方便,但是现在这个程序不行了,请大家帮我看看是什么问题呢?程序如下:
Title "Type your title here";
Parameters p51,p52,p53,p54,p55,p61,p62,p63,p64,p65,p71,p72,p73,p74,p75;
Constant R=83.1451,Tc=647.096;
Variable V,pLG,VG,VL,T,y;
Function y=R*T*(1/V-(1/10)*(60*ln(pLG*VL/(R*T))*R*T*VL^6*VG^6-60*ln(pLG*VG/(R*T))*R*T*VL^6*VG^6
+10*VL^5*R*T*(p61+p62*T/Tc+p63*T^2/Tc^2+p64*Tc/T+p65*Tc^2/T^2)*VG^2
+10*VL^5*VG*R*T*(p71+p72*T/Tc+p73*T^2/Tc^2+p74*Tc/T+p75*Tc^2/T^2)
+10*VL^5*R*T*(p51+p52*T/Tc+p53*T^2/Tc^2+p54*Tc/T+p55*Tc^2/T^2)*VG^3
-10*VL^2*VG^4*R*T*(p71+p72*T/Tc+p73*T^2/Tc^2+p74*Tc/T+p75*Tc^2/T^2)
-10*VL^3*VG^4*R*T*(p61+p62*T/Tc+p63*T^2/Tc^2+p64*Tc/T+p65*Tc^2/T^2)
+10*VL^4*VG^3*R*T*(p61+p62*T/Tc+p63*T^2/Tc^2+p64*Tc/T+p65*Tc^2/T^2)
+10*VL^4*VG^2*R*T*(p71+p72*T/Tc+p73*T^2/Tc^2+p74*Tc/T+p75*Tc^2/T^2)
-10*VL*VG^5*R*T*(p71+p72*T/Tc+p73*T^2/Tc^2+p74*Tc/T+p75*Tc^2/T^2)
-10*VL^3*VG^5*R*T*(p51+p52*T/Tc+p53*T^2/Tc^2+p54*Tc/T+p55*Tc^2/T^2)
-10*VL^2*VG^5*R*T*(p61+p62*T/Tc+p63*T^2/Tc^2+p64*Tc/T+p65*Tc^2/T^2)
-10*VL^5*pLG*VG^8-80*pLG*VL^7*VG^6+80*pLG*VG^7*VL^6+10*VL^8*VG^5*pLG
+10*VL^9*VG^4*pLG-10*VL^4*pLG*VG^9+10*VL^5*R*T*VG^7-10*VL^7*R*T*VG^5
-10*VL^8*R*T*VG^4+10*VL^4*R*T*VG^8-(8*(p61+p62*T/Tc+p63*T^2/Tc^2
+p64*Tc/T+p65*Tc^2/T^2))*R*T*VL^6*VG-(5*(p51+p52*T/Tc+p53*T^2/Tc^2
+p54*Tc/T+p55*Tc^2/T^2))*R*T*VL^6*VG^2+(8*(p61+p62*T/Tc+p63*T^2/Tc^2
+p64*Tc/T+p65*Tc^2/T^2))*R*T*VL*VG^6+(5*(p51+p52*T/Tc+p53*T^2/Tc^2
+p54*Tc/T+p55*Tc^2/T^2))*R*T*VL^2*VG^6+(10*(p71+p72*T/Tc+p73*T^2/Tc^2
+p74*Tc/T+p75*Tc^2/T^2))*R*T*VG^6-(10*(p71+p72*T/Tc+p73*T^2/Tc^2+p74*Tc/T
+p75*Tc^2/T^2))*R*T*VL^6)/(VG^4*VL^4*R*T*(-VL^3+3*VL^2*VG+VG^3-3*VL*VG^2)*V^2)
+(1/10)*(60*ln(pLG*VL/(R*T))*R*T*VL^6*VG^7-60*ln(pLG*VG/(R*T))*R*T*VL^6*VG^7
+60*ln(pLG*VL/(R*T))*R*T*VL^7*VG^6-60*ln(pLG*VG/(R*T))*R*T*VL^7*VG^6
-10*VG*(p71+p72*T/Tc+p73*T^2/Tc^2+p74*Tc/T+p75*Tc^2/T^2)*R*T*VL^6
+40*VL^5*R*T*(p61+p62*T/Tc+p63*T^2/Tc^2+p64*Tc/T+p65*Tc^2/T^2)*VG^3
+40*VL^5*VG^2*R*T*(p71+p72*T/Tc+p73*T^2/Tc^2+p74*Tc/T+p75*Tc^2/T^2)
+40*VL^5*R*T*(p51+p52*T/Tc+p53*T^2/Tc^2+p54*Tc/T+p55*Tc^2/T^2)*VG^4
+10*VL*VG^6*R*T*(p71+p72*T/Tc+p73*T^2/Tc^2+p74*Tc/T+p75*Tc^2/T^2)
+5*VL^3*VG^6*R*T*(p51+p52*T/Tc+p53*T^2/Tc^2+p54*Tc/T+p55*Tc^2/T^2)
+8*VL^2*VG^6*R*T*(p61+p62*T/Tc+p63*T^2/Tc^2+p64*Tc/T+p65*Tc^2/T^2)
-40*VL^2*VG^5*R*T*(p71+p72*T/Tc+p73*T^2/Tc^2+p74*Tc/T+p75*Tc^2/T^2)
-40*VL^4*VG^5*R*T*(p51+p52*T/Tc+p53*T^2/Tc^2+p54*Tc/T+p55*Tc^2/T^2)
-40*VL^3*VG^5*R*T*(p61+p62*T/Tc+p63*T^2/Tc^2+p64*Tc/T+p65*Tc^2/T^2)
+80*pLG*VG^8*VL^6-80*VL^8*VG^6*pLG+40*VL^9*VG^5*pLG-40*VL^5*pLG*VG^9
+40*VL^5*R*T*VG^8-40*VL^8*R*T*VG^5+(8*(p61+p62*T/Tc+p63*T^2/Tc^2+p64*Tc/T
+p65*Tc^2/T^2))*R*T*VL*VG^7-(8*(p61+p62*T/Tc+p63*T^2/Tc^2
+p64*Tc/T+p65*Tc^2/T^2))*R*T*VL^6*VG^2-(5*(p51+p52*T/Tc+p53*T^2/Tc^2
+p54*Tc/T+p55*Tc^2/T^2))*R*T*VL^6*VG^3+(10*(p71+p72*T/Tc+p73*T^2/Tc^2
+p74*Tc/T+p75*Tc^2/T^2))*R*T*VG^7-(8*(p61+p62*T/Tc+p63*T^2/Tc^2
+p64*Tc/T+p65*Tc^2/T^2))*R*T*VL^7*VG-(5*(p51+p52*T/Tc+p53*T^2/Tc^2
+p54*Tc/T+p55*Tc^2/T^2))*R*T*VL^7*VG^2+(5*(p51+p52*T/Tc+p53*T^2/Tc^2
+p54*Tc/T+p55*Tc^2/T^2))*R*T*VL^2*VG^7-(10*(p71+p72*T/Tc+p73*T^2/Tc^2
+p74*Tc/T+p75*Tc^2/T^2))*R*T*VL^7)/(VL^4*VG^4*R*T*(-VL^3+3*VL^2*VG+VG^3
-3*VL*VG^2)*V^3)-(1/10)*((15*(p51+p52*T/Tc+p53*T^2/Tc^2+p54*Tc/T
+p55*Tc^2/T^2))*R*T*VL^2*VG^6+(18*(p61+p62*T/Tc+p63*T^2/Tc^2+p64*Tc/T
+p65*Tc^2/T^2))*R*T*VL*VG^6+(20*(p71+p72*T/Tc+p73*T^2/Tc^2+p74*Tc/T
+p75*Tc^2/T^2))*R*T*VG^6-(15*(p51+p52*T/Tc+p53*T^2/Tc^2+p54*Tc/T
+p55*Tc^2/T^2))*R*T*VL^6*VG^2-(18*(p61+p62*T/Tc+p63*T^2/Tc^2+p64*Tc/T
+p65*Tc^2/T^2))*R*T*VL^6*VG-(20*(p71+p72*T/Tc+p73*T^2/Tc^2+p74*Tc/T
+p75*Tc^2/T^2))*R*T*VL^6+60*ln(pLG*VL/(R*T))*R*T*VL^6*VG^6
-60*ln(pLG*VG/(R*T))*R*T*VL^6*VG^6-30*VL^5*pLG*VG^8-90*pLG*VL^7*VG^6
+90*pLG*VG^7*VL^6+30*VL^5*R*T*(p61+p62*T/Tc+p63*T^2/Tc^2+p64*Tc/T
+p65*Tc^2/T^2)*VG^2+30*VL^5*VG*R*T*(p71+p72*T/Tc+p73*T^2/Tc^2+p74*Tc/T
+p75*Tc^2/T^2)+30*VL^5*R*T*(p51+p52*T/Tc+p53*T^2/Tc^2+p54*Tc/T
+p55*Tc^2/T^2)*VG^3-30*VL*VG^5*R*T*(p71+p72*T/Tc+p73*T^2/Tc^2
+p74*Tc/T+p75*Tc^2/T^2)-30*VL^3*VG^5*R*T*(p51+p52*T/Tc+p53*T^2/Tc^2
+p54*Tc/T+p55*Tc^2/T^2)-30*VL^2*VG^5*R*T*(p61+p62*T/Tc+p63*T^2/Tc^2
+p64*Tc/T+p65*Tc^2/T^2)+30*VL^5*R*T*VG^7-30*VL^7*R*T*VG^5
+30*VL^8*VG^5*pLG)/(R*T*VL^3*VG^3*(-VL^3+3*VL^2*VG+VG^3-3*VL*VG^2)*V^4)
+(p51+p52*T/Tc+p53*T^2/Tc^2+p54*Tc/T+p55*Tc^2/T^2)/V^5+(p61+p62*T/Tc
+p63*T^2/Tc^2+p64*Tc/T+p65*Tc^2/T^2)/V^6+(p71+p72*T/Tc+p73*T^2/Tc^2
+p74*Tc/T+p75*Tc^2/T^2)/V^7); 方程项很多
Data;
V,pLG,VG,VL,T,y=
227113342.7 0.006116548 3710987.122 18.01900658 273.16 0.0001
113555658.6 0.006116548 3710987.122 18.01900658 273.16 0.0002
75703097.02 0.006116548 3710987.122 18.01900658 273.16 0.0003
56776816.07 0.006116548 3710987.122 18.01900658 273.16 0.0004
45421047.33 0.006116548 3710987.122 18.01900658 273.16 0.0005
37850534.66 0.006116548 3710987.122 18.01900658 273.16 0.0006
32443025.44 0.006116548 3710987.122 18.01900658 273.16 0.0007
28387393.35 0.006116548 3710987.122 18.01900658 273.16 0.0008
25233012.66 0.006116548 3710987.122 18.01900658 273.16 0.0009
22709507.93 0.006116548 3710987.122 18.01900658 273.16 0.001
20644822.06 0.006116548 3710987.122 18.01900658 273.16 0.0011
18924250.32 0.006116548 3710987.122 18.01900658 273.16 0.0012
(数据很多,这里就略了)
一运行,就提示“读数据出错”,都急死了,不知道到底是什么问题,初步估计是函数的问题,可是折腾了半天,还是没有找到点上。
请高人帮我找找,到底是什么问题。非常感谢!!!!!!!!!!
[ 本帖最后由 wlxydgq 于 2010-1-5 18:03 编辑 ] |