function [ yp] = shoot(t,y)
% Diffyhtaloryhma, kuten odexx sen haluaa.
%   
  %y1'=y2
  %y2'=y1^2 -1
yp=[y(2);y(1).^2-1];  
end

