ニュートン法で高次方程式を解く

数値解法で方程式を解くときにpythonでニュートン法のライブラリが使える。プログラムの例は以下だ。

from scipy.optimize import newton

salary = 30
total = 30000
year = 40
def f(x): #方程式を記述
    return salary*12*(x+1)**year - salary*12*(x+1) - total*x
root = newton(f, 0.1) #0.1は探索の初期値
print(root) #答

必要なscipyモジュールをmacにインストールするには以下のコマンドを使う。依存関係で何か足りないと言われたら足りないものも同じコマンドでインストールすればよい。

python3 -m pip install --upgrade scipy