程序如下:
Title "Type your title here";
Parameters p51,p52,p53,p54,p55,p56,p57,p58,p59,p510,p511,p61,p62,p63,p64,p65,p66,p67,p68,p69,p610,p611,p71,p72,p73,p74,p75,p76,p77,p78,p79,p710,p711;
Constant R=83.1451,Tc=647.096,p50=-10097698.7020424,p60=242005341.091914,p70=-1926309306.61686;
Variable V,pLG,VG,VL,T,y;
Function y=R*T*(1/V-(1/10)*((10*(p70+p71*(Tc/T-1)+p72*(Tc^2/T^2-1)+p73*(Tc^3/T^3-1)+p74*(Tc^4/T^4-1)+p75*(Tc^5/T^5-1)+p76*(Tc^6/T^6-1)+p77*(Tc^7/T^7-1)))*R*T*VG^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-(10*(p70+p71*(Tc/T-1)+p72*(Tc^2/T^2-1)+p73*(Tc^3/T^3-1)+p74*(Tc^4/T^4-1)+p75*(Tc^5/T^5-1)+
p76*(Tc^6/T^6-1)+p77*(Tc^7/T^7-1)))*R*T*VL^6-80*pLG*VL^7*VG^6+80*pLG*VG^7*VL^6+10*VL^5*R*T*VG^7-10*VL^7*R*T*VG^5-10*VL^5*pLG*VG^8+10*VL^8*VG^5*pLG-
10*VL^8*R*T*VG^4+10*VL^4*R*T*VG^8+10*VL^9*VG^4*pLG-10*VL^4*pLG*VG^9+10*VL^5*R*T*(p60+p61*(Tc/T-1)+p62*(Tc^2/T^2-1)+p63*(Tc^3/T^3-1)+p64*(Tc^4/T^4-1)+p65*(Tc^5/T^5-1)+
p66*(Tc^6/T^6-1)+p67*(Tc^7/T^7-1))*VG^2+10*VL^5*VG*R*T*(p70+p71*(Tc/T-1)+p72*(Tc^2/T^2-1)+p73*(Tc^3/T^3-1)+p74*(Tc^4/T^4-1)+p75*(Tc^5/T^5-1)+p76*(Tc^6/T^6-1)+
p77*(Tc^7/T^7-1))+10*VL^5*R*T*(p50+p51*(Tc/T-1)+p52*(Tc^2/T^2-1)+p53*(Tc^3/T^3-1)+p54*(Tc^4/T^4-1)+p55*(Tc^5/T^5-1)+p56*(Tc^6/T^6-1)+p57*(Tc^7/T^7-1))*VG^3-
10*VL^2*VG^4*R*T*(p70+p71*(Tc/T-1)+p72*(Tc^2/T^2-1)+p73*(Tc^3/T^3-1)+p74*(Tc^4/T^4-1)+p75*(Tc^5/T^5-1)+p76*(Tc^6/T^6-1)+p77*(Tc^7/T^7-1))-10*VL^3*VG^4*R*T*(p60+
p61*(Tc/T-1)+p62*(Tc^2/T^2-1)+p63*(Tc^3/T^3-1)+p64*(Tc^4/T^4-1)+p65*(Tc^5/T^5-1)+p66*(Tc^6/T^6-1)+p67*(Tc^7/T^7-1))+10*VL^4*VG^3*R*T*(p60+p61*(Tc/T-1)+p62*(Tc^2/T^2-1)+
p63*(Tc^3/T^3-1)+p64*(Tc^4/T^4-1)+p65*(Tc^5/T^5-1)+p66*(Tc^6/T^6-1)+p67*(Tc^7/T^7-1))+10*VL^4*VG^2*R*T*(p70+p71*(Tc/T-1)+p72*(Tc^2/T^2-1)+p73*(Tc^3/T^3-1)+p74*(Tc^4/T^4-1)+
p75*(Tc^5/T^5-1)+p76*(Tc^6/T^6-1)+p77*(Tc^7/T^7-1))-10*VL*VG^5*R*T*(p70+p71*(Tc/T-1)+p72*(Tc^2/T^2-1)+p73*(Tc^3/T^3-1)+p74*(Tc^4/T^4-1)+p75*(Tc^5/T^5-1)+p76*(Tc^6/T^6-1)+
p77*(Tc^7/T^7-1))-10*VL^3*VG^5*R*T*(p50+p51*(Tc/T-1)+p52*(Tc^2/T^2-1)+p53*(Tc^3/T^3-1)+p54*(Tc^4/T^4-1)+p55*(Tc^5/T^5-1)+p56*(Tc^6/T^6-1)+p57*(Tc^7/T^7-1))-
10*VL^2*VG^5*R*T*(p60+p61*(Tc/T-1)+p62*(Tc^2/T^2-1)+p63*(Tc^3/T^3-1)+p64*(Tc^4/T^4-1)+p65*(Tc^5/T^5-1)+p66*(Tc^6/T^6-1)+p67*(Tc^7/T^7-1))+(5*(p50+p51*(Tc/T-1)+p52*(Tc^2/T^2-
1)+p53*(Tc^3/T^3-1)+p54*(Tc^4/T^4-1)+p55*(Tc^5/T^5-1)+p56*(Tc^6/T^6-1)+p57*(Tc^7/T^7-1)))*R*T*VL^2*VG^6-(8*(p60+p61*(Tc/T-1)+p62*(Tc^2/T^2-1)+p63*(Tc^3/T^3-1)+p64*(Tc^4/T^4-
1)+p65*(Tc^5/T^5-1)+p66*(Tc^6/T^6-1)+p67*(Tc^7/T^7-1)))*R*T*VL^6*VG-(5*(p50+p51*(Tc/T-1)+p52*(Tc^2/T^2-1)+p53*(Tc^3/T^3-1)+p54*(Tc^4/T^4-1)+p55*(Tc^5/T^5-1)+p56*(Tc^6/T^6-
1)+p57*(Tc^7/T^7-1)))*R*T*VL^6*VG^2+(8*(p60+p61*(Tc/T-1)+p62*(Tc^2/T^2-1)+p63*(Tc^3/T^3-1)+p64*(Tc^4/T^4-1)+p65*(Tc^5/T^5-1)+p66*(Tc^6/T^6-1)+
p67*(Tc^7/T^7-1)))*R*T*VL*VG^6)/(VG^4*VL^4*R*T*(-VL^3+3*VL^2*VG+VG^3-3*VL*VG^2)*V^2)+(1/10)*(-(8*(p60+p61*(Tc/T-1)+p62*(Tc^2/T^2-1)+p63*(Tc^3/T^3-1)+p64*(Tc^4/T^4-1)+
p65*(Tc^5/T^5-1)+p66*(Tc^6/T^6-1)+p67*(Tc^7/T^7-1)))*R*T*VL^6*VG^2+40*VL^5*R*T*VG^8-40*VL^8*R*T*VG^5+80*pLG*VG^8*VL^6-40*VL^5*pLG*VG^9-80*VL^8*VG^6*pLG+40*VL^9*VG^5*pLG+
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*(p70+p71*(Tc/T-1)+p72*(Tc^2/T^2-
1)+p73*(Tc^3/T^3-1)+p74*(Tc^4/T^4-1)+p75*(Tc^5/T^5-1)+p76*(Tc^6/T^6-1)+p77*(Tc^7/T^7-1)))*R*T*VG^7-(10*(p70+p71*(Tc/T-1)+p72*(Tc^2/T^2-1)+p73*(Tc^3/T^3-1)+p74*(Tc^4/T^4-1)+
p75*(Tc^5/T^5-1)+p76*(Tc^6/T^6-1)+p77*(Tc^7/T^7-1)))*R*T*VL^7+(5*(p50+p51*(Tc/T-1)+p52*(Tc^2/T^2-1)+p53*(Tc^3/T^3-1)+p54*(Tc^4/T^4-1)+p55*(Tc^5/T^5-1)+p56*(Tc^6/T^6-1)+
p57*(Tc^7/T^7-1)))*R*T*VL^2*VG^7-(5*(p50+p51*(Tc/T-1)+p52*(Tc^2/T^2-1)+p53*(Tc^3/T^3-1)+p54*(Tc^4/T^4-1)+p55*(Tc^5/T^5-1)+p56*(Tc^6/T^6-1)+p57*(Tc^7/T^7-1)))*R*T*VL^7*VG^2-
10*VG*(p70+p71*(Tc/T-1)+p72*(Tc^2/T^2-1)+p73*(Tc^3/T^3-1)+p74*(Tc^4/T^4-1)+p75*(Tc^5/T^5-1)+p76*(Tc^6/T^6-1)+p77*(Tc^7/T^7-1))*R*T*VL^6+40*VL^5*R*T*(p60+p61*(Tc/T-1)+
p62*(Tc^2/T^2-1)+p63*(Tc^3/T^3-1)+p64*(Tc^4/T^4-1)+p65*(Tc^5/T^5-1)+p66*(Tc^6/T^6-1)+p67*(Tc^7/T^7-1))*VG^3+40*VL^5*VG^2*R*T*(p70+p71*(Tc/T-1)+p72*(Tc^2/T^2-1)+p73*(Tc^3/T^3-
1)+p74*(Tc^4/T^4-1)+p75*(Tc^5/T^5-1)+p76*(Tc^6/T^6-1)+p77*(Tc^7/T^7-1))+40*VL^5*R*T*(p50+p51*(Tc/T-1)+p52*(Tc^2/T^2-1)+p53*(Tc^3/T^3-1)+p54*(Tc^4/T^4-1)+p55*(Tc^5/T^5-1)+
p56*(Tc^6/T^6-1)+p57*(Tc^7/T^7-1))*VG^4+10*VL*VG^6*R*T*(p70+p71*(Tc/T-1)+p72*(Tc^2/T^2-1)+p73*(Tc^3/T^3-1)+p74*(Tc^4/T^4-1)+p75*(Tc^5/T^5-1)+p76*(Tc^6/T^6-1)+p77*(Tc^7/T^7-
1))+5*VL^3*VG^6*R*T*(p50+p51*(Tc/T-1)+p52*(Tc^2/T^2-1)+p53*(Tc^3/T^3-1)+p54*(Tc^4/T^4-1)+p55*(Tc^5/T^5-1)+p56*(Tc^6/T^6-1)+p57*(Tc^7/T^7-1))+8*VL^2*VG^6*R*T*(p60+p61*(Tc/T-
1)+p62*(Tc^2/T^2-1)+p63*(Tc^3/T^3-1)+p64*(Tc^4/T^4-1)+p65*(Tc^5/T^5-1)+p66*(Tc^6/T^6-1)+p67*(Tc^7/T^7-1))-40*VL^2*VG^5*R*T*(p70+p71*(Tc/T-1)+p72*(Tc^2/T^2-1)+
p73*(Tc^3/T^3-1)+p74*(Tc^4/T^4-1)+p75*(Tc^5/T^5-1)+p76*(Tc^6/T^6-1)+p77*(Tc^7/T^7-1))-40*VL^4*VG^5*R*T*(p50+p51*(Tc/T-1)+p52*(Tc^2/T^2-1)+p53*(Tc^3/T^3-1)+p54*(Tc^4/T^4-1)+
p55*(Tc^5/T^5-1)+p56*(Tc^6/T^6-1)+p57*(Tc^7/T^7-1))-40*VL^3*VG^5*R*T*(p60+p61*(Tc/T-1)+p62*(Tc^2/T^2-1)+p63*(Tc^3/T^3-1)+p64*(Tc^4/T^4-1)+p65*(Tc^5/T^5-1)+p66*(Tc^6/T^6-1)+
p67*(Tc^7/T^7-1))-(5*(p50+p51*(Tc/T-1)+p52*(Tc^2/T^2-1)+p53*(Tc^3/T^3-1)+p54*(Tc^4/T^4-1)+p55*(Tc^5/T^5-1)+p56*(Tc^6/T^6-1)+p57*(Tc^7/T^7-1)))*R*T*VL^6*VG^3+(8*(p60+
p61*(Tc/T-1)+p62*(Tc^2/T^2-1)+p63*(Tc^3/T^3-1)+p64*(Tc^4/T^4-1)+p65*(Tc^5/T^5-1)+p66*(Tc^6/T^6-1)+p67*(Tc^7/T^7-1)))*R*T*VL*VG^7-(8*(p60+p61*(Tc/T-1)+p62*(Tc^2/T^2-1)+
p63*(Tc^3/T^3-1)+p64*(Tc^4/T^4-1)+p65*(Tc^5/T^5-1)+p66*(Tc^6/T^6-1)+p67*(Tc^7/T^7-1)))*R*T*VL^7*VG)/(VL^4*VG^4*R*T*(-VL^3+3*VL^2*VG+VG^3-3*VL*VG^2)*V^3)-(1/10)*((15*(p50+
p51*(Tc/T-1)+p52*(Tc^2/T^2-1)+p53*(Tc^3/T^3-1)+p54*(Tc^4/T^4-1)+p55*(Tc^5/T^5-1)+p56*(Tc^6/T^6-1)+p57*(Tc^7/T^7-1)))*R*T*VL^2*VG^6+(18*(p60+p61*(Tc/T-1)+p62*(Tc^2/T^2-1)+
p63*(Tc^3/T^3-1)+p64*(Tc^4/T^4-1)+p65*(Tc^5/T^5-1)+p66*(Tc^6/T^6-1)+p67*(Tc^7/T^7-1)))*R*T*VL*VG^6+(20*(p70+p71*(Tc/T-1)+p72*(Tc^2/T^2-1)+p73*(Tc^3/T^3-1)+p74*(Tc^4/T^4-1)+
p75*(Tc^5/T^5-1)+p76*(Tc^6/T^6-1)+p77*(Tc^7/T^7-1)))*R*T*VG^6-(15*(p50+p51*(Tc/T-1)+p52*(Tc^2/T^2-1)+p53*(Tc^3/T^3-1)+p54*(Tc^4/T^4-1)+p55*(Tc^5/T^5-1)+p56*(Tc^6/T^6-1)+
p57*(Tc^7/T^7-1)))*R*T*VL^6*VG^2-(18*(p60+p61*(Tc/T-1)+p62*(Tc^2/T^2-1)+p63*(Tc^3/T^3-1)+p64*(Tc^4/T^4-1)+p65*(Tc^5/T^5-1)+p66*(Tc^6/T^6-1)+p67*(Tc^7/T^7-1)))*R*T*VL^6*VG-
(20*(p70+p71*(Tc/T-1)+p72*(Tc^2/T^2-1)+p73*(Tc^3/T^3-1)+p74*(Tc^4/T^4-1)+p75*(Tc^5/T^5-1)+p76*(Tc^6/T^6-1)+p77*(Tc^7/T^7-1)))*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-90*pLG*VL^7*VG^6+90*pLG*VG^7*VL^6+30*VL^5*R*T*VG^7-30*VL^7*R*T*VG^5-30*VL^5*pLG*VG^8+30*VL^8*VG^5*pLG+30*VL^5*R*T*(p60+p61*(Tc/T-1)+
p62*(Tc^2/T^2-1)+p63*(Tc^3/T^3-1)+p64*(Tc^4/T^4-1)+p65*(Tc^5/T^5-1)+p66*(Tc^6/T^6-1)+p67*(Tc^7/T^7-1))*VG^2+30*VL^5*VG*R*T*(p70+p71*(Tc/T-1)+p72*(Tc^2/T^2-1)+p73*(Tc^3/T^3-
1)+p74*(Tc^4/T^4-1)+p75*(Tc^5/T^5-1)+p76*(Tc^6/T^6-1)+p77*(Tc^7/T^7-1))+30*VL^5*R*T*(p50+p51*(Tc/T-1)+p52*(Tc^2/T^2-1)+p53*(Tc^3/T^3-1)+p54*(Tc^4/T^4-1)+p55*(Tc^5/T^5-1)+
p56*(Tc^6/T^6-1)+p57*(Tc^7/T^7-1))*VG^3-30*VL*VG^5*R*T*(p70+p71*(Tc/T-1)+p72*(Tc^2/T^2-1)+p73*(Tc^3/T^3-1)+p74*(Tc^4/T^4-1)+p75*(Tc^5/T^5-1)+p76*(Tc^6/T^6-1)+p77*(Tc^7/T^7-
1))-30*VL^3*VG^5*R*T*(p50+p51*(Tc/T-1)+p52*(Tc^2/T^2-1)+p53*(Tc^3/T^3-1)+p54*(Tc^4/T^4-1)+p55*(Tc^5/T^5-1)+p56*(Tc^6/T^6-1)+p57*(Tc^7/T^7-1))-30*VL^2*VG^5*R*T*(p60+
p61*(Tc/T-1)+p62*(Tc^2/T^2-1)+p63*(Tc^3/T^3-1)+p64*(Tc^4/T^4-1)+p65*(Tc^5/T^5-1)+p66*(Tc^6/T^6-1)+p67*(Tc^7/T^7-1)))/(R*T*VL^3*VG^3*(-VL^3+3*VL^2*VG+VG^3-3*VL*VG^2)*V^4)+
(p50+p51*(Tc/T-1)+p52*(Tc^2/T^2-1)+p53*(Tc^3/T^3-1)+p54*(Tc^4/T^4-1)+p55*(Tc^5/T^5-1)+p56*(Tc^6/T^6-1)+p57*(Tc^7/T^7-1))/V^5+(p60+p61*(Tc/T-1)+p62*(Tc^2/T^2-1)+
p63*(Tc^3/T^3-1)+p64*(Tc^4/T^4-1)+p65*(Tc^5/T^5-1)+p66*(Tc^6/T^6-1)+p67*(Tc^7/T^7-1))/V^6+(p70+p71*(Tc/T-1)+p72*(Tc^2/T^2-1)+p73*(Tc^3/T^3-1)+p74*(Tc^4/T^4-1)+
p75*(Tc^5/T^5-1)+p76*(Tc^6/T^6-1)+p77*(Tc^7/T^7-1))/V^7);
Data;
V,pLG,VG,VL,T,y=
45421047.33 0.006116548 3710987.122 18.01900658 273.16 0.0005
22709507.93 0.006116548 3710987.122 18.01900658 273.16 0.001
15138991.54 0.006116548 3710987.122 18.01900658 273.16 0.0015
11353729.76 0.006116548 3710987.122 18.01900658 273.16 0.002
9082568.923 0.006116548 3710987.122 18.01900658 273.16 0.0025
7568457.791 0.006116548 3710987.122 18.01900658 273.16 0.003
6486945.84 0.006116548 3710987.122 18.01900658 273.16 0.0035
5675807.806 0.006116548 3710987.122 18.01900658 273.16 0.004
5044918.542 0.006116548 3710987.122 18.01900658 273.16 0.0045
4540202.959 0.006116548 3710987.122 18.01900658 273.16 0.005
4127249.636 0.006116548 3710987.122 18.01900658 273.16 0.0055
3783117.628 0.006116548 3710987.122 18.01900658 273.16 0.006
。。。。。。。。。。。。。。。。。。。 |