#include <stdio.h>
#include <math.h>
#define Qnet 18359.0
#define tlk 30
#define ilk 193
#define q3 0.5
#define q4 1.5
#define q5 0.66
#define igr 450
#define igr 3332.4
#define ibh 1109.5
#define tgs 115
#define igs 485.93
#define D 36.1

main()
{
	float Tpy,Ipy,Bj,B,q2,Qky,qky,y;
	float ggl,Qgl,trk,irk;

	float Qlk,Qa,Ta,Tl2,Il2,Vc,krl,kul,kl,kpsl,T2ls,I2ls,Ql,qf;

    float Tn1,In1,Tn2,In2,Qny,Qnc,Xne;

    float Tgg1,Igg1,tgg1,Tgg2,Igg2,Xgge,Xggez,igg1,Qgr2;

	float Tdg1,Idg1,Tdg2,Idg2,idg1,tdg1,idg2,tdg2,Qgr1,XD,Xdge;

	float Qsmz;

    float tss1,iss1,tss2,iss2,Tss1,Iss1,Tss2,Iss2,Qss,Xsse;

	float tsk1,isk1,tsk2,isk2,Tsk1,Isk1,Tsk2,Isk2,Qsk,Xske;

	float txs1,ixs1,txs2,ixs2,Txs1,Ixs1,Txs2,Ixs2,Qxs,Xxse;

	float txk1,ixk1,txk2,ixk2,Txk1,Ixk1,Txk2,Ixk2,Qxk,Xxke;

	Tn2=1008;                                     /*nz假定*/
	tgg1=320;                                     /*gg假定*/
    XD=2.4;                                       /*假定*/
	tss1=186;                                     /*假定*/
	iss1=789.77;                                  /*臨時計算*/
	tsk1=166;                                     /*假定*/
    isk1=643+0.01*(tsk1-100)*(1294-643);          /*臨時計算*/
    txs2=186;                                       /*假定*/
	ixs2=794.22;                                    /*臨時計算*/
    
	/*************熱平衡及燃料燃燒計算****************/
	{
		Tpy=155;                         /*暫時假定*/
		Ipy=967.0+(Tpy-100)*9.91;
		trk=340;                         /*暫時假定*/
		irk=1958.0+(2634.0-1958.0)*0.01*(trk-300.0);
        q2=(Ipy-1.33*ilk)*(100.0-q4)/Qnet;
	    ggl=100.0-(q4+q3+q2+q5);
	    Qgl=36.1*(igr-igs)+0.02*36.1*(ibh-igs);
	    B=100.0*Qgl/(ggl*Qnet);
	    Bj=B*(1.0-q4*0.01);
	    Qky=1.06*(irk-ilk);
	    qky=Qky*100.0/Qnet;
	    y=1.0/(1.0+q5/(ggl+qky));
		printf("/*************熱平衡及燃料燃燒計算****************/\n");
	    printf("Tpy=%4.1f,Ipy=%4.1f,q2=%2.2f\n",Tpy,Ipy,q2);
	    printf("ggl=%2.2f,Qgl=%6.1f,B=%1.3f\n",ggl,Qgl,B);
	    printf("Bj=%1.4f,Qky=%4.1f,qky=%2.2f,y=%1.4f\n\n",Bj,Qky,qky,y);

	}
	/*************爐膛傳熱計算***********/
	{   
	float ah,a1,x,qv;
	float x1,x2;
	Qlk=1.06*irk+0.14*ilk;
	Qa=Qlk+Qnet*(100.0-q3-q4)/(100.0-q4);
	Ta=1900.0+100.0*(Qa-20337.0)/(21518.0-20337.0);
	Tl2=1090.0;               /*假定*/
    Il2=9926.0+0.01*(Tl2-1000)*(11028-9926);
	Vc=(Qa-Il2)/(Ta-Tl2);
	krl=(1.0-0.37*(Tl2+273)*0.001)*0.95/(1.0-0.37*1.323);
	kul=2.52*(pow(1323,0.6667))/(pow((Tl2+273),0.6667));
	kl=krl+kul+0.5;
	kpsl=kl*0.45;
	ah=1.0-(exp(-kpsl));
    a1=ah/(ah+(1-ah)*0.4397);
	x1=5.67*0.000001*0.4397*465.75*a1*(pow(Ta+273,3));
	x2=0.9936*100000*Bj*11.6;
    x=x1/x2;
	T2ls=-273+(Ta+273)/(1.0+0.4575*(pow(x,0.6)));
	I2ls=9926.0+0.01*(T2ls-1000.0)*(11028-9926);
    Ql=0.9936*(Qa-I2ls);
	qf=Bj*Ql/455;
	printf("/*************爐膛傳熱計算***********/\n");
    printf("Qlk=%4.1f,Qa=%6.1f,Ta=%4.1f,Tl2=%4.1f,T2ls=%4.1f\n",Qlk,Qa,Ta,Tl2,T2ls);
	printf("x=%6.6f,I2ls=%6.1f,I2l=%6.1f,Vc=%2.1f,krl=%1.2f,kul=%1.2f\n",x,I2ls,Il2,Vc,krl,kul);
	printf("kl=%1.2f,kpsl=%1.4f,ah=%1.4f,a1=%1.4f,Ql=%4.1f,qf=%3.1f\n\n",kl,kpsl,ah,a1,Ql,qf);
}
	/*********************凝渣管傳熱計算***************/
	{
    float Tpj,wy,ad,kr,ku,kps,a,t,th,af,a1,k,xt;
    Tn1=T2ls;
	In1=I2ls;
	t=255;
	th=t+80;
	In2=9926+0.01*(Tn2-1000)*(11028-9926);
	Qny=0.9936*(In1-In2);
	Tpj=(Tn1+Tn2)/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;
	Qnc=k*xt*89.6/(1000*Bj);
	Xne=100*(Qny-Qnc)/Qny;
	printf("/*********************凝渣管傳熱計算***************/\n");
	printf("Tn2=%4.1f,In2=%6.1f,Qny=%3.1f,Tpj=%4.1f,wy=%1.2f,ad=%f\n",Tn2,In2,Qny,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,Qnc=%3.1f,Xne=%1.1f\n\n",af,a1,k,xt,Qnc,Xne);
	}
	/**************************二級過熱器傳熱計算****************/
	{
    float h1,h2,Qq,Q1s,Q3s,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,id,v1,v2,v3;
	float Qf,tgg2,igg2,a0,a0c,a0s,adc,ads,a0q,ad;
	Qf=qf*0.68*11.75/Bj;
    Tgg1=Tn2;
	Igg1=In2;
	tgg2=450;
	igg2=3332.4;
	h1=3014.3+0.1*(tgg1-320)*(-3014.3+3040.5);
	h2=2984.0+0.1*(tgg1-320)*(-2984.0+3012.6);
	igg1=h1+0.1*(h2-h1);
	Qq=D*(igg2-igg1)/Bj;
	Q1s=55;                                /*假定*/
	Q3s=230;                               /*假定*/
    Q13=Q1s+Q3s;
	Qgr2=Qq-Qf+Q13;
	Igg2=Igg1-Qgr2/0.9936+0.015*193;
	if (6769<=Igg2<=7838)
		Tgg2=700+100*(Igg2-6769)/(7838-6769);
	else if (7838<=Igg2<=8925)
		Tgg2=800+100*(Igg2-7838)/(8925-7838);
	Tpj=(Tgg1+Tgg2)/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=(tgg1+tgg2)/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((Tgg1+273)*0.001,0.25)*pow((0.665/1.11),0.07));
	a1=ad+aff;
    K=0.65*a1*a2/(a1+a2);
	td=Tgg1-tgg2;
	tx=Tgg2-tgg1;
	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;
	Xgge=100*(Qgr2-Qe)/Qgr2;
	Xggez=100*(Qgr2-Q1-Q3-Q)/(Qgr2-Q1-Q3);
	id=2799+Q3*Bj/(36.1-2.4);           /*需重新計算*/
	printf("/**************************二級過熱器傳熱計算****************/\n");
	printf("igg1=%f,Qq=%f,Q13=%f,Igg2=%f,Tgg2=%f\n",igg1,Qq,Q13,Igg2,Tgg2);
	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,Qgr2=%f\n",Q1,Q3,Qe,Qgr2);
	printf("Xgge=%1.2f,Xggez=%4.4f,id=%f\n\n",Xgge,Xggez,id);
	}
	/**********************一級過熱器傳熱計算************************/
	{
	float Tpj,wy,ad,a0,tpj,v,wq,a0q,a2,kr;
	float ku,kps,a,th,af,aff,a1,K,td,tx,Xt,Q;
	float h1,h2,v1,v2;
	
	tdg1=257;
	idg1=2808.2;
	Tdg1=Tgg2;
	Idg1=Igg2;
	idg2=(igg1*36.1-XD*1109.5)/(36.1-XD);       /*(108706.847-2662.08)/(36.1-2.4)=3146.7*/
	h1=3116.3-0.1*(3116.3-3093.7);
	h2=3140.8-0.1*(3140.8-3119.6);
	tdg2=10*(idg2-h1)/(h2-h1)+360;
	Qgr1=(D-XD)*(idg2-idg1)/Bj;
	Idg2=Idg1-Qgr1/0.9936+0.015*193;
    Tdg2=600+100*(Idg2-5787)/(6840-5787);
	Tpj=(Tdg1+Tdg2)/2;
	wy=Bj*6.377*(Tpj+273)/(273*15.63);
	a0=76;             /*臨時計算*/
	ad=a0*1*1*0.96;
	tpj=(tdg1+tdg2)/2;
    v1=0.058821-0.1*(0.058821-0.045301);
	v2=0.060422-0.1*(0.060422-0.046723);
	v=v1+(tpj-300)*0.1*(v2-v1);
	wq=(D-XD)*v/0.095;
	a0q=1415;           /*臨時計算*/
	a2=a0q*0.99;
    kr=(1.0-0.37*(Tpj+273)*0.001)*6.19/(1.0-0.37*0.961);
    ku=3.06*(pow(961,0.6667))/(pow((Tpj+273),0.6667));
	kps=(kr+ku)*0.0199;
	a=1.0-(exp(-kps));
	th=tpj+Qgr1*1000*Bj*(0.0048+1/a2)/470.7;
	af=118.6*a;
	aff=af*(1+0.4*pow((Tpj+273)*0.001,0.25)*pow((0.955/1.164),0.07));
    a1=ad+aff;
	K=0.65*a1*a2/(a1+a2);
	td=Tdg1-tdg2;
	tx=Tdg2-tdg1;
	Xt=(td-tx)/log(td/tx);
	Q=K*Xt*470.7/(1000*Bj);
	Xdge=100*(Qgr1-Q)/Qgr1;
	printf("/**********************一級過熱器傳熱計算************************/\n");
	printf("Qgr1=%f,idg2=%f,tdg2=%f,a0q=%f\n",Qgr1,idg2,tdg2,a0q);
	printf("Idg2=%f,Tdg2=%f,Tpj=%f,ad=%f,tpj=%f,wy=%f,wq=%f\n",Idg2,Tdg2,Tpj,ad,tpj,wy,wq);
	printf("ad=%f,v=%f,wq=%f,a2=%f,kr=%f,ku=%f,kps=%f,a=%f\n",ad,v,wq,a2,kr,ku,kps,a);
	printf("th=%f,af=%f,aff=%f,a1=%f,K=%f,Xt=%f,Q=%f,Xdge=%f\n\n",th,af,aff,a1,K,Xt,Q,Xdge);
	}
	 /***********************受熱面熱量分配*************************/
	{
    float Qlz,Qnzz,Qg2z,Qg1z,Qkyz,Qky;
	Qlz=Bj*Ql;
	Qnzz=Bj*Qny;
	Qg2z=Bj*Qgr2;                            /*未整理完*/
	Qg1z=Bj*Qgr1;
	Qsmz=Qgl-(Qlz+Qnzz+Qg2z+Qg1z);
    Qky=(1.06+0.03)*(irk-ilk);
	Qkyz=Bj*Qky;
    Ipy=Idg2-(Qsmz+Qkyz)/(Bj*0.9936)+0.04*ilk+0.06*(irk+ilk)/2;
	Tpy=100+100*(Ipy-967)/991;
	printf("/***********************受熱面熱量分配*************************/\n");
	printf("Qlz=%f,Qnzz=%f,Qg2z=%f,Qg1z=%f,Qsmz=%f,Qky=%f,Qkyz=%f,Ipy=%f,Tpy=%f\n\n",Qlz,Qnzz,Qg2z,Qg1z,Qsmz,Qky,Qkyz,Ipy,Tpy);
	}
	/******************上級省煤器傳熱熱力計算******************/
	{
    float Tpj,tpj,wy,td,tx,tb;
	float a0,ad,kr,ku,kps,a,tv,af,aff,K,Xt,Q;
	Tss1=Tdg2;
	Iss1=Idg2;
	iss2=igs+(XD*1688+Qsmz)/(1.02*D);
    tss2=234+(iss2-1008.8)/(1013.5-1008.8);      /*臨時計算*/
	Qss=1.02*D*(iss2-iss1)/Bj;
	Iss2=Iss1-Qss/0.9936+0.02*193;
	if (3811<=Iss2<=4831)
		Tss2=400+100*(Iss2-3811)/(4831-3811);
	else if (4831<Iss2<=5868)
		Tss2=500+100*(Iss2-4831)/(5868-4831);
    Tpj=(Tss1+Tss2)/2;
	tpj=(tss1+tss2)/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=Tss1-tss2;
	tx=Tss2-tss1;
	Xt=(td-tx)/(log(td/tx));
	Q=K*Xt*297.5/(1000*Bj);
	Xsse=100*(Qss-Q)/Qss;
	printf("/******************上級省煤器傳熱熱力計算******************/\n\n");
	printf("iss2=%f,tss2=%f,Qss=%f,Iss2=%f,Tss2=%f,Tpj=%f\n",iss2,tss2,Qss,Iss2,Tss2,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,Xsse=%f,tss1=%f,iss1=%f,\n",Xt,Q,Xsse,tss1,iss1);
	}
	/***********************上級空氣預熱器傳熱計算***************/
	{
	float tpj,ipj,Tpj,wy,a0,l;
	float ay,wk,a0q,ak,K,xtx,xtd,td,tx,Xtnl,P,R,Xt,Q;
    Tsk1=Tss2;
	Isk1=Iss2;
	tsk2=335;
	isk2=1958+0.01*(tsk2-300)*(2634-1958);
	tpj=(tsk1+tsk2)/2;
	ipj=(isk1+isk2)/2;
	Qsk=(1.06+0.015)*(isk2-isk1);
	Isk2=Isk1-Qsk/0.9936+0.03*ipj;
	Tsk2=300+100*(Isk2-2876)/(3890-2876);
	Tpj=(Tsk1+Tsk2)/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=Tsk1-tsk2;
	xtd=Tsk2-tsk1;
	Xtnl=(xtd-xtx)/(log(xtd/xtx));
	td=tsk2-tsk1;
	tx=Tsk1-Tsk2;
	P=tx/(Tsk1-tsk1);
	R=td/tx;
	l=0.925;                /*臨時計算*/ 
	Xt=l*Xtnl;
	Q=K*Xt*1767.8/(1000*Bj);
	Xske=100*(Qsk-Q)/Qsk;
	printf("/***********************上級空氣預熱器傳熱計算***************/\n\n");
	printf("tsk2=%f,isk2=%f\n",tsk2,isk2);
	printf("tsk1=%f,isk1=%f,tpj=%f,ipj=%f,Qk=%f,Isk2=%f,Tsk2=%f,Tpj=%f\n",tsk1,isk1,tpj,ipj,Qsk,Isk2,Tsk2,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,Xske=%f\n\n",P,R,l,Xt,Q,Xske);
	}
	/************************下級省煤器傳熱計算***********************/
	{
	float Tpjj,Wy,tpj,tb,kr,Qw,ku,kps,e,ad,a1,ea1;
    float af,aff,K,a0,a,td;
	float h1,h2,Tpj,tx,tlog,Xt,Q;
	Txs1=Tsk2;
	Ixs1=Isk2;
	Qw=(XD*1688)/(1.02*D);
	ixs1=igs+Qw;
    txs1=142+(ixs1-597.8)/(602.1-597.8);
    Qxs=1.02*D*(ixs2-ixs1)/Bj;
    Ixs2=Ixs1-Qxs/0.9936+0.02*ilk;
	h1=1919.0;
	h2=2915.0;
    Txs2=100.0*(Ixs2-h1)/(h2-h1)+200.0;
	Tpj=(Txs1+Txs2)/2;
	Tpjj=Tpj+273.0;
    Wy=Bj*6.709*Tpjj/(11*273.0);
    a0=81;                        /*臨時計算*/
	ad=a0*1*1*1;
	tpj=(txs1+txs2)/2;
    tb=tpj+25.0;
    kr=8.96*(1-370.0*Tpjj)/(1.0-370.0*590.2);
    ku=4.0*(pow(Tpjj,-0.6667))/(pow(590.2,-0.6667));
    kps=(kr+ku)*0.0133;
    a=1.0-(exp(-kps));
	af=33.7*(1.0-exp(-kps));
	aff=af*(1+0.4*pow(Tpjj/1000,0.25)*pow(1/1.8,0.07));
    e=0.0027;
    a1=ad+aff;
	ea1=1/(1+e*a1);
	K=a1*ea1;
	td=Txs1-txs2;
	tx=Txs2-txs1;
    tlog=log(td/tx);
	Xt=(td-tx)/tlog;
	Q=K*697.7*Xt/(1000.0*Bj);
	Xxse=100.0*(Qxs-Q)/Qxs;
	printf("/************************下級省煤器傳熱計算***********************/\n");
	printf("a=%f,af=%f,aff=%f,a1=%f,txs1=%f,txs2=%f\n",a,af,aff,a1,txs1,txs2);
	printf("Qw=%f,ixs1=%f,ixs2=%6.1f,Ixs2=%f,Txs2=%f,K=%4.1f\n",Qw,ixs1,ixs2,Ixs2,Txs2,K);
	printf("Wy=%f,Qxs=%f,Tpj=%f,tpj=%f,tb=%f,kr=%f,ku=%f,kps=%f\n",Wy,Qxs,Tpj,tpj,tb,kr,ku,kps);
	printf("Txs1=%f,Ixs1=%f,ad=%f,td=%3.1f,tx=%3.1f,tlog=%2.6f\n",Txs1,Ixs1,ad,td,tx,tlog);
	printf("txs1=%4.1f,Qxs=%5.2f,Xt=%4.1f,Q=%5.2f,Xxse=%3.2f\n\n",txs1,Qxs,Xt,Q,Xxse);
	}
	/*******************************下級空氣預熱器傳熱計算**********************/
	{
    float Xt1,Xt2,Xt,Xtt,td,tx,P,R,Q,l;
    float tpj,ipj,Tpj,wy,a0,ay,wk,a0q,ak,k,kk;
	Txk1=Txs2;
	Ixk1=Ixs2;
	txk1=tlk;
	ixk1=ilk;
	txk2=tsk1;
	ixk2=643+0.01*(txk2-100)*(1294-643);
    tpj=(txk1+txk2)/2;
	ipj=(ixk1+ixk2)/2;
	Qxk=1.105*(ixk2-ixk1);
	Ixk2=Ixk1-Qxk/0.9936+0.03*ipj;
	Txk2=100+100*(Ixk2-967)/991;
	Tpj=(Txk1+Txk2)/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=Txk1-txk2;
	Xt2=Txk2-txk1;
	Xt=(Xt2-Xt1)/(log(Xt2/Xt1));
	td=txk2-txk1;
	tx=Txk1-Txk2;
	P=tx/(Txk1-txk1);
	R=td/tx;
	l=0.89;               /*臨時計算*/
    Xtt=l*Xt;
	Q=k*Xtt*3913/(1000*Bj);
	Xxke=100.0*(Qxk-Q)/Qxk;
	printf("/*******************************下級空氣預熱器傳熱計算**********************/\n");
	printf("Txk1=%f,Ixk1=%f\n",Txk1,Ixk1);
	printf("tpj=%f,ipj=%f,ak=%f,ay=%f,ixk2=%f\n",tpj,ipj,ak,ay,ixk2);
	printf("Ixk2=%f,Tpj=%f,Qxk=%f,Txk2=%f,wy=%f,wk=%f,k=%f\n",Ixk2,Tpj,Qxk,Txk2,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,Xxke=%f\n",Xt1,Xt2,Xt,td,tx,P,R,Xtt,Q,Xxke);
	}
	printf("/*********綜合變化***********/\n");
	printf("y=%1.4f,Bj=%f\n",y,Bj);
	printf("Tl2=%f,T2ls=%f\n",Tl2,T2ls);
	printf("Xne=%f\n",Xne);
	printf("Xgge=%f,Xggez=%f\n",Xgge,Xggez);
	printf("Xdge=%f\n",Xdge);
	printf("Xsse=%f\n",Xsse);
	printf("Xske=%f\n",Xske);
	printf("Xxse=%f\n",Xxse);
	printf("Xxke=%f\n",Xxke);
	}