1.cpp

用simpson公式计算三重积分（采用C++）

#include<iostream>
#include<iomanip>
#include<string>
using namespace std;
#define pi 3.14159
 int m,n,p; 
double a(double theta)//theta的下标
 { double y;
   y=0;
 return y;
 }
 double b(double theta) //theta的上标
 {double  y;
y=pi;
 return  y;
 }
 double c(double fai) // 的下标
 {double z;
z=0;
 return z;
 }
 double d(double fai) // 的上标
 {
	  double z;
	 z=2*pi;
	 return z;
 }
 double f(double r,double theta,double fai) //  
 {
	 double   I;
	 I=r*sin(theta)*sin(theta)*cos(fai);
	 return I;
 }
 double g(double r,double theta) //用Simpson公式计算函数 的值
{double fai;
	double subscript; //subscript为下标，下同
subscript=c(fai);
double supscript;//subscript为上标，下同
supscript=d(fai); 
double sum,step;
step=(supscript-subscript)/(2*p); // step为步长，下同
sum=f(r,theta,subscript)+f(r,theta,supscript);
for(int i=1;i<=(2*p-1);i++)
{ double z;
z=subscript+i*step;
if(i%2)
	sum+=4*f(r,theta, fai);
else
	sum+=2*f(r,theta, fai);
}
return sum*step/3;
}
 
double h(double r) //用simpson公式计算函数 的值
{ 
	double subscript;
subscript=a(r);
double  supscript;
supscript=b(r);
double sum,step;
step=(supscript-subscript)/(2*n);
sum=g(r,subscript)+g(r,supscript);
for(int j=1;j<=(2*n-1);j++)
{double theta;
theta=subscript+j*step;
if(j%2)
	sum+=4*g(r,theta);
else
	sum+=2*g(r,theta);
}
return sum*step/3;
} 
double Simpson(double subscript, double supscript) //用Simpson公式计算函数 
{
	double sum,step;
step=(supscript-subscript)/(2*m);
sum=h(subscript)+h(supscript);
for(int k=1;k<=(2*m-1);k++)
{
	double r;
r=subscript+k*step;
if(k%2)
	sum+=4*h(r);
else
	sum+=2*h(r);
}
return sum*step/3;
}

int main()
 {int a,b;
	 cout<<"请输入区间a，b和等分数m，n，p："<<endl;
 cin>>a>>b>>m>>n>>p;
 double I =Simpson(a,b);
 cout<<"请输出三重积分值:"<<I <<endl;
 return 0;
 }