def Runge_Kutta_4(f, a, b, N, y0, A, B, alpha, beta):

    h = n((b-a)/N)          # paso
    nodos = [n(a)]          # Inicializamos una lista nodos con el primer nodo
    aprox = [n(y0)]         # Inicializamos una lista aprox con la condición inicial

    for k in range(N):

        K1 = n(f(nodos[k], aprox[k]))
        K2 = n(f(nodos[k]+alpha*h,aprox[k]+beta*K1*h))
        K3 = n(f(nodos[k]+alpha*h,aprox[k]+beta*K2*h))
        K4 = n(f(nodos[k]*h,aprox[k]*K3*h))
       nodos.append(n(nodos[k]*h))
       yRK4=n(aprox[k]+h*(K1+2*K2+2*K3+K4)/6)
       aprox.append(n(yRK4))

       return matrix(zip(nodos, aprox))