how to converge equation

1 次查看（过去 30 天）

shiv gaur 2022-2-26

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1658855-how-to-converge-equation

评论： shiv gaur 2022-2-26

xnew=1;
k0=2*pi/632;
t3=1000;
t2=2000;
n1=1.512;n2=1.521;n3=4.1-1i*0.211;
n4=1;
m=0;
m=0;
k1=k0*sqrt(n1^2-xnew^2);
k2=k0*sqrt(n2^2-xnew^2);
k3=k0*sqrt(n3^2-xnew^2);
k4=k0*sqrt(n4^2-xnew^2);
n=1;
while d>0.0001 & n< 1000
    n=n+1;
    xold=xnew
    k1=k0*sqrt(n1^2-xnew^2);
    k2=k0*sqrt(n2^2-xnew^2);
    k3=k0*sqrt(n3^2-xnew^2);
    k4=k0*sqrt(n4^2-xnew^2);
    
    y=(k2)*t2-atan(k1/1i*k2)-atan((k3/k2)*tan(atan(k4/1i*k2)-k3*t3))+m*pi;
    d=abs(xnew-xold);
end
disp([n,xnew,xold])
 pl show the value  by varying the t3  1 t0 1000
plot(t3,xnew)