Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- import random
- import math
- sigma = 0.00000005670373
- SUN_TEMP_K = 5778
- SUN_RADIUS_M = 695700 * 1000
- SUN_LUM_WATTS = 1368
- m = round(random.uniform(0.1,1.8), 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
- Lw = l * SUN_LUM_WATTS
- print "Luminosity (watts):\t%s" % Lw
- R = m ** 0.85
- print "Radius: %s" % R
- R_km = R * SUN_RADIUS_M
- print "Radius (meters): %s" % R_km
- T = (Lw / (4.0 * math.pi * sigma * (R_km**2.0))) ** (1.0/4.0)
- print "Surface temp: %s" % T
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement