对分发,迭代法,牛顿法RK,SRK方程.docxVIP

对分法SRK方程#includestdio.h#includemath.h# define R 8.314float fun(float A,float B, float Z){floatfz;fz=pow(Z,3)-pow(Z,2)+(A-B-B*B)*Z-A*B;returnfz;}main(){inti=1; float p0,t0,w,p1,t1,h,Z0,Z1,Z2;float a,b,m,A,B,t2,at,y,y1,y2;a=b=0.0;printf(请输入临界参数:\n p0= );scanf(%f,p0);printf(请输入临界参数:\n t0= );scanf(%f,t0);printf(请输入偏心因子:\n w= );scanf(%f,w);printf(请输入实际温度:\n t1=);scanf(%f,t1);printf(请输入实际压力:\n p1=);scanf(%f,p1); m=0.480+1.574*w-0.176*w*w;t2=t1/t0;at=(1+m*(1-sqrt(t2)))*(1+m*(1-sqrt(t2))); a=0.42748*R*R*pow(t0,2)/p0*at; b=0.08664*R*t0/p0; A=a*p1/(R*R*pow(t1,2)); B=b*p1/(R*t1); Z1