import random
import math

sigma = 0.00000005670373

m = round(random.uniform(0.8,2.0), 2)
print "Working mass:\t%s" % m

if m < 0.43:
    l = 0.23*pow(m,2.3)
elif 0.43 < m < 2:
    l = pow(m,4)
elif 2.0 < m < 20:
    l = 1.5*pow(m,3.5)
else:
    l = 3200*m
print "Luminosity:\t%s" % l

r = pow(m,0.85)
print "Radius:\t\t%s" % r

L = l
b = 4 * math.pi * sigma * r**2

T = (l/b)**0.25
print T