问题描述

组合辛普森公式求函数f(x)=2+sin(2*sqrt(x))的积分近似值。

输入形式

在屏幕上依次输入积分上限、下限和等距子区间个数。

输出形式

输出使用组合辛普森公式求得的积分近似值。

样例输入

1 6 5

样例输出

8.18301549

样例说明

输入:积分上限a为1、下限b为6和等距子区间个数m为5。输出:积分近似值(保留小数点后8位有效数字)

代码

# 组合辛普森公式
import mathdef Input():a, b, m = input().split()a = float(a)b = float(b)m = int(m)return a, b, mdef f(x):return 2 + math.sin(2 * math.sqrt(x))def simprl(a, b, m):h = (b - a) / (2 * m)s1 = 0s2 = 0for k in range(1, m + 1):x = a + h * (2 * k - 1)s1 += f(x)for k in range(1, m):x = a + h * 2 * ks2 += f(x)s = h * (f(a) + f(b) + 4 * s1 + 2 * s2) / 3return sdef out(x):print(x)def main():a, b, m = Input()ans = simprl(a, b, m)out('%.8f' % ans)if __name__ == '__main__':main()

组合辛普森公式(数值积分)相关推荐

  1. 【计算方法数值分析】复化梯形公式、复化辛普森公式和龙贝格数值积分

    [计算方法数值分析]复化梯形公式.复化辛普森公式和龙贝格数值积分 1. 复化梯形公式 %复化梯形公式 function t=agui_trapz(fname,d2fname,a,b,e) %fname ...

  2. 【数值分析实验MATLAB】数值积分:梯形公式、辛普森公式、复化梯形公式、复化辛普森公式、龙贝格算法、自适应求积方法

    (1)Matlab代码: f1: if x==0y=1; else%y=exp(-x^2);y=(sin(x))/x; end 梯形公式: a=0; b=1; tx=(b-a)/2*(f1(a)+f1 ...

  3. 数值积分-求积公式余项,牛顿-柯特斯公式,辛普森公式,复合梯形公式,复合辛普森公式

    文章目录 1.求积公式余项 1.1 定义 1.2 Python实现求积公式余项 2.牛顿-柯特斯公式 2.1 定义 2.2 Python实现牛顿-柯特斯公式 3.复合梯形公式 3.1 定义 3.2 P ...

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

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

  5. 数值分析c语言编程辛普森公式,数值计算方法:矩形法、梯形法与辛普森公式...

    在数值分析中,数值积分是计算定积分数值的方法和理论.在数学分析中,给定函数的定积分的计算不总是可行的.许多定积分不能用已知的积分公式得到精确值.数值积分是利用黎曼积分和积分中值等数学定义和定理,用数值 ...

  6. Python实现梯形公式 、辛普森公式、复化梯形公式、复化辛普森公式

    数值分析:梯形公式 .辛普森公式.复化梯形公式.复化辛普森公式 Python实现梯形公式 .辛普森公式.复化梯形公式.复化辛普森公式 数值求积公式概念 梯形公式与辛普森公式 梯形公式与辛普森公式的余项 ...

  7. 自适应复化辛普森公式求积算法(C语言实现)

    自适应复化辛普森公式求积算法(C语言实现) 利用复化辛普森公式求积分自适应步骤 基于C语言实现的代码 利用复化辛普森公式求积分自适应步骤 h为步长,a为积分下限,b为积分上限,f为积分函数,n为划分的 ...

  8. 一张图快速get浅层辛普森公式

    一张图快速get浅层辛普森公式

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

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

最新文章

  1. 为了更高效的开发代码,这里列出了一些webstorm的快捷键和zencoding
  2. 【原创】ucos信号量的操作及原理
  3. 值得一用的Windows软件
  4. Dapper.Common基于Dapper的开源LINQ超轻量扩展
  5. springmvc log4j配置
  6. 【转】unity地形插件T4M使用帮助
  7. Java高并发编程详解系列-不可变对象设计模式
  8. ES6_解构赋值_note
  9. IReport无法打开的原因
  10. [面试] C/C++ 语法(六)—— RTTI(运行时类型信息)
  11. office2007安装时,提示找不到Office.zh-cn下的OfficeMUI.msi解决方法
  12. dev c 编程语言,devc如何自定义头文件并使用
  13. 大华网络摄像头ip搜索工具_【技术篇】NVR4.0接第三方摄像头,安排!
  14. HDU_4585_Shaolin
  15. python 监听键盘输入并收集数据进行分析
  16. shiro介绍(一)
  17. 小白如何在阿里云快速建站(pc站+手机站+公众号+小程序)
  18. 2.5元组tuple
  19. 江理工单片机实训(自动窗帘模拟控制系统)
  20. 【漏洞通告】Gitlab存在多个高危漏洞

热门文章

  1. html网页底部弹窗,【HTML】底部弹窗插件代码
  2. 避免全局变量污染的一些措施!
  3. 你使用计算机来做过什么,我来交你使用软做计算机日常维护!.doc
  4. react 自定义日历 手把手教你
  5. Springboot 报错 template might not exist or might not be accessible by any of the configured Template
  6. De Bruijn 序列生成的本原多项式方法
  7. linux 键盘记录,Linux键盘记录sh2log keylogger notes
  8. 验证码(正则表达式)
  9. 基于粒子群优化的灰色预测方法(Matlab)——附代码
  10. android 最简洁的日历,简洁日历 - 日历与 ToDo 的完美结合 #Android