组合辛普森公式(数值积分)
问题描述
组合辛普森公式求函数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. 复化梯形公式 %复化梯形公式 function t=agui_trapz(fname,d2fname,a,b,e) %fname ...
- 【数值分析实验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 ...
- 数值积分-求积公式余项,牛顿-柯特斯公式,辛普森公式,复合梯形公式,复合辛普森公式
文章目录 1.求积公式余项 1.1 定义 1.2 Python实现求积公式余项 2.牛顿-柯特斯公式 2.1 定义 2.2 Python实现牛顿-柯特斯公式 3.复合梯形公式 3.1 定义 3.2 P ...
- matlab用辛普森公式求积分_标准正态分布概率密度函数的定积分计算方法及Python实现代码...
最近利用碎片时间在读Allen B.Downey的<贝叶斯思维:统计建模的Python学习法>,顺便用手机上的Pythonista写实例.因为Pythonista没有scipy科学计算包, ...
- 数值分析c语言编程辛普森公式,数值计算方法:矩形法、梯形法与辛普森公式...
在数值分析中,数值积分是计算定积分数值的方法和理论.在数学分析中,给定函数的定积分的计算不总是可行的.许多定积分不能用已知的积分公式得到精确值.数值积分是利用黎曼积分和积分中值等数学定义和定理,用数值 ...
- Python实现梯形公式 、辛普森公式、复化梯形公式、复化辛普森公式
数值分析:梯形公式 .辛普森公式.复化梯形公式.复化辛普森公式 Python实现梯形公式 .辛普森公式.复化梯形公式.复化辛普森公式 数值求积公式概念 梯形公式与辛普森公式 梯形公式与辛普森公式的余项 ...
- 自适应复化辛普森公式求积算法(C语言实现)
自适应复化辛普森公式求积算法(C语言实现) 利用复化辛普森公式求积分自适应步骤 基于C语言实现的代码 利用复化辛普森公式求积分自适应步骤 h为步长,a为积分下限,b为积分上限,f为积分函数,n为划分的 ...
- 一张图快速get浅层辛普森公式
一张图快速get浅层辛普森公式
- 2017杭电ACM集训队单人排位赛 - 2 -1002 地狱飞龙 (辛普森公式求积分)(模板)
题干: 最近clover迷上了皇室战争,他抽到了一种地狱飞龙,很开心.假设地域飞龙会对距离为d的敌人每秒造成k/d2k/d^2伤害.假设地域飞龙位于坐标轴原点,以每秒v1的速度向y轴正方向移动,敌人在 ...
最新文章
- 为了更高效的开发代码,这里列出了一些webstorm的快捷键和zencoding
- 【原创】ucos信号量的操作及原理
- 值得一用的Windows软件
- Dapper.Common基于Dapper的开源LINQ超轻量扩展
- springmvc log4j配置
- 【转】unity地形插件T4M使用帮助
- Java高并发编程详解系列-不可变对象设计模式
- ES6_解构赋值_note
- IReport无法打开的原因
- [面试] C/C++ 语法(六)—— RTTI(运行时类型信息)
- office2007安装时,提示找不到Office.zh-cn下的OfficeMUI.msi解决方法
- dev c 编程语言,devc如何自定义头文件并使用
- 大华网络摄像头ip搜索工具_【技术篇】NVR4.0接第三方摄像头,安排!
- HDU_4585_Shaolin
- python 监听键盘输入并收集数据进行分析
- shiro介绍(一)
- 小白如何在阿里云快速建站(pc站+手机站+公众号+小程序)
- 2.5元组tuple
- 江理工单片机实训(自动窗帘模拟控制系统)
- 【漏洞通告】Gitlab存在多个高危漏洞
热门文章
- html网页底部弹窗,【HTML】底部弹窗插件代码
- 避免全局变量污染的一些措施!
- 你使用计算机来做过什么,我来交你使用软做计算机日常维护!.doc
- react 自定义日历 手把手教你
- Springboot 报错 template might not exist or might not be accessible by any of the configured Template
- De Bruijn 序列生成的本原多项式方法
- linux 键盘记录,Linux键盘记录sh2log keylogger notes
- 验证码(正则表达式)
- 基于粒子群优化的灰色预测方法(Matlab)——附代码
- android 最简洁的日历,简洁日历 - 日历与 ToDo 的完美结合 #Android