ほぼ
NAGのサンプルのまま．
ライブラリとして使いやすいようサブルーチン化しただけ．
subroutine eig_sym(
a, m, w, z, info)
a [in]
real(8) array of shape (n,n) (n>0)
対称行列．
m [out]
integer
見つかった固有値の数．0 <= m <= n.
w [out]
real(8) array of shape (n)
固有値．
1から
m番目までの要素に固有値が入っている．
z [out]
real(8) array of shape (n,n)
固有ベクトル．
1から
m番目の列に固有ベクトルが入っている．
info [out]
integer
= 0: 成功
< 0:
info=-iのとき，
i番目の引数が不正（よく分からん）
> 0: 内部エラー
A = \begin{pmatrix}
a & b \\
c & d
\end{pmatrix}
結果