本文实例讲述了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程序设计有所帮助。
CentOS中使用virtualenv搭建python3环境
问题描述环境:CentOS6.5想在此环境下使用python3进行开发,但CentOS6.5默认的python环境是2.6.6版本。之前的做法是直接从源码安装python3,替换掉现有的开发
Python判断字符串与大小写转换
判断字符串s.isalnum()#所有字符都是数字或者字母s.isalpha()#所有字符都是字母s.isdigit()#所有字符都是数字s.islower()#所有字符都是小写s.isupper()#所有字符都
Python文件右键找不到IDLE打开项解决办法
经常会碰到,双击.py文件运行不了,或右键没有IDLE编辑的项,在WIN7系统中比较常见.双击*.py文件运行不了解决办法:右键点击-打开方式-选择默认程序-浏