%
% Metodo di Collatz o Eulero Migliorato per sistemi di equazioni differenziali
% 
% Ingresso
%  f      = funzione per il sistema y'=f(x,y)
%  x0, y0 = dato iniziale del problema y(x0) = y0
%  xf     = estremo superiore dell'intervallo di calcolo [x0,xf]
%  N      = numero di intervalli in cui suddividere  [x0,xf]
% 
% Uscita
%  x(:)  = vettore soluzione delle 'x'
%  y(:)  = vettore soluzione delle 'y'
%
function [x,y]=Eulero_esplicito(f,x0,y0,X,N)
 x(1)   = x0;
 y(1,:) = y0;
 h      = X/N ;
 for k=1:N
   % griglia spazio
   x(k+1) = x(k) + h;
   % soluzione approssimata col metodo di Eulero
   y(k+1,:) = y(k,:) + h*feval(f,x(k),y(k,:));
 end;
end