%% 近月点、远月点的速度求解。
clc;clear;
% 万有引力常量G=6.672x10^-11N*m^2/kg^2
G=6.672*10^-11;
%月球质量 M=7.3477*10^22
M=7.3477*10^22;
%% 利用机械能守恒去计算
%椭圆的长半轴a、短半轴b、焦点c均通过小计算器算得到。
c=42.5*10^3;
a=1794.5*10^3;
b=1793.99665*10^3;
Va=(a+c)/b*sqrt(G*M/a) %近日点速度
Vb=(a-c)/b*sqrt(G*M/a) %远日点速度


