使用Runge-Kutta求解四个数值耦合微分方程

问题描述

如何在Python中以数值方式求解以下四个耦合方程式？

我需要源代码!!!!!

y'(1) = (-wa+wp)*y(2)+g*y(4)
y'(2) = (-wa+2wp)*y(1)-g*y(3)
y'(3) = (wa-wp)*y(4)+g*y(2)
y'(4) =-(wa+wp)*y(3)-g*y(1)

解决方法

签出Integration and ODEs (scipy.integrate)。究竟使用哪种求解器取决于您所拥有的边界条件：特别是对于初值问题，Solve_ivp；对于边界值问题，solve_bvp。