复化辛普森公式

/*
 * 复化梯形公式
 * 复化辛普森公式
 */
#include <iostream>
#include <iomanip>
#include <cmath>
using namespace std;
/*
 * 复化梯形公式
 * a表示区间起点
 * b表示区间终点
 * n表示区间长度
 */
double Trapezoid(double(*fun)(double x), double a, double b, int n)
{
   
	double result = 0.0;
	if(a > b)
		swap(a, b);
	result = (fun(a) + fun(b))/2;
	double h = (1.0/n