我是编程新手,目前我正在做一个需要使用Fortran 95的项目。是否有任何用于求解线性方程的子程序,例如在A*x=b
中查找矩阵x
,其中A
是 2*2 矩阵。
如果您给我任何有用的链接来帮助我解决这个问题,我将不胜感激。感谢您的帮助
对于 2x2 方程组,您应该编写 Cramer 规则,因为 det(A) 的表达式相当简单(也适用于 3x3)。https://en.wikipedia.org/wiki/Cramer%27s_rule
您可以使用许多库。事实上的标准是 LAPACK 包,您可以选择更多算法。有许多相同子程序的免费和商业实施,例如 OpenBLAS、Intel MKL 或 Sun Performance Library。
如果您的系统非常大,则需要寻找一些迭代求解器。其中有很多,只需搜索GMRES,BiCGSTAB或类似的方法及其实现(例如,http://people.sc.fsu.edu/~%20jburkardt/f_src/templates/templates.html)。