石黒/日誌/2013-05-10
をテンプレートにして作成
[
トップ
] [
新規
|
一覧
|
単語検索
|
最終更新
|
ヘルプ
|
ログイン
]
開始行:
*1変数のニュートン法のコード [#f5680ced]
この関数だとループ2回で終わる。前回の勾配法よりはるかに速い
from sympy import *
f = lambda x: (x-5)**2 + 3
df = lambda x: 2*(x-5)
ddf = lambda x: 2
def newton(df, ddf):
#1.
x = 1
#2.
while True:
_x = x
x = _x - df(_x) / float(ddf(_x))
#3.
if abs(x - _x) < 0.000001: #epsilon
break
return x
print newton(df, ddf)
終了行:
*1変数のニュートン法のコード [#f5680ced]
この関数だとループ2回で終わる。前回の勾配法よりはるかに速い
from sympy import *
f = lambda x: (x-5)**2 + 3
df = lambda x: 2*(x-5)
ddf = lambda x: 2
def newton(df, ddf):
#1.
x = 1
#2.
while True:
_x = x
x = _x - df(_x) / float(ddf(_x))
#3.
if abs(x - _x) < 0.000001: #epsilon
break
return x
print newton(df, ddf)
ページ名: