coding Shift_JIS;	// 文字コードの明示（文字化け予防）
import Math;	// 数学関数ライブラリの読み込み
import tool.Graph2D;	// グラフ描画ライブラリの読み込み


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

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

// ---------------------------------------------
// 以下、台形法（台形近似）による数値積分の処理
// ---------------------------------------------

double delta = (B - A) / N;	// 台形の幅 Δx
double value = 0.0;	// 計算結果の積分値(区間全体の面積)

for(int i=0; i<N; i++) {
	
	// i番目の台形の左下X値
	double x = A + i * delta;
	
	// i 番目の地点までの積分値をコンソールに出力
	println(x, value); // ※ N が大きい時は重いのでコメントアウト
	
	// 積分値に i 番目の台形の面積を足す
	value += ( f(x) + f(x+delta) ) * delta / 2.0;
}

// 最終的な積分結果(区間全体の面積)をコンソールに出力
println(B, value);

// コンソール内容をグラフにプロット
// ※ N が大きい時は重いのでコメントアウト
setGraph2DData(newGraph2D(), load());