Matlab实现列选主元高斯消去解方程组

function[x,res]=myPLU(A,b)
AA=A;bb=b;[m,n]=size(A);
A=[A,b];y=zeros(n,1);x=zeros(n,1);
for k=1:n-1
    [m,p]=max(abs(A(k:n,k)));
    p=p+k-1;
    c=A(p,:);
    A(p,:)=A(k,:);
    A(k,:)=c;
    if A(k,k)~=0
        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);
    else
        stop
    end
end

b=A(:,n+1);
l=tril(A(:,1:n));
u=triu(A(:,1:n));
for i=1:n
       l(i,i)=1;
end

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);