自己编的多重插补算法，可能对大家有用

Xe = input('输入终点横坐标 X\nXe = ');
Ye = input('输入终点纵坐标 Y\nYe = ');
h = input('输入步长\nh = ');
Tstep= (abs(Xe)+abs(Ye))/h;
Step=0;
F=0;
Xi=0;
Yi=0;
x = 0;
y = 0;
if(Xe>0&Ye>0) Xi=Xi+h; x = [x Xi]; y = [y Yi];
end
if(Xe<0&Ye>0) Yi=Yi+h; x = [x Xi]; y = [y Yi];
end
if(Xe<0&Ye<0) Xi=Xi-h; x = [x Xi]; y = [y Yi];
end
if(Xe>0&Ye<0) Yi=Yi-h; x = [x Xi]; y = [y Yi];
end
Step=1;
while (Step<Tstep)
if(Xi>0&Yi>=0) Dir = 1;
end
if(Xi<=0&Yi>0) Dir = 2;
end
if(Xi<0&Yi<=0) Dir = 3;
end
if(Xi>=0&Yi<0) Dir = 4;
end
F=Xe*Yi-Xi*Ye; 
if(Dir==1)