beta=10^(-9);
l=.25:.25:5.5;
for j=1:22
w=[.8 3 5 10 20 50 100 150 200 250 300 400 500 600 700 800 1000 1500 2000 2500];
w1=(w.^1.5).*beta;
alfa=[1.23 2.38 2.95 3.96 5.09 6.37 7.15 7.64 8.37 9.22 10.01 11.31 12.62 13.75 14.7 15.66  17.4 21.5 23.88 26.36];
alfa1=alfa./(-10);alfa2=alfa1.*l(j);
hc=10.^alfa2;
in=log2(1+hc./w1);
s=0;
for i=1:length(w)-1
s=s+(w(i+1)-w(i))*(in(i)+in(i+1));
end
c(j)=1000*s;
end
plot (l,c)
 %s-a omis factorul .5 la cal cularea ariilor trapezelor!
% Se va lua
%c.*.5 adica 1/2 din valoarea afisata