T0 = 300
sigma = 5.67E-08
def temp(T):
return (715 - 0.4*sigma*(T**4-T0**4) + 32*T0)/32
print temp(301)
322.266821424
T = 301
while T-temp(T) != 0:
T = temp(T)
print T
322.266821424 320.440012888 320.611882213 320.595837164 320.59733616 320.597196127 320.597209208 320.597207986 320.597208101 320.59720809 320.597208091 320.597208091 320.597208091 320.597208091 320.597208091
def rk_2_step(T):
return ?