泡点温度c语言编程,PR泡点温度计算小程序

#include

#include

main()

{

FILE *fp;

Double Pc[4]={3390000,5050000,7390000,22090000},Tc[4]={126,154,304,647

.3},Vc[4]={0.0000899,0.0000734,0.000094,0.000056},w[4]={0.040,0.021,0.225,0.344},P=100000,

M[4],qz,K1,K2,K3,K4,H,T,

x[4]={0.7,0.05,0.15,0.1}，y1=0,y2=0,y3=0,y4=0,y5,E,a,z,f[4],A,B，

q1,q2,q3,q4，aT[4],b[4],aT0=0,b0=0,aTc,aTr,Tr,k0,v;

int i,j;

double k[4][4];

for(T=100;T<900;T=T+0.01)

{

for(i=0;i<4;i++)

{

for(j=i+1;j<4;j++)

{

y5=pow(Vc[i],1.0/3)+pow(Vc[j],1.0/3);

k[i][j]=1-8*sqrt(Vc[i]*Vc[j])/pow(y5,3);

k[j][i]=k[i][j];

}

k[i][i]=0;

}

for(i=0;i<4;i++)