用Python代码实现龙格库塔方法

好的，以下是一个简单的Python代码实现龙格库塔方法：

def runge_kutta(f, x0, y0, h, n):
    """
    执行龙格库塔方法来计算ODE的解
    :param f: 函数f(x, y)
    :param x0: 初始点x0
    :param y0: 初始值y0
    :param h: 步长h
    :param n: 步数
    :return: 迭代解的列表
    """
    res = [y0]
    for i in range(n):
        k1 = f(x0, res[-1])
        k2 = f(x0 + h / 2, res[-1] + h / 2 * k1)
        k3 = f(x0 + h / 2, res[-1] + h / 2 * k2)
        k4 = f(x0 + h, res[-1] + h * k3)
        y_next = res[-1] + h / 6 * (k1 + 2 * k2 + 2 * k3 + k4)
        res.append(y_next)
        x0 += h
    return res

其中，f是ODE的函数，x0是初始点，y0是初始值，h是步长，n是步数，函数返回一个迭代解的列表。您可以调用这个函数并使用您自己的ODE函数来测试它。