示例代碼:
def quad_custom(func, a, b, N=1000):
dx = (b - a) / N
integral = 0
for i in range(N):
x1 = a + i*dx
x2 = a + (i+1)*dx
integral += func(x1) * dx
return integral
# 定義要積分的函數(shù)
def f(x):
return x**2
# 使用自定義的積分策略進(jìn)行積分
result = quad_custom(f, 0, 1)
print(result)
在上面的代碼中,quad_custom
函數(shù)接受一個函數(shù) func
、積分區(qū)間 [a, b]
和可選的積分精度 N
,并返回對函數(shù) func
在 [a, b]
區(qū)間上的積分值。這個函數(shù)使用矩形法計(jì)算積分,將積分區(qū)間分成 N
個小矩形,并計(jì)算每個小矩形的面積,然后將它們相加得到總的積分值。
在示例中,我們定義了一個函數(shù) f(x) = x**2
,然后使用自定義的積分策略 quad_custom
對這個函數(shù)在區(qū)間 [0, 1]
上進(jìn)行積分。最后打印出積分結(jié)果。