你好,下面是我需要解决P的函数,但是当我使用下面的代码时,它给了我错误的P,我有什么办法可以修改它吗?
函数[L, U, P] = lu_pp(A)
n =长度(一个);
U = 0 (n, n);
L = 0 (n, n);
P = 0 (n, n);
对于k = 1:n
U (k, k) = (k, k);
对于I = k+1:n
L (i (k) = (i (k)。/ U (k, k);
U (k, i) = (k, i);
结束
对于I = k+1:n
对于j = k+1:n
一个(i, j) = (i, j) - l (i (k) * U (k, j);
结束
结束
结束
P = L *你*发票(一个);
结束