解线性方程组就就是解一个矩阵方程Ax=b.A为系数矩阵，x为未知的列向量，b为一个常数列向量。
手工解矩阵方程一般都是进行初等变换，转化成最简行向量，然后得出一个基础解系。
请问如何在c++中解一个矩阵方程。是不是也进行初等变换，如果是这样，那变换的终止条件是什么，如何实现初等变换。
（此处解得矩阵方程是有解的，行列式D不为0的）
本人在编一个实现拟合函数的代码中遇到了这个问题，还希望有大虾来拉一把。