# -*- coding:utf-8 -*-
import math
def simpr(f, a, b, n):
"""simpr函式為用復合辛普森公式求積分
f是被積函式
a,b分別為積分的上下限
n是子區間的個數
s是梯形總面積，即所求積分數值"""
h = (b - a) / (2 * n)
s1 = 0
s2 = 0
for k in range(1,n+1):
x = a + h * (2 * k -1)
s1 = s1 + f(x)
for k in range(1,n):
x = a + h * 2 * k
s2 = s2 + f(x)
s = h * (f(a) + f(b) + 4 * s1 + 2 * s2) / 3
return s
def f(g):
'''f函式為求積函式，此函式返回相應的函式值'''
return math.exp(g) + 10*g - 1
if __name__ == '__main__':
a = 0.0
b = 1.0
n = 5
sum = simpr(f,a, b, n)
print "用復合simpson公式求得積分值為：",sum
print "精確的積分值為：", 3+math.e
辛普森積分（自適應辛普森公式求積分）
自適應辛普森公式求積分 第一回接觸辛普森積分，至於這個辛普森是幹嘛的呢，在這裡就有必要好好地講一講了。應用 立體幾何中用來求擬柱體體積的公式。這裡就不詳細說辛普森公式了，有需要的朋友可以看這裡 接下來我們好好的說說自適應辛普森公式求積分，自適應辛普森公式求積分是很重要的乙個知識點，弄懂了自適應辛普森...
自適應辛普森公式求積分
假設我們求以下積分 b af x dx 比較特殊的情況，就是可以推導出來最後的形式。但是比較一般的情況是，我們只能大致得到乙個xy 座標系裡的曲線，我們求的就是曲線和 x 軸所圍成的面積。因此我們有自適應辛普森公式，他會根據實際情況來自動的調整精度。它的大致過程就是，給定乙個要求達到的精度eps，演...
python復合型別 python復合型別
鏈結 list python的列表的特點 可以隨時新增和刪除元素 可以做的操作有 切片，新增，刪除，查詢，修改，排序 列表元素可以是不同型別，在型別不同的情況下就無法排序 列表是有序的 定義小馬列表 li twilight pinkiepie rainbow dash 插入 insert self ...