を解きたいシチュエーションに出会った。これは一般にn次方程式と呼ばれている気がする。(1変数の代数方程式ともいうのか) これは5次以上のものについては代数的な解がないことが知られている。 とはいっても5次についてはいろんな特殊関数を使えば表現できるようだし、近似解については数値計算もできる。 4次以上についてはそもそも数値計算を用いるのが普通のようだ。 Matlabなんかではrootsという代数方程式の数値計算を行う関数が用意されているようだが、pythonでは特にライブラリが見当たらなかった。というわけで実装してみることにした。 アルゴリズム 解を求めるとなるとNewton法が最初に思いつくように思われるが、n次方程式の解を求めるNewton法に近いアルゴリズムとしてDKA法というものが知られているようだ。ただこのアルゴリズムを実装するのは面倒だし、Matlabの実装をみると、固有値計算
![pythonでn次方程式の解を求める - Qiita](https://cdn-ak-scissors.b.st-hatena.com/image/square/19f0141ab18ff17db1c4d73e67c6a444ffac22b1/height=288;version=1;width=512/https%3A%2F%2Fqiita-user-contents.imgix.net%2Fhttps%253A%252F%252Fcdn.qiita.com%252Fassets%252Fpublic%252Farticle-ogp-background-9f5428127621718a910c8b63951390ad.png%3Fixlib%3Drb-4.0.0%26w%3D1200%26mark64%3DaHR0cHM6Ly9xaWl0YS11c2VyLWNvbnRlbnRzLmltZ2l4Lm5ldC9-dGV4dD9peGxpYj1yYi00LjAuMCZ3PTkxNiZoPTMzNiZ0eHQ9cHl0aG9uJUUzJTgxJUE3biVFNiVBQyVBMSVFNiU5NiVCOSVFNyVBOCU4QiVFNSVCQyU4RiVFMyU4MSVBRSVFOCVBNyVBMyVFMyU4MiU5MiVFNiVCMSU4MiVFMyU4MiU4MSVFMyU4MiU4QiZ0eHQtY29sb3I9JTIzMjEyMTIxJnR4dC1mb250PUhpcmFnaW5vJTIwU2FucyUyMFc2JnR4dC1zaXplPTU2JnR4dC1jbGlwPWVsbGlwc2lzJnR4dC1hbGlnbj1sZWZ0JTJDdG9wJnM9MmFmZmRiODU4ZDBjYmVlNjc3NjFhZGFjZTRmMDNhNzU%26mark-x%3D142%26mark-y%3D112%26blend64%3DaHR0cHM6Ly9xaWl0YS11c2VyLWNvbnRlbnRzLmltZ2l4Lm5ldC9-dGV4dD9peGxpYj1yYi00LjAuMCZ3PTYxNiZ0eHQ9JTQwY2FwcmVzdCZ0eHQtY29sb3I9JTIzMjEyMTIxJnR4dC1mb250PUhpcmFnaW5vJTIwU2FucyUyMFc2JnR4dC1zaXplPTM2JnR4dC1hbGlnbj1sZWZ0JTJDdG9wJnM9Y2M1MDllNzhlMDYwZDRmY2FlMTE1OTVhZGZhMGFkMzk%26blend-x%3D142%26blend-y%3D491%26blend-mode%3Dnormal%26s%3D6a11d166eec127bdcf5a399c227d39fe)