Nordsieck Methods with an Off-Step Point

Abstract

Hybrid methods, incorporating one or more off-step points, are difficult to implement in a variable stepsize situation using the standard representation of input and output data in each step. However, instead of representing this data in terms of solution values and derivative values at a sequence of step points, it is possible to reformulate the method so that it operates on a Nordsieck vector. This has the consequence of reducing stepsize adjustments to nothing more than rescaling the components of the Nordsieck vector. This paper shows how to derive methods in both formulations and considers some implementation details. It is also possible to derive a new type of hybrid method using the Norsieck representation as the starting point and this is also discussed in the paper. The new method is found to have comparable accuracy for corresponding work expended as for standard methods.