>>385
0〜θ の面積分率は

-π/2 ≦ θ <π/2 のとき
　{arctan(√(1-ee)･sinθ/(e+cosθ)) - e√(1-ee)･sinθ/(1+e･cosθ)}/2π,
π/2 < θ < 3π/2 のとき
　1/2 + {arctan(√(1-ee)･sinθ/(e+cosθ)) - e√(1-ee)･sinθ/(1+e･cosθ)}/2π,
3π/2 < θ < 5π /2のとき
　1 + {arctan(√(1-ee)･sinθ/(e+cosθ)) - e√(1-ee)･sinθ/(1+e･cosθ)}/2π,

・地球と月の距離r
　1/r = (1+e･cosθ)/L,　　　　　　L = a(1-ee).