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
杜利特尔(直接三角分解法)的Matlab
最新推荐文章于 2022-03-23 20:05:17 发布
![](https://img-home.csdnimg.cn/images/20240711042549.png)