【Matlab】Matlab实现方阵LU分解——Doolittle方法

1、消元过程：求解L/U三角方阵，并在原方阵A中存储

2、回代过程

A = [5 7 9 10;
     6 8 10 9;
     7 10 8 7;
     5 7 6 5];
b = [1;1; 1; 1];
xc = A\b

n = 4

%求解三角方阵
for r = 1 : n
    for j = r : n
        if r > 1;
            for k = 1 : r-1
                A(r,j) = A(r,j) - A(r,k)*A(k,j);
            end
        end
    end
    for i = r+1 : n
        if r < n;
            for k = 1 : r-1
                A(i,r) = A(i,r) -A(i,k)*A(k,r);
            end
            A(i,r) = A(i,r)/A(r,r);
        end
    end
    disp(['A(', r, '): ']); A
end

%回代过程（略）
%disp(['The maximum of the HRIR is',num2str(123456)])