質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

新規登録して質問してみよう
ただいま回答率
85.48%
Python 3.x

Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

Q&A

1回答

772閲覧

Pythonによるニュートン法の実装

yield

総合スコア12

Python 3.x

Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

0グッド

0クリップ

投稿2019/10/07 12:24

発生している問題・エラーメッセージ

x**2-2をニュートン法で解くプログラムを作成しましたがうまく動作しません。 実行結果 <0>: f(2.0) = 4.0 <1>: f(2.0) = 4.0 <2>: f(2.0) = 4.0 <3>: f(2.0) = 4.0 <4>: f(2.0) = 4.0 <5>: f(2.0) = 4.0 <6>: f(2.0) = 4.0 <7>: f(2.0) = 4.0 <8>: f(2.0) = 4.0 <9>: f(2.0) = 4.0 <10>: f(2.0) = 4.0

python

1import numpy 2#--------------------------------------- 3__以下コード__ 4#--------------------------------------- 5#--------------------------------------- 6# ここを書きかえます 7#--------------------------------------- 8def f(x): 9 # 次の行を正しく書きかえてください 10 return x**2-2 11 12#--------------------------------------- 13# f'(x) のこと.ここは書きかえません! 14#--------------------------------------- 15def df(x): 16 dx = 1.0e-10 17 return (f(x+dx)-f(x)) / dx 18 19#--------------------------------------- 20# ここを書きかえます 21#--------------------------------------- 22def g(x): 23 # 次の行を正しく書きかえてください 24 return 2*x 25 26#--------------------------------------- 27# Newton法の途中経過を表示する関数 28#--------------------------------------- 29def show(x,k): 30 print("<" + str(k) + ">: f(" + str(x) + ") = " + str(f(x))) 31 32#--------------------------------------- 33# Newton法.ここは書きかえません! 34#--------------------------------------- 35def Newton_itr(x0): 36 K_MAX = 100 37 EPS = 1.0e-10 38 x = x0 39 k = 0 40 show(x0,k) 41 if (abs(f(x)) <= EPS): 42 return (x0,k) 43 while (k <= K_MAX): 44 x = g(x) 45 k = k + 1 46 show(x,k) 47 if (abs(f(x)) <= EPS): 48 break 49 return (x,k) 50 51#--------------------------------------- 52# main.ここは書きかえません! 53#--------------------------------------- 54if __name__ == '__main__': 55 x0 = 2.0 56 x, k = Newton_itr(x0)

前提・実現したいこと

ニュートン法の正常な動作。

気になる質問をクリップする

クリップした質問は、後からいつでもMYページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

guest

回答1

0

ニュートン法ならg(x)は x-f(x)/f'(x)にしないと正しく動かないと思います。

投稿2019/10/07 13:34

gnbrganchan

総合スコア438

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

まだベストアンサーが選ばれていません

会員登録して回答してみよう

アカウントをお持ちの方は

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
85.48%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問