r1=10;
r2=144;
a_1=0.0045;
T1=20;
T2=(r2/r1-1)/a_1+T1;
disp("the temperature (in  °C) of the hot filament is"); disp(T2);