已知一个方阵A，我想要U*A=L，其中U是一个上三角矩阵，L是一个下三角矩阵，并且都要求出来。用python实现

可以使用scipy库中的lu函数来实现LU分解，代码如下：

import numpy as np
from scipy.linalg import lu

# 定义一个方阵A
A = np.array([[1, 2, 3], [4, 5, 6], [7, 8, 9]])

# 对A进行LU分解
P, L, U = lu(A)

# 输出U和L
print("U:")
print(U)
print("L:")
print(L)

输出结果为：

U:
[[7.         8.         9.        ]
 [0.         0.85714286 1.71428571]
 [0.         0.         0.        ]]
L:
[[1.         0.         0.        ]
 [0.14285714 1.         0.        ]
 [0.57142857 0.5        1.        ]]

其中，P是置换矩阵，由于原矩阵A没有经过行交换，所以P为单位矩阵，可以忽略不计。