“…In [33], we employed a combination of two multi-step formulae which gives a rigorous strategy and produces satisfactory results. However, we have to solve two linear algebraic systems at each time level which leads to a higher computational cost than the approach described in [30,31].…”