MATLAB 计算 画图

炮弹飞行过程中，存在空气阻力，水平方向空气阻力与速度之间的关系是F=kv^2，试求炮弹的轨迹方程，并利用Matlab进行数字模拟，画出炮弹轨迹图 k为常数自行设定，查阅资料获得炮弹的初速度v0，炮弹初速度与水平方向的夹角自行设定

m=100;  %设置质量单位/kg
v0=600; %设置初速度单位m/s
k=0.01;  %设置空气阻力系数
g=9.8; %重力加速度
pi=3.1415926;
seta=pi/3;%设置初始发射角度
h=1;  %设置初始炮台高度
y=h;
detat=0.0001; %设置时间步长
vx=v0*cos(seta);
disp(vx)
vy=v0*sin(seta);
x=0;
plotjz=[];
i=1;
while y>0.001
   ax=k*vx^2/m;
   ay=g;
   x=x+vx*detat;
   y=y+vy*detat;
   vx=vx-ax*detat;
   vy=vy-ay*detat;
   plotjz(1,i)=x;
   plotjz(2,i)=y;
   i=i+1;
end

plot(plotjz(1,:),plotjz(2,:));

自己改一下初始参数吧

里面的常数都是自行设定吗，没有给出设定值吗