fortran - 在 Fortran 95 中求解线性方程

Question

我是编程新手，目前我正在做一个需要使用 Fortran 95 的项目。是否有任何求解线性方程的子程序，例如x在2*2 矩阵中查找矩阵A*x=b。A

如果您给我任何有用的链接来帮助我解决这个问题，我将不胜感激。感谢您的帮助

Answer 1

对于 2x2 方程组，您应该编写 Cramer 规则，因为 det(A) 的表达式相当简单（也适用于 3x3）。 https://en.wikipedia.org/wiki/Cramer%27s_rule

Answer 2

您可以使用许多库。事实上的标准是LAPACK包，您可以选择更多算法。相同子例程有许多免费和商业实现，例如 OpenBLAS、Intel MKL 或 Sun Performance Library。

如果您的系统非常大，您会寻找一些迭代求解器。其中有很多，只需搜索 GMRES、BiCGSTAB 或类似方法及其实现（例如，http ://people.sc.fsu.edu/~%20jburkardt/f_src/templates/templates.html ）。