Python中的Simpson规则

Python中的Simpson规则

您可能还忘记了

x

在第二个循环之前进行初始化，而且启动条件和迭代次数均已关闭。这是正确的方法：

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)

您的错误与循环不变性的概念有关。不必过多讨论细节，通常更容易理解和调试在循环结束时而不是在开始时进行的循环，这里我将

x += 2 *h

行移到末尾，这使验证合计从何处开始变得容易。在您的实现中，有必要

x = a - h

为第一个循环分配一个怪异的东西，只是将其添加

2 *h

为循环中的第一行。