C言語 Simpson法 積分

Question

プログラムについての質問です。 シンプソン法を用いて区間[0,1]でexpxを積分した値を求めるプログラムを作ろうとしているのですがエラーが出てしまっていてどの部分が間違っているのかわからないです。そしてこのプログラムについてですがEPSは分割数Nの場合と分割数2Nの場合の定積分値を比較しその差の絶対値がEPSよりも小さいときに定積分値の計算が終了しそれまでは分割数Nを二倍ずつ増加させる再帰的なプログラムとするためのものです。

#include<stdio.h>
 #include <math.h> 

double funk(double x);
double simpson(double a, double b, int n, double (*f)());

int main(void)
{
double *g;
double EPS=*g;
printf("sei no seisuu wo nyuuryoku ->",EPS); 
scanf("%d",EPS);

printf("exp(x) wo [0,1] de sekibun simasu. bunkatu suu ha %d desu \n",2*EPS); 
printf("result is %20.15f\n",simpson(0, 1.0, EPS, funk));
return 0;
}

double simpson(double a, double b, int n, double (*f)()) {

double S,h; 
int i;

h=(b-a)/(2.0*n);

S=((*f)(a)+(*f)(b));

for (i=1;i<n;i++)
{
S += 4.0*(*f)(a+(2.0*i-1.0)*h)+2.0*(*f)(a+2.0*i*h); }

S += 4.0*(*f)(a+(2.0*n-1.0)*h); S *= h/3.0; return S;
}

double EPS(double g)
{
int n;

*g = (simpson(0, 1.0, 2.0*n, funk)- simpson(0, 1.0, n, funk));

if(*g < EPS)

return simpson(double a, double b, int 2.0*n, double (*f)()); 

else return 1;
}

エラーメッセージ(コピペできなかったため写真を添付します)

よろしくお願いします

Comments

[y_waiwai]

このままではコードが読みづらいので、質問を編集し、<code>ボタンを押し、出てくる’’’の枠の中にコードを貼り付けてください

[y_waiwai]

エラーが出たなら、エラーメッセージを提示しましょう
エラーメッセージは、いらぬ翻訳省略しないで、出たそのママをコピペで提示してください

[KNTN]

teratailを初めて使用したため、使い方を十分に把握できていませんでした。申し訳ございません。

[episteme]

関数:EPS() がどこからも呼び出されていません。
関数:EPS() は再帰していません。

[KNTN]

指摘ありがとうございます。
EPS()の再帰にする方法を教えていただきたいです。

[episteme]

再帰する必要あるんですか? 分割数を増やしながら繰り返すだけちゃいますの?
そもそも呼ばれもしない関数定義して何の得があるんです?

[KNTN]

なるほど、確かに指示は分割数を増やすだけで、ESPの値は固定です。上記のコードもESPは再帰ではないですよね？関数:EPS() がどこからも呼び出されていません。の詳細を詳しくお聞かせ願います。

[episteme]

いや、main() の中に EPS() を呼んでるところがひとつもない。

[kazuma-s]

> 区間[0,1]でexpxを積分した値を求める

expx とは何ですか？
コード中の funk とは何ですか？
もしかすると、double funk(double x) { return exp(x); } ですか？

Answer 1

エラーメッセージの意味が分からないのですか？

39行目の *g の g は、35行目の double g ですから、型は double です。
単項演算子 * をポインタに適用して、演算結果がそのポインタの指すオブジェクト
になるべきなのに、「g がポインタではない」というエラーです。

41行目も同様。

43行目は関数simpson を呼び出して、返された値を return するつもりなのでしょうが、
関数呼出しの引数に書くのは、定数や変数や式などの値でなければなりません。
double a というのは値ではなく、宣言です。

次に、コードを main から見ていきます。

double *g; でポインタ変数 g を用意しましたが、初期化はしていないので
このポインタはどこを指しているのかまだ分かりません。

double EPS=*g; で doubleの変数 EPS を用意して、それをポインタ g が指す先の
double の値で初期化しようとしています。しかし、g はどこを指しているのか
分からないので、無意味な初期化です。

printf の第1引数の書式には "%d" や "%s" などの変換指定が無いので、
第2引数の EPS は無用です。

scanf の第1引数の書式の "%d" という変換指定は、第2引数に int * を要求します。
しかし、EPS は「int へのポインタ」ではなく、double です。

C のポインタが、全く理解できていないようです。
関数の宣言、関数の定義、関数の呼び出しも理解できていないようです。

追記
再帰ではなく、ループでやるとこんな感じでしょうか？

#include <stdio.h>  // printf
#include <math.h>   // fabs, exp

#define EPS  1e-6

double simpson(double a, double b, int n, double f(double))
{
	double h = (b - a) / (2 * n);
	double s = f(a) + f(b);
	for (int i = 1; i < n; i++)
		s += 4 * f(a + (2 * i - 1) * h) + 2 * f(a + 2 * i * h);
	s += 4 * f(a + (2 * n - 1) * h);
	s *= h / 3;
	return s;
}

int main(void)
{
	int n = 1;
	double y0 = simpson(0, 1, n, exp), y;
	for (n = 2; ; n *= 2) {
		//printf("n = %d, y0 = %f\n", n, y0);  // 途中経過を見る
		y = simpson(0, 1, n, exp);
		if (fabs(y - y0) < EPS) break;
		y0 = y;
	}
	printf("n = %d, y = %f\n", n, y);
}

コメントをお願いします。

Comments

[kazuma-s]

解決済みになっていますが、どのようなコードに修正したのですか？