# READ

(*
restart :
read("D:/OneDrive/docs/Maple/observing.mpl") :
*)

# PACKAGES

with(Units) :

# ALIASES

alias(log100 = log[100], pi = Pi, S_ellipsoid = ellipsoid, u = Unit) :

# UNITS

_km := u(km) :
_lm := u(lm) :
_m := u(m) :

# FUNCTIONS

E100 := x -> 100^x :
S_ellipse := (r_e, r_p) -> pi*r_e*r_p :
S_spheroid := (r_e, r_p) -> S_ellipsoid(r_e, r_e, r_p) :

S_circle := r -> S_ellipse(r, r) :
S_sphere := r -> S_spheroid(r, r) :

# ASSIGNMENTS

E_lx := u(lx) :
L_Sol := 3.75E28*_lm :
m1 := -14.18 :
n_filter := 100000. :
theta_rev2 := u(semicircle) :

# TERRA

a_Terra := 149598023.*_km :
delta_Terra := theta_rev2 :
r_e_Terra := 6378137.0*_m :

# LUNA

a_Luna := 384399.*_km :
A_Luna := 0.136 :
r_Luna := 1737.4*_km :

Delta_Luna := a_Luna - r_e_Terra :
J_Luna := A_Luna*L_Sol/S_sphere(a_Terra) :
sigma_Luna := S_circle(r_Luna) :

Phi_Luna := J_Luna*sigma_Luna/2 :

E_Luna := Phi_Luna/S_sphere(Delta_Luna) :

m_Luna := m1 - 5*log100(E_Luna/E_lx) : # -10.88046678

# SOL

Delta_Sol := a_Terra - r_e_Terra :

E_Sol := L_Sol/S_sphere(Delta_Sol) :

m_Sol := m1 - 5*log100(E_Sol/E_lx) : # -26.99251682140882611192725773641007

m_filtered_Sol := -5*log100(E100(-m_Sol/5)/n_filter) : # -14.49251682