求解一阶微分方程和二阶微分方程组（非线性）

Question

问题

我目前有一个由四个方程组成的系统。 两个是二阶微分方程，两个是一阶微分方程：

初始条件是：

x = 0 |
y = 0.3 | 
f(x) = 2.05 |
f(y) = 0.55 |

所有角度都以度为单位。

我试过的

我曾尝试使用 Google Colab，并与 SciPy 和 NumPy 一起工作。 不幸的是，我无法弄清楚如何对其进行编程，因为这些方程是非线性的。 有人可以给我一些有关要使用的任何其他模块的提示吗？

Answer 1

Scipy 有一个Runge Kutta 求解器。 首先，您必须将 ODE 转换为一阶系统（您始终可以通过设置 z=y' 来实现），然后尝试 RK 求解器。