博客

下面是两种办法

（一）函数表达式已知

（二）读入一组数据构成一条曲线，求积分

#include <iostream>

using namespace std;
//方法1，曲线的函数可以用表达式表示的情况
double func1(double x) {
    return x * x;
}

double integral(double(*f)(double), double min, double max) {
    double result = 0;
    const int N = 1000;//N是分割的份数
    double delta = (max - min) / N;//min是积分的下限，max是积分的上限
    for(double i = min + delta; i < max; i += delta) {
        result += f(min + i) * delta;
    }
    return result;
}

//方法2，读入一组数据构成一条曲线用梯形法进行计算
double Area(double *y, double *t, int n) {
    int i;
    double sum = 0, a;
    for(i = 0; i < n - 1; i++) {
        a = (y[i] + y[i + 1]) * (t[i + 1] - t[i]) / 2.0;
        sum = sum + a;
    }
    return sum;
}

int main() {
    double result = 0;
    //下面调用方法1
    result = integral(func1, 0, 2.0);
    cout << result << endl;
    //下面调用方法2
    double area = 0;
    double x[] = {0, 1, 2, 3, 4, 5};
    double y[] = {0, 3, 7, 2, 1, 9};
    area = Area(y, x, 2);  //参数 y，x是数组，n是点数。
    cout << area << endl;
    return 0;
}