In Fig. , further comparisons for the conditions described for Fig. are performed with other numerical methods (explicit Euler, Heun, Ralston, and Runge-Kutta of 3 rd and 5 th order), also showing excellent agreement. For the relative difference defined by Eq. (68), a maximum of 𝒪(10 0? )% was found for 𝑇 * and