clear all;
 close all;
 L1=(sqrt(2)/(pi*10000))*8
 C1=((1/(2*pi*5000))^2)/(L1)
 s2=L1*C1
 (1/(2*pi*5*10^3))^2
 L2=(sqrt(2)/(pi*10000))*8
 C2=L2/(2*pi*5*1000)

 for f=1:10000
    s=(1j*2*pi*f);
    x(f)=f;
    HD(f)=1/(((1j*2*pi*f)^2*(L1*C1))+(1j*2*pi*f*(L1/8))+1);
    HW(f)=(s^2*L2*C2)/((s^2*(C2*L2))+(s*(L2/8))+1);
 end
 plot(x,abs(HD));

 plot(x,abs(HW));
 H=HD+HW;
 plot(x,abs(H));