neff=3.5;
k0=500;
L=1e-2;
pe=221e-9;?
wl=1546e-9;
l=L/500;?
z=zeros（1501）;
B=zeros（22500）;
C=zeros（1100）;
z=0:2e-5:1e-2;?
for?x=1:2000
????wl=wl+0.001e-9;
????cc=2*pi*neff/wl;
????dc=cc-pi/pe;
????
for?t=1:500
????k=k0*exp（-10*（z（t）/L）^2）;
????sm=sqrt（（k^2-dc^2））;
????if?t==150
????T1=[cosh（sm*l）+（i*dc/sm）*sinh（sm*l）]*exp（-i*dc*l）*exp（i*pi）;
????T2=i*（k/sm）*sinh（sm*l）*exp（-i*dc*（z（t）+z（t+1）））;
????T3=-（i*k/sm）*sinh（sm*l）*exp（i*dc*（z（t）+z（t+1）））;
????T4=[cosh（sm*l）-（i*dc/sm）*sinh（sm*l）]*exp（i*dc*l）*exp（-i*pi）;
????else
????T1=[cosh（sm*l）+（i*dc/s