time_elapsed_s_list = # list with all time stamps
volume_water_l_list = # list with corresponding volume
volume_water_m3_list = # convert into cubic meters
m_rocket_list = # list with total mass of rocket (empty rocket + remaining water)
...
# lists for rocket's acceleration, velocity & altitude. Set initial values
a_rocket_list = [0] # acceleration of rocket in m/s^2
v_rocket_list = [0] # velocity of rocket in m/s
h_rocket_list = [0] # altitude of rocket in m
for i in range(len(time_elapsed_s_list) - 1):
dt = time_elapsed_s_list[i+1] - time_elapsed_s_list[i] # delta time [s]
# volume and mass of water expelled
# velocity and momentum of water expelled
# forces acting on rocket
# acceleration
# velocity, height (integrate step wise)
sca