clc; 
close all; 
clear all;

delta_t=0.01;
t0=0;
v0=0;
a0=0;
tfin=30;
N=3000;
delta_t=(tfin-t0)/N;

%INICIALIZACIÓN 
M=zeros(N+1,3);
M(1,1)=t0;
M(1,2)=v0;
M(1,3)=a0;

t=t0;
i=1;
tol_ac=1e-5;
while(t<tfin) 

M(i+1,1)=M(i,1)+delta_t;
M(i+1,2)=M(i,2)+M(i,3)*delta_t;
M(i+1,3)=aaereo(M(i,2),t);

t=M(i+1,1);
i=i+1;
cadenainfo=sprintf('iteraccion %u tiempo=%f velocidad=%f aceleracion=%f', i,t, M(i+1,2), M(i+1,3));
disp(cadenainfo);


end