杜利托尔分解

function [L,U]=doolittle_decomposition(A)
    %{
        杜利托尔分解
        param A：二维矩阵（方阵）
        return L：下三角矩阵
        return U：上三角矩阵
    %}
    
    rows=size(A,1);
    cols=size(A,2);
    L=eye(rows,cols);
    U=A;

    for ii=1:rows
        for jj=ii+1:rows%从当前行开始消元
            temp=-U(jj,ii)/U(ii,ii);
            U(jj,:)=U(jj,:)+temp*U(ii,:);
            
            %对单位阵做另一行变换
            L(ii,:)=L(ii,:)-temp*L(jj,:);
        end
    end

    L=L';%取转置
end