clc //initialisation si=5.735*10^-8//j m^-2 sec ^-1 deg^-4 t=1227+273//k r=0.003//m //CALCULATIONS e=3.14*r*r*si*t^4*60/4.2 //results printf(' \n energy= % 1f cal',e)