Skip to content

Instantly share code, notes, and snippets.

@komasaru
Created November 27, 2017 04:42
Show Gist options
  • Save komasaru/06220ef30c9d447153f633ce24b4edbc to your computer and use it in GitHub Desktop.
Save komasaru/06220ef30c9d447153f633ce24b4edbc to your computer and use it in GitHub Desktop.
C++ source code to compute definite integral by trapizoidal rule.
/*********************************************
* 台形則による定積分
*********************************************/
#include <iostream> // for cout
#include <math.h> // for sqrt
#include <stdio.h> // for printf()
# define f(x) sqrt(4-x*x) // 被積分関数
using namespace std;
/*
* 計算クラス
*/
class Calc
{
// 各種定数
static const int m = 100; // 積分区間分割数
// 各種変数
double h; // 1区間の幅
double x; // X 値
double s; // 面積
int k; // ループインデックス
public:
// 定積分計算
void calcIntegral(double a, double b);
};
/*
* 定積分計算
*/
void Calc::calcIntegral(double a, double b)
{
// 1区間の幅
h = (b - a) / m;
// 初期化
x = a; // X 値を a で初期化
s = 0; // 面積初期化
// 計算
for (k = 1; k <= m-1; k++) {
x = x + h;
s = s + f(x);
}
s = h * ((f(a) + f(b)) / 2 + s);
// 結果表示
printf(" /%f\n",b);
printf(" | f(x)dx = %f\n", s);
printf(" /%f\n",a);
}
/*
* メイン処理
*/
int main()
{
// 積分区間
double a, b;
try
{
// データ入力
cout << "積分区間 A, B :";
scanf("%lf %lf", &a, &b);
// 計算クラスインスタンス化
Calc objCalc;
// 定積分計算
objCalc.calcIntegral(a, b);
}
catch (...) {
cout << "例外発生!" << endl;
return -1;
}
// 正常終了
return 0;
}
@FelixOyoo
Copy link

thank you

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment