function Sistema_Masa_Resorte(Metodo)
    
    % Numero de puntos
    N = 1000; 
    
    % Valores Iniciales
    tini = 0;
    tfin = 40;
    d_ini = 2;
    v_ini = 0;
    
    h = (tfin-tini)/N;
    
    [t,x] = Metodo(tini, [d_ini, v_ini], h, N, @Sistema_MR);
        
    title('Sistema Masa Resorte');
  
    subplot(2,2,1)
    plot(t, x(:,1));
    xlabel('tiempo (seg.)');
    ylabel('Desplazamiento (m.)');
    
    subplot(2,2,2)
    plot(t, x(:,2));
    xlabel('tiempo (seg.)');
    ylabel('Velocidad (m/s.)');

    subplot(2,2,3)
    plot(x(:,1), x(:,2));
    xlabel('Desplazamient (m.)');
    ylabel('Velocidad (m/s.)');


    % Ecuación diferencial del circuito RL

    function D = Sistema_MR(t, r)
        M = 10.0;  % Masa  en Kg 
        C = 1;   % Coeficiente de Amortiguamiento Dinnamica
        K = 5;   % Constante del Resorte
        F = 0;    % Fuerza Aplicada
        d = r(1);
        v = r(2);
        
        D(1) = v;   % derivada del desplazamiento
        D(2) = (F - C*v - K*d)/M; % derivada de la velocidad
    end
end
