%
% Esempio di metodo basato sullo sviluppo di Taylor
%
% calcolo della soluzione approssimata
%
function [x,y]=test_1_taylor(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 basato su Taylor
    yp  = test_1_yp(x(k),y(k)) ;
    ypp = test_1_ypp(x(k),y(k)) ;
    y(k+1) = y(k) + h*yp + (h*h/2)*ypp ;
  end;
end