public class IntegralForPlot {

	// =======================================================================
	// 以下、積分内容の設定
	// =======================================================================

	// 被積分関数 (積分対象の関数)
	static final double f(double x) {
		return x * Math.cos(x);
	}

	static final double A = 0.0;   // 積分区間の下端 (x軸上の左端)
	static final double B = 40.0;  // 積分区間の上端 (y軸上の左端)
	static final int N = 10000;    // 数値積分における区間分割数 (刻み数)


	// =======================================================================
	// 以下、データ出力関連の設定
	// =======================================================================

	static final int PRINT_POINTS = 1000; // 出力する座標点の数 (グラフのプロット点数)
	static final int PRINT_INTERVAL = N / PRINT_POINTS; // [自動設定] ループ内での出力周期


	// =======================================================================
	// 以下、積分の値を数値的に求める処理
	// =======================================================================

	public static void main(String[] args) {
		double delta = (B - A) / N;  // 足し上げる微小区間の幅 (刻み幅)
		double value = 0.0;		  // この変数に積分値を足し上げていく

		// 積分区間をN分割した微小区間を、左(下端)から右(上端)へ順に辿るループ
		for(int i=0; i<N; ++i) {
		
			// i 番目の微小区間における、左端の x 座標値を求める
			double x = A + i * delta;

			// 設定周期ごとに、グラフにプロットするための座標値を出力する
			if (i % PRINT_INTERVAL == 0) {
				System.out.println(x + "\t" + value);
			}

			// -------------------------------------------------------------------
			// 以下、i 番目の微小区間の面積を求め、積分値に加算 (どれか 1 つ選択)
			// -------------------------------------------------------------------

			// 矩形法 (長方形で近似した微小領域の面積を足す)
			// value += f(x) * delta;

			// 台形法 (台形で近似した微小領域の面積を足す)
			// value += ( f(x) + f(x+delta) ) * delta / 2.0;

			// シンプソン法 (二次関数で近似した微小領域の面積を足す)
			value += ( f(x) + f(x+delta) + 4.0 * f(x+delta/2.0) ) * delta / 6.0;
		}

		// 最後の一点を出力
		System.out.println(B + "\t" + value);
	}
}