我有3x3和4x4矩阵A，我需要找到
的矩阵X解决方案A*X=B
其中B是矩阵，但A的行列式约为0.0！。 A的大小是3x3（或4x4），B的大小是3xr（或4xr）。由于A的行列式接近0.0，我不想（或者不能）使用A的精确矩阵逆来找到X或Cramer算法（因为除以矩阵的行列式）。
请注意
A*x(:,s) = B(:,s)
其中B（：，s）是矩阵B的列号s。所以，我的问题是解决一个带有几个右边项的线性系统。
我使用MKL库在fortran中运行我的代码（用于解决大型稀疏矩阵），因此我可以使用该库中的某些内容。
解决这个问题有什么好主意？也许LU分解？