当前位置: 代码迷 >> C语言 >> 定积分的通用函数
  详细解决方案

定积分的通用函数

热度:151   发布时间:2007-04-22 15:07:06.0
定积分的通用函数

/* 题目:编写一个求定积分的通用函数,用它分别求以下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
*************************************

搜索更多相关的解决方案: 函数  积分  

----------------解决方案--------------------------------------------------------
  相关解决方案