化工热力学公式 C语言实现 （一） 迭代

化工热力学公式 C语言实现 （一） 迭代

#include <stdio.h>
#include <math.h>
int main()
{
	double Vm,B,C,R,T,P,Tc,Pc,w,a,b,d,e,m;
	R=8.314;
	float V[500];
	int i=0,c,yn;
	printf ("请选择您计算的方程：1 维里方程 2 R-K方程 3 SRK方程");
	scanf ("%d",&c);
	if (c==1)
	{
	 printf("请依次输入 B C T P值");
	 scanf("%lf%lf%lf%lf",&B,&C,&T,&P) ;
	 V[0]=R*T/P; 
	 do{
		i++;
 		V[i]=R*T/P*(1+B/V[i-1]+C/pow(V[i-1],2));
 	    printf("第%d次迭代的值是%lf\n:",i,V[i]);
 	   } while((V[i]-V[i-1])<10e-4);
	}
	if (c==2)
	{
	 printf("请依次输入 Tc Pc T P值");
	 scanf("%lf%lf%lf%lf",&Tc,&Pc,&T,&P) ;
	 V[0]=R*T/P;
	 printf("V0是%f",V[0]) ;
	 a=0.42748*pow(R,2)*pow(Tc,2.5)/Pc;
	 b=0.08664*pow(R,1)*pow(Tc,1)/Pc;
	 printf("修正因子a是%f b是%f",a,b);
	 printf("是否采用牛顿法 :1 yes 2 no") ;
	 scanf("%d",&yn);
        if(yn==2)