da=[.02 .05 .1];k2=2.33;
u=10.^(-da/5);
w=ones(1,3);
v=(w-u)/k2;
 m=-5*log10(v);