#include <stdio.h>
#include <math.h>
#define Bj 6.1115
#define T1 1094.4
#define I1 10966.6
main()
{
	float T2,I2,Qy,Tpj,wy,ad,kr,ku,kps,a,t,th,af,a1,k,xt,Qc,Xe;
	T2=1008;              /*假定*/
	t=255;
	th=t+80;
	I2=9926+0.01*(T2-1000)*(11028-9926);
	Qy=0.9936*(I1-I2);
	Tpj=(T1+T2)/2;
	wy=Bj*6.264*(Tpj+273)/(28.1*273);
	ad=55.75*0.88*0.95*0.93;
	kr=(1.0-0.37*(Tpj+273)*0.001)*2.022/(1.0-0.37*1.2927);
	ku=2.576*(pow(1323,0.6667))/(pow((Tpj+273),0.6667));
	kps=(kr+ku)*0.117;
    a=1.0-exp(-kps);
    af=a*197.7;
	a1=ad+af;
	k=a1*0.65;
	xt=Tpj-255;
	Qc=k*xt*89.6/(1000*Bj);
	Xe=100*(Qy-Qc)/Qy;
	printf("T2=%4.1f,I2=%6.1f,Qy=%3.1f,Tpj=%4.1f,wy=%1.2f,ad=%f\n",T2,I2,Qy,Tpj,wy,ad);
	printf("kr=%1.3f,ku=%1.3f,kps=%1.3f,a=%1.3f,th=%3.1f\n",kr,ku,kps,a,th);
	printf("af=%2.2f,a1=%2.2f,k=%2.2f,xt=%3.1f,Qc=%3.1f,Xe=%1.1f\n",af,a1,k,xt,Qc,Xe);
}