desolve_system error "unable to make sense of Maxima expression"
t = var('t')
y = function('y', t)
w = function('w', t)
p = function('p', t)
de1 = diff(y,t) == - 0.5*y - 0.9*p + 0.5*w
de2 = diff(w,t) == 1.0*p + 0.7*w
de3 = diff(p,t) == + 0.5*y - 0.3*w
desolve_system([de1, de2, de3], [y,p,w])
The function desolve_system stops with the error:
TypeError: unable to make sense of Maxima expression 'y(t)=ilt((200y(0)?g1090327^2+(-140y(0)+100w(0)-180p(0))?g1090327+60y(0)+54w(0)+226p(0))/(200?g1090327^3-40?g1090327^2+80?g1090327-83),?g1090327,t)' in Sage
Why? Thanks in advance!
Same holds also with ivar and ics arguments in desolve_system.
This is already http://trac.sagemath.org/sage_trac/ticket/9825.
Thanks for the answer. It is a little disappointing. Is there a workaround? Maybe I could solve it in another way with Sage (using the e-function for matrizes ...) or by using directly the Maxima functions.. Without being able to solve linear differential equations I would have to switch to other math programms.
Those both sound like good workarounds. As noted in the ticket, this is really an upstream problem to some extent; we should not see the `?g123` type symbols.
It seems like the worst problem here may have been solved upstream, see https://sourceforge.net/p/maxima/bugs/2774/