function y=fsys(t,y)
L=0.1; C=1.e-03; R=10; e=5;
yy(1) = y(2);
y(2) = -(L+R*C^2)/(L*C*R)*y(2)-2/(L*C)*y(1)+e/(L*C);
y(1)=yy(1);
return