%
%
%
clear;
w0=-50;
w1=50;
w=w0:0.01:w1;
t=0:0.01:6;
L=input('input the value of L:');
C=input('input the value of C:');
wc=1/sqrt(L*C);
H=sqrt(1./(((1-(w/wc).^2).^2+(w/wc).^2)));
phi=-57.3*atan((w/wc)./((1-(w/wc).^2)+eps));
h=2*wc/sqrt(3)*exp(-0.5*wc*t).*sin(0.866*wc*t);

clf;
subplot(3,1,1)
plot(w,H);grid on;title('the magnitude response');
subplot(3,1,2)
plot(w,phi);grid on;title('the phase response');
subplot(3,1,3)
plot(t,h);grid on;title('the impulse response')