The RiemannPsols routine returns either a list of solutions of a second order linear differential equation having three regular singular points or an equation of which the solutions have been multiplied by an exponential integral.

There are two input forms. The first has as the first argument a linear differential equation in diff or D form and as the second argument the variable in the differential equation.

A second input sequence accepts for the first argument a list of coefficients of the linear ode, and for the second argument the independent variable of the lode. This input sequence is useful for programming with the RiemannPsols 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.

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