なんとかボナッチ行列の固有値をより早く、より正確に

 
 


m=12
p=zeros(m,m)
for n=1:m
a1=ones(1,n)
a2=eye(n-1,n-1)
a3=zeros(n-1,1)
A=[a1;a2,a3]
x=poly(0,'x')
q=(x^(m-n))*det(x*eye(A)-A)
p(:,n)=roots(q)
re=real(p)
im=imag(p)
end

(mを変えるだけでいくらでも：だいたい500くらいまでだったら有限時間内に終わるんじゃないすかね)

scilabの力を借りると、こんなにもあっさりと解けてしまうんですねえ
フリー言語なのにすごいなぁ
解析的な代数方程式をなかば仮想的にでも定義できるし
それを複素数の範囲で解くこともできるなんて。


まあ、超越方程式の複素根とかはさすがに無茶ブリですよね＾＾；