El objetivo de este método es resolver un sistema de ecuaciones lineales con la forma Ax=b creando una matriz del tipo A=LU, donde A es una matriz de coeficientes, siendo x el vector desconocido y b los términos independientes. U es una matriz triangular superior y L una matriz triangular inferior y se crean aplicando operaciones elementales de fila a A.
Con L se calcula una nueva matriz con la expresión Lz=b aplicando la sustitución progresiva, y finalmente se calcula el vector solución x con la fórmula Ux=z aplicando la sustitución regresiva.
Este método también puede incluir un pivote, en este caso un pivote parcial o permutaciones de filas con la forma PA=LU, donde P es una matriz de permutación que al multiplicar a A reordena las filas del coeficiente.
Deberá proporcionar al código una matriz de coeficientes A (una matriz al cuadrado), los términos independientes b como un vector y una selección de pivote o no.
https://drive.google.com/file/d/17t-BONtgXT3rHOcCdeDeINuHI9m4v1Mw/view?usp=sharing
https://docs.google.com/document/d/1pDJBuscQWi4S6nAa1KxDEWEJf-CF45wCTK26-z7BIKM/edit?usp=sharing