微分方程表示

Question

我正在尝试使用sympy处理微分方程系统：

from sympy import symbols, Function
t = symbols('t')

x_2 = Function('x_2')
x_3 = Function('x_3')

eq = x_3(t).diff(t) + x_2(t).diff(t)
eq1 = eq.subs(x_2(t), x_3(t) + x_3(t).diff(t))

答案是：

但我需要以下形式的结果：

我尝试使用

eq1.simplify()

但是结果是一样的。

我怎么能得到这个？ 谢谢。

Answer 1

您可以使用doit：

In [2]: eq1                                                                                                                       
Out[2]: 
d ⎛        d        ⎞   d        
──⎜x₃(t) + ──(x₃(t))⎟ + ──(x₃(t))
dt⎝        dt       ⎠   dt       

In [3]: eq1.doit()                                                                                                                
Out[3]: 
                2       
  d            d        
2⋅──(x₃(t)) + ───(x₃(t))
  dt            2       
              dt