“…The resulting method is of order 10, ensuring outstanding accuracy, and its implicit nature ensures its robustness and stability. This implicit Runge-Kutta algorithm has the following form Deriving the method requires sophisticated numerical analysis techniques, such as continuous collocation and interpolation, as well as the properties of Gauss-Legendre quadrature formulas and implicit Runge-Kutta methods (Laurie, 2001;Press et al, 2007;Kulikov, 2009, Sahu & Saha Ray, 2016, Mohamed et al, 2020. Gauss-Legendre methods are founded on the quadrature formula with the highest possible order, 𝑝 = 2𝑆.…”