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

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

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

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

Q&A

解決済

1回答

1001閲覧

【添削のお願い】二分法のプログラムってこれでいいと思います?

OHASI

総合スコア96

C

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

0グッド

0クリップ

投稿2021/10/16 12:27

前提・実現したいこと

【C言語】2x^2 - 12sin(x) = 0の解のうち,0より大きいものを二分法で精度10-6で求めるプログラムを書きました。

普通、検算結果は限りなく0に近づくはずですが、

実際やってみたら、検算結果がf(xm) = -1.172601 になるのでコードが当たっているのか心配です。

このコードであってると思いますか?

###前置き
はじめにf(x) = 2x^2 - 12sin(x)グラフはエクセルで作ったところ以下のようになったので、

初期値はそれぞれxa = 2,xb = 3としました。

グラフ

###ソースコード

c

1#include <stdio.h> 2#include <math.h> 3#define EPS 1.0e-6 4float f(float x); 5int main(void) 6{ 7 float xa,xb,xm; 8 printf("xa="); 9 scanf("%f",& xa); 10 printf("xb="); 11 scanf("%f",& xb); 12 printf("f(x)="); 13 printf("%.6f\n", f(xa)); 14 15 do{ 16 xm = (xa+xb)/2; 17 printf("%.6f\n", xm); 18 if( f(xa)*f(xb) < 0 ){ 19 xb = xm; 20 }else{ 21 xa = xm; 22 } 23 }while( xb-xa>EPS ); 24 25 printf(“xm = %.6f, xm); 26 printf("検算の結果\nf(xm)=%.6f\n", f(xm)); 27 printf("よって解はx=%.6f\n", xm); 28} 29 30float f(float x){ 31 float y; 32 y = 2*x*x-12*sin(x); 33 return(y); 34} 35

###結果
結果

以上、よろしくお願いします。

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

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

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

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

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

guest

回答1

0

ベストアンサー

c

1 if( f(xa)*f(xb) < 0 ){ 2 xb = xm; 3 }else{ 4 xa = xm; 5 }

範囲を狭める条件はこれで正しいのでしょうか?f(xm)を見ずに決めてしまうのが不適当と思われます。

投稿2021/10/16 12:37

maisumakun

総合スコア145208

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

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

OHASI

2021/10/16 17:37

分かりました。  ありがとうございました。
guest

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

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

ただいまの回答率
85.46%

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

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

質問する

関連した質問