#include <stdio.h>
#include <math.h>
#define Bj 6.1115
#define Qf 167.6
#define T1 1008
#define I1 10014.2
#define t2 450
#define i2 3332.4
#define D 36.1
main()
{
	float t1,h1,h2,i1,T2,I2,Qq,Q1s,Q3s,Qy,Tpj,wy,tpj,v,wq,kr,ku,kps,a,th,af,aff,a1,a2,K,td,tx,Xtnl,Xt,Q;
	float Q1,Q3,Q13,Xt1,Xt3,Qe,Xe,Xez,id,v1,v2,v3;
	float a0,a0c,a0s,adc,ads,a0q,ad;
	t1=320;                              /*假定*/
	h1=3014.3+0.1*(t1-320)*(-3014.3+3040.5);
	h2=2984.0+0.1*(t1-320)*(-2984.0+3012.6);
	i1=h1+0.1*(h2-h1);
	Qq=D*(i2-i1)/Bj;
	Q1s=55;                                /*假定*/
	Q3s=230;                               /*假定*/
    Q13=Q1s+Q3s;
	Qy=Qq-Qf+Q13;
	I2=I1-Qy/0.9936+0.015*193;
	if (6769<=I2<=7838)
		T2=700+100*(I2-6769)/(7838-6769);
	else if (7838<=I2<=8925)
		T2=800+100*(I2-7838)/(8925-7838);
	Tpj=(T1+T2)/2;
	wy=Bj*6.303*(Tpj+273)/(19.2*273);
	a0c=77;                   /*臨時計算*/
	a0s=65;                   /*臨時計算*/
    adc=a0c*0.91*0.94*0.96;
	ads=a0s*0.96*0.95*1;
	ad=(adc*74+ads*208.6)/282.6;
	tpj=(t1+t2)/2;
	v1=0.069277+(0.069277-0.054342)*0.1;
	v2=0.070668+(0.070668-0.055514)*0.1;
	v3=0.072042+(0.072042-0.056667)*0.1;
	if (370<=tpj<=380)
		v=v1+(tpj-370)*0.1*(v2-v1);
	else if (380<tpj<=390)
		v=v2+(tpj-380)*0.1*(v3-v2);
	wq=D*v/0.119;
	a0q=1250;                   /*臨時計算*/
	a2=a0q*0.96;
    kr=(1.0-0.37*(Tpj+273)*0.001)*4/(1.0-0.37*1.1477);
    ku=2.75*(pow(1147.7,0.6667))/(pow((Tpj+273),0.6667));
	kps=(kr+ku)*0.0377;
	a=1.0-(exp(-kps));
	th=tpj+Qq*1000*Bj*(0.0048+1/1228)/282.6;
	af=204*a;
	aff=af*(1+0.4*pow((1004+273)*0.001,0.25)*pow((0.665/1.11),0.07));
	a1=ad+aff;
    K=0.65*a1*a2/(a1+a2);
	td=T1-t2;
	tx=T2-t1;
	Xtnl=(td-tx)/log(td/tx);
	Xt=Xtnl*0.985;
	Q=K*Xt*282.6*0.001/Bj;
	Xt1=Tpj-255;
    Q1=K*Xt1*30.27*0.001/Bj;
	Xt3=Xt1;
	Q3=K*Xt3*7.22*0.001/Bj;
	Qe=Q+Q1+Q3;
	Xe=100*(Qy-Qe)/Qy;
	Xez=100*(Qy-Q1-Q3-Q)/(Qy-Q1-Q3);
	id=2799+Q3*Bj/(36.1-2.4);           /*需重新計算*/
	printf("i1=%f,Qq=%f,Q13=%f,I2=%f,T2=%f\n",i1,Qq,Q13,I2,T2);
	printf("Tpj=%f,wy=%f,a0=%f,adc=%f,ads=%f\n",Tpj,wy,a0,adc,ads);
	printf("ad=%f,tpj=%f,v=%f,wq=%f,a0q=%f\n",ad,tpj,v,wq,a0q);
	printf("a2=%f,kr=%f,ku=%f,kps=%f,a=%f\n",a2,kr,ku,kps,a);
	printf("th=%f,af=%f,aff=%f,Xtnl=%f,Xt1=%f",th,af,aff,Xtnl,Xt1);
	printf("a1=%f,K=%f,Xt=%f,Q=%3.1f\n",a1,K,Xt,Q);
	printf("Q1=%f,Q3=%f,Qe=%f,Qy=%f\n",Q1,Q3,Qe,Qy);
	printf("Xe=%1.2f,Xez=%4.4f,id=%f\n",Xe,Xez,id);
}