R= 1000;
C=60.*10.^-9
L1=0.159;
L2=0.053;
H=tf([-C*L1*L2 0 0 0],[C*L1*L2 4*R*L1*L2 L1 R]);
bode(H)
%step(H)