%ROCKET_CAR Summary of this function goes here
%   Detailed explanation goes here
initial_vel=0;
seconds = 20;
drag=0.001;
mass = .14;
thrust=20;
[t, vel_numerical]=ode23('g',0,seconds,initial_vel);
acc=143-0.0071*vel_numerical.^2;
subplot(2,1,1), plot(t,vel_numerical), title('velocity'), ylabel('m/s')
subplot(2,1,2),plot(t,acc), title('Acceleration'), xlabel ('Time,s'),ylabel('m/s^2')