In this paper we consider the stability of variable step-size Runge-Kutta methods approximating bounded, stable, and time-dependent solutions of ordinary differential equation initial value problems. We use Lyapunov exponent theory to determine conditions on the maximum allowable step-size that guarantees the numerical solution of an asymptotically decaying time-dependent linear problem also de...

The equations governing the flow of an electrically conducting, incompressible viscous fluid over an infinite flat plate in the presence of a magnetic field are investigated using the homotopy perturbation method (HPM) with Padé approximants (PA) and 4 order Runge–Kutta method (4RKM). Approximate analytical and numerical solutions for the velocity field and heat transfer are obtained and compar...

In recent years, implicit stochastic Runge–Kutta (SRK) methods have been developed both for strong and weak approximations. For these methods, the stage values are only given implicitly. However, in practice these implicit equations are solved by iterative schemes such as simple iteration, modified Newton iteration or full Newton iteration. We employ a unifying approach for the construction of ...

This paper considers the asymptotic stability analysis of both exact and numerical solutions of the following neutral delay differential equation with pantograph delay. ⎧⎨ ⎩ x′(t) +Bx(t) + Cx′(qt) +Dx(qt) = 0, t > 0, x(0) = x0, where B,C,D ∈ Cd×d, q ∈ (0, 1), and B is regular. After transforming the above equation to non-automatic neutral equation with constant delay, we determine sufficient co...

Runge-Kutta and Adams methods are the most popular codes to solve numerically nonstiff ODEs. The Adams methods are useful to reduce the number of function calls, but they usually require more CPU time than the Runge-Kutta methods. In this work we develop a numerical study of a variable step length Adams implementation, which can only take preassigned step-size ratios. Our aim is the reduction o...

This work generalizes the additively partitioned Runge-Kutta methods by allowing for different stage values as arguments of different components of the right hand side. An order conditions theory is developed for the new family of generalized additive methods, and stability and monotonicity investigations are carried out. The paper discusses the construction and properties of implicit-explicit ...

In this work we introduce a new approach to Dynamical Monte Carlo methods to simulate markovian processes. We apply this approach to formulate and study an epidemic generalized SIRS model. The results are in excellent agreement with the fourth order Runge-Kutta method in a region of deterministic solution. Introducing local stochastic interactions, the Runge-Kutta method is no longer applicable...

We consider the preservation of weak solution invariants in the time integration of ordinary diier-ential equations (ODEs). Recent research has concentrated on obtaining symplectic discretizations of Hamiltonian systems and schemes that preserve certain rst integrals (i.e. strong invariants). In this article, we examine the connection between constrained systems and ODEs with weak invariants fo...

We will consider the efficient implementation of a fourth order two stage implicit Runge-Kutta method to solve periodic second order initial value problems. To solve the resulting systems, we will use the factorization of the discretized operator. Such proposed factorization involves both complex and real arithmetic. The latter case is considered here. The resulting system will be efficient and...

In this paper the performance of a parallel iterated Runge-Kutta method is compared versus those of the serial fouth order Runge-Kutta and Dormand-Prince methods. It was found that, typically, the runtime for the parallel method is comparable to that of the serial versions, thought it uses considerably more computational resources. A new algorithm is proposed where full parallelization is used ...