oderf.spad line 688 [edit on github]
Solution of linear ordinary differential equations, constant coefficient case.
constDsolve(op, g, x) returns [f, [y1, ..., ym]] where f is a particular solution of the equation op y = g, and the yi's form a basis for the solutions of op y = 0.