% Metodo de Euler para integrar EDOs numericamente
% ENTRADA:
% 1) Una funcion anonima f, de dos variables (importante)
% 2) Un vector con las coordenadas x de los puntos en que aproximar la solucion
% 3) un valor inicial y0 correspondiente a x0 (que es x(1))

% SALIDA:
% un vector de valores y(i) que aproximan la solucion en x(i)
function [y] = euler(f, x, y0)
    % se creea la lista y con la misma longitud que x
    y = zeros(size(x));
    % se almacena la condicion inicial en el primer punto
    y(1) = y0;

    % Bucle de Euler
    for s = 2:length(x)
         y(s) = y(s-1) + f(x(s-1), y(s-1)).*(x(s) - x(s-1));
    end
end