scholarly journals Runge–Kutta Pairs of Orders 5(4) Trained to Best Address Keplerian Type Orbits

Mathematics ◽  
2021 ◽  
Vol 9 (19) ◽  
pp. 2400
Author(s):  
Vladislav N. Kovalnogov ◽  
Ruslan V. Fedorov ◽  
Tamara V. Karpukhina ◽  
Theodore E. Simos ◽  
Charalampos Tsitouras

The derivation of Runge–Kutta pairs of orders five and four that effectively uses six stages per step is considered. The coefficients provided by such a method are 27 and have to satisfy a system of 25 nonlinear equations. Traditionally, various solutions have been tried. Each of these solutions makes use of some simplified assumptions and offers different families of methods. Here, we make use of the most celebrated family to appear in the literature, where we may use as the last stage the first function evaluation from the next step (FSAL property). The family under consideration has the advantage of being solved explicitly. Actually, we arrive at a subsystem where all the coefficients are found with respect to five free parameters. These free parameters are adjusted (trained) in order to deliver a pair that outperforms other similar pairs of orders 5(4) in Keplerian type orbits, e.g., Kepler, perturbed Kepler, Arenstorf orbit or Pleiades. The training uses differential evolution technique. The finally proposed pair has a remarkable performance and offers on average more than a digit of accuracy in a variety of orbits.

Mathematics ◽  
2021 ◽  
Vol 9 (23) ◽  
pp. 3071
Author(s):  
Vladislav N. Kovalnogov ◽  
Ruslan V. Fedorov ◽  
Andrey V. Chukalin ◽  
Theodore E. Simos ◽  
Charalampos Tsitouras

The family of Numerov-type methods that effectively uses seven stages per step is considered. All the coefficients of the methods belonging to this family can be expressed analytically with respect to four free parameters. These coefficients are trained through a differential evolution technique in order to perform best in a wide range of Keplerian-type orbits. Then it is observed with extended numerical tests that a certain method behaves extremely well in a variety of orbits (e.g., Kepler, perturbed Kepler, Arenstorf, Pleiades) for various steplengths used by the methods and for various intervals of integration.


Mathematics ◽  
2021 ◽  
Vol 9 (18) ◽  
pp. 2306
Author(s):  
Vladislav N. Kovalnogov ◽  
Ruslan V. Fedorov ◽  
Andrey V. Chukalin ◽  
Theodore E. Simos ◽  
Charalampos Tsitouras

The purpose of the present work is to construct a new Runge–Kutta pair of orders five and four to outperform the state-of-the-art in these kind of methods when addressing problems with periodic solutions. We consider the family of such pairs that the celebrated Dormand–Prince pair also belongs. The chosen family comes with coefficients that all depend on five free parameters. These latter parameters are tuned in a way to furnish a new method that performs best on a couple of oscillators. Then, we observe that this trained pair outperforms other well known methods in the relevant literature in a standard set of problems with periodic solutions. This is remarkable since no special property holds such as high phase-lag order or an extended interval of periodicity.


Mathematics ◽  
2021 ◽  
Vol 9 (16) ◽  
pp. 1842
Author(s):  
Vladislav N. Kovalnogov ◽  
Ruslan V. Fedorov ◽  
Yuri A. Khakhalev ◽  
Theodore E. Simos ◽  
Charalampos Tsitouras

We consider the scalar autonomous initial value problem as solved by an explicit Runge-Kutta pair of orders 6 and 5. We focus on an efficient family of such pairs, which were studied extensively in previous decades. This family comes with 5 coefficients that one is able to select arbitrarily. We set, as a fitness function, a certain measure, which is evaluated after running the pair in a couple of relevant problems. Thus, we may adjust the coefficients of the pair, minimizing this fitness function using the differential evolution technique. We conclude with a method (i.e. a Runge-Kutta pair) which outperforms other pairs of the same two orders in a variety of scalar autonomous problems.


Algorithms ◽  
2021 ◽  
Vol 14 (4) ◽  
pp. 101
Author(s):  
Alicia Cordero ◽  
Marlon Moscoso-Martínez ◽  
Juan R. Torregrosa

In this paper, we present a new parametric family of three-step iterative for solving nonlinear equations. First, we design a fourth-order triparametric family that, by holding only one of its parameters, we get to accelerate its convergence and finally obtain a sixth-order uniparametric family. With this last family, we study its convergence, its complex dynamics (stability), and its numerical behavior. The parameter spaces and dynamical planes are presented showing the complexity of the family. From the parameter spaces, we have been able to determine different members of the family that have bad convergence properties, as attracting periodic orbits and attracting strange fixed points appear in their dynamical planes. Moreover, this same study has allowed us to detect family members with especially stable behavior and suitable for solving practical problems. Several numerical tests are performed to illustrate the efficiency and stability of the presented family.


Author(s):  
Janak Raj Sharma ◽  
Sunil Kumar ◽  
Ioannis K. Argyros

In this paper, a class of efficient iterative methods with increasing order of convergence for solving systems of nonlinear equations is developed and analyzed. The methodology uses well-known third-order Potra–Pták iteration in the first step and Newton-like iterations in the subsequent steps. Novelty of the methods is the increase in convergence order by an amount three per step at the cost of only one additional function evaluation. In addition, the algorithm uses a single inverse operator in each iteration, which makes it computationally more efficient and attractive. Local convergence is studied in the more general setting of a Banach space under suitable assumptions. Theoretical results of convergence and computational efficiency are verified through numerical experimentation. Comparison of numerical results indicates that the developed algorithms outperform the other similar algorithms available in the literature, particularly when applied to solve the large systems of equations. The basins of attraction of some of the existing methods along with the proposed method are given to exhibit their performance.


2019 ◽  
Vol 17 (1) ◽  
pp. 1567-1598
Author(s):  
Tianbao Liu ◽  
Xiwen Qin ◽  
Qiuyue Li

Abstract In this paper, we derive and analyze a new one-parameter family of modified Cauchy method free from second derivative for obtaining simple roots of nonlinear equations by using Padé approximant. The convergence analysis of the family is also considered, and the methods have convergence order three. Based on the family of third-order method, in order to increase the order of the convergence, a new optimal fourth-order family of modified Cauchy methods is obtained by using weight function. We also perform some numerical tests and the comparison with existing optimal fourth-order methods to show the high computational efficiency of the proposed scheme, which confirm our theoretical results. The basins of attraction of this optimal fourth-order family and existing fourth-order methods are presented and compared to illustrate some elements of the proposed family have equal or better stable behavior in many aspects. Furthermore, from the fractal graphics, with the increase of the value m of the series in iterative methods, the chaotic behaviors of the methods become more and more complex, which also reflected in some existing fourth-order methods.


2018 ◽  
Vol 24 (1) ◽  
pp. 3
Author(s):  
Himani Arora ◽  
Juan Torregrosa ◽  
Alicia Cordero

In this study, an iterative scheme of sixth order of convergence for solving systems of nonlinear equations is presented. The scheme is composed of three steps, of which the first two steps are that of third order Potra-Pták method and last is weighted-Newton step. Furthermore, we generalize our work to derive a family of multi-step iterative methods with order of convergence 3 r + 6 , r = 0 , 1 , 2 , … . The sixth order method is the special case of this multi-step scheme for r = 0 . The family gives a four-step ninth order method for r = 1 . As much higher order methods are not used in practice, so we study sixth and ninth order methods in detail. Numerical examples are included to confirm theoretical results and to compare the methods with some existing ones. Different numerical tests, containing academical functions and systems resulting from the discretization of boundary problems, are introduced to show the efficiency and reliability of the proposed methods.


2000 ◽  
Vol 4 (2) ◽  
pp. 143-150 ◽  
Author(s):  
Lin-Yi Chou ◽  
P. W. Sharp

Order five symplectic explicit Runge-Kutta Nyström methods of five stages are known to exist. However, these methods do not have free parameters with which to minimise the principal error coefficients. By adding one derivative evaluation per step, to give either a six-stage non-FSAL family or a seven-stage FSAL family of methods, two free parameters become available for the minimisation. This raises the possibility of improving the efficiency of order five methods despite the extra cost of taking a step.We perform a minimisation of the two families to obtain an optimal method and then compare its numerical performance with published methods of orders four to seven. These comparisons along with those based on the principal error coefficients show the new method is significantly more efficient than the five-stage, order five methods. The numerical comparisons also suggest the new methods can be more efficient than published methods of other orders.


Sign in / Sign up

Export Citation Format

Share Document