w=[48.7 100.9  231.2 489 974 1460 2128  4672 6959 11020];
at=3.3
 for i=1:10
ae(i)=20*log10(abs((600+w(i))/(600-w(i))));
end
ad=ae+2*at;
 disp(ad)
