%求解兩層、1階第2個模式的色散曲線   @@@(程序4）

%運用割線法

clear all;

l=1;
wl1=500e-9;
wlN=694.453188e-9;
deltawl=5e-9;
wl=wl1:deltawl:wlN;

N=1+fix((wlN-wl1)./deltawl);
neff=zeros(1,N);
neff1=1.45294294183074;
neff(1)=neff1;



for i=2:1:N
   u=secant(neff(i-1),neff(i-1)-0.001,wl(i),l);
    neff(i)=u;
end
f=3.*10.^8./wl;

r1=4e-6;
nc=1.458;
n=1.45;
k=2.*pi./wl;
beta=k.*neff;
b=(beta./k-n)./(nc-n);%歸一化
v=r1.*k.*sqrt(nc.^2-n.^2);
plot(v,b,'k')
grid on;

xlabel('歸一化頻率')
ylabel('歸一化有效折射率')


hold on;