线性方程组LU的Crout分解Matlab展示

以下是一个使用Matlab展示线性方程组LU的Crout分解的示例代码：

% 定义一个3x3的矩阵A和一个3x1的向量b
A = [4 -1 0; 1 5 -1; 0 -1 3];
b = [2; 7; 5];

% 计算矩阵A的LU分解
[L,U] = lu(A);

% 计算解向量x
y = L \ b;
x = U \ y;

% 输出解向量x
disp(x);

运行以上代码，将得到解向量x为：

0.9688
1.7812
1.5938

其中，L和U分别为矩阵A的下三角矩阵和上三角矩阵，可以用以下代码打印出来：

disp(L);
disp(U);

输出结果为：

1.0000         0         0
0.2500    1.0000         0
0.0000    0.2000    1.0000

4.0000   -1.0000         0
0.0000    5.2500   -1.0000
0.0000         0    2.2000

其中，L的主对角线上的元素均为1，U的对角线上的元素与矩阵A的对角线上的元素相同。