function Circuito_RL()
    
    % Numero de puntos
    N = 100; 
    
    % Valores Iniciales
    tini = 0;
    tfin = 1;
    i_ini = 0;
    h = (tfin-tini)/N;
    

    [t0, i0] = Solucion_Real(tini, i_ini, h, N);
    [t1, i1] = Euler(tini, i_ini, h, N, @Circuito_RL);
    
    [t2, i2] = RK(tini, i_ini, h, N, @Circuito_RL);
    [t3, i3] = RK3(tini, i_ini, h, N, @Circuito_RL);
    [t4, i4] = RK4(tini, i_ini, h, N, @Circuito_RL);
   
    [t5, i5] = Heun(tini, i_ini, h, N, @Circuito_RL);
 

    
    plot(t0, i0, 'ko-', t1, i1, 'b .-', ...
         t2, i2, 'r*-', t3, i3, 'g--',  ...
         t4, i4, 'y+-', t5, i5, 'c');

     xlabel('tiempo (seg.)');
     ylabel('Corriente (Amp.)');
     title('Solucion Circuito RL');

    % Solucion exacta de la ecuacion Diferencial
    
    function [t, ic] =Solucion_Real(t0, ir, h, N) 
       %Parametros

        R = 2.0;  % Ohms
        L = 0.3;  % Henries
        V = 10.0; % Volts
        
        t = [t0:h:N*h]';
        Tau = L/R;
        ic = 5 *(1 - exp(-t/Tau));   
    end

    % Ecuación diferencial del circuito RL

    function di = Circuito_RL(t, i)
        %Parametros

        R = 2.0;  % Ohms
        L = 0.3;  % Henries
        V = 10.0; % Volts
        di = (V - R *i)/L;
    end
end
