Matlab实现高斯消去解方程组

function [x,res]=MyLU(A,b)
AA=A;[m,n]=size(A);
x=zeros(n,1);y=zeros(n,1);
for k=1:n-1
    A(k+1:n,k)=A(k+1:n,k)/A(k,k);
    A(k+1:n,k+1:n)=A(k+1:n,k+1:n)-A(k+1:n,k)*A(k,k+1:n);
end
l=tril(A);
for i=1:n
    l(i,i)=1;
end
u=triu(A);
for i=1:n
    sum=0;
    for j=1:i-1
        sum=sum+l(i,j)*y(j,1);
    end
    y(i,1)=(b(i,1)-sum)/l(i,i);
end
    for i=n:-1:1
        sum=0;
        for j=i+1:n
            sum=sum+u(i,j)*x(j,1);
        end
        x(i,1)=(y(i,1)-sum)/u(i,i);
    end
    res=norm(AA*x-bb);