nclu.m

数值线性代数的Matlab应用程序包 共13个程序函数

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