我新手在編程和此刻的我工作的一個項目,我需要使用的Fortran 95有求解線性方程,例如在發現矩陣x
任何子程序A*x=b
其中A
是2 * 2矩陣。求解的Fortran線性方程組95
如果你給我任何有用的鏈接,幫我解決這個問題,我將不勝感激。 感謝您的幫助
我新手在編程和此刻的我工作的一個項目,我需要使用的Fortran 95有求解線性方程,例如在發現矩陣x
任何子程序A*x=b
其中A
是2 * 2矩陣。求解的Fortran線性方程組95
如果你給我任何有用的鏈接,幫我解決這個問題,我將不勝感激。 感謝您的幫助
有您可以使用許多庫。事實上的標準是LAPACK包含更多算法,您可以選擇。相同子例程有許多免費和商業實現,例如OpenBLAS,Intel MKL或Sun Performance Library。
如果你的系統是非常大的,你會尋找一些迭代求解。其中有很多,只要搜索GMRES,BiCGSTAB或類似的方法及其實現(例如,http://people.sc.fsu.edu/~%20jburkardt/f_src/templates/templates.html)。
對於方程你應該在代碼的克萊姆法則的2x2的系統中,因爲DET(A)的表達式是相當簡單的(也爲3×3)。 https://en.wikipedia.org/wiki/Cramer%27s_rule
我不同意這個答案,但我認爲對於一個2 * 2矩陣,我只是寫下Fortran方程式的解決方案並繼續。 –
是的,但我認爲OP不限於2 * 2,只是用它來表示他通過線性系統表達的意思。 –