l=600:1000:19600;beta=10^(-9);l0=18000;d=18000000;
a=-1.158;f=1250:2500:d+1250;
for i=1:20
b=a*l(i)/l0*10^(-1.5);
u(i)=exp(b.*f.^0.5);beta1=beta*10^(-4.5);
v(i)=beta1.*f.^1.5;
g(i)=log2(1+u./v);
%c=quad('g',0.1,1000);
z(i)=2500*trapz(g(i));
%for r=1:max(size(f))
%s=0;
%s=s+log2(1+u(r)/v(r));
%end
%zz(i)=2500*s;
end
plot(l,z)
hold on
plot(l,zz)
