function [x,y] = Punto_Medio(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

        xm = x(k) + h/2;
        ym = y(k,:) + f(x(k), y(k,:))*h/2;

        x(k+1) = x(k) + h;
        y(k+1,:) = y(k,:) + f(xm, ym)*h;        
    end;
end

