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程序员
Jun 12 Python
Python中动态检测编码chardet的使用教程
Jul 06 Python
python实现在IDLE中输入多行的方法
Apr 19 Python
python使用PIL给图片添加文字生成海报示例
Aug 17 Python
python实现对指定字符串补足固定长度倍数截断输出的方法
Nov 15 Python
Python提取频域特征知识点浅析
Mar 04 Python
python如何制作缩略图
Apr 30 Python
pyinstaller 3.6版本通过pip安装失败的解决办法(推荐)
Jan 18 Python
jupyter notebook 实现matplotlib图动态刷新
Apr 22 Python
python求numpy中array按列非零元素的平均值案例
Jun 08 Python
python中pyplot基础图标函数整理
Nov 10 Python
filter使用python3代码进行迭代元素的实例详解
Dec 03 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
【星际争霸1】人族1v7家ZBath
2020/03/04 星际争霸
php中\r \r\n \t的区别示例介绍
2014/02/08 PHP
php去除字符串换行符示例分享
2014/02/13 PHP
laravel 5 实现模板主题功能(续)
2015/03/02 PHP
php动态函数调用方法
2015/05/21 PHP
Zend Framework前端控制器用法示例
2016/12/11 PHP
记录几个javascript有关的小细节
2007/04/02 Javascript
javascript之函数直接量(function(){})()
2007/06/29 Javascript
javascript轻松实现当鼠标移开时已弹出子菜单自动消失
2013/12/29 Javascript
BOOTSTRAP时间控件显示在模态框下面的bug修复
2015/02/05 Javascript
Jquery使用val方法读写value值
2015/05/18 Javascript
JavaScript必知必会(九)function 说起 闭包问题
2016/06/08 Javascript
简单三步实现报表页面集成天气
2016/12/15 Javascript
canvas学习之API整理笔记(一)
2016/12/29 Javascript
jq.ajax+php+mysql实现关键字模糊查询(示例讲解)
2018/01/02 Javascript
jQuery内容过滤选择器与子元素过滤选择器用法实例分析
2019/02/20 jQuery
详解vue中axios的使用与封装
2019/03/20 Javascript
jQuery实现鼠标放置名字上显示详细内容气泡提示框效果的方法分析
2020/04/04 jQuery
vue中组件通信详解(父子组件, 爷孙组件, 兄弟组件)
2020/07/27 Javascript
vue 动态给每个页面添加title、关键词和描述的方法
2020/08/28 Javascript
python 读写、创建 文件的方法(必看)
2016/09/12 Python
Python实现将蓝底照片转化为白底照片功能完整实例
2019/12/13 Python
is_file和file_exists效率比较
2021/03/14 PHP
中外合拍动画首获奥斯卡提名,“上海出品”《飞奔去月球》能否拿下最终大奖?
2021/03/16 国漫
英国评分最高的女性剃须刀订阅盒:FFS Beauty
2018/01/25 全球购物
Made in Design德国:设计师家具、灯具和装饰
2019/10/31 全球购物
巴西葡萄酒商店:Divvino
2020/02/22 全球购物
大专应届生个人简历的自我评价
2013/10/15 职场文书
《蝙蝠和雷达》教学反思
2014/04/23 职场文书
节水标语大全
2014/06/11 职场文书
2014年最新学校运动会广播稿
2014/09/17 职场文书
2015年建党94周年演讲稿
2015/03/19 职场文书
不同意离婚答辩状
2015/05/22 职场文书
志愿者服务宣传标语口号
2015/12/26 职场文书
Redis Stream类型的使用详解
2021/11/11 Redis
Mysql的Table doesn't exist问题及解决
2022/12/24 MySQL