nclu.m

来自「数值线性代数的Matlab应用程序包 共13个程序函数」· M 代码 · 共 11 行

function [L,U]=NCLU(A)
[m,n]=size(A);
L=eye(m);
U=A;
for k=1:m-1
    for j=k+1:m
        L(j,k)=U(j,k)/U(k,k);
        U(j,k:m)=U(j,k:m)-L(j,k)*U(k,k:m);
    end
end