function ecuacion_simple() 
    % Numero de puntos
    N = 8; 
    
    % Valores Iniciales
    xini = 0;
    xfin = 4;
    y_ini = 1;
    h = (xfin-xini)/N;
    x = [xini:h:xfin]';
    

    [x1,y1] = Euler(xini, y_ini, h, N, @ecuacion);
    [x2,y2] = Heun(xini, y_ini, h, N, @ecuacion);
    [x3,y3] = Punto_Medio(xini, y_ini, h, N, @ecuacion);
    [x4,y4] = RK(xini, y_ini, h, N, @ecuacion);
    [x5,y5] = RK3(xini, y_ini, h, N, @ecuacion);
    [x6,y6] = RK4(xini, y_ini, h, N, @ecuacion);
    y2
  
    [x, solucion_real(x), y1, y2, y3, y4, y5, y6]

  
      [x, abs(solucion_real(x)-y1),...
          abs(solucion_real(x)-y2),...
          abs(solucion_real(x)-y3),...
          abs(solucion_real(x)-y4),... 
          abs(solucion_real(x)-y5),...
          abs(solucion_real(x)-y6)]
    
    function dy = ecuacion(x, y)
        dy = -2*x^3 + 12*x^2- 20*x + 8.5;
    end

    function y = solucion_real(x)
        y = -0.5*x.^4 + 4*x.^3 - 10*x.^2 + 8.5*x +1;
    end
end