from scipy import constants as phys
from astropy import constants as astro
from astropy import units as u
MT = astro.M_earth
ML = MT/81
RT = astro.R_earth
RL = RT*0.2728
la presion en la luna sera
print((RT/RL)**4*ML/MT)
2.2291408838837024
tau luna sera
tau = (RT/RL)**2*0.87
print(tau)
11.690430938846417
temp superficial
tsup = (1+3/4*tau)**0.25*279
print(tsup)
493.2347322774534
Ms = astro.M_sun
Mj = astro.M_jup
Rj = astro.R_jup
Sig = astro.sigma_sb
Pi = 3.1415926
G = astro.G
variacion del radio en metros por segundo
dRdt = 20*Pi*Rj**4*Sig.decompose()/Mj**2/G/3*(124**4 - 110**4)*(u.K)**4
en metros por segundo
print(dRdt)
1.1612945088959213e-11 m / s
en cms por año
print(dRdt.value*100*3600*24*365.25)
0.03664766759393393