Problem

Change the adaptive solver so that it uses the stability factor $S_c(T)$ of the problem. Apply your new adaptive solver on different IVP, e.g. the Akzo-Nobel problem.