# IntegrateSolutions(F, L)ΒΆ

F: Join(Field, CharacteristicZero, RetractableTo Integer, RetractableTo Fraction Integer)

L: UnivariateSkewPolynomialCategory F with

adjoint: % -> %

IntegrateSolutions implements a method to compute integral of solution of differential operator.

- integrate_sols: (L, (L, F) -> Record(particular: Union(F, failed), basis: List F)) -> Record(ltilde: L, r: Union(L, failed))
`integrate_sols(l, rat_solve)`

integrates the solutions of an operator`l`

given rat_solve(op,`g`

) that returns`["failed", []]`

if the equation`op y = g`

has no rational solution. Otherwise, the rat_solve function should return`[f, [y1, ..., ym]]`

where`f`

is a particular rational solution and the`yi`

`'s`

form a basis for the rational solutions of the homogeneous equation.