我在离散一维网格上定义了两个函数 (f,g)，它们来自 pde 的离散化：
我想在给定初始条件 X0、Y0 和周期性边界条件的情况下求解它们
x[0] = x[-1] y[0] = y[-1]
如何使用 python 包作为 scipy.integrate 或类似的包来做到这一点？我已经实现了一个 RK 方案，但我猜它很慢并且提供低程度的灵活性。