数值分析（六）：C++实现非奇异矩阵A的LU分解

矩阵的LU分解是指非奇异矩阵A能够分解成单位下三角矩阵L和上三角矩阵U的乘的形式，这部分内容是高斯消元法章节讲解的内容，因为最近开始复习数值分析这门功课，所以又回头来开始学习，于是就碰到了矩阵A的LU分解，看到这个分解就是要死磕一下，就想着用C++实现一下，于是就有了这篇博客，O(∩_∩)O哈哈哈~。比较简单，所以原计划是只安排半小时实现的，但是最后磨磨蹭蹭到了两个小时，中间磕磕绊绊出现了一些错误，最后还是慢慢发现并解决，调试运行得到了正确的答案。
同样，这次也不再繁琐地介绍LU分解的数学原理了，只是简单地带过一下，非奇异矩阵A的LU分解，形式如下图：

然后是怎么根据A矩阵得到LU的呢？是怎么根据A中的元素得到LU中的元素的呢？最基本的答案就是这是直接根据矩阵的乘法，即L的行乘以U的列计算的规则得到的（感觉自己像智障一样，这样回答问题真的是找打，，，，233333）当然最基本的还是矩阵的乘法计算规则，但是实际计算过程中当然不会这么傻乎乎地行乘以列，而是采用一种比较直接的，也比较适用于手算，且更适用于C++实现的简便运算。
先上代码，然后稍微简单地介绍一下这种简便方法。

#include<iostream>

using namespace std;

int main()
{
	int n;
	cout << "输入方阵的阶数n:" <<