#include <stdio.h>
#include <math.h>
#define T1 631.5
#define I1 6118.9
#define Bj 6.1115
#define igs 485.93
#define D 36.1
#define XD 2.4
#define Qsmz 15121.3
#define Ilk 193
main()
{
	float i2,t2,t1,i1,Qs,I2,T2,Tpj,tpj,wy,td,tx,tb;
	float a0,ad,kr,ku,kps,a,tv,af,aff,K,Xt,Q,Xe;
	i2=igs+(XD*1688+Qsmz)/(1.02*D);
    t2=234+(i2-1008.8)/(1013.5-1008.8);
	t1=186;                /*假定*/
	i1=789.77;            /*臨時計算*/
	Qs=1.02*D*(i2-i1)/Bj;
	I2=I1-Qs/0.9936+0.02*193;
	if (3811<=I2<=4831)
		T2=400+100*(I2-3811)/(4831-3811);
	else if (4831<I2<=5868)
		T2=500+100*(I2-4831)/(5868-4831);
    Tpj=(T1+T2)/2;
	tpj=(t1+t2)/2;
	wy=Bj*6.462*(Tpj+273)/(11.88*273);
	a0=95;                /*臨時計算*/
	ad=a0*1*0.97*1;
	kr=(1.0-0.37*(Tpj+273)*0.001)*7.5/(1.0-0.37*0.813);
    ku=3.377*(pow(813,0.6667))/(pow((Tpj+273),0.6667));
	kps=(kr+ku)*0.0158;
	a=1.0-(exp(-kps));
	tb=tpj+60;
	af=75.6*a;
	aff=af*(1+0.4*pow((Tpj+273)*0.001,0.25)*pow((5000/900),0.07));
	K=(ad+aff)/(1+0.004*(ad+aff));
    td=T1-t2;
	tx=T2-t1;
	Xt=(td-tx)/(log(td/tx));
	Q=K*Xt*297.5/(1000*Bj);
	Xe=100*(Qs-Q)/Qs;
	printf("i2=%f,t2=%f,Qs=%f,I2=%f,T2=%f,Tpj=%f\n",i2,t2,Qs,I2,T2,Tpj);
	printf("tpj=%f,wy=%f,a0=%f,ad=%f,kr=%f,ku=%f\n",tpj,wy,a0,ad,kr,ku);
	printf("kps=%f,a=%f,tb=%f,af=%f,aff=%f,K=%f\n",kps,a,tb,af,aff,K);
	printf("Xt=%f,Q=%f,Xe=%f,t1=%f,i1=%f,\n",Xt,Q,Xe,t1,i1);
}