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


// 計算条件のパラメータ定数
#define A 0.0	// 積分区間の下端
#define B 1.0	// 積分区間の上端
#define N 10	// 微小区間の刻み数（値を大きく＝細かく刻むほど高精度）

// 積分したい関数 f(x) = cos(x)
double f(double x) {
	return cos(x);
}


int main(void) {
	
	// --------------------------------------------------
	// 以下、シンプソン法による数値積分の処理
	// --------------------------------------------------
	
	double delta = (B - A) / N;	// 台形の幅 Δx
	double value = 0.0;	// 計算結果の積分値(区間全体の面積)
	int i;
	
	for(i=0; i<N; i++) {
		
		// i番目の台形の左下X値
		double x = A + i * delta;
		
		// i 番目の地点までの積分値を出力
		printf("%f\t%f\n", x, value);
		
		// 積分値に i 番目の微小区間の面積を足す
		value += ( f(x) + f(x+delta) + 4.0 * f(x+delta/2.0) ) * delta / 6.0;
	}
	
	// 最終的な積分結果(区間全体の面積)をコンソールに出力
	printf("%f\t%f\n", B, value);
	
	return 0;
}