n阶方阵的克劳特（crout）分解matlab实现

代码说明：

应用背景n阶方阵的克劳特（crout）分解matlab实现关键技术在矩阵分析课上学习了n阶方阵的三角分解或叫克劳特(crout)分解,便尝试着用MATLAB编程实现。给一个例子如下： 原矩阵a=[2 1 -5 1;1 -3 0 -6;0 2 -1 2;1 4 -7 6] 执行 [l,u]=crout(a); 分解后得： l=[2 0 0 0;1 -3.5 0 0;0 2 0.4286 0;1 3.5 -2 -9] u=[1 0.2 2.5 0.5;0 1 -0.7143 1.8571;0 0 1 -4;0 0 0 1]