Simpson复合求积分.doc

数值分析实验报告六

一、 实验目的
理解Simpson复合求积分，掌握Simpson复合求积分算法设计

二、 实验内容
使用vc++编程，实现该方法，即Simpson复合求积分

三、 实验步骤
程序如下：
#include <math.h>
#include <math.h>
double f(double);
double mp(double a,double b,double eps)

{
int n,k;
double h,t1,t2,s1,s2,ep,p,x;
n=1;h=b-a;
t1=h*(f(a)+f(b))/2;
s1=t1;
ep=eps+1;

while(ep>=eps)
{
p=0;
for(k=0;k<=n-1;k++)
{
x=a+(k+0.5)*h;
p=p+f(x);
}
t2=(t1+h*p)/2;
s2=(4*t2-t1)/3;
ep=fabs(s2-s1);
t1=t2;s1=s2;n=n+n;h=h/2;
}
return(s2);
}

#include <iostream.h>
//#include <f.c>
void main()
{
double a,b,eps,t;
a=0;b=1;eps=0.000001;
t=mp(a,b,eps);
cout<<"t="<<t<<endl;
}
double f(double x)

{
double y;
y=log(1+x)/(1+x*x);
return(y);
}

四、 运行结果