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

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

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

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

for

for文は、様々なプログラミング言語で使われている制御構造です。for文に定義している条件から外れるまで、for文内の命令文を繰り返し実行します。

ループ

ループとは、プログラミングにおいて、条件に合致している間、複数回繰り返し実行される箇所や、その制御構造を指します

関数

関数(ファンクション・メソッド・サブルーチンとも呼ばれる)は、はプログラムのコードの一部であり、ある特定のタスクを処理するように設計されたものです。

関数型プログラミング

関数型プログラミングとは、関数を用いて演算子を構築し、算出し、コンピュータプログラムを構成する枠組みです。

Q&A

解決済

1回答

949閲覧

forループで求めた和が反映されない

Ryuuse

総合スコア27

C

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

for

for文は、様々なプログラミング言語で使われている制御構造です。for文に定義している条件から外れるまで、for文内の命令文を繰り返し実行します。

ループ

ループとは、プログラミングにおいて、条件に合致している間、複数回繰り返し実行される箇所や、その制御構造を指します

関数

関数(ファンクション・メソッド・サブルーチンとも呼ばれる)は、はプログラムのコードの一部であり、ある特定のタスクを処理するように設計されたものです。

関数型プログラミング

関数型プログラミングとは、関数を用いて演算子を構築し、算出し、コンピュータプログラムを構成する枠組みです。

0グッド

0クリップ

投稿2019/12/18 07:42

編集2019/12/18 07:45

関数f(x)=4/x^2+1を区間0<=x<=1で積分したときの値をシンプソンの公式を用いて計算するプログラムを作っています。
シンプソンの公式を使えば、台形公式よりも高い精密度で計算することができます。公式については添付された画像を参照して下さい。
この積分値はπであり、プログラムを実行すればπに近い値になるはずです。

プログラムを実際に作成して実行したところ、偶数部分の和だけ反映されず、πとは全く異なる値になってしまいました。
この原因が分かれば教えて下さい。よろしくお願いします。
プログラムは画像の下側にあります。

イメージ説明イメージ説明イメージ説明

/* 被積分関数 */ double f(double x) { /return 4/(1+(x*x)); /* 関数f(x)の値を計算して返す */ } double integ(double a, double b, int n) { /* * ★ 関数f(x)の区間[a,b]における積分を *  分割数nのシンプソン公式で求め、その値を返す。 */ double I = 0; int m; double h = (b-a)/n; I += h*1/3*f(a)+h*1/3*f(b); for( m=1; m<=n-1; m += 2){ I += (4/3)*h*f(a+m*h); } for( m=2; m<=n-2; m += 2){ I += 2/3*h*f(a+m*h); } return I; } int main(void) { double a, b; int n; a = 0; b = 1; /* nを4から20まで4刻みで変えながら繰り返す */ for (n=4; n<=20; n += 4) { printf("%d %10.7f\n", n, integ(a, b, n)); } return 0; }

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

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

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

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

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

guest

回答1

0

ベストアンサー

C

1I += 2/3*h*f(a+m*h);

この式の右辺は左結合となりますが、2/3整数除算なので0となります。0に何をかけようが0です。

投稿2019/12/18 07:45

maisumakun

総合スコア145123

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

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

Ryuuse

2019/12/18 07:53

回答ありがとうございます。「I += 2*h*f(a+m*h)/3;」としたところ、正しい 値を得ることができました。奇数部分についても同様でした。
maisumakun

2019/12/18 07:57

よく使う技として、「2.0/3*(後略)」のように数値リテラルを浮動小数点数にしてしまう方法があります。
guest

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

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

ただいまの回答率
85.50%

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

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

質問する

関連した質問