function Presa_Depredador(Metodo)
    
    % Numero de puntos
    N = 1000; 
    
    % Valores Iniciales
    tini = 0;
    tfin = 20;
    
    presas_ini = 1;
    depredador_ini = 1;
    
    h = (tfin-tini)/N;
    
    [t,x] = Metodo(tini, [presas_ini, depredador_ini], h, N, @predprey);
        
    title('Solucion Ecuaciones Presa depredador');
  
    plot(t, x(:,1),  'k', t, x(:,2));
    xlabel('tiempo (seg.)');   


    % Ecuación diferencial 

    function yp = predprey(t, y) 
        yp = [
            1.2*y(1) - 0.6*y(1)*y(2), ...
           -0.8*y(2) + 0.3*y(1)*y(2)];
    end
end
