ニュートン法の実装 # (解を求める方程式, 初期値, 微分で用いる微小量, 許容(絶対)誤差, 反復の最大数) def newton_method (f, x0, h= 1e4, tol= 148e8, maxiter= 50) for cnt in range (1, maxiter 1) # 中心差分で微分を求める df = (f(x0 h) f(x0 h)) / (2 * h) next_x = x0 f(x0) / df # 反復回数, 解の近似値x, f(x)を表示 print ('{2dニュートン法(ニュートン・ラフソン法) とは︖ 非線形方程式f(x) = 0 を数値的に解く方法の1つ 微分可能な方程式であれば、たとえ微分しなくても解が求まる 繰り返し計算により解に近づく 初期値を変えて何回か解く方がよい 1 ocaml 微分 – OCamlでニュートン法(2) 僕の所属する研究室では毎年Software Foundationsの輪講をやっているんですが,Brzozowski derivativeに基づいて正規表現エンジンを実装し,あまつさえCoqで検証してしまおうといった練習問題が追加されていて大変興味深 ニュートン法の例題 ルート3の近似値を与える方法 理数アラカルト 微分 ニュートン法