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