求教高人1stopt的拟合问题，总是提示出错

wlxydgq

学习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 编辑 ]

dingd

把“Data;”前的汉字“方程项很多”去掉，运行一下没问题。

wlxydgq

:@L 怎么会呢？那我运行就说：读 数据有问题。

wlxydgq

问题解决了，感谢dingd！！！！！！！！！！！！！！！