Journal of Mathematics

Volume 2018, Article ID 2960237, 7 pages

https://doi.org/10.1155/2018/2960237

## Solving Oscillatory Delay Differential Equations Using Block Hybrid Methods

^{1}Department of Mathematics, Faculty of Science, Universiti Putra Malaysia, 43400 Serdang, Selangor, Malaysia^{2}Institute for Mathematical Research, Universiti Putra Malaysia, 43400 Serdang, Selangor, Malaysia

Correspondence should be addressed to Fudziah Ismail; ym.moc.oohay@i_haizduf

Received 17 April 2018; Accepted 3 September 2018; Published 1 October 2018

Guest Editor: Muhammad Ozair

Copyright © 2018 Sufia Zulfa Ahmad et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

A set of order condition for block explicit hybrid method up to order five is presented and, based on the order conditions, two-point block explicit hybrid method of order five for the approximation of special second order delay differential equations is derived. The method is then trigonometrically fitted and used to integrate second-order delay differential equations with oscillatory solutions. The efficiency curves based on the log of maximum errors versus the CPU time taken to do the integration are plotted, which clearly demonstrated the superiority of the trigonometrically fitted block hybrid method.

#### 1. Introduction

Differential equations with a time delay are used to model the process which does not only depend on the current state of a system but also the past states. This type of equation is called delay differential equations (DDEs) in which the derivative at any time depends on the solution at prior times. The special second order DDE can be written in the form ofwhere is the delay term and the first derivative does not appear explicitly. It is a more realistic model which includes some of the past history of the system to determine the future behavior. DDEs have become an important criteria to investigate the complexities of the real-world problems concerning infectious diseases, biotic population, neuronal networks, and population dynamics.

Methods such as Runge-Kutta (RK), Runge-Kutta Nyström (RKN), hybrid, and multistep are widely used for solving DDEs. Ismail et al. [1] used RK method and Hermite interpolation to solve first-order DDEs. Taiwo and Odetunde [2] worked on decomposition method as an integrator for delay differential equations. Some authors also derived block linear multistep method (LMM) to solve DDEs; and such work can be seen in [3–6]. Hoo et al. [7] constructed Adams-Moulton Method for directly solving second-order DDEs. Mechee et al. [8] in their paper has adapted RKN for directly solving second-order DDEs.

However, all the studies previously mentioned have not been applied for solving DDE problems with oscillatory properties. Hence, in this paper we derived order condition for block explicit hybrid method up to order five using computer algebra system as proposed in Gander and Gruntz [9]. The reason to derive block hybrid method is that a faster numerical solutions can be obtained since the method approximates the solution at more than one point per step. From the order conditions we constructed a two-point three-stage fifth-order block explicit hybrid method which is then trigonometrically fitted so that it is suitable for solving second-order DDEs which are oscillatory in nature. Finally, we tested the new methods using DDEs test problems to indicate that it is superior and more efficient for solving oscillatory second order DDEs.

#### 2. Derivation of Order Condition for Block Explicit Hybrid Method

The general formula of two-step explicit hybrid method for solving the special second-order ordinary differential equations is given aswhere , and . The method coefficients of , , and can be represented in Butcher tableau as follows:where , , and .

In this section, we derived the order condition for block explicit hybrid method (BEHM). The Taylor series expansion for and is as follows:where is step sizes and is the number of step forward. Adding (5) and (6), we obtainorwhich can be expressed asThe increment function isThen, consider the general formula for block explicit hybrid method in the form ofEquation (11) can be simplified assuch that the Taylor series increment may be written aswhere the first few elementary differential are , , .

The local truncation errors of the solution is obtain by subtracting (12) from (9). This gives usTherefore, from (14), we obtain the order condition for block explicit hybrid method up to order 5 as follows:For the first point (, the order condition is the same as the order conditions for explicit hybrid method given in Coleman [10], while, for the second point (, the order conditions are given as follows: The general formula for block explicit hybrid method (BEHM) for can be written as

#### 3. Construction of Trigonometrically Fitted Block Explicit Hybrid Method

In order to derive the two-point BEHM, we used the algebraic coefficients of the original three-stage explicit hybrid method in Franco [11] as the first point (). The block method can be written as follows:First and second point of the BEHM share the same values of and . By solving (16)-(19) simultaneously, we obtain a unique solution of By checking the fifth-order conditions for the first point, we noticed that the method at the first point ( satisfies , , .

Hence, the method in Franco [11] is order five.

And by checking fifth-order conditions for the second point, we noticed that the method at the second point satisfies Hence, the block explicit hybrid method that we have derived is three-stage and order five, denoted as BEHM3.

To trigonometrically fit the method, we require (23), (24), and (25) to integrate exactly the linear combination of the functions for . Hence, the following equations are obtained:where , is step size, and is the fitted frequency of the problem ( depends on the problems).

By solving (29), (30) with , , and , we obtain the remaining values in terms of as follows:To find values, we solve (31), (32) with choice of coefficients , , and , we obtained the following:Then, for the values, we solve (33), (34) with choice of coefficients and letting and , we obtained in terms of asHence, we denote the new method as three-stage fifth-order trigonometrically fitted block explicit hybrid method (TF-BEHM3). The method is shown as follows:The coefficients can be written in Taylor expansion, to avoid heavy cancellation in the implementation of the method.

#### 4. Problems Tested and Numerical Results

In this section, the new methods, BEHM3 and TF-BEHM3, are used to solve oscillatory delay differential equations problems. The delay terms are evaluated using Newton divided different interpolation. A measure of the accuracy is examined using absolute error which is defined by where is the exact solution and is the computed solution.

The test problems are listed as follows.

*Problem 1 (source: Schmidt [12]). *

*Problem 2 (source: Schmidt [12]). *

*Problem 3 (source: Ladas and Stavroulakis [13]). *

*Problem 4 (source: Bhagat Singh[14]). *

The following notations are used in Figures 1–4: **TF-BEHM3(****5)**: A three-stage fifth-order trigonometrically fitted block explicit hybrid method derived in this paper. **BEHM3(****5)**: A three-stage fifth-order new block explicit hybrid method derived in this paper. **SIHM4(****5)**: A four-stage fifth-order semi-implicit hybrid method in Ahmad et al.[15] **MPAFRKN4(****4)**: A modified phase-fitted and amplification-fitted RKN method of four-stage fourth-order by Papadopoulos et al. [16]. **PFRKN4(****4)**: A phase-fitted RKN method of four-stage fourth-order by Papadopoulos et al. [17]. **DIRKN4(****4)**: A four-stage fourth-order diagonally implicit RKN method by Senu et al. [18]. **EHM4(****5)**: A four-stage fifth-order phase-fitted hybrid method by Franco [11].