function f=eqsys(t,y)
%
% Equazioni per il sistema differenziale
%
% x' = y
% y' = -x
%
f(1) = y(2);
f(2) = -y(1);
end