L=0.01;
neff=1.447;
c=1*1e-7/（2*neff）;
deltaneff=0.0004;
N=100;
M=2000;
lamda1=1548;
lamda2=1552;
lamda=linspace（lamda1*1e-9lamda2*1e-9M）;
deltalamda=（lamda2*1e-9-lamda1*1e-9）/（M-1）;
for?k=1:M
????F=[10;01];
????for?n=1:N

deltaneff_z=deltaneff;
??%g=4;
???%?deltaneff_z=deltaneff*exp（-g*g*（（2*n-1-N）/（2*N））^2）;
??%??H=0.5;
??%deltaneff_z=deltaneff*（1+H*cos（pi*（2*n-1-N）/N））/（1+H）;
????lamdaD=1550*1e-9-L*neff*c+2*neff*c*（L*n/N-0.5*L/N）;
????sigma=2*pi*neff*（1/lamda（k）-1/lamdaD）+2*pi*deltaneff_z/lamda（k）+（4*pi*neff）*c*（2*（neff+deltaneff））*（-L/2+n*L/N）/lamdaD^2;
????kac=pi*deltaneff_z/lamda（k）;
????RB=sqrt（kac^2-sigma^2）+eps;
????F=F*[cosh（RB*L/N）-（sigma/RB）*sinh（RB*