#include <stdio.h>
#include <math.h>

#define	dt		0.01				// ��‚Ý•�[sec]
#define pi		3.1415926			// ‰~Žü—¦
#define Loop	5000000				// ŒvŽZ‰ñ�”
#define STEP	10000

#define To		300					// �‰ŠúƒAƒ‹ƒ~‰·“x[K]
#define Ts		3					// ‰F’ˆ”wŒi‰·“x[K]

#define xx		0.25				// �¡–@x[m]
#define yy		0.25				// �¡–@y[m]
#define zz		0.005				// �¡–@z[m] Œú‚Ý

#define sig		5.67E-8				// ½Ã̧Ý�EÎÞÙÂÏÝ’è�”
#define A		(xx*yy)				// ”–Ê�Ï[m2]
#define V		(A*zz)				// ƒAƒ‹ƒ~”‚̑Ì�Ï[m3]
#define ep		0.8					// •úŽË—¦ ±ÙЕ\–Ê
#define rho		2700				// ƒAƒ‹ƒ~‚Ì–§“x[kg/m3]
#define C		881					// ƒAƒ‹ƒ~‚Ì”ä”M[J/kg�EK]
#define Cal		(C*V*rho)			// ƒAƒ‹ƒ~‚Ì”M—e—Ê[J/K]26.0111

#define Qh		17.6				// ƒq�[ƒ^�o—Í[W]

static double 	Tn;					// ƒAƒ‹ƒ~‰·“x[deg]–¢—ˆ
static double 	Tq;					// ƒAƒ‹ƒ~‰·“x[deg]Œ»�Ý
static double 	Qr;					// •úŽË

int main(void)
{
	int i;
	double Qr = 0;
	double Tq = Tn = To;
	
	for(i=0;i<Loop;i++)
	{
		// Qr‚ÌŒvŽZ
		Qr = sig*A*ep*(Tq*Tq*Tq*Tq - (double)Ts*Ts*Ts*Ts);
		
		// ”M•û’öŽ®‚̉ð–@
		Tn = Tq + (-Qr*2 + Qh)*dt/Cal;
		
		// ŒvŽZŒ‹‰Ê‚Ì•\Ž¦
		if (i % STEP == 0) printf("%f\t%f\t%f\n",i*dt,Tq,Qr);
		
		// ’l�X�V
		Tq	= Tn;
	}
	
	return 0;
}