Question

C语言，用梯形法编程求定积分x^3+x/2+1的值

用梯形法编程求定积分x^3+x/2+1从0积到10的值，要求使用指针
h=（b-a）/2,函数f在（a,b）区间的定积分的公式为s=h/2(f(a)+f(b))+h(i=1,n-1)f(xi)求和，取n为1000,a为0,b为10

Answer 1

这个其实很简单，给你做了一个很简单的通用的积分函数，
已经做过优化，自己一调用就可以了，效率很高。

#include <stdio.h>

double f(double x) //返回函数的值
{ double y;
y=x*x+x*2+1.0;
return y;
}

double Integrate(double (*f)(double),double a,double b,int n) //计算函数的定积分
{ double dx,y,sy,x;
if(n<=1)return 0.0;
dx=(b-a)/n;
x=a;
sy=0.0;
for(int i=0;i<=n;i++,x+=dx)
sy+=f(x);
sy-=f(a)/2.0;
sy-=f(b)/2.0;
return sy*dx;
}

void main()
{ double s;
s=Integrate(f, 1, 3, 100000);
printf("s=%lf\n",s);
}