用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<stdio.h>

#include<局卖math.h>//pow()指数函数头文件

#define e 2.71828//定义自然常数e

int main ()

{

  double a,b,n1,n2,n3,n4,n5,t

  scanf("%lf%lf",&a,&b)

  if(a<b)//保证上限比下限大；

  {

  t=a

  a=b

b=t

}

  n1=a-b+0.5*(a*a-b*b)//a和b是未知数，带袭积分运算后，原式=a和b混合运算的代数式；

  n2=a*a-b*b+3*a-3*b

  n3=e*0.5*(a*a-b*b)+a-b

  n4=(pow(a+1,3)-pow(b+1,3))/3

  n5=(pow(a,4)-pow(b,4))/4

  printf("%.2lf %.2lf %.2lf %.2lf %.2lf",n1,n2,n3,n4,n5)

  return 0

}

图1是运行结果，图2是原题的积分结果；通过数学的积分蠢腊兄结果可以看出，积分符号的不在了；只剩下上下限a，b组成的代数式；

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

以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)

}

---