function [f] = sho(y,t)
   f = y;
   f(1) = y(2);
   f(2) = -y(1);
end