From the research carried out in publications of works related to the creation of R-K methods for solving Ordinary Differential Equations it was observed that the description of the creation of these methods was done in a general way and for all classes of normal differential equations. The approach to create this method was not simple but was partly complex. So it was decided to create and propose a process for creating a R-K method that will be simple, understandable and applicable.

A system of ordinary differential equations of the form

y ′ = f ( x , y ) , y ( x 0 ) = y 0 (1)

with x 0 ∈ ℝ , y , y ′ ∈ ℝ m and f : ℝ × ℝ m → ℝ m , is called Initial Value Problem (IVP).

Runge-Kutta methods are commonly used numerical methods for addressing (1). They usually presented in a so-called Butcher table (Table 1) [2] [3]:

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 T , c ∈ ℝ m while A ∈ ℝ m × m . Then, the method shares m stages and in case that c₁ = 0 and A is strictly lower triangular, it is evaluated explicitly.

The solution of a differential equation is a continuous curve y(x) that passes through the point (x₀,y₀) and satisfies y ′ ( x ) = f ( x , y ) . Numerical solution of a differential equation is a distinct set of values of y(x) which is an approach to the continuous solution of the y(x) curve.

Carl David Tolmé Runge [4] and Martin Wilhelm Kutta [5] introduced the methods bearing their names almost in the turning of the 19th century. Runge and Kutta observed that the derivation of high-order derivatives that appear in the Taylor method can be avoided. In this method we place the problem with indeterminate parameters and make the result at the highest order using calculations of f(x,y) inside (x_n,y_n) and (x_n+1,y_n+1)intervals. The derivatives in the Taylor form are replaced by calculating f(x,y)at a number of points inside (x_n,y_n) and (x_n+1,y_n+1) intervals.

Runge was the first to present a 2nd order R-K method by combining a sequence of Euler formulas [4]. Some years later, Kutta managed to construct a 4 stages 4th order method [5]. Nyström showed a method (5,6) of 5^th order and 6 stages [6]. Fehlberg [7], Shanks [8] and Lawson [9] showed 5^th order methods of 6 stages too. 6^th order methods have been presented by Butcher [2], Fehlberg [7], Shanks [8] and Lawson [9]. Huta’s 6^th order method of 8 stages is the most popular [10]. Higher order R-K methods have been presented by Shanks [8], Felhberg [7], Feagin [11] [12] Hairer [13], Butcher [14] [15], Curtis [16], Famelis [17], Papakostas [17] [18], Tsitouras [17] [18] [19] and others.

Some problems that could be solved in this paper:

· We want with analytic way to derive the RK(7,9) method and we introduce our method for that.

· We give first arbitrary variables with values of the existing table of RK(7,9) method in order to compare our method of solving the non-linear system.

· We suggest and some others arbitrary variables which lead to desired Tables of the RK(7,9) method, because not all the arbitrary values lead to desired results of the method.

Firstly, we present the Introduction of the paper with historical references, and then in section 2 we give an analytic approach of the Runge-Kutta method (7,9) 7^th order with 9 stages method. Finally, we give the conclusions of our work of a certain set of values of the parameters of the method.

The reduction formula of R-K methods for an ordinary differential equation is given by the relation y n + 1 = y n + ∑ i = 1 ν w i K i (2), with w_i acting as coefficients of weight, ν the number of steps and K i = h f ( x n + c i h , y n + ∑ j = 1 i − 1 α i j K j ) ,   c 1 = 0 ,   i = 2 , 3 , ⋯ , ν (3) withh the step of the method. The parameters w_i, c_i and α_ij must be specified. In every R-K method the relations ∑ i = 1 ν w i = 1 (4) και c i = ∑ j = 1 i − 1 α i j ,   i = 2 , 3 , ⋯ , ν (5) must be valid.

Runge-Kutta (7,9) method is a method of 7^th order and 9 stages and we use the coefficients obtained for r = 1 , 2 , ⋯ , 7 (6), where r is the order of ODE, from Butcher’s Table[1] from whom the equations of the nonlinear 85 × 53 system result.

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 differential equation.

The equations of the system are numbered from (8), (9), ···, (92) and introducing the abbreviation:

P κ λ = α κ 2 c 2 λ + α κ 3 c 3 λ + ⋯ + α κ κ − 1 c κ − 1 λ     with   κ = 3 , 4 , ⋯ , 9     and     λ = 1 , 2 , 3 , 4 , 5 (7)

[7] the following system is obtained:

∑ κ = 1 9 w κ = 1 (8)

∑ κ = 2 9 w κ c κ = 1 2 (9)

∑ κ = 3 9 w κ P κ 1 = 1 6 (10)

∑ κ = 2 9 w κ c κ 2 = 1 3 (11)

∑ κ = 4 9 w κ ( ∑ λ = 3 κ − 1 α κ λ P λ 1 ) = 1 24 (12)

∑ κ = 3 9 w κ P κ 2 = 1 12 (13)

∑ κ = 3 9 w κ c κ P κ 1 = 1 8 (14)

∑ κ = 2 9 w κ c κ 3 = 1 4 (15)

∑ κ = 5 9 w κ [ ∑ λ = 4 κ − 1 α κ λ ( ∑ μ = 3 λ − 1 α λ μ P μ 1 ) ] = 1 120 (16)

∑ κ = 4 9 w κ ( ∑ λ = 3 κ − 1 α κ λ P λ 2 ) = 1 60 (17)

∑ κ = 4 9 w κ ( ∑ λ = 3 κ − 1 α κ λ c λ P λ 1 ) = 1 40 (18)

∑ κ = 3 9 w κ P κ 3 = 1 20 (19)

∑ κ = 4 9 w κ c κ ( ∑ λ = 3 κ − 1 α κ λ P λ 1 ) = 1 30 (20)

∑ κ = 3 9 w κ c κ P κ 2 = 1 15 (21)

∑ κ = 3 9 w κ P κ 1 2 = 1 20 (22)

∑ κ = 3 9 w κ c κ 2 P κ 1 = 1 10 (23)

∑ κ = 2 9 w κ c κ 4 = 1 5 (24)

∑ κ = 6 9 w κ { ∑ λ = 5 κ − 1 α κ λ [ ∑ μ = 4 λ − 1 α λ μ ( ∑ ν = 3 μ − 1 α μ ν P ν 1 ) ] } = 1 720 (25)

∑ κ = 5 9 w κ [ ∑ λ = 4 κ − 1 α κ λ ( ∑ μ = 3 λ − 1 α λ μ P μ 2 ) ] = 1 360 (26)

∑ κ = 5 9 w κ [ ∑ λ = 4 κ − 1 α κ λ ( ∑ μ = 3 λ − 1 α λ μ c μ P μ 1 ) ] = 1 240 (27)

∑ κ = 4 9 w κ ( ∑ λ = 3 κ − 1 α κ λ P λ 1 ) = 1 120 (28)

∑ κ = 5 9 w κ [ ∑ λ = 4 κ − 1 α κ λ c λ ( ∑ μ = 3 λ − 1 α λ μ P μ 1 ) ] = 1 180 (29)

∑ κ = 4 9 w κ ( ∑ λ = 3 κ − 1 α κ λ c λ P λ 2 ) = 1 90 (30)

∑ κ = 4 9 w κ ( ∑ λ = 3 κ − 1 α κ λ P λ 1 2 ) = 1 120 (31)

∑ κ = 4 9 w κ ( ∑ λ = 3 κ − 1 α κ λ c λ 2 P λ 1 ) = 1 60 (32)

∑ κ = 3 9 w κ P κ 4 = 1 30 (33)

∑ κ = 5 9 w κ c κ [ ∑ λ = 4 κ − 1 α κ λ ( ∑ μ = 3 λ − 1 α λ μ P μ 1 ) ] = 1 144 (34)

∑ κ = 4 9 w κ c κ ( ∑ λ = 3 κ − 1 α κ λ P λ 2 ) = 1 72 (35)

∑ κ = 4 9 w κ c κ ( ∑ λ = 3 κ − 1 α κ λ c λ P λ 1 ) = 1 48 (36)

∑ κ = 3 9 w κ c κ P κ 3 = 1 24 (37)

∑ κ = 4 9 w κ P κ 1 ( ∑ λ = 3 κ − 1 α κ λ P λ 1 ) = 1 72 (38)

∑ κ = 3 9 w κ P κ 1 P κ 2 = 1 36 (39)

∑ κ = 4 9 w κ c κ 2 ( ∑ λ = 3 κ − 1 α κ λ P λ 1 ) = 1 36 (40)

∑ κ = 3 9 w κ c κ 2 P κ 2 = 1 18 (41)

∑ κ = 3 9 w κ c κ P κ 1 2 = 1 24 (42)

∑ κ = 3 9 w κ c κ 3 P κ 1 = 1 12 (43)

∑ κ = 2 9 w κ c κ 5 = 1 6 (44)

∑ κ = 7 9 w κ 〈 ∑ λ = 6 κ − 1 α κ λ { ∑ μ = 5 λ − 1 α λ μ [ ∑ ν = 4 μ − 1 α μ ν ( ∑ ρ = 3 ν − 1 α ν ρ P ρ 1 ) ] } 〉 = 1 5040 (45)

∑ κ = 6 9 w κ { ∑ λ = 5 κ − 1 α κ λ [ ∑ μ = 4 λ − 1 α λ μ ( ∑ ν = 3 μ − 1 α μ ν P ν 2 ) ] } = 1 2520 (46)

∑ κ = 6 9 w κ { ∑ λ = 5 κ − 1 α κ λ [ ∑ μ = 4 λ − 1 α λ μ ( ∑ ν = 3 μ − 1 α μ ν c ν P ν 1 ) ] } = 1 1680 (47)

∑ κ = 5 9 w κ [ ∑ λ = 4 κ − 1 α κ λ ( ∑ μ = 3 λ − 1 α λ μ P μ 3 ) ] = 1 840 (48)

∑ κ = 6 9 w κ { ∑ λ = 5 κ − 1 α κ λ [ ∑ μ = 4 λ − 1 α λ μ c μ ( ∑ ν = 3 μ − 1 α μ ν P ν 1 ) ] } = 1 1260 (49)

∑ κ = 5 9 w κ [ ∑ λ = 4 κ − 1 α κ λ ( ∑ μ = 3 λ − 1 α λ μ c μ P μ 2 ) ] = 1 630 (50)

∑ κ = 5 9 w κ [ ∑ λ = 4 κ − 1 α κ λ ( ∑ μ = 3 λ − 1 α λ μ P μ 1 2 ) ] = 1 840 (51)

∑ κ = 5 9 w κ [ ∑ λ = 4 κ − 1 α κ λ ( ∑ μ = 3 λ − 1 α λ μ c μ 2 P μ 1 ) ] = 1 420 (52)

∑ κ = 4 9 w κ ( ∑ λ = 3 κ − 1 α κ λ P λ 4 ) = 1 210 (53)

∑ κ = 6 9 w κ { ∑ λ = 5 κ − 1 α κ λ c λ [ ∑ μ = 4 λ − 1 α λ μ ( ∑ ν = 3 μ − 1 α μ ν P ν 1 ) ] } = 1 1008 (54)

∑ κ = 5 9 w κ [ ∑ λ = 4 κ − 1 α κ λ c λ ( ∑ μ = 3 λ − 1 α λ μ P μ 2 ) ] = 1 504 (55)

∑ κ = 5 9 w κ [ ∑ λ = 4 κ − 1 α κ λ c λ ( ∑ μ = 3 λ − 1 α λ μ c μ P μ 1 ) ] = 1 336 (56)

∑ κ = 4 9 w κ ( ∑ λ = 3 κ − 1 α κ λ c λ P λ 3 ) = 1 168 (57)

∑ κ = 5 9 w κ [ ∑ λ = 4 κ − 1 α κ λ P λ 1 ( ∑ μ = 3 λ − 1 α λ μ P μ 1 ) ] = 1 504 (58)

∑ κ = 4 9 w κ ( ∑ λ = 3 κ − 1 α κ λ P λ 1 P λ 3 ) = 1 252 (59)

∑ κ = 5 9 w κ [ ∑ λ = 4 κ − 1 α κ λ c λ 2 ( ∑ μ = 3 λ − 1 α λ μ P μ 1 ) ] = 1 252 (60)

∑ κ = 4 9 w κ ( ∑ λ = 3 κ − 1 α κ λ c λ 2 P λ 2 ) = 1 126 (61)

∑ κ = 4 9 w κ ( ∑ λ = 3 κ − 1 α κ λ c λ P λ 1 2 ) = 1 168 (62)

∑ κ = 4 9 w κ ( ∑ λ = 3 κ − 1 α κ λ c λ 3 P λ 1 ) = 1 84 (63)

∑ κ = 3 9 w κ P κ 5 = 1 42 (64)

∑ κ = 6 9 w κ c κ { ∑ λ = 5 κ − 1 α κ λ [ ∑ μ = 4 λ − 1 α λ μ ( ∑ ν = 3 μ − 1 α μ ν P ν 1 ) ] } = 1 840 (65)

∑ κ = 5 9 w κ c κ [ ∑ λ = 4 κ − 1 α κ λ ( ∑ μ = 3 λ − 1 α λ μ P μ 2 ) ] = 1 420 (66)

∑ κ = 5 9 w κ c κ [ ∑ λ = 4 κ − 1 α κ λ ( ∑ μ = 3 λ − 1 α λ μ c μ P μ 1 ) ] = 1 280 (67)

∑ κ = 4 9 w κ c κ ( ∑ λ = 3 κ − 1 α κ λ P λ 3 ) = 1 140 (68)

∑ κ = 5 9 w κ c κ [ ∑ λ = 4 κ − 1 α κ λ c λ ( ∑ μ = 3 λ − 1 α λ μ P μ 1 ) ] = 1 210 (69)

∑ κ = 4 9 w κ c κ ( ∑ λ = 3 κ − 1 α κ λ c λ P λ 2 ) = 1 105 (70)

∑ κ = 4 9 w κ c κ ( ∑ λ = 3 κ − 1 α κ λ P λ 1 2 ) = 1 140 (71)

∑ κ = 4 9 w κ c κ ( ∑ λ = 3 κ − 1 α κ λ c λ 2 P λ 1 ) = 1 70 (72)

∑ κ = 3 9 w κ c κ P κ 4 = 1 35 (73)

∑ κ = 5 9 w κ P κ 1 [ ∑ λ = 4 κ − 1 α κ λ ( ∑ μ = 3 λ − 1 α λ μ P μ 1 ) ] = 1 336 (74)

∑ κ = 4 9 w κ P κ 1 ( ∑ λ = 3 κ − 1 α κ λ P λ 2 ) = 1 168 (75)

∑ κ = 4 9 w κ P κ 1 ( ∑ λ = 3 κ − 1 α κ λ c λ P λ 1 ) = 1 112 (76)

∑ κ = 3 9 w κ P κ 1 P κ 3 = 1 56 (77)

∑ κ = 5 9 w κ c κ 2 [ ∑ λ = 4 κ − 1 α κ λ ( ∑ μ = 3 λ − 1 α λ μ P μ 1 ) ] = 1 168 (78)

∑ κ = 4 9 w κ c k 2 ( ∑ λ = 3 κ − 1 α κ λ P λ 2 ) = 1 84 (79)

∑ κ = 4 9 w κ c k 2 ( ∑ λ = 3 κ − 1 α κ λ c λ P λ 1 ) = 1 56 (80)

∑ κ = 3 9 w κ c κ 2 P κ 3 = 1 28 (81)

∑ κ = 4 9 w κ ( ∑ λ = 3 κ − 1 α κ λ P λ 1 ) 2 = 1 252 (82)

∑ κ = 4 9 w κ P κ 2 ( ∑ λ = 3 κ − 1 α κ λ P λ 1 ) = 1 126 (83)

∑ κ = 3 9 w κ P κ 2 2 = 1 63 (84)

∑ κ = 4 9 w κ c κ P κ 1 ( ∑ λ = 3 κ − 1 α κ λ P λ 1 ) = 1 84 (85)

∑ κ = 3 9 w κ c κ P κ 1 P κ 2 = 1 42 (86)

∑ κ = 4 9 w κ c κ 3 ( ∑ λ = 3 κ − 1 α κ λ P λ 1 ) = 1 42 (87)

∑ κ = 3 9 w κ c κ 3 P κ 2 = 1 21 (88)

∑ κ = 3 9 w κ P κ 1 3 = 1 56 (89)

∑ κ = 3 9 w κ c κ 2 P κ 1 2 = 1 28 (90)

∑ κ = 3 9 w κ c κ 4 P κ 1 = 1 14 (91)

∑ κ = 2 9 w κ c κ 6 = 1 7 (92)

In the system of (8), (9), (11), (15), (24), (44) and (92) equations we set as c₂ = c₃ = 1/12, c₄ = 1/6, c₅ = 2/6, c₆ = 3/6,c₇ = 4/6, c₈ = 5/6, c₉ = 6/6. Τhe values of c₂, c₃, ··· , c₉ are chosen to be in ascending order and as small and different from each other as possible. We set in addition w 2 = w 3 = 0 (93) and the resulting solution is:

w 1 = w 9 = 41 840 (94)

w 4 = w 8 = 216 840 (95)

w 5 = w 7 = 27 840 (96)

and

w 6 = 272 840 (97)

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

P 41 = α 42 c 2 + α 43 c 3 = α 42 + α 43 12 = S 4 12 (98)

P 51 = α 52 c 2 + α 53 c 3 + α 54 c 4 = α 52 + α 53 + 2   α 54 12 = S 5 12 (99)

P 61 = α 62 + α 63 + 2   α 44 + 4   α 65 12 = S 6 12 (100)

P 71 = α 72 + α 73 + 2   α 74 + 4   α 75 + 6   α 76 12 = S 7 12 (101)

P 81 = α 82 + α 83 + 2   α 84 + 4   α 85 + 6   α 86 + 8   α 87 12 = S 8 12 (102)

P 91 = α 92 + α 93 + 2   α 94 + 4   α 95 + 6   α 96 + 8   α 97 + 10   α 98 12 = S 9 12 (103)

P 42 = α 42 c 2 2 + α 43 c 3 2 = α 42 + α 43 144 = S 4 144 (104)

P 52 = α 52 c 2 2 + α 53 c 3 2 + α 54 c 4 2 = α 52 + α 53 + 4   α 54 144 = S 5 + 2   α 54 144 (105)

P 62 = α 62 c 2 2 + α 63 c 3 2 + α 64 c 4 2 + α 65 c 5 2 = S 6 + 2   α 64 + 12   α 65 144 (106)

P 72 = α 72 c 2 2 + α 73 c 3 2 + α 74 c 4 2 + α 75 c 5 2 + α 76 c 6 2 = S 7 + 2   α 74 + 12   α 75 + 30   a 76 144 (107)

P 82 = S 8 + 2   α 84 + 12   α 85 + 30   a 86 + 56   a 87 144 (108)

P 92 = S 9 + 2   α 94 + 12   α 95 + 30   a 96 + 56   a 97 + 90   a 98 144 (109)

P 43 = a 42 + a 43 1728 = S 4 1728 (110)

P 53 = a 52 + a 53 + 8   a 54 1728 = S 5 + 6   α 54 1728 (111)

P 63 = S 6 + 6   α 64 + 60   α 65 1728 (112)

P 73 = S 7 + 6   α 74 + 60   α 75 + 210   a 76 1728 (113)

P 83 = S 8 + 6   α 84 + 60   α 85 + 210   a 86 + 504   a 87 1728 (114)

P 93 = S 9 + 6   α 94 + 60   α 95 + 210   a 96 + 504   a 97 + 990   a 98 1728 (115)

P 44 = S 4 20736 (116)

P 54 = S 5 + 14   α 54 20736 (117)

P 64 = S 6 + 14   a 64 + 252   α 65 20736 (118)

P 74 = S 7 + 14   a 74 + 252   α 75 + 1290   a 76 20736 (119)

P 84 = S 8 + 14   a 84 + 252   α 85 + 1290   a 86 + 4088   α 87 20736 (120)

P 94 = S 9 + 14   a 94 + 252   α 95 + 1290   a 96 + 4088   α 97 + 9990   α 98 20736 (121)

P 45 = S 4 248832 (122)

P 55 = S 5 + 30   α 54 248832 (123)

P 65 = S 6 + 30   α 64 + 1020   α 65 248832 (124)

P 75 = S 7 + 30   α 74 + 1020   α 75 + 7770   α 76 248832 (125)

P 85 = S 8 + 30   α 84 + 1020   α 85 + 7770   α 86 + 32760   α 87 248832 (126)

P 95 = S 9 + 30   α 94 + 1020   α 95 + 7770   α 96 + 32760   α 97 + 99990   α 98 248832 (127)

Then we substitute the defined abbreviations in the original system, as well as the found values of c₂, c₃, c₄, c₅, c₆, c₇, c₈, c₉, w₂, w₃, w₄, w₅, w₆, w₇, w₈, w₉, and as a result the system is simplified.

In the system of (10), (14), (23), (43) and (91) we express S₅, S₆, S₇, S₈ and S₉ as a function of S₄ and by substituting them in (22) we find that: S 4 = 1 6 (128) and S 5 = 4 6 (129), S 6 = 9 6 (130), S 7 = 16 6 (131), S 8 = 25 6 (132), S 9 = 36 6 = 6 (133).

To continue we set α 42 = 0 (134). From the abbreviation S 4 = α 42 + α 43 and from the relation α 41 + α 42 + α 43 = c 4 we obtain that α 43 = 2 12 (135) and α 41 = 0 (136).

In the system of equations (13), (18), (21), (32), (33), (36), (37), (39), (41), (62), (63), (64), (72), (73), (76), (77), (80), (81), (86) and (88) we substitute the values found above, omitting the equations which are a linear combination of equations of the system and also considering α₉₄, α₉₅, α₉₆, α₉₇, α₉₈ as parameters, the 10 × 15 linear system (A) is obtained:

27   α 54 = − 822 + 41 ( α 94 + 6   α 95 + 15   α 96 + 28   α 97 + 45   α 98 ) (137)

272   α 64 = 102 − 41 ( 4   α 94 + 18   α 95 + 30   α 96 + 28   α 97 ) (138)

272   α 65 = 653 − 41 ( α 95 + 5   α 96 + 14   α 97 + 30   α 98 ) (139)

27   α 74 = 102 + 41 ( 6   α 94 + 18   α 95 + 15   α 96 ) (140)

27   α 75 = − 208 + 41 ( 3   α 95 + 10   α 96 + 14   α 97 ) SYSTEM (A) (141)

270   α 76 = − 2560 + 41 ( 10   α 96 + 56   α 97 + 180   α 98 ) (142)

216   α 84 = − 822 − 41 ( 4   α 94 + 6   α 95 ) (143)

216   α 85 = 653 − 41 ( 3   α 95 + 5   α 96 ) (144)

216 ( 5   α 86 ) = − 1280 − 41 ( 10   α 96 + 28   α 97 ) (145)

216 ( 7   α 87 ) = 1683 − 41 ( 7   α 97 + 45   α 98 ) (146)

From (20) and (85) equations we obtain:

27 ( α 54 ) + 272 [ 4 ( α 64 + 4   α 65 ) ] + 27 [ 10 ( α 74 + 4   α 75 + 9   α 76 ) ]   + 216 [ 20   ( α 84 + 4   α 85 + 9   α 86 + 16   α 87 ) ]   + 41 [ 35 ( α 94 + 4   α 95 + 9   α 96 + 16   α 97 + 25   α 98 ) ] = 49824 (147)

which along with the equations of system (A) and after setting: α 65 = α 75 = α 84 = α 86 = 0 (148) results that:

α 54 = 160 9 (149)

α 64 = − 297 68 (150)

α 65 = 0 (151)

α 74 = 157 3 (152)

α 75 = 0 (153)

α 76 = 157 45 (154)

α 84 = 0 (155)

α 85 = − 65 72 (156)

α 86 = 0 (157)

α 87 = 29 90 (158)

α 94 = − 1211 326 (159)

α 95 = − 205 237 (160)

α 96 = 419 90 (161)

α 96 = − 25 9 (162)

α 98 = 27 25 (163)

We found above α₄₃ = 1/6 and setting α 32 = 11 12 (164) ( α 31 = − 10 12 (165)), from the system of (30), (34), (58), (60) equations, implies that: α 53 = 4 9 (166) α 63 = 18 5 (167) α 73 = − 23 2 (168) α 83 = − 3 14 (169) and from Equation (12) result α 93 = − 119 6 (170).

From the abbreviations: S₅ = 4/6, S₆ = 9/6, S₇ = 16/6, S₈ = 25/6 and S₉ = 36/6 = 6 results that: α 52 = − 106 3 (171) α 62 = 199 30 (172) α 72 = − 683 6 (173) α 82 = 379 70 (174) and α 92 = 283 14 (175).

From relations ∑ j = 1 i − 1 α i j = c i , i = 2 , 3 , 4 , 5 , 6 , 7 , 8 , 9 we obtain: α 51 = 157 9 (176) α 61 = − 161 30 (177) α 71 = 3158 45 (178) α 81 = − 53 14 (179) and α 91 = 56 25 (180).

According to the so-called Butcher’s Table(Table2) the (7,9) R-K method is given as below:

therefore

K 1 = h f ( x n , y n ) (181)

K 2 = h f ( x n + h 12 , y n + K 1 12 ) (182)

K 3 = h f ( x n + h 12 , y n + − 10 K 1 + 11 K 2 12 ) (183)

K 4 = h f ( x n + 2 h 12 , y n + 2 K 3 12 ) (184)

K 5 = h f ( x n + 4 h 12 , y n + 157 K 1 − 318 K 2 + 4 K 3 + 160 K 4 9 ) (185)

K 6 = h f ( x n + 6 h 12 , y n + − 322 K 1 + 199 K 2 + 108 K 3 − 131 K 5 30 ) (186)

K 7 = h f ( x n + 8 h 12 , y n + 3158 K 1 45 − 638 K 2 6 − 23 K 3 2 + 157 K 4 3 + 157 K 6 45 ) (187)

K 8 = h f ( x n + 10 h 12 , y n − 53 K 1 14 + 38 K 2 7 − 3 K 3 14 − 65 K 5 72 + 29 K 7 90 ) (188)

K 9 = h f ( x n + h , y n + 56 K 1 25 + 283 k 2 14 − 119 K 3 6 − 26 K 4 7 − 13 K 5 15                 + 149 K 6 32 − 25 K 7 9 + 27 K 8 25 ) (189)

and the reduction formula for the solution of the Differential Equation is:

y n + 1 = y n + 41 K 1 + 216 K 4 + 27 K 5 + 272 K 6 + 27 K 7 + 216 K 8 + 41 K 9 840 (190)

This paper is concerned with training the coefficients of a 7^th order and 9 stages Runge-Kutta method for addressing initial value problems. As the presented method is 9 stages, we use a set of 9 free parameters. After optimizing the free parameters (coefficients), we concluded to a certain set of values of them. This set of values was found to outperform other representatives in a wide range of relevant problems.

The authors declare no conflicts of interest regarding the publication of this paper.

Trikkaliotis, G.D. and Gousidou-Koutita, M.Ch. (2022) Production of the Reduction Formula of Seventh Order Runge-Kutta Method with Step Size Control of an Ordinary Differential Equation. Applied Mathematics, 13, 325-337. https://doi.org/10.4236/am.2022.134023