# PureAlgebraicLODE(F, UP, UPUP, R)¶

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.