import math

GRAVITATIONAL_CONSTANT = 6.6743e-11
MARS_RADIUS = 3389.5e3 #m
MARS_MASS = 6.39e23    #kg
MARS_PRESSURE = 1e3    #0,7-1,155 Pa

EARTH_PRESSURE = 101325 #Pa 

MOON_RADIUS = 1737.1e3 #m
MOON_MASS = 7.3477e22  #kg
MOON_PRESSURE = 10e-9  #Pa

TIMESPAN = 100*365

def calc_atmospere_mass(radius, mass, pressure):
    area = 4*math.pi*radius**2
    g = GRAVITATIONAL_CONSTANT*mass/radius**2
    atmospere_mass = pressure *area / g
    return atmospere_mass

if __name__ == "__main__":
    mars_atm = calc_atmospere_mass(MARS_RADIUS, MARS_MASS, MARS_PRESSURE)
    print(f"Поточна атмосфера Марса важить {mars_atm:5.2g}")
    mars_teraform_atm = calc_atmospere_mass(MARS_RADIUS, MARS_MASS, EARTH_PRESSURE)
    print(f"Атмосфера Марса після тераформації {mars_teraform_atm:5.2g}")
    print(f"Атмосфера Марса важитиме більше на {mars_teraform_atm - mars_atm:5.2g}")
    print(f"Треба додавати {(mars_teraform_atm - mars_atm)/TIMESPAN:5.2g} на день 100 років")

    moon_atm = calc_atmospere_mass(MOON_RADIUS, MOON_MASS, MOON_PRESSURE)
    print(f"Поточна атмосфера Місяця важить {moon_atm:5.2g}")
    moon_teraform_atm = calc_atmospere_mass(MOON_RADIUS, MOON_MASS, EARTH_PRESSURE)
    print(f"Атмосфера Місяця після тераформації {moon_teraform_atm:5.2g}")
    print(f"Атмосфера Місяця важитиме більше на {moon_teraform_atm - moon_atm:5.2g}")
    print(f"Треба додавати {(moon_teraform_atm - moon_atm)/TIMESPAN:5.2g} на день 100 років")