LU分解法，列主元三角分解法MATLAB代码按层次

clear;augm=input（‘Please input the augumental matric:‘）;
[nt]=size（augm）;
for i=1:n             %d控制列三角主元的层数
    b=augm（i:）;p=i;  %4到10行是选列主元并交换
    for w=i+1:n
        if b（11）<=augm（wi）
            b=augm（w:）;p=w;
        end
    end
    c=augm（i:）; augm（i:）=b;augm（p:）=c;
    if p~=i               %只有发生了换行才将这种效果输出
        augm
    end

    for j=i:t         %首先变换与augm（ii）同一行的元素，其列指标从i到t
        s=0;
        for k=1:i-1
            s=s+