B.2.A.sagews

from sympy.physics.units import *
from sympy import *

G = 6.674 * 10**(-11) *m**3/kg/s**2
g = 9.80665 *m/s**2

r = 6371 *km
m1 = 5.97 * 10**(24) *kg

print("\n--- a: -------------------------------")
m2 = 0.1 *kg
F2 = G*m1*m2/r**2
pprint("\n")
pprint("F2:")
pprint(F2)
pprint(N(F2,2))

print("\n--- b: -------------------------------")
pprint("\n")
pprint("F2:")
m2 = 1 *kg
F2 = G*m1*m2/r**2
pprint(N(F2,2))

print("\n--- c: -------------------------------")
m2=var("m2")
a2= G*m1/r**2

print("\n--- d: -------------------------------")
pprint("e:")
e = abs(g-a2)
pprint(N(e,6))

--- a: ------------------------------- F2: 0.981624350902734*kg*m ---------------------- 2 s 0.98*kg*m --------- 2 s --- b: ------------------------------- F2: 9.8*kg*m -------- 2 s --- c: ------------------------------- --- d: ------------------------------- e: 0.00959351*m ------------ 2 s