clear;clf
h=0.001;%时间步长
k=1;
t(1)=0;x(1)=0;y(1)=0;%初始值
while x<=30
x(k+1)=x(k)+400*h/sqrt(1+((90*t(k)-y(k))/(30-x(k)))^2)+eps;
y(k+1)=y(k)+400*h/sqrt(1+((30-x(k))/(90*t(k)-y(k)))^2)+eps;
t(k+1)=h*k;
k=k+1;
endfor i=1:30000
      plot(x(1:i),y(1:i),'r-',[-1 32],[0 0],'k',[30 30],[0 y(i)],'b.-');
      text(4,7,'导弹追踪'); axis equal; hold on;
      m(i)=getframe;
      hold off;
     end
movie(m);
   导弹追踪问题   运行时
>> clear
>> missile
Warning: Divide by zero.
> In missile at 8
??? Index exceeds matrix dimensions.

Error in ==> missile at 15
      plot(x(1:i),y(1:i),'r-',[-1 32],[0 0],'k',[30 30],[0 y(i)],'b.-');
请问 这是哪里出错了  还有运行时 坐标会动

> size(x)

ans =

     1    80
i=30000>80
i超过了x的维数，取i=80就可以了

提示很清楚

谢谢  再问一下  
我方战舰发现敌舰位于正东30km，并以90km/h的速度向正北方向行驶，我方迅即向敌舰发射导弹，导弹速度是400km/h，且方向始终对着敌舰。
  这题可以得到解析解吗  
dsolve('(Dx)^2+(Dy)^2=400^2','Dy=Dx*(90*t-y)/(30-x)','x(0)=0','y(0)=0')
  这样不对   哪里出错了

怎么好像是个圆？另外应该用不到dsolve