L = 0.5e9 T = 26.7468*1.e6 ; Boltzmann constant - k = 1.38e-16 erg K^-1 Ko = 1.e-6 ;max # needed to add to Lcond+Lrad+Udot to make it =0 mx = 1.1280883e+028 ;Lcond where mx is greatest lc = 9.18456e+025 nlc = lc + mx top = nlc*L bottom = 4*Ko*T^(7/2) nA = top/bottom print, na end