Untitled

m=0.1 # масса квадратного метра листа
C=1000 # теплоемкость дж на килограм кельвин
dt=0.0001 # шаг по времени
T1=273+90 # начальная температура
SB=5.67e-8 # постоянная стефана больцмана
A=0.8 # степень черноты

# падение температуры за одну секунуд
def cool(T):
	t = 0.0
	while t < 1.0:
		dE=A*SB*T*T*T*T*dt # закон стефана больцмана
		dT=dE/m/C
		T=T-dT
		t=t+dt
	return T

T = T1
print(f'T1={T1}')
# таблица температуры на 50 сек
for i in range(50):
	T = cool(T)
	print(f'{i+1}: {T}')

print('ok')