1.cpp

用C++写的辛普生算法

#include<math.h>
#include<iostream.h>
#define f(x) exp(x)
#define n 100
double simpson(double a,double b)
{
	double t=0.0,h,tn,j;
	h=(b-a)/n;
	for(int i=1;i<=n-1;i++)
	{
	j=a+h*i;
	t+=f(j);
	}
    t*=2;
	for( i=0;i<=n-1;i++)
	{
	j=a+h*(i+1/2);
    t+=4*f(j);
	}
	t+=f(a)+f(b);
	tn=h*t/6;
	return tn;
}
 void main()
 {
	double M,N;
    cout<<"Enter 下界M,上界N"<<endl;
    cin>>M>>N;
    double Sn=simpson(M,N);
    cout<<"Sn="<<Sn<<endl;
 }