clc;
Cd=xlsread('trial',1,'b2:b6');
rho = xlsread('trial',1,'a1:a6');
S=3.2365;
m=1000;
gama=10;
H=input('enter altitude in meter (m)');
V=input('enter the velocity (m/s)');
dt=1;
t=0:dt:10^3;
n=(10^3/dt+1);
Ux=zeros(1,n);
Ux = V*cos(gama);
Uy = V*sin(gama);
for i=1:n
D=0.5.*rho*V^2*S.*Cd;
Ux(1,i+1)=Ux(1,i)-(dt)*(D/m*cos(gama));
Uy(1,i+1)=Uy(1,i)+(dt)*(g-(D/m*sin(gama)));
V(1,i+1)=((Ux(1,i+1))^2+(Uy(1,i+1)^2)^0.5;
i=i+1;
end
Best Answer