Help:=proc(): print(``):end: SMdata:=proc(): dS:=149600000; diS:=1391400; deS:=1.408; #tidal force: c*M/d^3= c*rho*(r/d)^3 dM:=384399; diM:=3500; evalf(dS/dM); evalf(diS/diM); deM:=3.34; end: