Revision history [back]

The answer returned consists of initial constants plus inverse Laplace transforms of fractions with quadratic numerators and cubic denominators. This is easier to see if you assign numeric values to all coefficients.

A fraction with a cubic denominator has a partial fraction decomposition that in general is three fractions with constant numerators and linear denominators. The constants that appear in the decomposition may be real or complex, depending on the values of the input coefficients. The inverse Laplace transform of each of these three fractions is an exponential or a circular function, depending on the values of the constants.

Your answer does consist of exponentials. It looks complicated because the 4x4 matrix of coefficients for the system, while relatively simple, leads to combinations of coefficients that are not simple, so that the inverse Laplace transforms cannot be evaluated in general form.

The answer returned consists of initial constants plus inverse Laplace transforms of fractions with quadratic numerators and cubic denominators. This is easier to see if you assign numeric values to all coefficients.coefficients. For example

eta=3; zeta=2; gamma=5
m_0=2; Qx=3; Qy=7; Qz=6; V=1
t = var('time')
w=function('w', t)
x=function('x', t)
y=function('y', t)
z=function('z', t)
dw = diff(w,t) - Qx/V*x + Qy/V*eta*w + Qz/V*eta*w == 0
dx = diff(x,t) - Qy/V*gamma*y - Qz/V*zeta*z + Qx/V*x == 0
dy = diff(y,t) - Qy/V*eta*w + Qy/V*gamma*y == 0
dz = diff(z,t) - Qz/V*eta*w + Qz/V*zeta*z == 0
eqs = [dw, dx, dy, dz]
vars = [w, x, y, z]
ics = [0, 0, m_0, 0, 0]
ans = desolve_system(eqs, vars, ics, ivar=t)
show( ans )

produces the answer

$$\left[w\left(\mathit{time}\right) = \mathcal{L}^{-1}\left(-\frac{2 {\left(10 g_{2272}^{2} + 407 g_{2272} + 2409\right)}}{161 {\left(g_{2272}^{3} + 89 g_{2272}^{2} + 2511 g_{2272} + 20286\right)}}, g_{2272}, \mathit{time}\right) + \frac{20}{161}, x\left(\mathit{time}\right) = \mathcal{L}^{-1}\left(\frac{2 {\left(31 g_{2272}^{2} + 2276 g_{2272} + 36303\right)}}{161 {\left(g_{2272}^{3} + 89 g_{2272}^{2} + 2511 g_{2272} + 20286\right)}}, g_{2272}, \mathit{time}\right) + \frac{260}{161}, y\left(\mathit{time}\right) = \mathcal{L}^{-1}\left(-\frac{6 {\left(2 g_{2272}^{2} + 178 g_{2272} + 1641\right)}}{161 {\left(g_{2272}^{3} + 89 g_{2272}^{2} + 2511 g_{2272} + 20286\right)}}, g_{2272}, \mathit{time}\right) + \frac{12}{161}, z\left(\mathit{time}\right) = \mathcal{L}^{-1}\left(-\frac{6 {\left(5 g_{2272}^{2} + 445 g_{2272} + 9657\right)}}{161 {\left(g_{2272}^{3} + 89 g_{2272}^{2} + 2511 g_{2272} + 20286\right)}}, g_{2272}, \mathit{time}\right) + \frac{30}{161}\right]$$

A fraction with a cubic denominator has a partial fraction decomposition that in general is three fractions with constant numerators and linear denominators. The constants that appear in the decomposition may be real or complex, depending on the values of the input coefficients. The inverse Laplace transform of each of these three fractions is an exponential or a circular function, depending on the values of the constants.

Your answer does consist of exponentials. It looks complicated because the 4x4 matrix of coefficients for the system, while relatively simple, leads to combinations of coefficients that are not simple, so that the inverse Laplace transforms cannot be evaluated in general form.