Introduction - If you have any usage issues, please Google them yourself
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