[MATLAB]Doolittle分解

[MATLAB代码]关于杜利尔特Doolittle分解法求线性方程组的数值解

A = [2 -1 4 -3 1;
    -1 1 2 1 3;
    4 2 3 3 -1;
    -3 1 3 2 4;
    1 3 -1 4 4];
b = [11;14;4;16;18];
c = [1 1 1 1 1 ];
L = zeros(5) + diag(c);
U = zeros(5);
y = zeros(5,1);
x = zeros(5,1);
%Doolittle分解
for i = 1:5
    U(1, i) = A(1, i);
    L(i, 1) = A(i, 1)/U(1, 1);
end
d = 0;
e = 0;
for i = 2:5
    for j = i:5
        for k = 1:(i-1)
            d = d + L(i, k)*U(k, j);
            e = e + L(j, k)*U(k, i);
        end
        U(i, j) = A(i, j) - d;
        L(j, i) = (A(j, i) - e)/U(i, i);
        d = 0;
        e = 0;
    end 
end
%Gauss消去法求解Ly=b
y(1) = b(1);
for i = 2:5
    for j = 1:(i-1)
        d = d + y(j)*L(i, j);
    end
    y(i) = (b(i)-d)/L(i,i);
    d = 0;
end
%Gauss消去法求解Ux=y，x即Ax=b的解
x(5) = y(5)/U(5,5);
for i = 4:-1:1
    for j = 5:-1:(i+1)
        d = d + x(j)*U(i, j);
    end
    x(i) = (y(i)-d)/U(i,i);
    d = 0;
end