#include <stdio.h>
#include <math.h>
#define V1 582.66
#define B 6.205
#define Bj 6.1115
#define Qnet 18359.0
#define trk 335.0
#define irk 2194.6
#define tlk 30.0
#define ilk 193.0
#define q3 0.5
#define q4 1.5
main()
{
	float Qk,Qa,Ta,T2,I2,Vc,kr,ku,k,kps,T2s,I2s,Q1,qf;
	float ah,a1,x,qv;
	float x1,x2;
	Qk=1.06*irk+0.14*ilk;
	Qa=Qk+Qnet*(100.0-q3-q4)/(100.0-q4);
	Ta=1900.0+100.0*(Qa-20337.0)/(21518.0-20337.0);
	T2=1090.0;               /*假定*/
    I2=9926.0+0.01*(T2-1000)*(11028-9926);
	Vc=(Qa-I2)/(Ta-T2);
	kr=(1.0-0.37*(T2+273)*0.001)*0.95/(1.0-0.37*1.323);
	ku=2.52*(pow(1323,0.6667))/(pow((T2+273),0.6667));
	k=kr+ku+0.5;
	kps=k*0.45;
	ah=1.0-(exp(-kps));
    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;
	T2s=-273+(Ta+273)/(1.0+0.4575*(pow(x,0.6)));
	I2s=9926.0+0.01*(T2s-1000.0)*(11028-9926);
    Q1=0.9936*(Qa-I2s);
	qf=Bj*Q1/455;
    printf("Qk=%4.1f,Qa=%6.1f,Ta=%4.1f,T2=%4.1f,T2s=%4.1f\n",Qk,Qa,Ta,T2,T2s);
	printf("x=%6.6f,I2s=%6.1f,I2=%6.1f,Vc=%2.1f,kr=%1.2f,ku=%1.2f\n",x,I2s,I2,Vc,kr,ku);
	printf("k=%1.2f,kps=%1.4f,ah=%1.4f,a1=%1.4f,Q1=%4.1f,qf=%3.1f\n",k,kps,ah,a1,Q1,qf);
}