Python ニュートン法を実装する - Qiita

ニュートン法とは 1変数では初期値 $x_0$ を設定し、以下の数列によって関数 $f$ の近似解を求める方法。 多変数では初期値 $\boldsymbol{x_0} \in \mathbb{R}^n$ を設定し 、関数 $f: \mathbb{R}^n \rightarrow \mathbb{R}^n$ の近似解を求める方法。 $$\boldsymbol{x_{n + 1}} = \boldsymbol{x_n} - H^{-1} f(\boldsymbol{x_n})$$ ここで $H$ はヘッセ行列。 1変数 こちらに置きました。 import sympy as sp import numpy as np import matplotlib.pyplot as plt def newton_method(f, x, init_value, num): df = sp.diff(f,