10下级空预.c

这个源码是从其它网站上转载来的

#include <stdio.h>
#include <math.h>
#define T1 252.8
#define I1 2444.7
#define t2 166              /*暂时假定*/
#define t1 30.0
#define i1 193.0
#define Bj 6.1115
main()
{
	float Qk,Xt1,Xt2,Xt,Xtt,td,tx,P,R,Q,Xe,l,i2;
    float tpj,T2,I2,ipj,Tpj,wy,a0,ay,wk,a0q,ak,k,kk;
	i2=643+0.01*(t2-100)*(1294-643);
    tpj=(t1+t2)/2;
	ipj=(i1+i2)/2;
	Qk=1.105*(i2-i1);
	I2=I1-Qk/0.9936+0.03*ipj;
	T2=100+100*(I2-967)/991;
	Tpj=(T1+T2)/2;
    wy=Bj*6.832*(Tpj+273)/(6.77*273);
    a0=34.9;              /*临时计算*/
	ay=a0*1*1.11;
	wk=5.76*Bj*(tpj+273)/(5.5761*(110.2+273));
	a0q=66;               /*临时计算*/
	ak=a0q*1.11*0.99*1;
	k=0.65*ak*ay/(ak+ay);
	Xt1=T1-t2;
	Xt2=T2-t1;
	Xt=(Xt2-Xt1)/(log(Xt2/Xt1));
	td=t2-t1;
	tx=T1-T2;
	P=tx/(T1-t1);
	R=td/tx;
	l=0.89;               /*临时计算*/
    Xtt=l*Xt;
	Q=k*Xtt*3913/(1000*Bj);
	Xe=100.0*(Qk-Q)/Qk;
	printf("tpj=%f,ipj=%f,ak=%f,ay=%f,i2=%f\n",tpj,ipj,ak,ay,i2);
	printf("I2=%f,Tpj=%f,Qk=%f,T2=%f,wy=%f,wk=%f,k=%f\n",I2,Tpj,Qk,T2,wy,wk,k);
	printf("ay=%f,ak=%f\n",ay,ak);
	printf("Xt1=%3.1f,Xt2=%3.1f,Xt=%3.1f,td=%3.1f,tx=%3.1f,P=%1.3f,R=%1.3f,Xtt=%3.1f,Q=%6.3f,Xe=%f\n",Xt1,Xt2,Xt,td,tx,P,R,Xtt,Q,Xe);
}