我正在寻找使用python在数学上对复杂系统建模的方法，但是，我对python还是很陌生。
我的系统由一组形式为耦合的微分方程
组成f1 = c1 * da（t）/ dt + c2 * db（t）/ dt + c3
f2 = c4 * dg（t）/ dt + c5 * dh（t）/ dt + c6
我试图用数值方法解决的。
我遵循了SciPy Cookbook，旨在将odeint用于我的程序，但现在我陷入了一个事实，即odeint需要dx / dt形式的输入，并且我无法将eq重写为该形式，即额外的差分db。 / dt不依赖于da / dt，并且不能被取消。
《 SciPy Cookbook》中的所有示例似乎都集中于以下形式的耦合微分方程
f = c1 * x'+ c2 * x''+ c3
或
f1 = c1 * x'
f2 = c2 * y'
是否有推荐的函数或其他接受形式为
的耦合微分方程的库f = c1 * da（t）/ dt + c2 * db（t）/ dt + c3
还是有办法让odeint以某种方式接受多个差异？
非常感谢您的任何建议！