scilabには複素行列などの構造的な変数が最初から入れられるため

for文を使う用事にあまり出くわさない気がする。


でもこれならどうだ

この解をピンポイントで求めて、nを上げてってn次方程式の解を次々表示するなんてのはどうだ！？


いまんとこscilabのグラフ機能はいまいちわかってないから、Excelにデータを移そうと思ってんだけど
その際にfor文を使う用事ができそうな気がする！



たとえば、そうだな。疑似言語で書くと

for(n++){

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


みたいなことをやりたい

reとimはベクトルじゃなく行列になるけど、領域の半分ぐらいが空になる感じ。

もう複素二分法のためのスツルム列とか考える意味なくなっちゃうよね・・・そこなんだよな、プログラミングのちょっと虚しいところ

せっかくこの図のソースになるExcelファイルを改良しようと思ってたのに、意欲が削がれる削がれる