LU紧凑格式 matlab

function[L,U,b,x]=LUb(A,b);
%L:下三角矩阵；Ux=b:消元法所得;x：函数的解
[m,n]=size(A); 
%A为方阵，所以也可以为m=size(A，2)
%先一行，再一列
for i=1:m
    for j=i:n
        if i>1
            for t=1:i-1
                A(i,j)=A(i,j)-A(i,t)*A(t,j);
            end
        end
end	
%上面为行矩阵元素求解，U相关元素
    if i>1 
        for t=1:i-1
            b(i)=b(i)-b(t)*A(i,t);
        end
end
%为b进行变换
    for z=i+1:m
        if i