2011-07-01から1ヶ月間の記事一覧

ベアストウ法

ベアストウ法を勉強したのでPythonで書いてみた。 ベアストウ法を使うと、代数方程式にしか使えないけど複素解も含むすべての解を求めることができる。無理やり2次式で因数分解して積の形で表すって感じになる。詳しくはググってください。コードを書く際、…

Pythonの再帰の深さについてあれこれ

この前、コードを書いていたら次のエラーが出た。 RuntimeError: maximum recursion depth exceededどうやら再帰の回数の限界を超えてしまったようだ。で、調べてみるとこれらが引っかかった。・http://hep.planet-koo.com/index.php?g=tech&sid=py_reclimit…

2変数のニュートン法

今度は2変数についてのニュートン法をやってみた。今度は2変数のテイラー展開を利用する。詳しくはググってください。 今回用いた式は某教科書に例として載っていたものをそのまま利用した。 import math import sys chk = 0.00000001 #誤差の範囲 count = 0…

ニュートン法

Pythonでニュートン法をやってみた。ニュートン法って何?って方はググってください。 def newton(x, num): #nextx = (2 * (x*x*x) + x*x - 3) / (3 * (x*x) + 2*x - 5) #x - (g / dg)をまとめた場合 g = f(x) dg = df(x) nextx = x - (g / dg) num.append(n…