マニュアルを引いたら、 大した苦労もなしに GVLSP なるルーチンが見つかった。 能書きに曰く、Purpose: Compute all of the eigenvalues of the generalized real symmetric eigenvalue problem A*z = w*B*z, with B symetric positive definite. これですね。 サンプル・プログラムも短いので、 IMSL の使い方の復習を兼ねて打ち込んで実行してみた。
コンパイルは
f77o sample393.f -limslib
で OK でした。実行して即 OK となりました。
うん、これで何とかなるでしょう。