杜利特尔（直接三角分解法）的Matlab

function [L,U,y,x] = Doolittle(A, b) 
[row_a, col_a] = size(A);
for j = 1:col_a    
 U(1,j) = A(1,j);
end     L(1,1) = 1;
for i = 2:row_a    
L(i,1) = A(i,1)/A(1,1); 
end 
for i = 2:row_a  
    for j = i:col_a 
        temp_sum = 0;
        for