clear
%lambda=[380 500 670 870 1020 1640]/1000;
%AOD=[0.4 0.3 0.22 0.15 0.06 0.02];
options = optimset('Display','iter','MaxIter',20);
N0=fminsearch('sizefun',[1000,1],options)
