請愛用lapack
若不會調用，有人已包好:
http://www.netlib.org/
直接用就是。
好比原po要的是實矩陣的特徵值。
那查 eigenvalue 這就可以。
討驗的是，裡面代碼絕大多數都是以fotran編寫
所以還須要個編譯器.
幸運的是，英特爾作曲家，非商業使用不用錢。
http://software.intel.com/en-us/articles/intel-compilers/
有本事的人就直接用作曲家裡面的數學內核庫吧。
那是(編程)優化過的lapack
矩陣問題，還是交給專業的lapack求解。
這也是為何fotran至今不會死的根本原因啊。
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 182.235.79.253
※ 編輯: DrStein 來自: 182.235.79.253 (10/02 22:01)