方程式メーカーベンチマーク＠scilab

やる気が出ないのでぼーっと、行列2重ループに頼らない方程式メーカーを考えてました。
2378次方程式の時点で、スタックが足りないとか言い出して計算できなくなりました。
ハードの問題なのか、それともコードが要改良なのか、知るか。

m=2577
x=poly(0,'x')
y=zeros(m+1,1);
for n=0:1:m
　y(n+1)=x^n;
end

z=y(1);
for n=1:1:m
　z=z+y(n+1);
end

z=z-2*x^m;

w=roots(z);
re=real(w);im=imag(w);
w2=[re,im];


出てきたw2をExcelのcomplex関数で複素数(文字型)c=a+ibにしてやって、imargument(c)でソートしてやれば
線のグラフにしても編み物にならずシンプルにグラフ化されますよっと