scilabで計算するフィボナッチ系501次方程式

a1=ones(1,500)
a2=eye(500,500)←2重ループなので多少時間かかります
a3=zeros(500,1)
A=[a1,1;a2,a3]←2重
x=poly(0,'x')
q=det(x*eye(A)-A)←ここで一番時間食いそうです
p=roots(q)
re=real(p)
im=imag(p)

これを慣れたExcelに張り付けて、
re
im
complex(re,im)
imargument←これでソートして、re対imでグラフ化


さすがにちょっと時間かかりましたね。パソコンが

ベンチマークとしていかがでしょうｗｗｗｗｗｗ


もうね、解が、ほぼほぼ線です。吹き出しみたくなってます

1001次は終わる気しなかったので止めました。


実数だけども、非対称行列だとspecとbdiagの結果が食い違いますねえ。
specのほうがもっともらしく見えるけど、
もしかしてどっちも合ってるなんてことあるかな？？
今度、実際に対角化して確かめてみましょうか。