Solution of Delay Differential Equations Using a Modified Power Series Method

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.


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.
In this paper, we consider a modified power series method for solving the delay differential equations of the form.
where n and ( ) i n i > denote the order of derivatives.

Description of the Method
According to the MPSM, the Nth degree approximate solution to the DDE (1)-( 2) is given by ( ) (3) where 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 ( ) and the derivatives of ( ) N y t up to order n.However, if the nonhomogeneous term or a coefficient of y or its derivative is not a polynomial but analytic at 0 t = , then replace it by its Taylor series expansion of degree N about 0 t = .
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 0 t up to N n t − ) together with the n equations generated by the associated given conditions to find j a , 0,1, , j N =  .
Step 5 Substitute the coefficients j a determined in Step 4 into Equation ( 3) to obtain an approximate solution of degree N to the delay differential Equation ( 1)-(2).

Illustrative Examples
Example 3.1 (see [8] [9] [15]) Consider the first-order nonlinear DDE ( ) Subject to ( ) Applying the proposed method illustrated in section 2 to this problem for the cases N = 3, 5, 7, 8, 11, we obtain as follows: with the proposed method illustrated in section 2, we obtain for the cases N = 5, 7, 9, 11 and 13 as follows: (