2

我想写一个四阶 Adams Bashforth 来解决这个系统。以下是我所拥有的：

系统在以下链接： 系统我们有

def AdamsBashforth4( f, x0, t ):
    """
    Fourth-order Adams-Bashforth method::

    u[n+1] = u[n] + dt/24.*(55.*f(u[n], t[n]) - 59*f(u[n-1], t[n-1]) +
                            37*f(u[n-2], t[n-2]) - 9*f(u[n-3], t[n-3]))

    for constant time step dt.

    RK2 is used as default solver for first steps.
    """
    n = len( t )
    x = numpy.array( [ x0 ] * n )


    for i in xrange( n - 1 ):
        h = t[i+1] - t[i]
        f0 = f( x[i], t[i] )
        k1 = h * f0
        k2 = h * f( x[i] + 0.5 * k1, t[i] + 0.5 * h )
        k3 = h * f( x[i] + 0.5 * k2, t[i] + 0.5 * h )
        k4 = h * f( x[i] + k3, t[i+1] )
        x[i+1] = x[i] + h * ( 55.0 * f0 - 59.0 * k1 + 37.0 * k2 - 9.0 * k3 ) / 24.0

    return x

我对吗？

当我执行它时，它给了我以下错误消息（PS：使用链接定义系统：链接）

>>>X = AdamsBashforth4(equation, init, t)

Traceback（最近一次调用最后一次）：
文件“ <pyshell#21>”，第 1 行，在<module> X = AdamsBashforth4(equation, init, t)
文件“ <pyshell#20>”中，第 15 行，在AdamsBashforth4 f0 = f( x[:-1], t[:-1] )
文件“ <pyshell#11>”中，第 2 行，在方程中x, y = X ValueError：要解压的值太多

4

0 回答 0