function res=sizefun(N)
lambda=[380 500 670 870 1020 1640]/1000;
AOD=[0.4 0.3 0.22 0.15 0.06 0.02];
tau=mie2lognormal(N,lambda);
res=sum((tau-AOD).^2);