部分主元法的LU分解

function [LUP]=partialLUDecomposition（A）  % A=[1 2 -3 4;4 8 12 -8;2 3 2 1;-3 -1 1 -4]
[mn]=size（A）;
if m~=n
    error（‘不是方阵!‘）;
end
L=zeros（m）;
P=eye（m）;
for i=1:m-1
    index=find（abs（A（i:mi））==max（abs（A（i:mi））））;
    k=index+i-1;
    if k~=i