Question

写了一个MATLAB程序，可是除了B=0以外，其余取值都是无限死循环T T急求帮助！

function [xp,yp,xm,ym,T]=euler6(h,B)
Vp=135;Vm=450;
xp(1)=0;
yp(1)=120;
xm(1)=0;
ym(1)=0;
d(1)=120;
T(1)=0;
A(1)=pi/2;
C(1)=A(1)-B;
k=1;
while d>0.1
k=k+1;
xp(k)=xp(k-1)+Vp*h*cos(C(k-1));
yp(k)=yp(k-1)+Vp*h*sin(C(k-1));
xm(k)=xm(k-1)+Vm*h*cos(A(k-1));
ym(k)=ym(k-1)+Vm*h*sin(A(k-1));
d(k)=sqrt((xp(k)-xm(k))^2+(yp(k)-ym(k))^2);
A(k)=asin((yp(k)-ym(k))/d(k));
C(k)=A(k)-B;
T(k)=T(k-1)+h;
end;
xp,yp,xm,ym,T

这个就是那个程序，当h=0.001，B=0时能得出正确答案，可是代入B的别的取值后就是无限死循环T T不知道到底是哪里的问题T T求帮助T T

Answer 1

粗略看了下，应该是while后面条件写错了，我估计你的意思应该是写成d(k)>0.1这样应该不会死循环了

追问

不是这个的问题的说。。。。。因为。。。h=0.0001，B=0的时候算出来的答案是对的。。。顺便可不可以问下。。。怎么改才能求出B在区间【0，pi】上那个值的时候，T最大？