There exists a variety of error control algorithms [1]}, [2]}, [3]}, [4]}, [5]}. In the numerical experiment that follows, we used four separate error control strategies. The standard time adaptivity I controller provides a prospective time step estimate entirely based on the current local error estimate
\(\beta _{n+1} = err_{n+1}^{-k_1 / p}.\) 
