aminit;

%Test of Bodes relations
%This is the weighting function in Bodes relation
w=logspace(-2,2,1000);
y=log(abs((w+1)./(w-1)));
I=sum(y)*log(w(2)/w(1))/pi;
%Calculate area for checks
I*2/pi
subplot(2,1,1)
pl=semilogx(w,y,'b-');
set(gca,'Ytick',[0 2 4 6]);
set(pl,'Linewidth',1.5);
%title('Weight');
axis([0.01 100 0 6]);
xlabel('x')
ylabel('y')

amprint('bodeweight.eps');
