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编写的com组件发生R6034错误的原因与解决办法
Apr 01 Python
python对html代码进行escape编码的方法
May 04 Python
python框架django基础指南
Sep 08 Python
Python实现在线音乐播放器
Mar 03 Python
python使用fork实现守护进程的方法
Nov 16 Python
浅谈django model postgres的json字段编码问题
Jan 05 Python
python查询文件夹下excel的sheet名代码实例
Apr 02 Python
Python将文字转成语音并读出来的实例详解
Jul 15 Python
Python爬虫实现的根据分类爬取豆瓣电影信息功能示例
Sep 15 Python
浅谈ROC曲线的最佳阈值如何选取
Feb 28 Python
Python try except finally资源回收的实现
Jan 25 Python
Python包argparse模块常用方法
Jun 04 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
Thinkphp中的volist标签用法简介
2014/06/18 PHP
C# WinForm中实现快捷键自定义设置实例
2015/01/23 PHP
win7安装php框架Yii的方法
2016/01/25 PHP
javascript实现列表滚动的方法
2015/07/30 Javascript
Angular2中select用法之设置默认值与事件详解
2017/05/07 Javascript
vue v-on监听事件详解
2017/05/17 Javascript
jQuery实现的简单获取索引功能示例
2018/06/04 jQuery
vue-cli3添加模式配置多环境变量的方法
2019/06/05 Javascript
python实现rest请求api示例
2014/04/22 Python
Python中列表、字典、元组数据结构的简单学习笔记
2016/03/20 Python
python基础while循环及if判断的实例讲解
2017/08/25 Python
Python 实现购物商城,含有用户入口和商家入口的示例
2017/09/15 Python
Python装饰器(decorator)定义与用法详解
2018/02/09 Python
解决python3捕获cx_oracle抛出的异常错误问题
2018/10/18 Python
网易有道2017内推编程题 洗牌(python)
2019/06/19 Python
用Python配平化学方程式的方法
2019/07/20 Python
Python 使用指定的网卡发送HTTP请求的实例
2019/08/21 Python
Python 实现毫秒级淘宝抢购脚本的示例代码
2019/09/16 Python
python词云库wordcloud的使用方法与实例详解
2020/02/17 Python
OpenCV图片漫画效果的实现示例
2020/08/18 Python
Python描述数据结构学习之哈夫曼树篇
2020/09/07 Python
Pytorch - TORCH.NN.INIT 参数初始化的操作
2021/02/27 Python
怎样实现H5+CSS3手指滑动切换图片的示例代码
2019/05/05 HTML / CSS
老海军美国官网:Old Navy
2016/09/05 全球购物
找到不普通的东西:Bonanza
2016/10/20 全球购物
大学生职业生涯规划方案
2014/01/03 职场文书
电力公司个人求职信范文
2014/02/04 职场文书
《小鹰学飞》教学反思
2014/04/23 职场文书
企业法人代表任命书
2014/06/06 职场文书
投标文件签署授权委托书范本
2014/10/12 职场文书
教师批评与自我批评剖析材料
2014/10/16 职场文书
中标通知书格式
2015/04/17 职场文书
美丽的大脚观后感
2015/06/03 职场文书
同学聚会感言一句话
2015/07/30 职场文书
一篇文章带你搞懂Python类的相关知识
2021/05/20 Python
mysql 生成连续日期及变量赋值
2022/03/20 MySQL