Strange constants in desolve_system result for first-order linear system

edit

At Sagecell (running Sage 10.2) and at my local Sage 10.5beta, result is

[x(t) == 1/3*((x(0) + y(0))*sin(3*t) + 3*cos(3*t)*x(0))*e^(-t),
 y(t) == -1/3*((10*x(0) + y(0))*sin(3*t) - 3*cos(3*t)*y(0))*e^(-t)]

I updated the test case. It seems that an earlier command had an effect.

It looks like a bug. Please report at https://github.com/sagemath/sage/issues

This may or may not be a bug :

sage: reset()
sage: t=var("t")
sage: x=function("x")(t)
sage: y=function("y")(t)
sage: des=[diff(x,t) == y, diff(y,t) == -10*x + -2*y]
sage: SolG=desolve_system([diff(x,t) == y, diff(y,t) == -10*x + -2*y], [x,y]) ;
....: SolG
[x(t) == 1/3*((x(0) + 1)*sin(3*t) + 3*cos(3*t)*x(0))*e^(-t),
 y(t) == -1/3*((10*x(0) + 1)*sin(3*t) - 3*cos(3*t))*e^(-t)]

This strange-looking solution does check the original equations :

sage: list(map(lambda c:bool(reduce(lambda a, b:a.substitute_function(b.lhs().operator(), b.rhs().function()), SolG, des[0])), des))
[True, True]

But this solution may be only part of the solution.

To be continued...

Follow-up.

Sympy offers a more complete solution :

sage: import sympy
sage: SolGS=[u._sage_() for u in sympy.dsolve(*map(sympy.sympify, ([diff(x,t) ==
....:  y, diff(y,t) == -10*x + -2*y], [x,y])))] ; SolGS
[x(t) == 1/10*(3*C1 - C2)*cos(3*t)*e^(-t) + 1/10*(C1 + 3*C2)*e^(-t)*sin(3*t),
 y(t) == C2*cos(3*t)*e^(-t) - C1*e^(-t)*sin(3*t)]

which does check the original differential equations:

sage: var("C1, C2")
(C1, C2)
sage: list(map(lambda c:bool(reduce(lambda a, b:a.substitute_function(b.lhs().op
....: erator(), b.rhs().function()), SolGS, des[0])), des))
[True, True]