Earth_position(t) = Earth_position(t - dt) + (rate_of_change_of__Earth_pos) * dt INIT Earth_position = initial_pos_Earth rate_of_change_of__Earth_pos = Earth_velocity Earth_velocity(t) = Earth_velocity(t - dt) + (rate_of_change_of_Earth_vel) * dt INIT Earth_velocity = 0 rate_of_change_of_Earth_vel = Earth_accel Mass_position_above_Earth(t) = Mass_position_above_Earth(t - dt) + (rate_of_change_of_mass_pos) * dt INIT Mass_position_above_Earth = initial_pos_mass rate_of_change_of_mass_pos = Mass_velocity Mass_velocity(t) = Mass_velocity(t - dt) + (rate_of_change_of_mass_vel) * dt INIT Mass_velocity = 0 rate_of_change_of_mass_vel = Mass_accel distance_between = ABS(6.38E6 + Mass_position_above_Earth-Earth_position) Earth_accel = net_force/Earth_mass Earth_mass = 5.98E24 initial_pos_Earth = 0 initial_pos_mass = 1500 {m} mass = 1 Mass_accel = -net_force/mass net_force = 6.67E-11*(Earth_mass*mass)/(distance_between^2)