clear
load thrust
load velocity

mass = 55;
t(1) = 0;
a(1) = thrust(1);
v(1) = 0;
d(1) = 0;
r = 32.174*.00273;

for ii = 2:3001
    t(ii) = (ii-1)/100;
    a(ii) = interp1(velocity,thrust,v(ii-1))/mass*32.174 - (.01*v(ii-1)^2 + .02*v(ii-1));
    v(ii) = v(ii-1) + (a(ii-1)+a(ii))/2*(t(ii)-t(ii-1));
    d(ii) = d(ii-1) + (v(ii-1)+v(ii))/2*(t(ii)-t(ii-1));
    if d(ii) > 170,break,end
end

plot(t,[a',v',d'])
legend('Acceleration','Velocity','Distance')
xlabel('Time (s)')
ylabel('^f^t/_s2    |    ^f^t/_s    |    ft')
figure(2)
plot(d,v)
xlabel('Groundroll Distance (ft)')
ylabel('Velocity (^f^t/_s)')