python实现数值积分的Simpson方法实例分析


Posted in Python onJune 05, 2015

本文实例讲述了python实现数值积分的Simpson方法。分享给大家供大家参考。具体如下:

#coding = utf-8
#simpson 法计算积分,数值积分,效果非常理想
from math import *
def func(x): 
 """
 定义被积分函数
 """
 return x*sin(x)
def Get_N(a,b,width):
 # width为步长
 N=int((b-a)/width + 1)
 if N%2 == 0:
  N=N+1
 return N
def GenerateData(a,b,n,width):
 datas = []
 r=a
 for i in range(0,n):
  datas.append(func(r))
  r = r+width
 return datas
def simpson_integral(datas,width,n):
 sum = datas[0]+datas[n-1]
 for i in range(2,n):
  if i%2== 0:
   sum = sum +4*datas[i-1]
  else:
   sum = sum +2*datas[i-1]
 return sum*width/3.0
if __name__ == "__main__":
 a=1.0 #积分上限
 b=3.0 #积分下限
 width=0.0625 #步长
 N=Get_N(a,b,width)
 datas = GenerateData(a,b,N,width)
 print simpson_integral(datas,width,N)

希望本文所述对大家的Python程序设计有所帮助。

Python 相关文章推荐
python和C语言混合编程实例
Jun 04 Python
Python实现抓取页面上链接的简单爬虫分享
Jan 21 Python
python中while循环语句用法简单实例
May 07 Python
Python3写入文件常用方法实例分析
May 22 Python
Python中对数组集进行按行打乱shuffle的方法
Nov 08 Python
浅析python的优势和不足之处
Nov 20 Python
Django集成搜索引擎Elasticserach的方法示例
Jun 04 Python
Python中的 is 和 == 以及字符串驻留机制详解
Jun 28 Python
Python字符串大小写转换拼接删除空白
Sep 19 Python
python中with用法讲解
Feb 07 Python
Anaconda+spyder+pycharm的pytorch配置详解(GPU)
Oct 18 Python
python实现一个简单RPC框架的示例
Oct 28 Python
python采集百度百科的方法
Jun 05 #Python
python中引用与复制用法实例分析
Jun 04 #Python
Python进程间通信用法实例
Jun 04 #Python
python中__slots__用法实例
Jun 04 #Python
python协程用法实例分析
Jun 04 #Python
Python获取系统默认字符编码的方法
Jun 04 #Python
python中__call__内置函数用法实例
Jun 04 #Python
You might like
PHP全概率运算函数(优化版) Webgame开发必备
2011/07/04 PHP
ThinkPHP多表联合查询的常用方法
2020/03/24 PHP
标准版Eclipse搭建PHP环境的详细步骤
2015/11/18 PHP
JavaScript语句可以不以;结尾的烦恼
2007/03/08 Javascript
什么是JavaScript
2009/08/13 Javascript
选择TreeView控件的树状数据节点的JS方法(jquery)
2010/02/06 Javascript
获取服务器传来的数据 用JS去空格的正则表达式
2012/03/26 Javascript
Javascript中Array.prototype.map()详解
2014/10/22 Javascript
angularjs学习笔记之完整的项目结构
2015/09/26 Javascript
jQuery 3.0 的 setter和getter 模式详解
2016/07/11 Javascript
Javascript 获取鼠标当前的位置实现方法
2016/10/27 Javascript
微信小程序 自定义Toast实例代码
2017/06/12 Javascript
vue中的计算属性的使用和vue实例的方法示例
2017/12/04 Javascript
vue使用echarts图表的详细方法
2018/10/22 Javascript
JavaScript实现Tab选项卡切换
2020/02/13 Javascript
python进阶教程之异常处理
2014/08/30 Python
python中os操作文件及文件路径实例汇总
2015/01/15 Python
使用Python编写一个简单的tic-tac-toe游戏的教程
2015/04/16 Python
Python基本语法经典教程
2016/03/11 Python
Python中一般处理中文的几种方法
2019/03/06 Python
python 自动轨迹绘制的实例代码
2019/07/05 Python
Python中变量的输入输出实例代码详解
2019/07/28 Python
Python图片的横坐标汉字实例
2019/12/04 Python
keras自定义损失函数并且模型加载的写法介绍
2020/06/15 Python
Python如何将模块打包并发布
2020/08/30 Python
Space NK美国站:英国高端美妆护肤商城
2017/05/22 全球购物
来自世界上最好大学的在线课程:edX
2018/10/16 全球购物
英国手机零售商:Metrofone
2019/03/18 全球购物
古驰英国官网:GUCCI英国
2020/03/07 全球购物
JMS中Topic和Queue有什么区别
2013/05/15 面试题
高中生期末评语大全
2014/01/28 职场文书
继承公证书
2014/04/09 职场文书
小学班主任经验交流材料
2014/12/16 职场文书
解放思想大讨论活动总结
2015/05/09 职场文书
大学校园招聘会感想
2015/08/10 职场文书
导游词之西安骊山
2019/12/03 职场文书