function sigma=rayleigh(wavel)
nu=10000000/wavel;

corr1=3.17e-60 + (7/180)*2.62e-60;
corr2=2.66e-60 + (7/180)*1.26e-60;

sigma=
