我有一个由 4 个微分方程组成的系统，它是一些电路网格分析的结果。 我的方程是：
i1 = Function('i1')
i2 = Function('i2')
i3 = Function('i3')
i4 = Function('i4')
t = symbols('t')
eq1 = diff(i1(t))+ integrate(i1(t)-i2(t), t)
eq5 = diff(eq1,t)
eq2 = -integrate(i1(t)-i2(t), t) + (i2(t) - i3(t)) + integrate(i2(t)-i3(t), t)
eq6 = diff(eq2,t)
eq3 = -integrate(i2(t)-i3(t), t) - (i2(t)-i3(t)) + diff(i3(t),t) +integrate(i3(t)-i4(t), t)
eq7 = diff(eq3, t)
eq4 = -integrate(i3(t)-i4(t), t) +i4(t)
eq8= diff(eq4,t)
dsolve 给我一个 MatrixError，使用时如下：
dsolve([eq5, eq6, eq7, eq8])
Jordan normal form is not implemented if the matrix have eigenvalues in CRootOf form
所以我想我必须尝试数字？