用c语言编程计算积分

#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))}

Matlab计算定积分的具体步骤如下：

1、以f(x)=e^2x+sin(x+π/3)  ，积分下限：a=0，积分上限：b=π/4  为例。首先，建立被积函数M文件；点击New,选择Function；

2、在编辑器中输入指令：function f=f(x)；f=exp(2*x).*sin(x+pi/3)（注意，用  .* )

3、保存函数文件到自己的工作路径，比如G:\matlab\work。

输入指令：cd G:\matlab\work，是该路径成为当前路径。

4、matlab中计算积分的两种指令：

1.F=quad('fname',a,b,tol,trace)        Simpson数值积分法

2.F=quad8('fname',a,b,tol,trace)     Newton-Cotes数值积分法

其中： fname是被积函数表达式或函数名，a,b分别是上下限，tol可以控制积分精度，省略则取0.001；trace=1则用图形表示积分过程，trace=0，没有图形。

两者 比较，quad8精度更高。

5、调用积分函数squad进行计算。输入指令： F=quad('f(x)',0,pi/4)，如图，回车后可得到计算结果1.8612。其他函数，只要修改函数文件中的表达式即可。

---