Numerical methods for the eigenvalue determination of second-order ordinary differential equations

An accurate method for the numerical solution of the eigenvalue problem of second-order ordinary differential equation using the shooting method is presented. The method has three steps. Firstly initial values for the eigenvalue and eigenfunction at both ends are obtained by using the discretized matrix eigenvalue method. Secondly the initial-value problem is solved using new, highly accurate formulas of the linear multistep method. Thirdly the eigenvalue is properly corrected at the matching point. The efficiency of the proposed methods is demonstrated by their applications to bound states for the one-dimensional harmonic oscillator, anharmonic oscillators, the Morse potential, and the modified Pöschl–Teller potential in quantum mechanics.     

Supplement: Efficient weak second order stochastic Runge-Kutta methods for non-commutative Stratonovich stochastic differential equations

This paper gives a modification of a class of stochastic Runge-Kutta methods proposed in a paper by Komori (2007). The slight modification can reduce the computational costs of the methods significantly.     

Diagonally implicit general linear methods for ordinary differential equations

We investigate some classes of general linear methods withs internal andr external approximations, with stage orderq and orderp, adjacent to the class withs=r=q=p considered by Butcher. We demonstrate that interesting methods exist also ifs+1=r=q, p=q orq+1,s=r+1=q, p=q orq+1, ands=r=q, p=q+1. Examples of such methods are constructed with stability function matching theA-acceptable generalized Padé approximations to the exponential function.The work of Z. Jackiewicz was partially supported by the National Science Foundation under grant NSF DMS-9208048.     

Weak second order S-ROCK methods for Stratonovich stochastic differential equations

It is well known that the numerical solution of stiff stochastic ordinary differential equations leads to a step size reduction when explicit methods are used. This has led to a plethora of implicit or semi-implicit methods with a wide variety of stability properties. However, for stiff stochastic problems in which the eigenvalues of a drift term lie near the negative real axis, such as those arising from stochastic partial differential equations, explicit methods with extended stability regions can be very effective. In the present paper our aim is to derive explicit Runge–Kutta schemes for non-commutative Stratonovich stochastic differential equations, which are of weak order two and which have large stability regions. This will be achieved by the use of a technique in Chebyshev methods for ordinary differential equations.     

A stability property ofA-stable natural Runge-Kutta methods for systems of delay differential equations

A natural Runge-Kutta method is a special type of Runge-Kutta method for delay differential equations (DDEs); it is known that any one-step collocation method is equivalent to one of such methods. In this paper, we consider a linear constant-coefficient system of DDEs with a constant delay, and discuss the application of natural Runge-Kutta methods to the system. We show that anA-stable method preserves the asymptotic stability property of the analytical solutions of the system.     

Effective computation of periodic orbits and bifurcation diagrams in delay equations

Summary By employing a numerical method which uses only rather classical tools of Numerical Analysis such as Newton's method and routines for ordinary differential equations, unstable periodic solutions of differential-difference equations can be computed. The method is applied to determine bifurcation diagrams with backward bifurcation.Dedicated to Professor Lothar Collatz on the occasion of his 70th birthdayThis paper has been read at the Conference of Numerical Mathematicians, Zeist, The Netherlands, October 12–15, 1979     

Numerical treatment of delay differential equations by Hermite Interpolation

Summary A class of numerical methods for the treatment of delay differential equations is developed. These methods are based on the wellknown Runge-Kutta-Fehlberg methods. The retarded argument is approximated by an appropriate multipoint Hermite Interpolation. The inherent jump discontinuities in the various derivatives of the solution are considered automatically.Problems with piecewise continuous right-hand side and initial function are treated too. Real-life problems are used for the numerical test and a comparison with other methods published in literature.     

Periodic orbits of delay differential equations under discretization

This paper deals with the long-time behaviour of numerical solutions of delay differential equations that have asymptotically stable periodic orbits. It is shown that Runge-Kutta discretizations of such equations have attractive invariant curves which approximate the periodic orbit with the order of the method. The research by this author has been made possible by a fellowship of the Royal Netherlands Academy of Arts and Sciences.     

Runge-Kutta methods for quadratic ordinary differential equations

Many systems of ordinary differential equations are quadratic: the derivative can be expressed as a quadratic function of the dependent variable. We demonstrate that this feature can be exploited in the numerical solution by Runge-Kutta methods, since the quadratic structure serves to decrease the number of order conditions. We discuss issues related to construction design and implementation and present a number of new methods of Runge-Kutta and Runge-Kutta-Nyström type that display superior behaviour when applied to quadratic ordinary differential equations.     

Stability of Runge-Kutta methods for linear delay differential equations

Summary. This paper investigates the stability of Runge-Kutta methods when they are applied to the complex linear scalar delay differential equation . This kind of stability is called stability. We give a characterization of stable Runge-Kutta methods and then we prove that implicit Euler method is stable. Received November 3, 1998 / Revised version received March 23, 1999 / Published online July 12, 2000     

Numerical approach to differential equations of fractional order

In this paper, the variational iteration method and the Adomian decomposition method are implemented to give approximate solutions for linear and nonlinear systems of differential equations of fractional order. The two methods in applied mathematics can be used as alternative methods for obtaining analytic and approximate solutions for different types of differential equations. In these schemes, the solution takes the form of a convergent series with easily computable components. This paper presents a numerical comparison between the two methods for solving systems of fractional differential equations. Numerical results show that the two approaches are easy to implement and accurate when applied to differential equations of fractional order.     

On the convergence of Runge–Kutta methods for stiff non linear differential equations

Summary. This paper studies the convergence properties of general Runge–Kutta methods when applied to the numerical solution of a special class of stiff non linear initial value problems. It is proved that under weaker assumptions on the coefficients of a Runge–Kutta method than in the standard theory of B-convergence, it is possible to ensure the convergence of the method for stiff non linear systems belonging to the above mentioned class. Thus, it is shown that some methods which are not algebraically stable, like the Lobatto IIIA or A-stable SIRK methods, are convergent for the class of stiff problems under consideration. Finally, some results on the existence and uniqueness of the Runge–Kutta solution are also presented. Received November 18, 1996 / Revised version received October 6, 1997     

Stability of the Radau IA and Lobatto IIIC methods for delay differential equations

Recently, we have proved that the Radau IA and Lobatto IIIC methods are P-stable, i.e., they have an analogous stability property to A-stability with respect to scalar delay differential equations (DDEs). In this paper, we study stability of those methods applied to multidimensional DDEs. We show that they have a similar property to P-stability with respect to multidimensional equations which satisfy certain conditions for asymptotic stability of the zero solutions. The conditions are closely related to stability criteria for DDEs considered in systems theory. Received October 8, 1996 / Revised version received February 21, 1997     

Instability of Runge-Kutta methods when applied to linear systems of delay differential equations

Summary. This paper investigates the stability of Runge-Kutta methods when they are applied to the complex linear system of delay differential equations , where . We prove that no Runge-Kutta method preserves asymptotic stability. Received January 24, 2000 / Revised version received July 19, 2000 / Published online June 7, 2001     

Asymptotic stability analysis of Runge-Kutta methods for nonlinear systems of delay differential equations

Summary. We consider systems of delay differential equations (DDEs) of the form with the initial condition . Recently, Torelli [10] introduced a concept of stability for numerical methods applied to dissipative nonlinear systems of DDEs (in some inner product norm), namely RN-stability, which is the straighforward generalization of the wellknown concept of BN-stability of numerical methods with respect to dissipative systems of ODEs. Dissipativity means that the solutions and corresponding to different initial functions and , respectively, satisfy the inequality , and is guaranteed by suitable conditions on the Lipschitz constants of the right-hand side function . A numerical method is said to be RN-stable if it preserves this contractivity property. After showing that, under slightly more stringent hypotheses on the Lipschitz constants and on the delay function , the solutions and are such that , in this paper we prove that RN-stable continuous Runge-Kutta methods preserve also this asymptotic stability property. Received March 29, 1996 / Revised version received August 12, 1996     

On dynamic iteration for delay differential equations

In this paper we study dynamic iteration techniques for systems of nonlinear delay differential equations. After pointing out a close connection to the truncated infinite embedding, as proposed by Feldstein, Iserles, and Levin, we give a proof of the superlinear convergence of the simple dynamic iteration scheme. Then we propose a more general scheme that in addition allows for a decoupling of the equations into disjoint subsystems, just like what we are used to from dynamic iteration schemes for ODEs. This scheme is also shown to converge superlinearly.     

Numerical solution of delay differential equations by uniform corrections to an implicit Runge-Kutta method

Summary In this paper we develop a class of numerical methods to approximate the solutions of delay differential equations. They are essentially based on a modified version, in a predictor-corrector mode, of the one-step collocation method atn Gaussian points. These methods, applied to ODE's, provide a continuous approximate solution which is accurate of order 2n at the nodes and of ordern+1 uniformly in the whole interval. In order to extend the methods to delay differential equations, the uniform accuracy is raised to the order 2n by some a posteriori corrections. Numerical tests and comparisons with other methods are made on real-life problems.This work was supported by CNR within the Progetto Finalizzato Informatica-Sottopr. P1-SOFMAT     

Rooted tree analysis of the order conditions of row-type scheme for stochastic differential equations

Numerical schemes for initial value problems of stochastic differential equations (SDEs) are considered so as to derive the order conditions of ROW-type schemes in the weak sense. Rooted tree analysis, the well-known useful technique for the counterpart of the ordinary differential equation case, is extended to be applicable to the SDE case. In our analysis, the roots are bi-colored corresponding to the ordinary and stochastic differential terms, whereas the vertices have four kinds of label corresponding to the terms derived from the ROW-schemes. The analysis brings a transparent way for the weak order conditions of the scheme. An example is given for illustration.     

Stability analysis of numerical methods for delay differential equations

Summary This paper deals with the stability analysis of step-by-step methods for the numerical solution of delay differential equations. We focus on the behaviour of such methods when they are applied to the linear testproblemU(t)=U(t)+U(t–) with >0 and , complex. A general theorem is presented which can be used to obtain complete characterizations of the stability regions of these methods.