#include <stdio.h>
#include <math.h>
#define I1 4807.9
#define T1 497.7
#define Bj 6.1115
main()
{
	float t1,i1,t2,i2,tpj,ipj,Qk,I2,T2,Tpj,wy,a0,l;
	float ay,wk,a0q,ak,K,xtx,xtd,td,tx,Xtnl,P,R,Xt,Q,Xe;
	t1=166;                 /*假定*/
    i1=643+0.01*(t1-100)*(1294-643);
	t2=335;
	i2=1958+0.01*(t2-300)*(2634-1958);
	tpj=(t1+t2)/2;
	ipj=(i1+i2)/2;
	Qk=(1.06+0.015)*(i2-i1);
	I2=I1-Qk/0.9936+0.03*ipj;
	T2=300+100*(I2-2876)/(3890-2876);
	Tpj=(T1+T2)/2;
	wy=Bj*6.585*(Tpj+273)/(6.54*273);
    a0=46;                /*臨時計算*/
	ay=a0*1*0.95;
	wk=8.7*Bj*(tpj+273)/(5.576*538.2);
	a0q=82;              /*臨時計算*/
    ak=a0q*1.11*0.91*1;
	K=0.85*ay*ak/(ay+ak);
	xtx=T1-t2;
	xtd=T2-t1;
	Xtnl=(xtd-xtx)/(log(xtd/xtx));
	td=t2-t1;
	tx=T1-T2;
	P=tx/(T1-t1);
	R=td/tx;
	l=0.925;                /*臨時計算*/ 
	Xt=l*Xtnl;
	Q=K*Xt*1767.8/(1000*Bj);
	Xe=100*(Qk-Q)/Qk;
	printf("t2=%f,i2=%f\n",t2,i2);
	printf("t1=%f,i1=%f,tpj=%f,ipj=%f,Qk=%f,I2=%f,T2=%f,Tpj=%f\n",t1,i1,tpj,ipj,Qk,I2,T2,Tpj);
    printf("wy=%f,a0=%f,ay=%f,wk=%f,a0q=%f,ak=%f\n",wy,a0,ay,wk,a0q,ak);
	printf("K=%f,xtx=%f,xtd=%f,Xtnl=%f,td=%f,tx=%f\n",K,xtx,xtd,Xtnl,td,tx);
	printf("P=%f,R=%f,l=%f,Xt=%f,Q=%f,Xe=%f\n",P,R,l,Xt,Q,Xe);
}