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

#include<stdio.h>

#include<math.h>

float

f1(float

x)

{

return(1.0+x)

}

float

f2(float

x)

{

return(2.0*x+3.0)

}

float

f3(float

x)

{

return(exp(x)+1)

}

float

f4(float

x)

{

return(pow(1+x,2))

}

float

f5(float

x)

{

return(pow(x,3))

}

float

fsimp(float

a,float

b,float

(*p)(float))

{

float

c,s

c=(a+b)/2

s=(b-a)/6*(p(a)+4*p(c)+p(b))

return

s

}

int

main()

{

float

a,b

printf("请输告橡入积分下限a的值："缓友毁)

scanf("%f",&a)

printf("请输入积分上限b的值：")

scanf("%f"扰备,&b)

printf("%f\n",fsimp(a,b,f1))

printf("%f\n",fsimp(a,b,f2))

printf("%f\n",fsimp(a,b,f3))

printf("%f\n",fsimp(a,b,f4))

printf("%f\n",fsimp(a,b,f5))

}

#include #include float f1(float x) { return(1.0+x)} float f2(float x) { return(2.0*x+3.0)} float f3(float x) { return(exp(x)+1)} float f4(float x) { return(pow(1+x,2))} float f5(float x) { return(pow(x,3))} float fsimp(float a,float b,float (*p)(float)) { float c,sc=(a+b)/2s=(b-a)/6*(p(a)+4*p(c)+p(b))return s} int main() { float a,bprintf("请输入积分下限a的腔并值：")scanf("%f",&a)printf("请答虚输入积分上限b的值清圆燃：")scanf("%f",&b)printf("%f\n",fsimp(a,b,f1))printf("%f\n",fsimp(a,b,f2))printf("%f\n",fsimp(a,b,f3))printf("%f\n",fsimp(a,b,f4))printf("%f\n",fsimp(a,b,f5))}

基本是这样的，用梯形发求定积分，对应于一个积分式就要有一段程序，不过你可以改变程序的一哗核顷小部分来改变你所要求的积分式。

以c为例：求f(x)=xsinx从1到2的积分

#include <math.h>

float integral(float(*fun)(float x),float a,float b,int,n)

{float s,h,y

int i

s=(fun(a)+fun(b))/2

h=(b-a)/n/*积分步长*/

for(i=1i<ni++)

s=s+fun(a+i*h)

y=s*h

return y/*返回积分值*/

}

float f(float x)

{return(x*sinx) /*修改此处可以改变被积函数*/

}

main()

{float y

y=integral(f,1.0,2.0,150)/*修改此氏圆处可以改变积分上下限和步长*/

printf("y=%f\乱陆n",y)

}

---