在python中求解具有非恒定系数的odes系统

Question

所以我有这种类型的耦合极点，我不能通过代数方法来“分离”。它们的格式为

dg/dt + df/dt*R(f,g) = T(f,g)
dg/dt + df/dt*S(f,g) = P(f,g)

其中R,T,S,P是变量f,g中的超越函数。 python中是否有直接函数来解决这类问题？