Rocket sim

t = 0:0.001:10;  %s
g = -9.81;   %m/s2
vi = 50; %m/s
height = vi .* t + 0.5*g*t.^2;
vf = vi + g .* t;

plot(t,height,t,vf)
xlabel('Time (s)')
ylabel('Height (m)')