运行不出来啊 郁闷 。 这是等误差查补程序 高手过目啊

梦想成真

x1=input('输入起点横坐标 X\n x1 = ');
y1=input('输入起点纵坐标 Y\n y1 = ');
x2=input('输入终点横坐标 X\n x2 = ');
y2=input('输入起点纵坐标 Y\n y2 = ');
xa=x1;
ya=y1;
x=x1;
y=y1;
xi=x1;
yi=y1;
while(xa<=x2)
       yp =ya-0.04;
while(yp-ya<=0.04)
m=(yp- ya)*(yp- ya)
xp=xa-sqrt(0.04*0.04-m);
k=sqrt(0.04*0.04/m -1);
xt=(k-1)/2;
yt =xt*xt+xt+1;
k1=(yt-yp)/(xt-xp);
if abs(k-k1)<=0.01
n=(1-k)*(1-k);
u=4*(1+k*xa-ya);
t=sqrt(n-u);
xi=k/2-0.5+t/2;
yi=-k*xa+k*xi+ya;
    xa=xi;
ya=yi;
xi=[x1,xi];
yi=[y1,yi];
yp = yp +0.01;
else
yp = yp +0.01;
end
end
end
plot(xi,yi);
    grid;

gaoshanliu水

不懂，高手快来。。。

数模鸿雁

我表示看不懂你的程序！能不能把你的原题拿出来看下！