对于积分区间[a, b],若

则成立

辛普森公式

辛普森公式可看作是改良的梯形公式。梯形公式是以直线逼近实际曲线,而辛普森公式则以二次曲线(即抛物线)逼近。

以二次曲线逼近实际曲线

根据辛普森公式可得各积分小区间内的积分值如下式所示:

其中

于是,全积分区间[a, b]的积分值为

为便于计算机计算,上式可写成

其中

C语言程序实现辛普森公式计算积分

  • 计算积分
  • 代码

//simpson.c

#include

#include

/*函数定义*/

double f(double x)

{

return (1.0-x)*exp(-1*x);

}

/*辛普森公式计算积分*/

double simpson(

double a, double b, int n)

{

double h = (b-a)/n;

double s0,s1=0.0,s2=0.0;

double S;

double x;

int i;

s0 = f(a) + f(b);

for( i = 1; i < n; i += 2 )

{

x = a + h * i;

s1 += f(x);

}

for( i = 2; i < n; i += 2 )

{

x = a + h * i;

s2 += f(x);

}

S = h*(s0+4*s1+2*s2)/3;

return S;

}

int main(void)

{

//积分区间[a, b]

double a = -1.0, b = 1.0;

//区间分割数

int N[5]={10,20,30,50,100};

int i = 0;

for( i = 0; i < 5; i++ )

printf("[%3d等分]积分近似值:

I = %.7lf

matlab用辛普森公式求积分_积分近似计算之辛普森公式相关推荐

  1. matlab用辛普森公式求积分_数值计算实验9 数值积分实验

    实验9 数值积分实验 成绩 实验类型:●验证性实验  ○综合性实验  ○设计性实验 实验目的:进一步熟练掌握变步长数值积分算法,提高编程能力和解决定积分问题的实践技能. 实验内容:用龙贝格积分算法计算 ...

  2. matlab用辛普森公式求积分_数值积分常用方法

    数值积分的基本思想 由积分中值定理可知,在积分区间 内存在一点 ,成立 式的几何意义即为:底为 而高为 的矩形的面积恰等于所求曲边梯形的面积 .因此,要想求出 式左端积分,我们只需要知道三个值: 即可 ...

  3. matlab用辛普森公式求积分_如何用Excel公式求最大值对应的行列序号

    微信公众号: Excel and Python 微信名搜索: 实用办公编程技能 如何用Excel公式求最大值对应的行列序号呢? 下面,我们来看看来自问题互动栏目的一个具体问题. 具体问题:求出哪一天哪 ...

  4. matlab用辛普森公式求积分_变限积分函数求导以及高阶导数求法的一些总结

    感谢 @聚创考研 的张帆老师,给我上了一堂生动的课.特此总结一下课上求导数的方法(怕自己忘了). 1.变限积分函数求导 变限积分函数求导简单的分为三类: 第一类(或者形如 这种)可以直接得到 ,第二. ...

  5. matlab用辛普森公式求积分_标准正态分布概率密度函数的定积分计算方法及Python实现代码...

    最近利用碎片时间在读Allen B.Downey的<贝叶斯思维:统计建模的Python学习法>,顺便用手机上的Pythonista写实例.因为Pythonista没有scipy科学计算包, ...

  6. matlab中几种求积分的方法

    matlab中几种求积分的方法 举例图形 Midpoint Rule Trapezoid Rule 13 Simpsons 利用integral 结果显示 matlab中几种求积分的方法 最近看了几天 ...

  7. 2017杭电ACM集训队单人排位赛 - 2 -1002 地狱飞龙 (辛普森公式求积分)(模板)

    题干: 最近clover迷上了皇室战争,他抽到了一种地狱飞龙,很开心.假设地域飞龙会对距离为d的敌人每秒造成k/d2k/d^2伤害.假设地域飞龙位于坐标轴原点,以每秒v1的速度向y轴正方向移动,敌人在 ...

  8. 【matlab初学】各种求积分的方法和分段函数求积分

    提示:文章写完后,目录可以自动生成,如何生成可参考右边的帮助文档 一.函数 1.int 符号积分 int(s) 符号表达式s的不定积分 int(s,x) 符号表达式s关于变量x的不定积分 int(s, ...

  9. matlab图像导数求积分_第二讲matlab求微分方程导数积分

    第二讲matlab求微分方程导数积分 第二讲 导数与微分方程 一.实验内容 1.实际引例 (牛顿冷却模型)警察上午9点钟发现一被谋杀者,并测得尸体温 度为32.4℃,一小时以后,尸体的温度变为31.7 ...

最新文章

  1. CSS3重新定义input中呆若木鸡的默认复选框CheckBox和单选框Radio样式
  2. 【系列】EOS开发1 开发环境搭建
  3. 排序算法二:归并排序(Merge sort)
  4. mybatis反射的核心类
  5. 【解题报告】Leecode 237. 删除链表中的节点——Leecode每日一题系列
  6. ZXing二维码扫描添加个人总结
  7. MSP430F5529 DriverLib 库函数学习笔记(十三)认识低功耗模式
  8. 只考虑用户估计的计算机时间,操作系统第四章进程调度和死锁习题及答案
  9. Litjson快速入门
  10. 使用MySql的yearweek()函数和通过java的Calendar获取的时间是当前年份第几周的数值不一致的问题
  11. 港科夜闻|香港科技大学(HKUST)及香港科大智能建造实验室(HKUST BIM Lab)荣获CIC建造数码化大奖组织类最高奖项...
  12. 使用电信光猫加路由器实现内网穿透,外网访问内网
  13. 静态分析android程序之阅读smali代码
  14. 深度学习之ISAR超分辨率成像
  15. PredRNN++: Towards A Resolution of the Deep-in-Time Dilemma in Spatiotemporal Predictive Learning 翻译
  16. MySQL原理与学习
  17. 7-15 福到了 c语言,7-54 福到了 (15 分)
  18. sparklines_Sparklines入门200年来功能强大的数据可视化形式
  19. 人行征信2.0对接服务:全业务种类数据,精细您的征信业务管理!
  20. backtrack 5 虚拟机 安装

热门文章

  1. 当我用ChatGPT面试项目经理,真有趣
  2. 【很全】各类型的电容是怎么失效的?
  3. RxJava2 Flowable blocking系列
  4. Jmeter做单接口测试-超详细步骤讲解
  5. 使用Python Pandas处理亿级数据
  6. 快速安装create-react-app脚手架
  7. Deformable Convolution Networks
  8. 我的第一篇课程设计(交通咨询系统)
  9. 数据结构与算法——课程设计报告(全国交通咨询模拟)
  10. 软件模块化设计的思路