[C++成长日记]使用梯形求积分的公式求积分

题目：已知用梯形法求积分的公式为： ，其中，h = (b-a) / n，n为积分区间的等分数，编程求如下积分的值。要求：把求积分公式编写成一个函数，并使用函数指针作为形式参数。调用该函数时，给出不同的被积函数作为实际参数求不同的积分。

例如：

结果分别为：3.14159        1.81009        1

by gzhu_fwm​​​​​​​

#include <iostream>
#include <cmath>
using namespace std;
#define n 10000//n为等分数
#define pi 3.141592657

typedef double fun(double);
fun fun1, fun2, fun3;//被积函数
double integral(fun* pfun, double a, double b);//求积分的函数

int main()
{
	cout << "求积分结果为：" << integral(fun1, 0, 1) << endl;
	cout << "求积分结果为：" << integral(fun2, 1, 2) << endl;
	cout << "求积分结果为：" << integral(fun3, 0, pi/2) << endl;
}
//被积函数
double fun1(double x)
{
	return 4 / (1 + pow(x, 2));
}
double fun2(double x)
{
	return sqrt(1 + pow(x, 2));
}
double fun3(double x)
{
	return sin(x);
}

double integral(fun* pfun, double a, double b)//求积分的函数,a，b为积分下限上限
{
	double h = (b - a) / n, result = 0;
	for (int i = 1; i <= n - 1; i++) result += pfun(a + i * h);
	result *= h;
	result += h * (pfun(a) + pfun(b)) / 2;
	return result;
}