n×nの実対称行列Aの固有値を大きい順にr<n個得たいです.lapackを使うと val, vec = eigen(Symmetric(A), n-r+1:n) でOKです.この時,eigenの第一引数はSymmetric型じゃないといけないので注意.(ところで,ここでeigenではどんなアルゴリズムが動いてるんですかね?全部の固有値を求める方法は調べると出てきますが,大きい順に,っていう時は何が裏で動いているんでしょう...?知っている人がいたら教えてほしいです.) またArpackを使って, using Arpack val, vec=eigs(A, nev=r, which=:LR); …