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

/***************** 表3.6の誤差を出力するプログラム **********************

	・表3.6の厳密な値のみの出力としましたが，＊行の下の条件文を
		取り除けば，その他のデータも取り出せます．

*************************************************************************/

#define Pi acos(-1.)		//πの値
#define rPi sqrt(Pi)		//ルートπの値
#define iniL 4			//Lの初期値(textに対応)
#define Lstep 1			//Lのステップ(初期値からの減算，正の整数)
#define nmax iniL*20		//区間[-L, L]の分割数の最大値(h=0.1まで)

double f(double w){
	double y;
	y = exp(-pow(sinh(w),2))*cosh(w);
	return y;
}

int main(void){
	int n, i;
	double I, L, h, xn;

	for(L=iniL; 1<L; L-=Lstep){
		for(n=1; n<nmax+1; n++){
			h = 2*L/n;
			I = 0.;
			for(i=0; i<n+1; i++){
				xn = -L+i*h;
				if(i==0 || i==n)
					I += .5*h*f(xn);
				else
					I += h*f(xn);
			}
//＊
			if( fmod(h*10, 1)==0 )		//hの10倍が整数のとき
				printf("L=%3lf, n=%3d, h=%5.3lf,							|I-rPi|=%11.5e\n",									L, n, h, fabs(I-rPi));
		}
	}
	return 0;
}