function y=ffocault(t,y)
%FSYS Termine forzante per il problema del
%   pendolo di Focault.
l=20; k2=9.8/l; psi=pi/4; omega=7.29*1.e-05;
yy=y;
y(1)=yy(3);
y(2)=yy(4);
y(3)=2*omega*sin(psi)*yy(4)-k2*yy(1);
y(4)=-2*omega*sin(psi)*yy(3)-k2*yy(2);

return
