The DFactorsols routine returns a basis of solutions of a linear ODE. Using bases of linear equations determined by DEtools[DFactor], it factorizes the equation into smaller parts.

There are two possible calling sequences. The first calling sequence accepts a linear differential equation in diff or  form as the first argument and the variable in the differential equation as the second argument.

The second form of the calling sequence accepts the list of coefficients of a linear ode as the first argument, the right-hand side of such an equation for the second argument, and the independent variable of the lode as the third argument. This form of calling sequence is convenient for programming with the DFactorsols routine.

In the second calling sequence, the list of coefficients is given in order from low differential order to high differential order and does not include the nonhomogeneous term.

In the case of a homogeneous equation, a basis is returned in the form of a list. In the nonhomogeneous case, the return value is a two-element list, with the first element a basis for the homogeneous case and the second element a particular solution.

This function is part of the DEtools package, and so it can be used in the form DFactorsols(..) only after executing the command with(DEtools). However, it can always be accessed through the long form of the command by using DEtools[DFactorsols](..).