标签: dsolve

我正在使用SymPy版本0.7.3,并在使用dsolve函数时遇到一些问题.当输入方程具有太多参数时,似乎dsolve有困难.

我试过解决以下等式:

from sympy import *
p = Function('p')
t, u1, u2, u3, u4, u5 = symbols('t u1 u2 u3 u4 u5')
eq = Eq(Derivative(p(t),t), -(u3 + u4)*p(t) + exp(-t*(u1 + u2)))
eq
     Out: Derivative(p(t), t) == (-u3 - u4)*p(t) + exp(-t*(u1 + u2))
%time dsolve(eq)

得到了:

CPU times: user 213.11 s, sys: 0.00 s, total: 213.11 s
Wall time: 213.12 s
       p(t) == (C1 + Piecewise((t*u1/(u1 + u2) + t*u2/(u1 + u2), u3 == u1 + u2 - u4), …

在 Sympy 中 ODE 的定义（和解决方案）中是否有充分的理由使用Derivative而不是？似乎做得很好：diffdiff