%
% Risolve la ODE
%  y'   = -a*y
%  y(0) = 1
%
% con il metodo
%  y(1)   = y(0)   + h f(x(0),y(0))
%  y(k+1) = y(k-1) + 2 h f(x(k),y(k))
%
function [x,y]=instabile(fun,xf,N)
  h      = xf / N ;
  x(1)   = 0 ;
  y(1,:) = 1 ;
  %% passo Eulero Esplicito
  x(2)   = x(1) + h ;
  f1     = feval(fun,x(1),y(1,:));
  y(2,:) = y(1,:) + h * f1 ;
  %% ciclo con il metodo multistep
  for k=2:N
	fk = feval(fun,x(k),y(k,:));
    x(k+1)   = x(k)+h ;
    y(k+1,:) = y(k-1,:)+2*h*fk ;
  end
end
