python实现数值积分的Simpson方法实例分析
作者:小萝莉 时间:2023-08-01 17:35:01
本文实例讲述了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,积分
![](/images/zang.png)
![](/images/jiucuo.png)
猜你喜欢
asp 类型转换函数大全第1/2页
2011-04-07 11:06:00
设计的技术含量
2009-01-12 18:20:00
PHP xpath()函数讲解
2023-06-04 07:03:38
![](https://img.aspxhome.com/file/2023/8/55398_0s.png)
Python pyecharts数据可视化实例详解
2023-06-28 01:35:30
![](https://img.aspxhome.com/file/2023/2/60712_0s.png)
存贮查询与运行时查询孰优孰劣?
2009-11-01 15:13:00
PHP禁止页面缓存的代码
2023-07-04 05:35:02
从事设计行业的十年
2008-04-01 09:44:00
解析Anaconda创建python虚拟环境的问题
2023-08-01 03:58:13
![](https://img.aspxhome.com/file/2023/6/61196_0s.png)
谈谈网页一屏有多大?
2007-12-21 12:28:00
![](https://img.aspxhome.com/file/UploadPic/200712/21/2007122112356981s.jpg)
可以举出一个最简单的计数器吗?
2009-11-01 15:37:00
Mysql Explain 详解
2010-12-03 16:09:00
用Dreamweaver MX实现网站批量更新
2009-09-13 18:39:00
![](https://img.aspxhome.com/file/UploadPic/20072/200723113154618s.jpg)
asp如何正确显示数据库里同时存在的GB码和BIG5编码?
2010-06-28 18:26:00
CSS hack浏览器兼容一览表
2007-08-14 10:35:00
![](https://img.aspxhome.com/file/UploadPic/20078/14/2007814103816867s.jpg)
用Flash来作浏览历史功能
2009-08-11 14:47:00
![](https://img.aspxhome.com/file/UploadPic/20098/11/y200988221313-19s.gif)
css设计的具有亲和力的表格
2008-03-16 19:28:00
![](https://img.aspxhome.com/file/UploadPic/20083/16/2008316193145866s.gif)
讲解使用SQL Server升级顾问的详细步骤
2009-01-04 14:14:00
PL/SQL数据类型及操作符
2009-02-26 11:17:00
insert select与select into 的用法使用说明
2012-01-05 18:47:58
一个滑动门菜单例子源码
2007-12-31 10:16:00