Embedded explicit Runge-Kutta type methods for directly solving special third order differential equations y'''=f(x, y)

Abstract In this paper three pairs of embedded Runge–Kutta type methods for directly solving special third order ordinary differential equations (ODEs) of the form y ‴ = f ( x , y ) denoted as RKD methods are presented. The first is the RKD4(3) pair which is third order embedded in fourth-order method has the property first same as last (FSAL) whereby the last row of the coefficient matrix is equal to the vector output. The second method is the RKD5(4) pair followed by the RKD6(5) pair. The methods are derived with the strategies such that the higher order methods are very accurate and the lower order methods will give the best error estimates. Variables stepsize codes are developed based on the methods and used to solve a set of special third order problems. Numerical results are compared with the existing embedded Runge–Kutta pairs which require the problems to be reduced into a system of first order ODEs. Numerical results have clearly shown the advantage and the efficiency of the new RKD pairs.

[1]  Ernst Hairer,et al.  Solving Ordinary Differential Equations I: Nonstiff Problems , 2009 .

[2]  Tim G. Myers,et al.  Thin Films with High Surface Tension , 1998, SIAM Rev..

[3]  Ke Chen,et al.  Applied Mathematics and Computation , 2022 .

[4]  Moawwad E. A. El-Mikkawy,et al.  Families of Runge-Kutta-Nystrom Formulae , 1987 .

[5]  J. Lambert.  Numerical Methods for Ordinary Differential Equations , 1991 .

[6]  J. Dormand,et al.  A family of embedded Runge-Kutta formulae , 1980 .

[7]  Walter Gander,et al.  Derivation of Numerical Methods Using Computer Algebra , 1999, SIAM Rev..

[8]  Brian Duffy,et al.  A third-order differential equation arising in thin-film flows and relevant to Tanner's Law , 1997 .

[9]  Xiong You,et al.  Direct integrators of Runge-Kutta type for special third-order ordinary differential equations , 2013 .

[10]  Samuel N. Jator,et al.  Solving second order initial value problems by a hybrid multistep method without predictors , 2010, Appl. Math. Comput..

[11]  Fudziah Ismail,et al.  A Three-Stage Fifth-Order Runge-Kutta Method for Directly Solving Special Third-Order Differential Equation with Application to Thin Film Flow Problem , 2013 .

[12]  Fudziah Ismail,et al.  An embedded explicit Runge–Kutta–Nyström method for solving oscillatory problems , 2009 .

[13]  Zanariah Abdul Majid,et al.  A new algorithm for solving higher order IVPs of ODEs , 2011 .

[14]  Fudziah Bt Ismail,et al.  Embedded singly diagonally implicit runge-kutta methods (4, 5) in (5, 6). for the integration of stiff systems of odes , 1998, Int. J. Comput. Math..

[15]  Ebrahim Momoniat.  Symmetries, first integrals and phase planes of a third-order ordinary differential equation from thin film flow , 2009, Math. Comput. Model..

[16]  E. Fehlberg,et al.  Low-order classical Runge-Kutta formulas with stepsize control and their application to some heat transfer problems , 1969 .

[17]  J. Dormand,et al.  High order embedded Runge-Kutta formulae , 1981 .

[18]  Moawwad E. A. El-Mikkawy,et al.  A new optimized non-FSAL embedded Runge-Kutta-Nystrom algorithm of orders 6 and 4 in six stages , 2003, Appl. Math. Comput..

[19]  E. Fehlberg.  Classical Fifth-, Sixth-, Seventh-, and Eighth-Order Runge-Kutta Formulas with Stepsize Control , 1968 .

[20]  D. O. Awoyemi,et al.  A class of hybrid collocation methods for third-order ordinary differential equations , 2005, Int. J. Comput. Math..