Two point block multistep methods with trigonometric−fitting for solving oscillatory problems

In this paper, we present the absolute stability of the existing 2-point implicit block multistep step methods of step number k = 3 and k = 5 and solving special second order ordinary differential equations (ODEs). The methods are then trigonometrically fitted so that they are suitable for solving highly oscillatory problems arising from the special second order ODEs. Their explicit counterparts are also trigonometrically fitted so that in the implementation the methods can act as a predictor-corrector pairs. The numerical results based on the integration over a large interval are given to show the performance of the proposed methods. From the numerical results we can conclude that the new trigonometrically-fitted methods are superior in terms of accuracy and execution time, compared to the existing methods in the scientific literature when used for solving problems which are oscillatory in nature

A two-step trigonometrically fitted semi-implicit hybrid method for solving special second order oscillatory differential equation

In this paper, we derived a semi-implicit hybrid method (SIHM) which is a two-step method to solve special second order ordinary differential equations (ODEs). The SIHM which is three-stage and fourth-order is then trigonometrically fitted and denoted by TF-SIHM3(4). The method is constructed using trigonometrically fitted properties instead of using phase-lag and amplification properties. Numerical integration show that TF-SIHM3(4) is more accurate in term of accuracy compared to the existing explicit and implicit methods of the same order

An eighth-order exponentially fitted two-step hybrid method of explicit type for solving orbital and oscillatory problems

The construction of an eighth-order exponentially fitted (EF) two-step hybrid method for the numerical integration of oscillatory second-order initial value problems (IVPs) is considered. The EF two-step hybrid methods integrate exactly differential systems whose solutions can be expressed as linear combinations of exponential or trigonometric functions and have variable coefficients depending on the frequency of each problem. Based on the order conditions and the EF conditions for this class of methods, we construct an explicit EF two-step hybrid method with symmetric nodes and algebraic order eight which only uses seven function evaluations per step. This new method has the highest algebraic order we know for the case of explicit EF two-step hybrid methods. The numerical experiments carried out with several orbital and oscillatory problems show that the new eighth-order EF scheme is more efficient than other standard and EF two-step hybrid codes recently proposed in the scientific literature

A Functionally-Fitted Block Numerov Method for Solving Second-Order Initial-Value Problems with Oscillatory Solutions

[EN] A functionally-fitted Numerov-type method is developed for the numerical solution of second-order initial-value problems with oscillatory solutions. The basis functions are considered among trigonometric and hyperbolic ones. The characteristics of the method are studied, particularly, it is shown that it has a third order of convergence for the general second-order ordinary differential equation, y′′=f(x,y,y′), it is a fourth order convergent method for the special second-order ordinary differential equation, y′′=f(x,y). Comparison with other methods in the literature, even of higher order, shows the good performance of the proposed method.Open Access funding provided thanks to the CRUE-CSIC agreement with Springer Nature.Publicación en abierto financiada por el Consorcio de Bibliotecas Universitarias de Castilla y León (BUCLE), con cargo al Programa Operativo 2014ES16RFOP009 FEDER 2014-2020 DE CASTILLA Y LEÓN, Actuación:20007-CL - Apoyo Consorcio BUCL

High Order Multistep Methods with Improved Phase-Lag Characteristics for the Integration of the Schr\"odinger Equation

In this work we introduce a new family of twelve-step linear multistep methods for the integration of the Schr\"odinger equation. The new methods are constructed by adopting a new methodology which improves the phase lag characteristics by vanishing both the phase lag function and its first derivatives at a specific frequency. This results in decreasing the sensitivity of the integration method on the estimated frequency of the problem. The efficiency of the new family of methods is proved via error analysis and numerical applications.Comment: 36 pages, 6 figure

A Trigonometrically Fitted Block Method for Solving Oscillatory Second-Order Initial Value Problems and Hamiltonian Systems

In this paper, we present a block hybrid trigonometrically fitted Runge-Kutta-Nyström method (BHTRKNM), whose coefficients are functions of the frequency and the step-size for directly solving general second-order initial value problems (IVPs), including Hamiltonian systems such as the energy conserving equations and systems arising from the semidiscretization of partial differential equations (PDEs). Four discrete hybrid formulas used to formulate the BHTRKNM are provided by a continuous one-step hybrid trigonometrically fitted method with an off-grid point. We implement BHTRKNM in a block-by-block fashion; in this way, the method does not suffer from the disadvantages of requiring starting values and predictors which are inherent in predictor-corrector methods. The stability property of the BHTRKNM is discussed and the performance of the method is demonstrated on some numerical examples to show accuracy and efficiency advantages

A class of explicit high-order exponentially-fitted two-step methods for solving oscillatory IVPs

The derivation of new exponentially fitted (EF) modified two-step hybrid (MTSH) methods for the numerical integration of oscillatory second-order IVPs is analyzed. These methods are modifications of classical two-step hybrid methods so that they integrate exactly differential systems whose solutions can be expressed as linear combinations of the set of functions {exp(¿t), exp(-¿t)}, ¿¿C, or equivalently {sin(¿t), cos(¿t)} when ¿=i¿, ¿¿R, where ¿ represents an approximation of the main frequency of the problem. The EF conditions and the conditions for this class of EF schemes to have algebraic order p (with p=8) are derived. With the help of these conditions we construct explicit EFMTSH methods with algebraic orders seven and eight which require five and six function evaluation per step, respectively. These new EFMTSH schemes are optimal among the two-step hybrid methods in the sense that they reach a certain order of accuracy with minimal computational cost per step. In order to show the efficiency of the new high order explicit EFMTSH methods in comparison to other EF and standard two-step hybrid codes from the literature some numerical experiments with several orbital and oscillatory problems are presented

Block hybrid method with trigonometric-fitting for solving oscillatory problems

In this paper, we develop algebraic order conditions for two-point block hybrid method up to order five using the approach of B-series. Based on the order conditions, we derive fifth order two-point block explicit hybrid method for solving special second order ordinary differential equations (ODEs), where the existing explicit hybrid method of order five is used to be the method at the first point. The method is then trigonometrically fitted so that it can be suitable for solving highly oscillatory problems arising from special second order ODEs. The new trigonometrically-fitted block method is tested using a set of oscillatory problems over a very large interval. Numerical results clearly showed the superiority of the method in terms of accuracy and execution time compared to other existing methods in the scientific literature