Computing Error Constants in High Precision

In Searching for Tendler-like formulas we presented new Tendler-like cyclic linear multistep formulas for the solution of stiff differential equations. Assume the cyclic linear multistep methods of the form $$ \sum_{j=-k+1}^\ell \left[ \alpha_{ij} y_{m\ell+j} - h \beta_{ij} \dot y_{m\ell+j} \right] = 0, \quad i=1,\ldots,\ell. $$ When testing these formulas numerically, we are particularly…