**Applied Mathematics**

Vol.06 No.04(2015), Article ID:55756,4 pages

10.4236/am.2015.64061

Solution of Delay Differential Equations Using a Modified Power Series Method

Oladotun Matthew Ogunlaran^{1}, Adeyemi Sunday Olagunju^{2}

^{1}Department of Mathematics and Statistics, Bowen University, Iwo, Nigeria

^{2}Department of Mathematics, Federal University, Lafia, Nigeria

Email: dothew2002@yahoo.com

Copyright © 2015 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

Received 29 January 2015; accepted 10 April 2015; published 17 April 2015

ABSTRACT

This paper presents a Modified Power Series Method (MPSM) for the solution of delay differential equations. Unlike the traditional power series method which is applied to solve only linear differential equations, this new approach is applicable to both linear and nonlinear problems. The method produces a system of algebraic equations which is solved to determine the coefficients in the trial solution. The method provides the solution in form of a rapid convergent series. The obtained results for numerical examples demonstrate the reliability and efficiency of the method.

**Keywords:**

Delay Differential Equations, Power Series, Taylor Series, Newton’s Method

1. Introduction

Although both ordinary differential equations (ODEs) and delay differential equations (DDEs) are used to describe physical phenomena, they are different. While in ODEs the derivatives of unknown functions are dependent on only the current value of the independent variable, in addition to this in DDEs, the derivatives of unknown functions are dependent on the values of the functions at previous time. This implies that the solution of DDEs requires the knowledge of the current state and the state at certain previous times.

Some of the application areas of delay differential equations are population dynamics, infectious disease, physiological and pharmaceutical kinetics, chemical kinetics, models of conveyor belts, urban traffic, heat exchangers, robotics, navigational control of ships and aircrafts, and more general control problems (see [1] -[4] ).

There are a few classes of nonlinear ODEs for which solutions can be easily found and despite the obvious similarities between ODEs and DDEs, solutions of DDE problems can differ from solutions for ODE problems in several striking and significant way [5] . Delay problems always lead to an infinite spectrum of frequencies. Therefore, they are solved by numerical methods, asymptotic solution, approximations and graphical approaches [6] . As a result of a dramatic increase in the application of delay models, many authors have investigated and proposed various methods for solving DDEs, Spline methods (see [7] -[9] ), Homotopy Analysis Method (HAM) [10] , Homotopy perturbation method (HPM) [11] , Adomain decomposition method (see [12] -[14] ), Iterative decomposition method [15] and Variation Iteration Method (VIM) [6] .

In this paper, we consider a modified power series method for solving the delay differential equations of the form.

(1)

(2)

where n and denote the order of derivatives.

2. Description of the Method

According to the MPSM, the Nth degree approximate solution to the DDE (1)-(2) is given by

(3)

where

(4)

The MPSM is described by the following five-step procedure:

Step 1

Rewrite Equation (1) such that only the nonhomogeneous term is on the right hand side of the equation

Step 2

On the left hand side of the nonhomogeneous differential equation, substitute

(5)

(6)

and the derivatives of up to order n. However, if the nonhomogeneous term or a coefficient of y or its derivative is not a polynomial but analytic at, then replace it by its Taylor series expansion of degree N about.

Step 3

Collect the power of t on the left hand side of the equation resulting from step 2 and set the coefficient of each power of t on the left hand side equal to the corresponding coefficient on the right hand side of the equation.

Step 4

Solve, using either the Newton’s method or forward substitution method, the first equations resulting from Step 3 (these are the equations obtained by equating the coefficients of the corresponding powers of t from up to) together with the n equations generated by the associated given conditions to find,.

Step 5

Substitute the coefficients determined in Step 4 into Equation (3) to obtain an approximate solution of degree N to the delay differential Equation (1)-(2).

3. Illustrative Examples

Example 3.1 (see [8] [9] [15] )

Consider the first-order nonlinear DDE

(7)

Subject to

(8)

Applying the proposed method illustrated in section 2 to this problem for the cases N = 3, 5, 7, 8, 11, we obtain as follows:

Obviously,

which is the exact solution to Example 3.1.

Example 3.2 (see [7] [8] [10] [15] )

Consider the second-order linear DDE

(9)

Subject to

(10)

Using the proposed method for, we obtain which is the exact solution to this example.

Example 3.3 (see [8] [9] [12] [15] )

Consider the third-order nonlinear DDE

(11)

Subject to the conditions

(12)

with the proposed method illustrated in section 2, we obtain for the cases N = 5, 7, 9, 11 and 13 as follows:

Clearly,

which is the exact solution to Example 3.3.

Example 3.4 (see [6] [10] [12] )

Consider the first-order linear DDE

(13)

Subject to

(14)

The exact solution is given by

For, the present method yields in turn

which converges to the exact solution to Example 3.4.

4. Conclusion

A simple and straight forward technique based on the power series method has been studied for the solution of delay differential equations. This new approach is implemented without using restrictive assumptions or adding perturbation term and it gives excellent performance compared with existing techniques for solving delay differential equations.

References

- Lv, C. and Yuan, Z. (2009) Stability Analysis of Delay Differentialm Equation Models of HIV-1 Therapy for Fighting a Virus with Anothervirus, Journal of Mathematical Analysis and Applications, 352, 672-683. http://dx.doi.org/10.1016/j.jmaa.2008.11.026
- Marchuk, G.I. (1994) Mathematical Modelling of Immune Response Infection Disease. Kluwer Academic Publishers, Dordrecht.
- Glass, L. and Mackey, J.C. (1979) Pathological Conditions Resulting from Instabilities in Physiological Control Systems. Annals of the New York Academy of Sciences, 316, 214-235. http://dx.doi.org/10.1111/j.1749-6632.1979.tb29471.x
- Bunsenberg, S. and Tang, B. (1994) Mathematical Models of the Early Embryonic Cell Cycle: The Role of MPF Activation and Cyclin Degradation. Journal of Mathematical Biology, 32, 573-596. http://dx.doi.org/10.1007/BF00573462
- Baker, C.T.H., Paul, C.A.H. and Wille, D.R. (1995) Issues in the Numericalsolution of Evolutionary Delay Differential Equations. Advances in Computational Mathematics, 3, 171-196. http://dx.doi.org/10.1007/BF03028370
- Rangkuti, Y.M. and Noorani, M.S.M. (2012) The Exact Solution of Delay Differential Equations Using Coupling Vari- ational Iteration with Taylor Series and Small Term. Bulletin of Mathematics, 4, 1-15.
- El-Safty, A. (1993) Approximate Solution of the Delay Differential Equation with Cubic Spline Functions. Bulletin of the Faculty Science, Assiut University, 22, 67-73.>http://html.scirp.org/file/4-7402632x39.png" class="200" /> with Cubic Spline Functions. Bulletin of the Faculty Science, Assiut University, 22, 67-73.
- EI-Safty, A. and Abo-Hasha, S.M. (1990) On the Application of Spline Functions to Initial Value Problem with Re- tarded Argument. International Journal of Computer Mathematics, 32, 137-179.
- Ibrahim, M.A.K., EI-Safty, A. and Abo-Hasha, S.M. (1995) 2h-Step Spline Method for the Solution of Delay Differential Equations. Computers & Mathematics with Applications, 29, 1-6. http://dx.doi.org/10.1016/0898-1221(95)00024-S
- Alomari, A.K., Noorani, M.S.M. and Nazar, R. (2009) Solution of Delay Differential Equation by Means of Homotopy Analysis Method. Acta Applicandae Mathematicae, 108, 395-412. http://dx.doi.org/10.1007/s10440-008-9318-z
- Shakeri, F. and Dehghan, M. (2008) Solution of Delay Differentialequations via a Homotopy Perturbation Method. Mathematical and Computer Modelling, 48, 486-498. http://dx.doi.org/10.1016/j.mcm.2007.09.016
- Evans, D.J. and Raslan, K.R. (2005) The Adomian Decomposition Methodfor Solving Delay Differential Equation. International Journal of Computer Mathematics, 82, 49-54. http://dx.doi.org/10.1080/00207160412331286815
- Adomian, G. and Rach, R. (1983) A Nonlinear Delay Differential Equations. Journal of Mathematical Analysis and Applications. 91, 301-304. http://dx.doi.org/10.1016/0022-247X(83)90152-X
- Adomian, G. and Rach, R. (1983) Nonlinear Stochastic Differential Delay Equations. Journal of Mathematical Analysis and Applications, 91, 94-101. http://dx.doi.org/10.1016/0022-247X(83)90094-X
- Taiwo, O.A. and Odetunde, O.S. (2010) On the Numerical Approximation of Delay Differential Equations by a Decomposition Method. Asian Journal of Mathematics and Statistics, 3, 237-243. http://dx.doi.org/10.3923/ajms.2010.237.243