odeef.spad line 252 [edit on github]
In-field solution of an linear ordinary differential equation, pure algebraic case.
algDsolve(op, lg) is a parametric version of algDsolve(op, g)
algDsolve(op, g) returns ["failed", []] if the equation op y = g has no solution in R. Otherwise, it returns [f, [y1, ..., ym]] where f is a particular rational solution and the y_i's form a basis for the solutions in R of the homogeneous equation.