299
2014-02-06 15:12:29
0
% Gravity calculation for MATLAB by ppaktion
km = 1000; % unit conversion
GM = 3.986e14;
R = 6378*km;
r_i = 10000*km; % initial position
r_f = R;
v_i = 0;
dt = 10;
r = r_i;
v = v_i;
R_t = [];
V_t = [];
t = 0;
T = [];
while (r > r_f)
t = t + dt;
dv = (-1*GM)/(r^2) * dt;
dr = dv * dt;
v = v + dv;
r = r + dr;
V_t = [V_t v];
R_t = [R_t r];
T = [T t];
end
plot(T, R_t/km)