function ss=sesantong(wl1,wlN,neff1,l)
%求解兩層色散曲線的通用程序 函數(shù) 
%運(yùn)用割線法
%wl1為最小波長
%wlN為此模式的截止波長
%neff1為wl1時的有效折射率
%l為模式的階數(shù)

deltawl=1e-9;
wl=wl1:deltawl:wlN;

N=1+fix((wlN-wl1)./deltawl);
neff=zeros(1,N);
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('歸一化有效折射率')
axis([0,12,0,1]);