Doolittle法矩阵LU分解求解方程组python编写

在网上找了很久都找不到python编写的，于是自己写了，在这里分享一下，代码调试通过，有详细注释。这里主要编写了一个自定义函数Doolittle（A,B）用于解AX=B的方程组，途中输出L、U矩阵和中间矩阵y和最终的解x。希望对大家有帮助！

import numpy as np

def Doolittle（Bb）:

    A = np.array（B）
    n = len（A）
    L = np.zeros（shape=（nn））
    U = np.zeros（shape=（nn））
    
    for k in range（n-1）:
        gauss_vector = A[:k]
        gauss_vector[k+1:] = gauss_vector[k+1:] / gauss_vector[k]
        gauss_vector[0:k+1] = np.zeros（k+1）
        L[:k] = gauss_vector
        L[k][k] = 1.0
        
        for l in range（k+1n）:
            B[l:] = B[l:] - gauss_vector[l] * B[k:]
        
        A = np.array（B）
    L[k+1][k+1] = 1.0
    U = A
    print（“U is“）
    print（U）
    print（“L is“）
    print（L）
    
    Y = np.zeros（n）
    X = np.zeros（n）
    Y[0]=b[0]
    for i in ra