formulas (SDGEBDFs) based on the second derivative linear multi-step formulas of Cash [1].
This class of second derivative linear multistep formulas is implemented as boundary value
methods on stiff problems. The order, error constant and the linear stability properties of the
new methods are discussed.