#include <iostream>
#include <iomanip>
#include <math.h>
using namespace std;
#define F(x) exp(x)*cos(x)
int main()
{
long double f,x,a,b,h;
int n;
cout <<"请输入积分区间： \n"<<"a= ";
cin >>a;
cout<<"b= ";
cin >>b;
cout <<"\n"<<"请输入平分后区间的个数：n= ";
cin >>n;
f=F(a)+F(b);
h=(b-a)/n;
cout <<"用复化辛普生公式求解:"<<endl;
x=a;
for ( int i=0;i<n;i++)
{
x+=h/2;
f+=4*F(x);
x+=h/2;
f+=2*F(x);
}
f*=h/6;
cout <<"f(x)在[ "<<a<<" , "<<b<<" ]区间内的积分结果为："<<f<< endl;
return 0;
}
评论0