辛普森积分方法-估计误差

Question

对于n的每个值，我希望对辛普森规则的误差有一个估计。 但我不确定如何做到。

代码是：

def simpson(f, a, b, n):
h=(b-a)/n
k=0.0
x=a + h
for i in range(1,n/2 + 1):
    k += 4*f(x)
    x += 2*h

x = a + 2*h
for i in range(1,n/2):
    k += 2*f(x)
    x += 2*h
return (h/3)*(f(a)+f(b)+k)