The ExtendRegularSolution(sol, deg) function returns the regular solution with initial terms of the Laurent series solution involved in sol extended to the specified formal degree deg.

The specified solution sol must be in the form returned by LinearFunctionalSystems[RegularSolution] or ExtendRegularSolution. In other words, sol must have an attribute of the special form as described in LinearFunctionalSystems[RegularSolution].

This function computes the additional terms of the involved series expansions using the invertible leading matrix of the matrix recurrence system corresponding to the linear functional system that was originally specified. These recurrences are part of the special structure stored in the attribute of the given solution.

The result of ExtendRegularSolution is returned in the same form as the result of LinearFunctionalSystems[RegularSolution].

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