Programar el m´etodo de Milne vectorial. M´as espec´ıficamente, crear en Python una funci´on Milne(F, a, b, Ua, n) cuyas variables de retorno sean T y U; T vector que contiene la discretizaci´on del intervalo [a, b], U es la matriz que contiene las aproximaciones del vector soluci´on. Las variables de entrada son F el campo vectorial asociado a la ecuaci´on diferencial vectorial, a y b extremos el intervalo, Ua es la condici´on inicial U(a) y n es el n´umero de pasos (n´umero de subintervalos en que se divide el intervalo [a, b]).

import numpy as np

def Milne(F, a, b, Ua, n):
  
  T = np.linspace(a, b, n)
  h = T[1] - T[0]
  
  U = np.zeros((T.size, Ua.size))
  U[0] = Ua
  
  for j in range(T.size - 1):
    u = U[j]
    U[j + 1] = u + h * F(u)
    
  return T, U