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

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

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

	・実行結果の一部を下に示します．
				　　・
				　　・
			L=15.000000, n= 30, h=1.000,							|I-rPi|=6.48682e-04
			L=15.000000, n= 50, h=0.600,							|I-rPi|=3.57797e-07
			L=15.000000, n= 60, h=0.500,							|I-rPi|=1.21538e-06
			L=15.000000, n= 75, h=0.400,							|I-rPi|=1.24012e-06
			L=15.000000, n=100, h=0.300,							|I-rPi|=1.23277e-06
			L=15.000000, n=150, h=0.200,							|I-rPi|=1.22769e-06
			L=15.000000, n=300, h=0.100,							|I-rPi|=1.22463e-06

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

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

double f(double w){
	double y;
	y = 1/cosh(w);
	return y;
}

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

	for(L=iniL; 14<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-Pi|=%11.5e\n",									L, n, h, fabs(I-Pi));
		}
	}
	return 0;
}