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

}

题主给出的积分函数，用matlab 求解，可以这样做

1、确定x在【-300，300】中，取若干个x(i)数据，如

x=-300:10:300

2、使用for循环语句，将x(i)代入积分函数I(x)中，使用int函数计算其积分值，同时将值赋值给Ix数组变量中

for n=1:600/t+1

syms x

A_1=int(cos((pi*(x^2))/2),x,a1(1,n),a2(1,n))

A_2=int(sin((pi*(x^2))/2),x,a1(1,n),a2(1,n))

Ix(n)=1000*(lambda/(2*L))*((A_1)^2+(A_2)^2)

end

3、最后用plot函数，绘出其x—I(x)的图形

运行程序，可以得到如下结果。

---