FBG反射谱，要用的人会非常需要，需要的人可以看看，很不错的-FBG id1=fopen( lambdaCMT.txt , wt )　

fid2=fopen( rhoCMT.txt , wt )　



f1=inline( -i.*delta.*y1-k.*y2 , t , lambda , wk , delta , k , y1 , y2 )　　　

f2=inline( -k.*y1+i.*delta.*y2 , t , lambda , wk , delta , k , y1 , y2 )　　　　



lambda=1547　

while lambda<1552

t0=-501426　　

h=2.0057e+003　

n=500　

y1=1　　　

y2=0