Georgios D. Trikkaliotis, Maria Ch. Gousidou-Koutita
Department of Mathematics, Aristotle University of Thessaloniki, Thessaloniki, Greece.
DOI: 10.4236/am.2022.134024   PDF    HTML   XML   333 Downloads   2,914 Views   Citations

Abstract

This paper is describing in detail the way we define the equations which give the formulas in the methods Runge-Kutta 6^th order 7 stages with the incorporated control step size in the numerical solution of Ordinary Differential Equations (ODE). The purpose of the present work is to construct a system of nonlinear equations and then by solving this system to find the values of all set parameters and finally the reduction formula of the Runge-Kutta (6,7) method (6^th order and 7 stages) for the solution of an Ordinary Differential Equation (ODE). Since the system of high order conditions required to be solved is complicated, all coefficients are found with respect to 7 free parameters. These free parameters, as well as some others in addition, are adjusted in such a way to furnish more efficient R-K methods. We use the MATLAB software to solve several of the created subsystems for the comparison of our results which have been solved analytically. Some examples for five different choices of the arbitrary values of the systems are presented in this paper.

Keywords

Initial Value Problem, Runge-Kutta Methods, Ordinary Differential Equations

Share and Cite:

1. Introduction

A system of ordinary differential equations of the form:

$y^{\prime }=f\left(x,y\right),y\left(x_0\right)=y_0$ (1)

with $x_0\in ℝ$, $y,y^{\prime }\in {ℝ}^{\text{m}}$ and $f:ℝ\times {ℝ}^m\to {ℝ}^m$ is called Initial Value Problem. Some related references are cited in the bibliography part [1] [2].

Carl David Tolmé Runge [3] and Martin Wilhelm Kutta [4] introduced the methods bearing their names almost in the turning of the 19th century. These methods are very practical and popular one-step methods for solving Ordinary Differential Equations.

The reduction formula [5] of the method Runge-Kutta for the ODE $y^{\prime }\left(x\right)=f\left(x,y\right)$, $y\left(x_0\right)=y_0$ is given by the relation $y_{n+1}=y_n+{\displaystyle {\sum }_{i=1}^{\nu }{w}_i}{K}_i$ (2) with w_ias weighting factors, ν is the number of stages and $K_i=hf\left(x_n+c_{i}h,y_n+{\displaystyle {\sum }_{j=1}^{i-1}{\alpha }_{ij}{K}_i}\right)\text{for}i=1,2,\cdots ,\nu \text{and}j=1,2,\cdots ,i-1$ (3) with h the step size of this method. The values of w_i and expressionsK_i will be determined by the values of the unknown variables α_ijand c_i.In every R-K method the relations ${\displaystyle {\sum }_{i=1}^{}{w}_i}=1\text{and}{c}_i={\displaystyle {\sum }_{j=1}^{i-1}{\alpha }_{ij}}\text{for}i=2,3,\cdots ,\nu$ (4) must be valid.

High order R-K methods have been presented by Butcher, Shanks, Lawson, Fehlberg, Feagin, Hairer, Curtis and others.

The motivation for this work is to create high order R-K (6,7) methods in a simple and practical way. In the present work after the Introduction follow the paragraph with the Presentation of the R-K method 6^th order and 7 stages, the choices of the arbitrary parameters (5 choices), the detailed presentation and the solutions of the respective systems that result, as well as the tables with the found values of the coefficients of the R-K method (6,7) and the summary expressions of K_i. Finally, the Conclusions follow.

2. Presentation of the Method

From BUTCHER’S TABLE I (Table 1) [6] the equations of the nonlinear system result, where the first column represent the order of the method, the second column the symbolism of the function f(x,y) and its derivatives, and the third column and so on the number of coefficients, for the study of Runge-Kutta methods [6], of the equations of the nonlinear system. The values of w_i, c_i and α_ij will be found by the solving this system as well as the K_i and the reduction formula for the solution of the ordinary differential equation.

The equations of the nonlinear system in their general form are: [7]

${\displaystyle {\sum }_{\kappa =1}^{\omega }{w}_{\kappa }}=1$ (5)

${\displaystyle {\sum }_{\kappa =2}^{\omega }{w}_{\kappa }}{c}_{\kappa }=\frac{1}{2}$ (6)

${\displaystyle {\sum }_{\kappa =3}^{\omega }{w}_{\kappa }}{P}_{\kappa 1}=\frac{1}{6}$ (7)

${\displaystyle {\sum }_{\kappa =2}^{\omega }{w}_{\kappa }}{c}_{\kappa }^2=\frac{1}{3}$ (8)

${\displaystyle {\sum }_{\kappa =4}^{\omega }{w}_{\kappa }}\left({\displaystyle {\sum }_{\lambda =3}^{\kappa -1}{\alpha }_{\kappa \lambda }}{P}_{\lambda 1}\right)=\frac{1}{24}$ (9)

${\displaystyle {\sum }_{\kappa =3}^{\omega }{w}_{\kappa }}{P}_{\kappa 2}=\frac{1}{12}$ (10)

${\displaystyle {\sum }_{\kappa =3}^{\omega }{w}_{\kappa }{c}_{\kappa }}{P}_{\kappa 1}=\frac{1}{8}$ (11)

${\displaystyle {\sum }_{\kappa =2}^{\omega }{w}_{\kappa }}{c}_{\kappa }^3=\frac{1}{4}$ (12)

${\displaystyle {\sum }_{\kappa =5}^{\omega }{w}_{\kappa }}\left[{\displaystyle {\sum }_{\lambda =4}^{\kappa -1}{\alpha }_{\kappa \lambda }}\left({\displaystyle {\sum }_{\mu =3}^{\lambda -1}{\alpha }_{\lambda \mu }}{P}_{\mu 1}\right)\right]=\frac{1}{120}$ (13)

${\displaystyle {\sum }_{\kappa =4}^{\omega }{w}_{\kappa }}\left({\displaystyle {\sum }_{\lambda =3}^{\omega }{\alpha }_{\kappa \lambda }}{P}_{\lambda 2}\right)=\frac{1}{60}$ (14)

${\displaystyle {\sum }_{\kappa =4}^{\omega }{w}_{\kappa }}\left({\displaystyle {\sum }_{\lambda =3}^{\kappa -1}{\alpha }_{\kappa \lambda }}{c}_{\lambda }{P}_{\lambda 1}\right)=\frac{1}{40}$ (15)

${\displaystyle {\sum }_{\kappa =3}^{\omega }{w}_{\kappa }}{P}_{\kappa 3}=\frac{1}{20}$ (16)

${\displaystyle {\sum }_{\kappa =4}^{\omega }{w}_{\kappa }{c}_{\kappa }}\left({\displaystyle {\sum }_{\lambda =3}^{\kappa -1}{\alpha }_{\kappa \lambda }}{P}_{\lambda 1}\right)=\frac{1}{30}$ (17)

${\displaystyle {\sum }_{\kappa =3}^{\omega }{w}_{\kappa }{c}_{\kappa }}{P}_{\kappa 2}=\frac{1}{15}$ (18)

${\displaystyle {\sum }_{\kappa =3}^{\omega }{w}_{\kappa }}{P}_{\kappa 1}^2=\frac{1}{20}$ (19)

${\displaystyle {\sum }_{\kappa =3}^{\omega }{w}_{\kappa }{c}_{\kappa }^2}{P}_{\kappa 1}=\frac{1}{10}$ (20)

${\displaystyle {\sum }_{\kappa =2}^{\omega }{w}_{\kappa }}{c}_{\kappa }^4=\frac{1}{5}$ (21)

${\displaystyle {\sum }_{\kappa =6}^{\omega }{w}_{\kappa }\left\{{\displaystyle {\sum }_{\lambda =5}^{\kappa -1}{\alpha }_{\kappa \lambda }\left[{\displaystyle {\sum }_{\mu =4}^{\lambda -1}{\alpha }_{\lambda \mu }\left({\displaystyle {\sum }_{\nu =3}^{\mu -1}{\alpha }_{\mu \nu }}{P}_{\nu 1}\right)}\right]}\right\}}=\frac{1}{720}$ (22)

${\displaystyle {\sum }_{\kappa =5}^{\omega }{w}_{\kappa }}\left[{\displaystyle {\sum }_{\lambda =4}^{\kappa -1}{\alpha }_{\kappa \lambda }}\left({\displaystyle {\sum }_{\mu =3}^{\lambda -1}{\alpha }_{\lambda \mu }}{P}_{\mu 2}\right)\right]=\frac{1}{360}$ (23)

${\displaystyle {\sum }_{\kappa =5}^{\omega }{w}_{\kappa }}\left[{\displaystyle {\sum }_{\lambda =4}^{\kappa -1}{\alpha }_{\kappa \lambda }}\left({\displaystyle {\sum }_{\mu =3}^{\lambda -1}{\alpha }_{\lambda \mu }}{c}_{\mu }{P}_{\mu 1}\right)\right]=\frac{1}{240}$ (24)

${\displaystyle {\sum }_{\kappa =4}^{\omega }{w}_{\kappa }}\left({\displaystyle {\sum }_{\lambda =3}^{\kappa -1}{\alpha }_{\kappa \lambda }}{P}_{\lambda 2}\right)=\frac{1}{120}$ (25)

${\displaystyle {\sum }_{\kappa =5}^{\omega }{w}_{\kappa }}\left[{\displaystyle {\sum }_{\lambda =4}^{\kappa -1}{\alpha }_{\kappa \lambda }{c}_{\lambda }}\left({\displaystyle {\sum }_{\mu =3}^{\lambda -1}{\alpha }_{\lambda \mu }}{P}_{\mu 1}\right)\right]=\frac{1}{180}$ (26)

${\displaystyle {\sum }_{\kappa =4}^{\omega }{w}_{\kappa }}\left({\displaystyle {\sum }_{\lambda =3}^{\kappa -1}{\alpha }_{\kappa \lambda }}{c}_{\lambda }{P}_{\lambda 2}\right)=\frac{1}{90}$ (27)

${\displaystyle {\sum }_{\kappa =4}^{\omega }{w}_{\kappa }}\left({\displaystyle {\sum }_{\lambda =3}^{\kappa -1}{\alpha }_{\kappa \lambda }}{P}_{\lambda 1}^2\right)=\frac{1}{120}$ (28)

${\displaystyle {\sum }_{\kappa =4}^{\omega }{w}_{\kappa }}\left({\displaystyle {\sum }_{\lambda =3}^{\kappa -1}{\alpha }_{\kappa \lambda }}{c}_{\lambda }^2{P}_{\lambda 1}\right)=\frac{1}{60}$ (29)

${\displaystyle {\sum }_{\kappa =3}^{\omega }{w}_{\kappa }}{P}_{\kappa 4}=\frac{1}{30}$ (30)

${\displaystyle {\sum }_{\kappa =5}^{\omega }{w}_{\kappa }{c}_{\kappa }}\left[{\displaystyle {\sum }_{\lambda =4}^{\kappa -1}{\alpha }_{\kappa \lambda }}\left({\displaystyle {\sum }_{\mu =3}^{\lambda -1}{\alpha }_{\lambda \mu }}{P}_{\mu 1}\right)\right]=\frac{1}{144}$ (31)

${\displaystyle {\sum }_{\kappa =4}^{\omega }{w}_{\kappa }{c}_{\kappa }}\left({\displaystyle {\sum }_{\lambda =3}^{\kappa -1}{\alpha }_{\kappa \lambda }}{P}_{\lambda 2}\right)=\frac{1}{72}$ (32)

${\displaystyle {\sum }_{\kappa =4}^{\omega }{w}_{\kappa }{c}_{\kappa }}\left({\displaystyle {\sum }_{\lambda =3}^{\kappa -1}{\alpha }_{\kappa \lambda }}{c}_{\lambda }{P}_{\lambda 1}\right)=\frac{1}{48}$ (33)

${\displaystyle {\sum }_{\kappa =3}^{\omega }{w}_{\kappa }{c}_{\kappa }}{P}_{\kappa 3}=\frac{1}{24}$ (34)

${\displaystyle {\sum }_{\kappa =4}^{\omega }{w}_{\kappa }{P}_{\kappa 1}}\left({\displaystyle {\sum }_{\lambda =3}^{\kappa -1}{\alpha }_{\kappa \lambda }}{P}_{\lambda 1}\right)=\frac{1}{72}$ (35)

${\displaystyle {\sum }_{\kappa =3}^{\omega }{w}_{\kappa }}{P}_{\kappa 1}{P}_{\kappa 2}=\frac{1}{36}$ (36)

${\displaystyle {\sum }_{\kappa =4}^{\omega }{w}_{\kappa }{c}_{\kappa }^2}\left({\displaystyle {\sum }_{\lambda =3}^{\kappa -1}{\alpha }_{\kappa \lambda }}{P}_{\lambda 1}\right)=\frac{1}{36}$ (37)

${\displaystyle {\sum }_{\kappa =3}^{\omega }{w}_{\kappa }{c}_{\kappa }^2}{P}_{\kappa 2}=\frac{1}{18}$ (38)

${\displaystyle {\sum }_{\kappa =3}^{\omega }{w}_{\kappa }{c}_{\kappa }}{P}_{\kappa 1}^2=\frac{1}{24}$ (39)

${\displaystyle {\sum }_{\kappa =3}^{\omega }{w}_{\kappa }{c}_{\kappa }^3}{P}_{\kappa 1}=\frac{1}{12}$ (40)

${\displaystyle {\sum }_{\kappa =2}^{\omega }{w}_{\kappa }}{c}_{\kappa }^5=\frac{1}{6}$ (41)

with

$\omega =7$,(42)

$P_{\kappa \lambda }={\alpha }_{\kappa 2}{c}_2^{\lambda }+{\alpha }_{\kappa 3}{c}_3^{\lambda }+\cdots +{\alpha }_{\kappa \kappa -1}{c}_{\kappa -1}^{\lambda }\text{with}\kappa =3,4,\cdots ,7\text{and}\lambda =1,2,3,4$ [5] (43)

2.1. 1^st Choice: c₂ = c₄ = 1/3, c₃ = 2/3, c₅ = c₆ = 1/2, c₇ = 1 [5]

From the system of (5), (6), (8), (12) και (41) results that:

$w_1=\frac{11}{120}$ (44) $w_2+w_4=\frac{27}{40}$ (45) $w_3=\frac{27}{40}$ (46) $w_5+w_6=-\frac{8}{15}$ (47) $w_7=\frac{11}{120}$ (48) and setting $w_2=0$ (49) and $w_5=w_6$ (50) we have $w_4=\frac{27}{40}$ (51) and $w_5=w_6=-\frac{4}{15}$ (52) while the Equation (21) verified.

Since the above equations become somewhat lengthy, we introduce the following abbreviations: [7]

${\alpha }_{42}{c}_2+a_{43}{c}_3=\frac{{\alpha }_{42}+2{\alpha }_{43}}{3}=\frac{S_4}{3}$ (53)

${\alpha }_{52}{c}_2+a_{53}{c}_3+{\alpha }_{54}{c}_4=\frac{{\alpha }_{52}+2{\alpha }_{53}+{\alpha }_{54}}{3}=\frac{S_5}{3}$ (54)

${\alpha }_{62}{c}_2+a_{63}{c}_3+{\alpha }_{64}{c}_4+{\alpha }_{65}{c}_5=\frac{2{\alpha }_{62}+4{\alpha }_{63}+2{\alpha }_{64}+3{\alpha }_{65}}{6}=\frac{S_6}{6}$ (55)

${\alpha }_{72}{c}_2+a_{73}{c}_3+{\alpha }_{74}{c}_4+{\alpha }_{75}{c}_5+{\alpha }_{76}{c}_6=\frac{2{\alpha }_{72}+4{\alpha }_{73}+2{\alpha }_{74}+3{\alpha }_{75}+3{\alpha }_{76}}{6}=\frac{S_7}{6}$ (56)

and

${\alpha }_{42}{c}_2^2+a_{43}{c}_3^2=\frac{{\alpha }_{42}+4{\alpha }_{43}}{9}=\frac{S_4+2{\alpha }_{43}}{9}$ (57)

${\alpha }_{52}{c}_2^2+a_{53}{c}_3^2+{\alpha }_{54}{c}_4^2=\frac{S_5+2{\alpha }_{53}}{9}$ (58)

${\alpha }_{62}{c}_2^2+a_{63}{c}_3^2+{\alpha }_{64}^{}{c}_4^2+{\alpha }_{65}{c}_5^2=\frac{2S_6+8{\alpha }_{63}+3{\alpha }_{65}}{36}$ (59)

${\alpha }_{72}{c}_2^2+a_{73}{c}_3^2+{\alpha }_{74}^{}{c}_4^2+{\alpha }_{75}{c}_5^2+a_{76}{c}_6^2=\frac{2S_7+8{\alpha }_{73}+3{\alpha }_{75}+3{\alpha }_{76}}{36}$ (60)

and

${\alpha }_{42}{c}_2^3+a_{43}{c}_3^3=\frac{S_4+6{\alpha }_{43}}{27}$ (61)

${\alpha }_{52}{c}_2^3+a_{53}{c}_3^3+a_{54}{c}_4^3=\frac{S_5+6{\alpha }_{53}}{27}$ (62)

${\alpha }_{62}{c}_2^3+a_{63}{c}_3^3+{\alpha }_{64}^{}{c}_4^3+{\alpha }_{65}{c}_5^3=\frac{4S_6+48{\alpha }_{63}+15{\alpha }_{65}}{216}$ (63)

${\alpha }_{72}{c}_2^3+a_{73}{c}_3^3+{\alpha }_{74}^{}{c}_4^3+{\alpha }_{75}{c}_5^3+a_{76}{c}_6^3=\frac{4S_7+48{\alpha }_{73}+15{\alpha }_{75}+15{\alpha }_{76}}{216}$ (64)

and

${\alpha }_{42}{c}_2^4+a_{43}{c}_3^4=\frac{S_4+14{\alpha }_{43}}{81}$ (65)

${\alpha }_{52}{c}_2^4+a_{53}{c}_3^4+{\alpha }_{54}^{}{c}_4^4=\frac{S_5+14{\alpha }_{53}}{81}$ (66)

${\alpha }_{62}{c}_2^4+a_{63}{c}_3^4+{\alpha }_{64}^{}{c}_4^4+{\alpha }_{65}{c}_5^4=\frac{8S_6+224{\alpha }_{63}+57{\alpha }_{65}}{1296}$ (67)

${\alpha }_{72}{c}_2^4+a_{73}{c}_3^4+{\alpha }_{74}^{}{c}_4^4+{\alpha }_{75}{c}_5^4+a_{76}{c}_6^4=\frac{8S_7+224{\alpha }_{73}+57{\alpha }_{75}+57{\alpha }_{76}}{1296}$ (68)

Then we substitute the defined abbreviations in the original system, as well as the found values of c₂, c₃, c₄, c₅, c₆, c₇, w₂, w₃, w₄, w₅, w₆, w₇, and as a result the system is simplified. Also omitting the equations which are a linear combination of equations of the system the 30 × 15 system is obtained:

$162{\alpha }_{32}+162S_4-64S_5-32S_6+11S_7=120$ (69)

$\begin{array}{l}162{\alpha }_{43}{\alpha }_{32}-64\left({\alpha }_{53}{\alpha }_{32}+{\alpha }_{54}{S}_4\right)-64\left({\alpha }_{63}{\alpha }_{32}+{\alpha }_{64}{S}_4+{\alpha }_{65}{S}_5\right)\\ +11\left(2{\alpha }_{73}{\alpha }_{32}+2{\alpha }_{74}{S}_4+2{\alpha }_{75}{S}_5+{\alpha }_{76}{S}_6\right)=30\end{array}$ (70)

$648{\alpha }_{43}-256{\alpha }_{53}-32\left(8{\alpha }_{63}+3{\alpha }_{65}\right)+11\left(8{\alpha }_{73}+3{\alpha }_{75}+3{\alpha }_{76}\right)=120$ (71)

$108{\alpha }_{32}+54S_4-32S_5-16S_6+11S_7=90$ (72)

$\begin{array}{l}-32{\alpha }_{54}{\alpha }_{43}{\alpha }_{32}-32\left[{\alpha }_{64}{\alpha }_{43}{\alpha }_{32}+{\alpha }_{65}\left({\alpha }_{53}{\alpha }_{32}+{\alpha }_{54}{S}_4\right)\right]\\ +11\left[{\alpha }_{74}{\alpha }_{43}{\alpha }_{32}+{\alpha }_{75}\left({\alpha }_{53}{\alpha }_{32}+{\alpha }_{54}{S}_4\right)+{\alpha }_{76}\left({\alpha }_{63}{\alpha }_{32}+{\alpha }_{64}{S}_4+{\alpha }_{65}{S}_5\right)\right]=3\end{array}$ (73)

$\begin{array}{l}256\left({\alpha }_{54}{\alpha }_{43}+{\alpha }_{64}{\alpha }_{43}+{\alpha }_{65}{\alpha }_{53}\right)-88\left({\alpha }_{74}{\alpha }_{43}+{\alpha }_{75}{\alpha }_{53}+{\alpha }_{76}{\alpha }_{63}\right)\\ -33{\alpha }_{76}{\alpha }_{65}=-12\end{array}$ (74)

$128{\alpha }_{54}{S}_4+128{\alpha }_{64}{S}_4+64{\alpha }_{65}{S}_5-44{\alpha }_{74}{S}_4-22{\alpha }_{75}{S}_5-11{\alpha }_{76}{S}_6=-12$ (75)

$32{\alpha }_{65}-11{\alpha }_{75}-11{\alpha }_{76}=32$ (76)

$54{\alpha }_{43}{\alpha }_{32}-16\left({\alpha }_{53}{\alpha }_{32}+{\alpha }_{54}{S}_4\right)-16\left({\alpha }_{63}{\alpha }_{32}+{\alpha }_{64}{S}_4+{\alpha }_{65}{S}_5\right)=3$ (77)

$108{\alpha }_{43}-32{\alpha }_{53}-32{\alpha }_{63}-12{\alpha }_{65}=3$ (78)

$324{\alpha }_{32}^2+324S_4^2-128S_5^2-32S_6^2+11S_7^2=216$ (79)

$72{\alpha }_{32}+18S_4-16S_5-8S_6+11S_7=72$ (80)

$64{\alpha }_{65}{\alpha }_{54}{\alpha }_{43}{\alpha }_{32}-22\left\{{\alpha }_{75}{\alpha }_{54}{\alpha }_{43}{\alpha }_{32}+{\alpha }_{76}\left[{\alpha }_{64}{\alpha }_{43}{\alpha }_{32}+{\alpha }_{65}\left({\alpha }_{53}{\alpha }_{32}+{\alpha }_{54}{S}_4\right)\right]\right\}=-1$

(81)

$32{\alpha }_{65}{\alpha }_{54}{\alpha }_{43}-11\left[{\alpha }_{75}{\alpha }_{54}{\alpha }_{43}+{\alpha }_{76}\left({\alpha }_{64}{\alpha }_{43}+{\alpha }_{65}{\alpha }_{53}\right)\right]=0$ (82)

$64{\alpha }_{65}{\alpha }_{54}{S}_4-22{\alpha }_{75}{\alpha }_{54}{S}_4-11{\alpha }_{76}\left(2{\alpha }_{64}{S}_4+{\alpha }_{65}{S}_5\right)=-3$ (83)

$11{\alpha }_{76}{\alpha }_{65}=-8$ (84)

$64{\alpha }_{65}{\alpha }_{53}{\alpha }_{32}-11\left[2{\alpha }_{75}{\alpha }_{53}{\alpha }_{32}+{\alpha }_{76}\left(2{\alpha }_{63}{\alpha }_{32}+{\alpha }_{65}{S}_5\right)\right]=-9$ (85)

$\begin{array}{l}648{\alpha }_{43}{\alpha }_{32}-256{\alpha }_{53}{\alpha }_{32}-128\left[2{\alpha }_{63}{\alpha }_{32}+{\alpha }_{65}\left(2{\alpha }_{53}+S_5\right)\right]\\ +11\left[8{\alpha }_{73}{\alpha }_{32}+4{\alpha }_{75}\left(2{\alpha }_{53}+S_5\right)+{\alpha }_{76}\left(8{\alpha }_{63}+3{\alpha }_{65}+2S_6\right)\right]=144\end{array}$ (86)

$\begin{array}{l}324{\alpha }_{43}{\alpha }_{32}^2-128\left({\alpha }_{53}{\alpha }_{32}^2+{\alpha }_{54}{S}_4^2\right)-128\left({\alpha }_{63}{\alpha }_{32}^2+{\alpha }_{64}{S}_4^2+{\alpha }_{65}{S}_5^2\right)\\ +11\left(4{\alpha }_{73}{\alpha }_{32}^2+4{\alpha }_{74}{S}_4^2+4{\alpha }_{75}{S}_5^2+{\alpha }_{76}{S}_6^2\right)=36\end{array}$ (87)

$\begin{array}{l}1944{\alpha }_{43}{\alpha }_{32}-768{\alpha }_{53}{\alpha }_{32}-768{\alpha }_{63}{\alpha }_{32}-320{\alpha }_{65}{S}_5+264a_{73}{a}_{32}\\ +110{\alpha }_{75}{S}_5+55{\alpha }_{76}{S}_6=312\end{array}$ (88)

$11\left[{\alpha }_{74}{\alpha }_{43}{\alpha }_{32}+{\alpha }_{75}\left({\alpha }_{53}{\alpha }_{32}+{\alpha }_{54}{S}_4\right)+{\alpha }_{76}\left({\alpha }_{63}{\alpha }_{32}+{\alpha }_{64}{S}_4+{\alpha }_{65}{S}_5\right)\right]=2$ (89)

$\begin{array}{l}-54{\alpha }_{43}{\alpha }_{32}+128{\alpha }_{54}{\alpha }_{43}+128\left({\alpha }_{64}{\alpha }_{43}+{\alpha }_{65}{\alpha }_{53}\right)\\ +11\left(2{\alpha }_{73}{\alpha }_{32}+2{\alpha }_{74}{S}_4+2{\alpha }_{75}{S}_5+2{\alpha }_{76}{S}_6\right)=18\end{array}$ (90)

$-216{\alpha }_{43}{\alpha }_{32}+88{\alpha }_{73}{\alpha }_{32}+44{\alpha }_{74}{S}_4+66{\alpha }_{75}{S}_5+33{\alpha }_{76}{S}_6=72$ (91)

$-54{\alpha }_{43}+22{\alpha }_{65}+22{\alpha }_{73}=47$ (92)

$\begin{array}{l}324S_4{\alpha }_{43}{\alpha }_{32}-128S_5\left({\alpha }_{53}{\alpha }_{32}+{\alpha }_{54}{S}_4\right)-64S_6\left({\alpha }_{63}{\alpha }_{32}+{\alpha }_{64}{S}_4+{\alpha }_{65}{S}_5\right)\\ +11S_7\left(2{\alpha }_{73}{\alpha }_{32}+2{\alpha }_{74}{S}_4+2{\alpha }_{75}{S}_5+{\alpha }_{76}{S}_6\right)=60\end{array}$ (93)

$\begin{array}{l}1296{\alpha }_{43}{S}_4-512{\alpha }_{53}{S}_5-\left(256{\alpha }_{63}+396{\alpha }_{65}\right)S_6\\ +\left(88{\alpha }_{73}+33{\alpha }_{75}+33{\alpha }_{76}\right)S_7=288\end{array}$ (94)

$18{\alpha }_{43}{\alpha }_{32}=-1$ (95)

$72{\alpha }_{43}-16{\alpha }_{53}-16{\alpha }_{63}-4{\alpha }_{65}=1$ (96)

$108{\alpha }_{32}^2-108S_4^2+11S_7^2=144$ (97)

$24{\alpha }_{32}-6S_4+11S_7=48$ (98)

From (69), (72), (80) και (98) result:

${\alpha }_{52}=\frac{2}{3}$ (99)

$S_4=\frac{1}{6}$ (100)

$S_7=3$ (101)

and

$4S_5+2S_6=3$ (102)

From (79) we find:

$32S_5^2+8S_6^2=9$ (103)

so:

$S_5=\frac{3}{8}$ (104)

$S_6=\frac{3}{4}$ (105)

From (71), (76), (78), (92) and (96) result:

$a_{43}=-\frac{1}{12}$ (106)

$a_{65}=\frac{1}{2}$ (107)

$a_{73}=\frac{63}{44}$ (108)

and

$a_{53}+a_{63}=-\frac{9}{16}$ (109)

$a_{75}+a_{76}=-\frac{16}{11}$ (110)

From (84):

$a_{76}=-\frac{16}{11}$ (111)

and from (110):

$a_{75}=0$ (112)

while the Equations (95) and (97) are verified.

By substitution the found values from (75), (77) and (109) we find

$a_{74}=\frac{18}{11}$ (113)

and

$a_{54}+a_{64}=-\frac{9}{8}$ (114)

from (81):

$4a_{53}+a_{54}=-\frac{9}{8}$ (115)

and from (90):

$192a_{53}-32\left(a_{54}+a_{64}\right)=0$ (116)

The solution of the system of Equations (109), (114), (115) and (116) is:

$a_{53}=-\frac{3}{16}$ (117)

$a_{54}=-\frac{6}{16}$ (118)

$a_{63}=-\frac{6}{16}$ (119)

and

$a_{64}=-\frac{12}{16}$ (120)

Using the defined abbreviations $S_4={\alpha }_{42}+2{\alpha }_{43}$, $S_5={\alpha }_{52}+2{\alpha }_{53}+{\alpha }_{54}$

$S_6=2{\alpha }_{62}+4{\alpha }_{63}+2{\alpha }_{64}+3{\alpha }_{65}$ and $S_7=2{\alpha }_{72}+4{\alpha }_{73}+2{\alpha }_{74}+3{\alpha }_{75}+3{\alpha }_{76}$ result:

$a_{42}=\frac{4}{12}$ (121)

$a_{52}=\frac{18}{16}$ (122)

$a_{62}=\frac{9}{8}$ (123)

$a_{72}=-\frac{9}{11}$ (124)

From the relations ${\displaystyle {\sum }_{j=1}^{i-1}{\alpha }_{ij}}=c_i$,for $i=2,3,4,5,6,7$ result:

$a_{21}=\frac{1}{3}$ (125)

$a_{41}=\frac{1}{12}$ (126)

$a_{51}=-\frac{1}{16}$ (127)

$a_{61}=0$ (128)

and

$a_{71}=\frac{9}{44}$ (129)

The remaining equations are verified.

Runge-Kutta methods usually presented in a so-called Butcher table (Table 2) [8].

The table contains on the 1^st column the coefficients c_i, the matrix A with the coefficients of a_ij, which appear in the Formulae of K_i, and w_i the coefficients in Formula of y_i+1.

In this type of table, we have $w^{\text{T}},c\in {ℝ}^{\text{m}}$ while $A\in {ℝ}^{m\times m}$. Then, the method shares m stages and in case that c₁ = 0 and A [5] is strictly lower triangular, it is evaluated explicitly.

Therefore we get Table 3:

with

$K_1=hf\left(x_n,y_n\right)$ (130)

$K_2=hf\left(x_n+\frac{h}{3},y_n+\frac{K_1}{3}\right)$ (131)

$K_3=hf\left(x_n+\frac{2h}{3},y_n+\frac{2K_2}{3}\right)$ (132)

$K_4=hf\left(x_n+\frac{h}{3},y_n+\frac{K_1+4K_2-K_3}{12}\right)$ (133)

$K_5=hf\left(x_n+\frac{h}{2},y_n+\frac{-K_1+18K_2-3K_3-6K_4}{16}\right)$ (134)

$K_6=hf\left(x_n+\frac{h}{2},y_n+\frac{9K_2-3K_3-6K_4+4K_5}{8}\right)$ (135)

$K_7=hf\left(x_n+h,y_n+\frac{9K_1-36K_2+63K_3+72K_4-64K_5}{44}\right)$ (136)

and

$y_{n+1}=y_n+\frac{11K_1+81K_3+81K_4-32K_5-32K_6+11K_7}{120}$ (137)

2.2. 2^nd Choice: c₂ = c₃ = 1/4, c₄ = 2/4, c₅ = c₆ = 3/4, c₇ = 4/4 [5]

From the system of (5), (6), (8), (12) και (41) results that:

$w_1=\frac{7}{90}$ (138)

$w_2+w_3=\frac{32}{90}$ (139)

$w_4=\frac{12}{90}$ (140)

$w_5+w_6=\frac{32}{90}$ (141)

$w_7=\frac{7}{90}$ (142) and setting $w_2=0$ (143) and $w_5=w_6$ (144) we get

$w_3=\frac{32}{90}$ (145) and $w_5=w_6=\frac{16}{90}$ (146) while the Equation (21) verified.

Since the above equations become somewhat lengthy, we introduce the following abbreviations

${\alpha }_{42}{c}_2+a_{43}{c}_3=\frac{{\alpha }_{42}+{\alpha }_{43}}{4}=\frac{S_4}{4}$ (147)

${\alpha }_{52}{c}_2^{}+a_{53}{c}_3^{}+{\alpha }_{54}{c}_4=\frac{{\alpha }_{52}+{\alpha }_{53}+2{\alpha }_{54}}{4}=\frac{S_5}{4}$ (148)

${\alpha }_{62}{c}_2^{}+a_{63}{c}_3^{}+{\alpha }_{64}^{}{c}_4^{}+{\alpha }_{65}{c}_5^{}=\frac{{\alpha }_{62}+{\alpha }_{63}+2{\alpha }_{64}+3{\alpha }_{65}}{4}=\frac{S_6}{4}$ (149)

${\alpha }_{72}{c}_2^{}+a_{73}{c}_3^{}+{\alpha }_{74}^{}{c}_4^{}+{\alpha }_{75}{c}_5^{}+{\alpha }_{76}{c}_6=\frac{{\alpha }_{72}+{\alpha }_{73}+2{\alpha }_{74}+3{\alpha }_{75}+3{\alpha }_{76}}{4}=\frac{S_7}{4}$ (150)

and

${\alpha }_{42}{c}_2^2+a_{43}{c}_3^2=\frac{{\alpha }_{42}+{\alpha }_{43}}{16}=\frac{S_4}{16}$ (151)

${\alpha }_{52}{c}_2^2+a_{53}{c}_3^2+{\alpha }_{54}{c}_4^2=\frac{S_5+2{\alpha }_{54}}{16}$ (152)

${\alpha }_{62}{c}_2^2+a_{63}{c}_3^2+{\alpha }_{64}^{}{c}_4^2+{\alpha }_{65}{c}_5^2=\frac{S_6+2{\alpha }_{64}+6{\alpha }_{65}}{16}$ (153)

${\alpha }_{72}{c}_2^2+a_{73}{c}_3^2+{\alpha }_{74}^{}{c}_4^2+{\alpha }_{75}{c}_5^2+a_{76}{c}_6^2=\frac{S_7+2{\alpha }_{74}+6{\alpha }_{75}+6{\alpha }_{76}}{16}$ (154)

and

${\alpha }_{42}{c}_2^3+a_{43}{c}_3^3=\frac{S_4}{64}$ (155)

${\alpha }_{52}{c}_2^3+a_{53}{c}_3^3+{\alpha }_{54}{c}_{{}_4}^3=\frac{S_5+6{\alpha }_{54}}{64}$ (156)

${\alpha }_{62}{c}_2^3+a_{63}{c}_3^3+{\alpha }_{64}^{}{c}_4^3+{\alpha }_{65}{c}_5^3=\frac{S_6+{\alpha }_{64}+24{\alpha }_{65}}{64}$ (157)

${\alpha }_{72}{c}_2^3+a_{73}{c}_3^3+{\alpha }_{74}^{}{c}_4^3+{\alpha }_{75}{c}_5^3+a_{76}{c}_6^3=\frac{S_7+6{\alpha }_{74}+24{\alpha }_{75}+24{\alpha }_{76}}{64}$ (158)

and

${\alpha }_{42}{c}_2^4+a_{43}{c}_3^4=\frac{S_4}{256}$ (159)

${\alpha }_{52}{c}_2^4+a_{53}{c}_3^4+{\alpha }_{54}^{}{c}_4^4=\frac{S_5+14{\alpha }_{54}}{256}$ (160)

${\alpha }_{62}{c}_2^4+a_{63}{c}_3^4+{\alpha }_{64}^{}{c}_4^4+{\alpha }_{65}{c}_5^4=\frac{S_6+4{\alpha }_{64}+78{\alpha }_{65}}{256}$ (161)

${\alpha }_{72}{c}_2^4+a_{73}{c}_3^4+{\alpha }_{74}^{}{c}_4^4+{\alpha }_{75}{c}_5^4+a_{76}{c}_6^4=\frac{S_7+14{\alpha }_{74}+78{\alpha }_{75}+78{\alpha }_{76}}{256}$ (162)

Then we substitute the defined abbreviations in the original system, as well as the found values of c₂, c₃, c₄, c₅, c₆, c₇, w₂, w₃, w₄, w₅, w₆, w₇, and as a result the system is simplified. Also omitting the equations which are a linear combination of equations of the system and the 30 × 15 system is obtained:

$32{\alpha }_{32}+12S_4+16S_5+16S_6+7S_7=60$ (163)

$\begin{array}{l}12{\alpha }_{43}{\alpha }_{32}+16\left({\alpha }_{53}{\alpha }_{32}+{\alpha }_{54}{S}_4\right)+16\left({\alpha }_{63}{\alpha }_{32}+{\alpha }_{64}{S}_4+{\alpha }_{65}{S}_5\right)\\ +7\left({\alpha }_{73}{\alpha }_{32}+{\alpha }_{74}{S}_4+{\alpha }_{75}{S}_5+{\alpha }_{76}{S}_6\right)=15\end{array}$ (164)

$16\left(2{\alpha }_{54}\right)+16\left(2{\alpha }_{64}+6{\alpha }_{65}\right)+7\left(2{\alpha }_{74}+6{\alpha }_{75}+6{\alpha }_{76}\right)=60$ (165)

$32{\alpha }_{32}+24S_4+48S_5+48S_6+28S_7=180$ (166)

$\begin{array}{l}16{\alpha }_{54}{\alpha }_{43}{\alpha }_{32}+16\left[{\alpha }_{64}{\alpha }_{43}{\alpha }_{32}+{\alpha }_{65}\left({\alpha }_{53}{\alpha }_{32}+{\alpha }_{54}{S}_4\right)\right]\\ +7\left[{\alpha }_{74}{\alpha }_{43}{\alpha }_{32}+{\alpha }_{75}\left({\alpha }_{53}{\alpha }_{32}+{\alpha }_{54}{S}_4\right)+{\alpha }_{76}\left({\alpha }_{63}{\alpha }_{32}+{\alpha }_{64}{S}_4+{\alpha }_{65}{S}_5\right)\right]=3\end{array}$ (167)

$16{\alpha }_{65}\left(2{\alpha }_{54}\right)+7\left[{\alpha }_{75}\left(2{\alpha }_{54}\right)+{\alpha }_{76}\left(2{\alpha }_{64}+6{\alpha }_{65}\right)\right]=9$ (168)

$16{\alpha }_{54}{S}_4+16\left({\alpha }_{64}{S}_4+2{\alpha }_{65}{S}_5\right)+7\left({\alpha }_{74}{S}_4+2{\alpha }_{75}{S}_5+2{\alpha }_{76}{S}_6\right)=21$ (169)

$16{\alpha }_{65}+7\left({\alpha }_{75}+{\alpha }_{76}\right)=8$ (170)

$6{\alpha }_{43}{\alpha }_{32}+4\left({\alpha }_{53}{\alpha }_{32}+{\alpha }_{54}{S}_4\right)+4\left({\alpha }_{63}{\alpha }_{32}+{\alpha }_{64}{S}_4+{\alpha }_{65}{S}_5\right)=3$ (171)

$7\left({\alpha }_{74}+3{\alpha }_{75}+3{\alpha }_{76}\right)=12$ (172)

$32{\alpha }_{32}^2+12S_4^2+16S_5^2+16S_6^2+7S_7^2=72$ (173)

$2{\alpha }_{32}+3S_4+9S_5+9S_6+7S_7=36$ (174)

$32{\alpha }_{65}{\alpha }_{54}{\alpha }_{43}{\alpha }_{32}+14\left\{{\alpha }_{75}{\alpha }_{54}{\alpha }_{43}{\alpha }_{32}+{\alpha }_{76}\left[{\alpha }_{64}{\alpha }_{43}{\alpha }_{32}+{\alpha }_{65}\left({\alpha }_{53}{\alpha }_{32}+{\alpha }_{54}{S}_4\right)\right]\right\}=1$

(175)

$7{\alpha }_{76}{\alpha }_{65}2{\alpha }_{54}=1$ (176)

$16{\alpha }_{65}{\alpha }_{54}{S}_4+7\left[{\alpha }_{75}{\alpha }_{54}{S}_4+{\alpha }_{76}\left({\alpha }_{64}{S}_4+2{\alpha }_{65}{S}_5\right)\right]=3$ (177)

$7{\alpha }_{76}{\alpha }_{65}=1$ (178)

$64{\alpha }_{65}{\alpha }_{53}{\alpha }_{32}+7{\alpha }_{75}{\alpha }_{53}{\alpha }_{32}+7\left[{\alpha }_{76}\left({\alpha }_{63}{\alpha }_{32}-{\alpha }_{65}{S}_5\right)\right]=-1$ (179)

$16\left(6{\alpha }_{65}\right){\alpha }_{54}+7\left(6{\alpha }_{75}\right){\alpha }_{54}+7\left(6{\alpha }_{76}\right)\left({\alpha }_{64}+3{\alpha }_{65}\right)=27$ (180)

$\begin{array}{l}12{\alpha }_{43}{\alpha }_{32}^2+16\left({\alpha }_{53}{\alpha }_{32}^2+{\alpha }_{54}{S}_4^2\right)+16\left({\alpha }_{63}{\alpha }_{32}^2+{\alpha }_{64}{S}_4^2+{\alpha }_{65}{S}_5^2\right)\\ +7\left({\alpha }_{73}{\alpha }_{32}^2+{\alpha }_{74}{S}_4^2+{\alpha }_{75}{S}_5^2+{\alpha }_{76}{S}_6^2\right)=12\end{array}$ (181)

$16{\alpha }_{54}{S}_4+16{\alpha }_{64}{S}_4+7{\alpha }_{74}{S}_4=3$ (182)

$7\left[{\alpha }_{74}{\alpha }_{43}{\alpha }_{32}+{\alpha }_{75}\left({\alpha }_{53}{\alpha }_{32}+{\alpha }_{54}{S}_4\right)+{\alpha }_{76}\left({\alpha }_{63}{\alpha }_{32}+{\alpha }_{64}{S}_4+{\alpha }_{65}{S}_5\right)\right]=1$ (183)

$8{\alpha }_{65}{\alpha }_{54}=1$ (184)

$7\left({\alpha }_{74}{S}_4+2{\alpha }_{75}{S}_5+2{\alpha }_{76}{S}_6\right)=9$ (185)

$64{\alpha }_{65}+7{\alpha }_{73}=20$ (186)

$\begin{array}{l}12S_4{\alpha }_{43}{\alpha }_{32}+16S_5\left({\alpha }_{53}{\alpha }_{32}+{\alpha }_{54}{S}_4\right)+16S_6\left({\alpha }_{63}{\alpha }_{32}+{\alpha }_{64}{S}_4+{\alpha }_{65}{S}_5\right)\\ +7S_7\left({\alpha }_{73}{\alpha }_{32}+{\alpha }_{74}{S}_4+{\alpha }_{75}{S}_5+{\alpha }_{76}{S}_6\right)=20\end{array}$ (187)

$16S_5{\alpha }_{54}+16S_6\left({\alpha }_{64}+3{\alpha }_{65}\right)+7S_7\left({\alpha }_{74}+3{\alpha }_{75}+3{\alpha }_{76}\right)=42$ (188)

$6{\alpha }_{43}{\alpha }_{32}=1$ (189)

$8{\alpha }_{54}+8\left({\alpha }_{64}+3{\alpha }_{65}\right)=9$ (190)

$8{\alpha }_{32}^2+6S_4^2+12S_5^2+12S_6^2+7S_7^2=60$ (191)

$2{\alpha }_{32}+6S_4+27S_5+27S_6+28S_7=120$ (192)

From (163), (166), (174) and (192) result:

${\alpha }_{32}=\frac{1}{8}$ (193)

$S_4=\frac{1}{2}$ (194)

$S_7=2$ (195)

and

$S_5+S_6=\frac{9}{4}$ (196)

From both Equations (173) and (191) we find

$S_5^2+S_6^2=\frac{81}{32}$ (197)

From the system $S_5^2+S_6^2=\frac{81}{32}$

and $S_5+S_6=\frac{9}{4}$ result:

$S_5=S_6=\frac{9}{8}$ (198)

From (176) and (178):

${\alpha }_{54}=\frac{1}{2}$ (199)

From (184): ${\alpha }_{65}=\frac{1}{4}$ (200) from (178): ${\alpha }_{76}=\frac{4}{7}$ (201)

From (189): ${\alpha }_{43}=\frac{4}{3}$ (202) from (170): ${\alpha }_{75}=0$ (203)

From (172): ${\alpha }_{74}=0$ (204) from (182): ${\alpha }_{64}=-\frac{1}{8}$ (205) from (186): ${\alpha }_{73}=\frac{4}{7}$ (206)

From (171) and (189): ${\alpha }_{53}=0$ (207) and ${\alpha }_{63}=\frac{1}{4}$ (208)

The remaining equations are verified.

From $S_4=\frac{1}{2}$ : ${\alpha }_{42}=-\frac{5}{6}$ (209) from $S_5=\frac{9}{8}$ : ${\alpha }_{52}=\frac{1}{8}$ (210)

From $S_6=\frac{9}{8}$ : ${\alpha }_{62}=\frac{3}{8}$ (211) from $S_7=2$ : ${\alpha }_{72}=-\frac{2}{7}$ (212)

From the relations ${\displaystyle {\sum }_{j=1}^{i-1}{\alpha }_{ij}}=c_i$, $i=2,3,4,5,6,7$ result:

${\alpha }_{21}=\frac{1}{4}$ (213)

${\alpha }_{31}=\frac{1}{8}$ (214)

${\alpha }_{41}=0$ (215)

${\alpha }_{51}=\frac{1}{8}$ (216)

${\alpha }_{61}=0$ (217)

and

${\alpha }_{71}=\frac{1}{7}$ (218)

Therefore we get Table 4:

with

$K_1=hf\left(x_n,y_n\right)$ (219)

$K_2=hf\left(x_n+\frac{h}{4},y_n+\frac{K_1}{4}\right)$ (220)

$K_3=hf\left(x_n+\frac{h}{4},y_n+\frac{K_1+K_2}{8}\right)$ (221)

$K_4=hf\left(x_n+\frac{2h}{4},y_n+\frac{-5K_2+8K_3}{6}\right)$ (222)

$K_5=hf\left(x_n+\frac{3h}{4},y_n+\frac{K_1+K_2+4K_4}{8}\right)$ (223)

$K_6=hf\left(x_n+\frac{3h}{4},y_n+\frac{3K_2+2K_3-K_4+2K_5}{8}\right)$ (224)

$K_7=hf\left(x_n+h,y_n+\frac{K_1-2K_2+4K_3+4K_6}{7}\right)$ (225)

and

$y_{n+1}=y_n+\frac{7K_1+32K_3+12K_4+16K_5+16K_6+7K_7}{90}$ (226)

Working in a similar way we are presented three other choices:

2.3. 3^rd Choice: c₂ =c₄= 4/12,c₃= 8/12, c₅= 3/12, c₆= 9/12, c₇= 12/12 [5]

For 3^rd choice we get Table 5:

with

$K_1=hf\left(x_n,y_n\right)$ (227)

$K_2=hf\left(x_n+\frac{h}{3},y_n+\frac{K_1}{3}\right)$ (228)

$K_3=hf\left(x_n+\frac{2h}{3},y_n+\frac{2K_2}{3}\right)$ (229)

$K_4=hf\left(x_n+\frac{h}{3},y_n+\frac{K_1+4K_2-K_3}{12}\right)$ (230)

$K_5=hf\left(x_n+\frac{h}{4},y_n+\frac{10K_1-9K_2+15K_4}{64}\right)$ (231)

$K_6=hf\left(x_n+\frac{3h}{4},y_n+\frac{-3K_1-27K_2+27K_3-45K_4+96K_5}{64}\right)$ (232)

$K_7=hf\left(x_n+h,y_n+\frac{\frac{275K_1}{4}-81K_2-\frac{459K_3}{4}+528K_4-384K_5+128K_6}{145}\right)$ (233)

and

$y_{n+1}=y_n+\frac{145K_1+243K_3+243K_4+512K_5+512K_6+145K_7}{1800}$ (234)

2.4. 4^th Choice: c₂ =c₃= 1/5,c₄=2/5,c₅ = 3/5, c₆ = 4/5και c₇ = 5/5 = 1 [5]

For 4th choice we get Table 6:

with

$K_1=hf\left(x_n,y_n\right)$ (235)

$K_2=hf\left(x_n+\frac{h}{5},y_n+\frac{K_1}{5}\right)$ (236)

$K_3=hf\left(x_n+\frac{h}{5},y_n+\frac{K_1+K_2}{10}\right)$ (237)

$K_4=hf\left(x_n+\frac{2h}{5},y_n+\frac{-9K_2+25K_3}{40}\right)$ (238)

$K_5=hf\left(x_n+\frac{3h}{5},y_n+\frac{4K_1-17K_2+21K_3+16K_4}{40}\right)$ (239)

$K_6=hf\left(x_n+\frac{4h}{5},y_n+\frac{-8K_1+21K_2+35K_3-40K_4+40K_5}{60}\right)$ (240)

$K_7=hf\left(x_n+h,y_n+\frac{66K_1-5K_2-165K_3+240K_4-120K_5+60K_6}{76}\right)$ (241)

and

$y_{n+1}=y_n+\frac{19K_1+75K_3+50K_4+50K_5+75K_6+19K_7}{288}$ (242)

2.5. 5^th Choice: c₂ = c₄ = 1/5, c₃ = 2/5, c₅ = 3/5, c₆ = 4/5και c₇ = 5/5 = 1 [5]

For 5th choice we get Table 7:

with

$K_1=hf\left(x_n,y_n\right)$ (243)

$K_2=hf\left(x_n+\frac{h}{5},y_n+\frac{K_1}{5}\right)$ (244)

$K_3=hf\left(x_n+\frac{2h}{5},y_n+\frac{2K_2}{5}\right)$ (245)

$K_4=hf\left(x_n+\frac{h}{5},y_n+\frac{11K_1-4K_2+5K_3}{60}\right)$ (246)

$K_5=hf\left(x_n+\frac{3h}{5},y_n+\frac{-3K_1+4K_2+3K_3+8K_4}{20}\right)$ (247)

$K_6=hf\left(x_n+\frac{4h}{5},y_n+\frac{-12K_2-8K_3+22K_4+10K_5}{15}\right)$ (248)

$K_7=hf\left(x_n+h,y_n+\frac{51K_1+140K_2+225K_3-280K_4-120K_5+60K_6}{76}\right)$ (249)

and

$y_{n+1}=y_n+\frac{19K_1+50K_3+75K_4+50K_5+75K_6+19K_7}{288}$ (250)

3. Conclusions

In the R-K methods it is c₁ = 0 and we choose c_ν = 1, that is, in method (6,7) presented here, it is c₇ = 1. For methods up to 4th order the number of steps is equal to the order of the method. For the higher order method, the number of steps exceeds its order. In the 6th order R-K method the steps are 7 and 8 or more. We look for the method with the fewest steps because increasing the steps also increases the number of parameters to be calculated. Because the system from which the parameters will be calculated is not linear, some “arbitrary” values [5] must be given to parameters to solve the system. We have given 5 choices as examples, to the values of the arbitrary constants and the solutions of the respective systems that result are made analytically.

The process of solving the system starts by giving values to c_i to calculate the w_i. The resulting system is a linear 6 × 7 system. To solve it we give the same value in two c_i eg c₂ = c₃ or c₂ = c₄ etc. while the values ofc_i we choose to be in ascending order. As the variablec₂does not occur often we prefer one of the two equals c_i to be c₂. In the choice eg c₂ = c₃ we consider w₂ + w₃ as one parameter and after it is found that w₂ + w₃ = k we set w₂ = 0 and w₃ = k. In choices 2.1 and 2.2 the same values were given to two pairs of parameters c_i and we omitted one equation to make the system 5 × 5 and whose solution must verify the omitted equation.