以前質問させていただき実際にプログラミングを作ったのですが、なぜか間違った答えが出てしまいます。
区分求積法・台形公式・シンプソンの公式を用いて、1/1+x*xを求めたいのですが、
1)台形公式の答えが区分求積法の答えより精度がが悪くなってしまう。
2)シンプソン公式が答えに収束しない。
となってしまいます。
以下がそのプログラム↓
#include <stdio.h>
#define FROM 0.0
#define TO 1.0
double func(double x)
{
double out;
out = 1.0 / ( 1.0 + x * x );
return (out);
}
double kubun(double start, double end, int num)
{
int i;
double h, s;
h = ( end - start ) / num;
s = 0.0;
for(i=0; i<num; i++) s += func( start + i * h + h / 2.0 );
return ( s * h );
}
double daikei(double start,double end,int num)
{
int i;
double h,s;
h = ( end - start ) / num;
s = 0.0;
for(i=1; i<num-1; i++) s += func( start + i * h );
return ((func(start) / 2.0 + s + func(end) / 2.0) * h );
}
double simpson(double start,double end,int num)
{
int i;
double h,s;
h = ( end - start ) / num;
s = 0.0;
for(i=1;i<num;i+=2)
s += 4.0 * func(start + h * i);
for(i=2;i<num;i+=2)
s += 2.0 * func(start + h * i);
return ( (func(start) + s + func(end))/ 3 );
}
int main()
{
double func(double);
double kubun(double, double, int);
double daikei(double, double, int);
double simpson(double, double, int);
printf("\n");
printf("### Square Integration\n");
printf(" ++ Partition = 10\t Answer = %10.6f\n", kubun(FROM, TO, 10));
printf(" ++ Partition = 50\t Answer = %10.6f\n", kubun(FROM, TO, 50));
printf("\n");
printf("### daikei Integration\n");
printf(" ++ Partition = 10\t Answer = %10.6f\n", daikei(FROM, TO, 10));
printf(" ++ Partition = 50\t Answer = %10.6f\n", daikei(FROM, TO, 50));
printf("\n");
printf("### simpson Integration\n");
printf(" ++ Partition = 10\t Answer = %10.6f\n", simpson(FROM, TO, 10));
printf(" ++ Partition = 50\t Answer = %10.6f\n", simpson(FROM, TO, 50));
return (0);
}
画面に表示する際に、それぞれ分割数を10と50にした際の値を表示するように作りました。
細かい点までご指摘いただけると幸いです。
よろしくお願いします。
No.1ベストアンサー
- 回答日時:
1)台形公式の答えが区分求積法の答えより精度がが悪くなってしまう。
for(i=1; i<num-1; i++) s += func( start + i * h );
↓
for(i=1; i<num; i++) s += func( start + i * h );
2)シンプソン公式が答えに収束しない。
return ( (func(start) + s + func(end))/ 3 );
↓
return ( (func(start) + s + func(end))/ 3 * h);
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- C言語・C++・C# プログラミングの授業の課題です 1 2023/01/17 22:15
- C言語・C++・C# 並列プログラミングのπ計算について 1 2022/07/16 22:30
- C言語・C++・C# c言語でユーザ関数を利用して複素数のべき乗と絶対値の数列を計算するプログラムが作りたいです。 3 2023/01/29 22:13
- C言語・C++・C# C 言語の Gauss Jordan 法について 2 2022/12/28 11:16
- C言語・C++・C# Cのdoubleの浮動小数点表示について 3 2023/04/17 13:14
- C言語・C++・C# 10個の実数に対する降順ソート結果を出力するプログラムを作りたいのですが、以下のプログラムをどう直せ 1 2022/07/09 22:16
- C言語・C++・C# C言語のマクローリン展開ローラン展開のコードについて 3 2022/12/15 14:45
- C言語・C++・C# バイナリファイルをコピーするのにかかる時間を測りたいのですが実行するとFatel error:gli 2 2022/11/03 01:10
- C言語・C++・C# LU分解法のピボッティングについて(C言語/gcc-9) 3 2022/07/11 23:10
- C言語・C++・C# (C言語)めちゃくちゃな値になってしまいます。 5 2022/08/13 11:55
関連するカテゴリからQ&Aを探す
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
プログラムでの数字につく”f”の...
-
C言語 関数プロトタイプ宣言の...
-
doubleの変数にintとintの割り...
-
(C言語)ニュートン法のプログ...
-
C言語のデバック 領域の二重解...
-
C++で割り算の結果を昇順に出力...
-
C言語のプログラムで質問です。
-
C言語で内積、、、わかりません。
-
関数におけるif文とreturn文に...
-
初めてCを触って、困っています
-
-1.#IND00と出てしまうのですが...
-
C言語の型による処理速度の違い
-
atan2 オーバーロード関数の呼...
-
float型とdouble型の変数の違い...
-
C言語のpow関数の不具合
-
c言語で、繰り返し文の中で、0....
-
C言語における複素数の四則演算...
-
C言語
-
#define _CRT_SECURE_NO_WARNIN...
-
【C++】関数ポインタの使い方
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
プログラムでの数字につく”f”の...
-
C言語を実行すると-infが出てき...
-
float型とdouble型の変数の違い...
-
doubleの変数にintとintの割り...
-
c言語で、繰り返し文の中で、0....
-
C言語 関数プロトタイプ宣言の...
-
関数におけるif文とreturn文に...
-
C 開放してるのにエラー(doubl...
-
C言語の型による処理速度の違い
-
至急です! マクロ定義で #defi...
-
doubleは常に%lfとするべきなのか
-
Cで3乗根を求める方法
-
-1.#IND00と出てしまうのですが...
-
C言語で-23乗を取り扱うには
-
C++で外積
-
2次方程式の解を求めるプログ...
-
方程式を2分法を用いて解くプロ...
-
ニュートン法
-
c言語のコンパイルエラー canno...
-
difftime()について
おすすめ情報