擅长:python、mysql、java
<p>我认为你可以扩展矩阵表达式。但是你得到的不是矩阵,而是两个符号矩阵(Matsymbols)的乘法。如果你把你的表达式变成一个矩阵,你就可以得到你想要的展开式。请参阅下面的附加行</p>
<pre><code>from sympy import init_session
init_session()
from sympy import *
A = MatrixSymbol('A', 3, 3)
x = MatrixSymbol('x', 3, 1)
dt = symbols('dt')
k1 = A*x
k2 = A*(x + S(1)/2*k1*dt)
k3 = A*(x + S(1)/2*k2*dt)
k4 = A*(x + k3*dt)
final = dt*S(1)/6*(k1 + k2 + k3 + k4)
Matrix(final).expand()
</code></pre>