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开启多个子进程并行运行的方法
Apr 18 Python
python绘制简单折线图代码示例
Dec 19 Python
python读取文件名称生成list的方法
Apr 27 Python
Python设置在shell脚本中自动补全功能的方法
Jun 25 Python
python 用下标截取字符串的实例
Dec 25 Python
对Python 多线程统计所有csv文件的行数方法详解
Feb 12 Python
python绘制地震散点图
Jun 18 Python
Python函数中参数是传递值还是引用详解
Jul 02 Python
python实现拼接图片
Mar 23 Python
Python Django搭建网站流程图解
Jun 13 Python
Python 开发工具通过 agent 代理使用的方法
Sep 27 Python
Python必备技巧之字符数据操作详解
Mar 23 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
Codeigniter框架实现获取分页数据和总条数的方法
2014/12/05 PHP
Laravel搭建后台登录系统步骤详解
2016/07/26 PHP
PHP封装返回Ajax字符串和JSON数组的方法
2017/02/17 PHP
PHP基于双向链表与排序操作实现的会员排名功能示例
2017/12/26 PHP
设为首页和收藏的Javascript代码(亲测兼容IE,Firefox,chrome等浏览器)
2013/11/18 Javascript
javascript中this指向详解
2016/04/23 Javascript
实例讲解使用原生JavaScript处理AJAX请求的方法
2016/05/10 Javascript
微信小程序 教程之模板
2016/10/18 Javascript
jQuery快速实现商品数量加减的方法
2017/02/06 Javascript
微信小程序本地缓存数据增删改查实例详解
2017/05/24 Javascript
基于AngularJS的拖拽文件上传的实例代码
2017/07/15 Javascript
jquery-file-upload 文件上传带进度条效果
2017/11/21 jQuery
Vue Router的懒加载路径的解决方法
2018/06/21 Javascript
JS实现匀速与减速缓慢运动的动画效果封装示例
2018/08/27 Javascript
详解基于React.js和Node.js的SSR实现方案
2019/03/21 Javascript
express框架中使用jwt实现验证的方法
2019/08/25 Javascript
jQuery实现轮播图效果
2019/11/26 jQuery
浅谈es6中的元编程
2020/12/01 Javascript
[01:16:01]VGJ.S vs Mski Supermajor小组赛C组 BO3 第一场 6.3
2018/06/04 DOTA
python的几种开发工具介绍
2007/03/07 Python
python实现sublime3的less编译插件示例
2014/04/27 Python
详解Python使用tensorflow入门指南
2018/02/09 Python
python网络应用开发知识点浅析
2019/05/28 Python
Python中 Global和Nonlocal的用法详解
2020/01/20 Python
python 回溯法模板详解
2020/02/26 Python
python实现井字棋小游戏
2020/03/04 Python
python编写实现抽奖器
2020/09/10 Python
完美解决torch.cuda.is_available()一直返回False的玄学方法
2021/02/06 Python
CSS3实现3D翻书效果
2016/06/20 HTML / CSS
香港化妆品经销商:我的公主
2016/08/05 全球购物
法国购买隐形眼镜和眼镜网站:Optical Center
2019/10/08 全球购物
半年思想汇报
2013/12/30 职场文书
经典导游欢迎词大全
2014/01/16 职场文书
公司联欢会策划方案
2014/05/19 职场文书
幼儿园毕业典礼园长致辞
2015/07/29 职场文书
2016年学校招生广告语
2016/01/28 职场文书