vc++6.0里怎么求定积分?

来源:百度知道 编辑:UC知道 时间:2024/06/27 19:08:34
请问vc++6.0里如何实现输入积分上限和积分下限之后求出给定函数定积分?

可以用辛普森积分公式,示例程序:

//---------------------------------------------------------------------------

#include <stdio.h>
#include <math.h>
#define PI 3.14159

double f(double x)
{
double y;
y=log(x);/*要求积分的函数*/
return y;
}
double simp(double a,double b,const double eps,double (*f)(double))/*积分计算函数,a和b是上下限,eps是精度*/
{
int n=1,k;
double h,t1,t2,s1,s2,ep,p,x;
if (a-b>0) {
x=a;
a=b;
b=x;
}
h=b-a;
t1=h*((*f)(a)+(*f)(b))/2.0;
s1=t1;
ep=eps+1.0;
while (ep>=eps)
{
p=0.0;
for (k=0; k<n; k++) {
x=a+(k+0.5)*h;
p+=(*f)(x);
}
t2=(t1+h*p)/2.0;
s2=(4.0*t2-t1)/3.0;
ep=fabs(s2-s1);
t1=t2;
s1=s2