[求助]关于配焦椭流线法的matlab编程实现问题 [复制链接]

求前辈指点迷津，感觉椭流线表达式有问题，又不知道哪里不对 ?8aPd"x  
2mVLR;s{_  
程序如下 BRGTCR  
9S$?2z".2  
clearclc u@$pOLI
 
h=1000;N=100; qD/FxR-!  
R=100;for i=2:101 /7-qb^V  
    theata(i)=asin(sqrt((i-1)/N));end TM1J1GU  
theata(1)=0;for i=1:101 lp%i%*EQ*  
%      r(i)=R*sqrt((i-1)/N);r(i)=sqrt(R^2/i) 1PdxoRa4=  
end% r(1)=0; :z}  
a=zeros(101,1);c=zeros(101,1); XqK\'8]\Mw  
x=zeros(101,1);y=zeros(101,1); }t9A#GOz  
a(1)=550;x(1)=-50; fE&wtw{gi  
y(1)=0; J t,7S4JL  
mB]Y;R<  
for i=2:100    k=2*a(i-1); PmT,*C`/X  
    F=[h,r(i-1)/2-r(i)/2];    XY=((k^2-norm(F)^2)/(2*k-2*norm(F)*cos(pi-theata(i))))*[cos(pi-theata(i)),sin(theata(i))]; 1q(o3%   
    x(i)=XY(1);    y(i)=XY(2); C*fSPdg?  
    FX=F-XY;    a(i)=(norm(FX)+norm(XY))/2;  gC}D0l[  
    i    x(i) *) ?Fo  
    y(i)end _I}rQfPJ  
plot(x，y）

请问楼主解决了吗