用辛普森公式求重积分在数值计算中，计算结果较为精确

二分法、牛顿迭代法、复合梯形公式、复合辛普森公式、改进欧拉公式、四阶龙格库塔公式matlab代码合集，带有一份数据分析word文档

C语言数值分析上机实践作业：用复合辛普森公式计算积分1/x*x+1

本程序使用java实现复合辛普森求积公式，使用者可自行修改求积函数

复化梯形公式和复化辛普森公式 以及matlab源代码

实现复合辛普森公式的简单C语言程序 #include #include #define f(x) (x)/(4+(x)*(x)) main() {int k,n; float a,b,h,s,x; scanf("%f,%f,%d",&a,&b,&n); h=(b-a)/(2*n); x=a; s=f(x)-f(b); k=1; do {x=x+h;s=s+4*f(x); x=x+h;s=s+2*f(x); k++; }while(k<=n); s=(s*h)/3; printf("s=%f\n",s); }

修正了yyhapy上传的算法，加入了一个控制精度e，使得积分精度可以控制，使得该算法更加实用。在此感谢yyhapy的源码

显示积分范围内的函数以及用于近似其下方区域的抛物线。 函数接受 f 是一个变量中的函数I 是一个 1x2 向量m 是用于创建 m 条抛物线的点数。 并输出 f 在 I 上的图形和将用于近似 f 下面积的 m 条抛物线。

% x = simpsonQuadrature(fun,tLow,tUpp,n) % % 此函数在 n 个均匀段中的每一个上使用辛普森正交来% 在区间 tLow <= t <= tUpp 上近似 fun(t) 的积分% % 输入： % fun = 函数句柄% f = 乐趣(t) % t = [1, nt] = [a,b] 上的时间查询点% f = [nx, nt] = 查询点的向量函数% tLow = 时间标量下限% tUpp = 时间标量上限% n = 间隔划分的段数% % 输出： % x = [nx, 1] = 沿每个维度的积分%

用C语言辛普森算法算积分 最终得到运行程序