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

The purpose of this paper is to study the numerical oscillations of Runge-Kutta methods for the solution of alternately advanced and retarded differential equations with piecewise constant arguments. The conditions of oscillations for the Runge-Kutta methods are obtained. It is proven that the Runge-Kutta methods preserve the oscillations of the analytic solution. In addition, the relationship ...

In this paper we consider the numerical solution of 2D systems of certain types of taxis-diiusion-reaction equations from mathematical biology. By spatial discretization these PDE systems are approximated by systems of positive, nonlinear ODEs (Method of Lines). The aim of this paper is to examine the numerical integration of these ODE systems for low to moderate accuracy by means of splitting ...

Implicit integration schemes for ODEs, such as Runge-Kutta and Runge-Kutta-Nyström methods, are widely used in mathematics and engineering to numerically solve ordinary differential equations. Every integration method requires one to choose a step-size, h, for the integration. If h is too large or too small the efficiency of an implicit scheme is relatively low. As every implicit integration sc...

This paper is concerned with time-stepping numerical methods for computing stiff semi-discrete systems of ordinary differential equations for transient hypersonic flows with thermo-chemical nonequilibrium. The stiffness of the equations is mainly caused by the viscous flux terms across the boundary layers and by the source terms modeling finite-rate thermo-chemical processes. Implicit methods a...

When one solves differential equations, modeling physical phenomena, it is of great importance to take physical constraints into account. More precisely, numerical schemes have to be designed such that discrete solutions satisfy the same constraints as exact solutions. Based on general theory for positivity, with an explicit third-order Runge-Kutta method (we will refer to it as RK3 method) pos...

Article history: Received 31 January 2008 Received in revised form 18 June 2008 Accepted 23 July 2008 Available online 19 August 2008

In this paper we present an a priori error estimate of the Runge-Kutta discontinuous Galerkin method for solving symmetrizable conservation laws, where the time is discretized with the third order explicit total variation diminishing Runge-Kutta method and the finite element space is made up of piecewise polynomials of degree k ≥ 2. Quasi-optimal error estimate is obtained by energy techniques,...