function [x,y] = Euler(x0, y0, h, N, f)
    D = length(y0);
    
    x = zeros(N+1, 1);
    y = zeros(N+1, D);
    
    x(1) = x0;
    y(1,:) = y0;
    
    for k = 1:N
        y(k+1,:) = y(k,:) + f(x(k), y(k,:))*h;
        x(k+1) = x(k) + h;
    end;
end