A \(C^ 1\) interpolant for codes based on backward differentiation formulae (Q1820562)

The objective of this work is to produce a \(C^ 1\) interpolant for a fixed-step backward differentiation formula with interpolatory step changing, with particular reference to a code which is released in the NAG Library. Although the method essentially consists of augmenting the usual Nordsieck vector, the author shows that it is unnecessary to form this directly, since the correction to the normal vector can be calculated in the interpolation routine. Numerical examples are presented to illustrate the continuity achieved in practice.