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

