/* 题目:编写一个求定积分的通用函数,用它分别求以下5个函数的定积分:
∫(1+x)dx,∫(2x+3)xdx,∫(exp(x)+1)dx(exp表示e为底x为指数),∫(1+x)(1+x)dx,∫x*x*xdx .上-下限为a和b.
*/
程序清单如下:
#include "stdio.h"
#include "conio.h"
#include "math.h"
float f1(float x)
{return (1+x);}
float f2(float x)
{return (3+2*x);}
float f3(float x)
{return (1+exp(x));}
float f4(float x)
{return ((1+x)*(1+x));}
float f5(float x)
{return (x*x*x);}
float Integral(float a,float b,float (*fun)(float))
{
float fsum=0.0;
long N=400000L,i;
for(i=0;i<=N;i++)
fsum+=(*fun)(a+i*(b-a)/N);
return (fsum*(b-a)/N);
}
void main()
{
printf("f1=1+x 时∫f1dx=%.3lf\n",Integral(2,3.5,f1));
printf("f2=3+x*x 时∫f2dx=%.3lf\n",Integral(-6,1,f3));
printf("f4=(1+x)*(1+x) 时∫f4dx=%.3lf\n",Integral(3,4,f4));
printf("f5=x*x*x 时∫f5dx=%.3lf\n",Integral(1,2,f5));
getch();
}
运行结果如下:
*************************************
f1=1+x 时∫f1dx=5.625
f2=3+x*x 时∫f2dx=9.716
f4=(1+x)*(1+x) 时∫f4dx=20.333
f5=x*x*x 时∫f5dx=3.750
*************************************