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

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

新規登録して質問してみよう
ただいま回答率
85.50%
C

C言語は、1972年にAT&Tベル研究所の、デニス・リッチーが主体となって作成したプログラミング言語です。 B言語の後継言語として開発されたことからC言語と命名。そのため、表記法などはB言語やALGOLに近いとされています。 Cの拡張版であるC++言語とともに、現在世界中でもっとも普及されているプログラミング言語です。

Q&A

解決済

4回答

1757閲覧

C言語 Linux

uv-

総合スコア26

C

C言語は、1972年にAT&Tベル研究所の、デニス・リッチーが主体となって作成したプログラミング言語です。 B言語の後継言語として開発されたことからC言語と命名。そのため、表記法などはB言語やALGOLに近いとされています。 Cの拡張版であるC++言語とともに、現在世界中でもっとも普及されているプログラミング言語です。

0グッド

0クリップ

投稿2017/01/06 07:24

f(x)=xInx-1=0の解をNewton法で求める問題です。
Gnome端末でコンパイルしています。下記のようにプログラミングした所、
x= nanと出てしまいました。正解は1.76です。間違いを指摘していただければありがたいです。

#include<stdio.h>
#include<math.h>
#define F(x) (x*log(x)-1.)
#define FD(x) (log(x)+1.)
#define EPS (1.E-5)

int main(void)
{

double x=0.0,fx,fdx,dx;
int n=0;

do{
fx=F(x);
fdx=FD(x);
dx=-fx/fdx;
n++;
x+=dx;
}
while(fabs(dx/x)>EPS);
printf("x=%4.2f\n",x);

}

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

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

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

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

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

guest

回答4

0

こんにちは。

log(x)はx=0の時、未定義です。
なので、最初にx=0で開始しているのがよろしくないです。いきなりlog(x)がinf(無限大)になってます。

投稿2017/01/06 07:37

Chironian

総合スコア23272

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

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

0

ベストアンサー

x の初期値が悪いだけだと思います。ニュートン法の場合、求めたい解のほうに向かってグラフが傾いている必要があります。この場合だと、 x > 1/e でないと解に向かって収束しません。特に x=0.0 は関数の値が計算できないところですので、初期値にはできません。

x=0.9 からスタートしてみたら、動きました。

投稿2017/01/06 07:41

mit0223

総合スコア3401

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

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

0

ニュートン法ってよく分からないけど、xの初期値が0なのが問題では? log(0)って有限の値じゃないですよね?

投稿2017/01/06 07:41

majiponi

総合スコア1720

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

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

0

xに0.0をセットして計算しているからだと思います。
0除算が発生して浮動小数点以外になっているため、結果が
nanとなっていると思われます。
初期値やロジックを見直してみてはどうでしょうか。

ちなみに変数nは意味を持たないですね。

投稿2017/01/06 07:39

Take-y

総合スコア91

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

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

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

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

ただいまの回答率
85.50%

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

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

質問する

関連した質問