python实现数值积分的Simpson方法实例分析
2015-06-05 00:00
1141 查看
本文实例讲述了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 tab 补全
- Python列表和元组
- Python 脚本传参
- Python time模块
- Python文件倒读,主要用于日志分析
- Python装饰器有趣实例探究
- Django1.8+Python3.4环境搭建填坑整合
- python学习笔记
- 基于python的新浪微博模拟登陆
- python中关闭文件
- python文本编辑器Sublime text
- Python笔记——01模块安装方法
- Python图像处理(7):利用轮廓分块处理
- Python面试题集答案(1)
- 使用Python扫描端口
- 使用Python扫描端口情况
- Python 代码调试技巧
- [笔记]python的StringIO与BytesIO模块
- [Python] 机器学习库资料汇总