55（数值积分）组合辛普森公式

组合辛普森公式（数值积分）
【问题描述】组合辛普森公式求f(x)=2+sin(2*sqrt(x))的积分近似值。

【输入形式】在屏幕上依次输入积分上限、下限和等距子区间个数。

【输出形式】输出使用组合辛普森公式求得的积分近似值，保留6位有效数字。

【样例1输入】

1 6 5

【样例1输出】

8.18302

【样例1说明】输入：积分上限a为1、下限b为6和等距子区间个数m为5。输出：积分近似值（保留6位有效数字）

【评分标准】根据输入得到的输出准确

#include<iostream>
#include<cmath>
using namespace std;
double a, b, m;//积分上限、下限、等距子区间个数  eg.1  6  5
double n;//积分近似值
double f(double x)//函数
{
	return 2 + sin(2 * sqrt(x));
}
double xinpu()
{
	double h = (b - a) / (2 * m);//h是每个子区间的宽度
	double n, n1 = 0, n2 = 0;//n是最后的近似值
	double x;
	//这里用的是第3个等价公式
	for (int i = 1; i < m+1; i++)//k是1~m
	{
		x = a + h * (2 * i - 1);//这是求点xk=a+h*(2*i-1)  奇
		n1 += f(x);
	}
	for (int j = 1; j < m; j++)//k是1~m-1
	{
		x = a + h * 2 * j;//这是求点xk=a+2*h*j  偶
		n2 += f(x);
	}
	n = h * (f(a) + f(b) + 4 * n1 + 2 * n2) / 3;
	return n;
}
int main()
{
	cin >> a >> b >> m;
	n=xinpu();
	cout << n;
	return 0;
}