c语言求定积分的通用函数 c++求积分函数

C程序设计：编制求解定积分的通用函数

#includestdio.h

网站建设哪家好，找创新互联建站！专注于网页设计、网站建设、微信开发、小程序制作、集团企业网站建设等服务项目。为回馈新老客户创新互联还提供了荷塘免费建站欢迎大家使用！

#includemath.h

float f1(float a)

{

float z;

z=a+1;

return(z);

}

float f2(float a)

{

float z;

z=2*a+3;

return(z);

}

void integral(float a,float b,float (*fun)(float ))

{

float i,trace,hight1,hight2,area,sum=0;

trace=(b-a)/1000;

for ( i=a;ib;i+=trace )

{

hight1=(*fun)(a);

a+=trace;

hight2=(*fun)(a);

area=(((hight1+hight2)*trace)/2);

sum+=area;

}

printf("sum=%f\n",sum);

}

void main()

{

float a,b,sum1;

printf("*************************\n");

printf("请输入该定积分的上下区间a和b：\n");

scanf("%f %f",a,b);

integral(a,b,f1(a));

integral(a,b,f2(a));

}

用C语言求定积分

实际问题描述：

求定积分近似值

程序代码如下：

#include

#include

void main()

{

int i,n=1000;

float a,b,h,t1,t2,s1,s2,x;

printf("请输入积分限a,b:");

scanf("%f,%f",a,b);

h=(b-a)/n;

for(s1=0,s2=0,i=1;i=n;i++)

{

x=a+(i-1)*h;

t1=(float)exp(-x*x/2);t2(float)=exp(-(x+h)*(x+h)/2);

s1=s1+t1*h;        /*矩形面积累加*/

s2=s2+(t1+t2)*h/2;        /*梯形面积累加*/

}

printf("矩形法算得积分值：%f.\n",s1);

printf("梯形法算得积分值：%f.\n",s2);

}

程序运行结果如下：

矩形法算得积分值：0.855821

梯形法算得积分值：0.855624

由上面的比较可知，梯形法的精度要高于矩形法。

用C语言编写一个求定积分的程序

这是辛普森积分法。

给你写了fun_1( ），fun_2(），请自己添加另外几个被积函数。

调用方法 t=fsimp(a,b,eps,fun_i);

a,b --上下限，eps -- 迭代精度要求。

#includestdio.h

#includestdlib.h

#include math.h

double fun_1(double x)

{

return 1.0 + x ;

}

double fun_2(double x)

{

return 2.0 * x + 3.0 ;

}

double fsimp(double a,double b,double eps, double (*P)(double))

{

int n,k;

double h,t1,t2,s1,s2,ep,p,x;

n=1; h=b-a;

t1=h*(P(a)+P(b))/2.0;

s1=t1;

ep=eps+1.0;

while (ep=eps)

{

p=0.0;

for (k=0;k=n-1;k++)

{

x=a+(k+0.5)*h;

p=p+P(x);

}

t2=(t1+h*p)/2.0;

s2=(4.0*t2-t1)/3.0;

ep=fabs(s2-s1);

t1=t2; s1=s2; n=n+n; h=h/2.0;

}

return(s2);

}

void main()

{

double a,b,eps,t;

a=0.0; b=3.141592653589793238; eps=0.0000001;

// a definite integral by Simpson Method.

t=fsimp(a,b,eps,fun_1);

printf("%g\n",t);

t=fsimp(a,b,eps,fun_2);

printf("%g\n",t);

// ...

printf("\n Press any key to quit...");

getch();

}