in the third stage of the explicit part of the method is arbitrary. In [1]}, the above value of \(\alpha \)  was chosen with the aim of maximizing the stability region of the method. However, if a stability and absolute monotonicity analysis is carried out, as discussed in detail in Appendix ,
it can be seen that different choices might be more advantageous, in order to improve the monotonicity of the method without compromising its stability. In particular, the value of \(\alpha =1/2 \)  appears to be a more appropriate choice, as also demonstrated by the numerical experiments reported in Section .
