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
??