**Applied Mathematics**

Vol.06 No.03(2015), Article ID:54971,8 pages

10.4236/am.2015.63053

Differential Transform Method for Some Delay Differential Equations

Baoqing Liu^{1}, Xiaojian Zhou^{2}, Qikui Du^{3}

^{1}School of Applied Mathematics, Jiangsu Provincial Key Laboratory for NSLSCS, Nanjing University of Finance and Economics, Nanjing, China

^{2}School of Science, Jiangsu Provincial Key Laboratory for NSLSCS, Nantong University, Nantong, China

^{3}School of Mathematical Sciences, Jiangsu Key Laboratory for NSLSCS, Nanjing Normal University, Nanjing, China

Email: lyberal@163.com, zxjntu@gmail.com, duqikui@njnu.edu.cn

Copyright © 2015 by authors and Scientific Research Publishing Inc.

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

Received 28 February 2015; accepted 18 March 2015; published 24 March 2015

ABSTRACT

This paper concentrates on the differential transform method (DTM) to solve some delay differential equations (DDEs). Based on the method of steps for DDEs and using the computer algebra system Mathematica, we successfully apply DTM to find the analytic solution to some DDEs, including a neural delay differential equation. The results confirm the feasibility and efficiency of DTM.

**Keywords:**

Differential Transform Method, Delay Differential Equation, Method of Steps, Analytic Solution, Approximate Solution

1. Introduction

The differential transform method (DTM) is a semi analytical-numerical technique depending on Taylor series for solving integral-differential equations (IDEs). The method was first introduced by Pukhov [1] for solving linear and nonlinear initial value problems in physical processes. Zhou, at the same time, had also introduced DTM to study electrical circuits [2] . Since the main advantage of this method is that it can be applied directly to nonlinear ordinary and partial differential equations without requiring linearization, discretization or perturba- tion, it has been studied and applied during the last two decades widely. DTM has been used to obtain numerical and analytical solutions of ordinary differential equations [3] , partial differential equations [4] , eigenvalue pro- blems [5] , differential algebraic equations [6] [7] , integral equations [8] and so on.

Delay differential equations (DDEs) arise in many applied fields, such as control technology, communication networks, and biological population management, and hence they have attracted considerable attention. There are many papers devoted to the problem of approximate solution of DDEs [9] - [15] . Recently, F. Karako and H. Bereketoğlu [13] extend the method of differential transformation for solving the following two types of DDEs:

(1)

and

(2)

with and and the constant.

It should be pointed out that the solution to DDEs (2) maybe be non-unique (see Section 2 in [16] ). So usually, researchers pay more attention to the following DDEs, instead of (2)

(3)

where is a given function, called initial function.

In this paper, we will apply DTM to find the analytic solution to DDEs (3) with the help of the computer algebra system Mathematica. Thus, in some sense, our work can be viewed as a supplement to [13] .

2. Differential Transform

The basic theory of differential transform can be found in [1] [2] , in this section we will state it in brief.

Consider a function be analytic in the time domain, and let. The function is then

represented by one series whose center is located at. The differential transform of the function is the form

(4)

where is the transformed function of the original function.

Differential inverse transformation of is defined as follows:

(5)

From (4) and (5), it is easy to see that the concept of the differential transformation is derived from the Taylor series expansion. By our assumption, is taken as zero, then the function is expressed by a finite series and (5) can be written as

In this study, we use the lower case letters to represent the original functions and upper case letters to stand for the transformed functions (T-functions). The fundamental mathematical operations performed by differential transform method are listed in Table 1.

3. DTM for DDEs (3)

There are many methods to deal with the delay differential Equation (3). For example, linear multistep (LM) methods, Runge-Kutta (RK) methods, waveform relaxation (WR) methods, etc. However, the basic idea to solve the DDE (3) is to solve the following system of ODEs step by step:

Table 1. Estimates of, , and from simulated data.

(6)

with. In brief, this idea is to shift the interval from to and extend the solution from to by using the component in the current interval. This procedure can, in

principle, be continued as far as desired. It is called, quite naturally, the method of steps [16] .

Using the basic idea of the method of steps, first, we apply the DTM to find the solution to the following ODEs:

with.

Suppose the approximate solution is given by

If, is the solution to (3). Otherwise, we should continue to find the solution in the interval. At this time, we should solve the following ODEs

Applying the DTM to the differential equation above again, we will obtain the following solution

Of course, we should go on if holds also. In generally, applying the DTM to ODEs (6), we can obtain the analytic solution

until for some,. In fact, after necessary steps, we have the following solution to (3)

Remark 1 If, we can conclude that is the analytic solution to (3) directly.

Remark 2 If for some integer, we can conclude that the analytic solution (3) is

Remark 3 If we want to improve the accuracy of the approximate solution in each interval, we can combine the above method with the multi-step method given by [17] .

Remark 4 In fact, the DTM based on the method of steps can also be applied to solve the following neutral delay differential equations

4. Numerical Experiments

In this section, four examples are given to show the performance of the DTM based on the method of steps. First, we want to solve the following simple but classical DDE to further illustrate the process of DTM.

Example 4.1 Consider the DDE [18]

(7)

First, since, we apply the DTM to obtain the solution in the interval. In this interval, (7) can be written as, and the initial condition is. Taking the differential transform, we have

It is easy to get

Thus we have the analytic solution of (7) defined on.

Second, we should continue to solve the following DDE:

or equivalently,

(8)

with initial condition.

From (8), we have the following differential transform

It is easy to get

Thus we have the analytic solution

Now, if we want to obtain the solution in the interval [2, 3], we should deal with the following DDE:

or equivalently,

with the initial condition.

Then we have the following differential transform

and get

Thus the analytic solution defined on [2, 3] is given by

The DTM can be proceed till the desire solution is obtained.

Example 4.2 Consider the nonlinear DDE of third-order [11] [13]

(9)

Since, according to the foregoing, we have the following ODE, defined in the interval [0, 0.3]

Thus, applying DTM to the equation above, we obtain

The initial conditions lead to, , and. It is easy to have. Thus we have the solution

Noting that, Remark 1 tells us is the analytic solution to (9) in the whole interval.

Remark 5 In [13] , F. Karako and H. Bereketoğlu also apply DTM to DDE (9) where the initial function, is omitted, i.e.

(10)

It’s worth pointing out that, using the method given in [13] , only approximate solution can be obtained. On the other hand, for DDEs (3), the initial function has the vital role. Without it, the DDEs may have un-unique

solution. In fact, Example 4.2 shows that is a solution to (10). Let’s suppose is

also a solution to (10), then should satisfy the following DDEs

(11)

Figure 1 shows the solution to (11), together with it’s first and second derivative value on, with the

initial function, which satisfies the initial conditions in (10) obviously. It can be seen, then is also the solution to (11). So, (11) has infinite solutions. Maybe, the authors “happen

to” get the approximate solution of.

Example 4.3 Consider a single delay equation with a stiffness parameter [10]

(12)

with, ,.

Similarly, we should solve the following DDE limited in the interval:

Thus applying DTM to this equation, we obtain

Figure 1. and of (11) with.

The initial values lead to. With the help of Mathematica, we have

Then, we obtain the solution to (12):

This is the analytic solution to (12). Particularly, if, then the above solution can be simplified to, which coincides with the definition on. So in this case, is the analytic solution to (12) in the whole interval.

As the last example, we apply the DTM based on the method of steps to solve a neutral delay differential equations.

Example 4.4 Consider the neutral delay differential equation

(13)

with the initial function,

According to the idea of the method of steps, DDE (13) becomes

Applying DTM to this equation, we have

From the initial function, we get, so. With the help of Mathematica, we have

Then, the solution to (13) is

5. Conclusion

Although the theory of differential transform method is not complete yet, it has been successfully applied to solve ordinary differential equations, partial differential equations, integral-differential equations, differential- algebraic equations and etc. In this paper, we apply DTM based on the method of steps to solve some delay differential equations, including neutral delay differential equations, successfully. Numerical experiments show that DTM is feasible and efficient for them. We believe that the operations of DTM presented in this paper also can be used to solve some partial delay differential equations (PDDEs), which is worth while studying in the future work.

Acknowledgements

This work is supported by the National Natural Science Foundation of China, contract/grant number 11371198 and 11401296, Jiangsu Provincial Natural Science Foundation of China, contact/grant no. BK20141008, Natural Science Fund for colleges and universities in Jiangsu Province contact/grant no. 14KJB110007, Jiangsu Provin- cial Key Laboratory for Numerical Simulation of Large Scale Complex Systems contract/grant number 201305 and 201401.

References

- Pukhov, G.E. (1986) Differential Transformations and Mathematical Modelling of Physical Processes. Naukova Dumka, Kiev.
- Zhou, J.K. (1986) Differential Transformation and Its Application for Electrical Circuits. Huazhong University Press, Wuhan. (In Chinese)
- Ayaz, F. (2004) Solutions of the System of Differential Equations by Differential Transform Method. Applied Mathematics and Computation, 147, 547-567. http://dx.doi.org/10.1016/S0096-3003(02)00794-4
- Kurnaz, A., Oturnaz, G. and Kiris, M.E. (2005) n-Dimensional Differential Transformation Method for Solving PDEs. International Journal of Computer Mathematics, 82, 369-380. http://dx.doi.org/10.1080/0020716042000301725
- Ayaz, F. (2004) Applications of Differential Transform Method to Differential-Algebraic Equations. Applied Mathematics and Computation, 152, 649-657. http://dx.doi.org/10.1016/S0096-3003(03)00581-2
- Liu, H. and Song, Y.Z. (2007) Differential Transform Method Applied to High Index Differential-Algebraic Equations. Applied Mathematics and Computation, 184, 748-753. http://dx.doi.org/10.1016/j.amc.2006.05.173
- Abdel-Halim Hassan, I.H. (2002) On Solving Same Eigenvalue Problems by Using a Differential Transformation. Applied Mathematics and Computation, 127, 1-22. http://dx.doi.org/10.1016/S0096-3003(00)00123-5
- Biazar, J., Eslami, M. and Islam, M.R. (2012) Differential Transform Method for Special Systems of Integral Equations. Journal of King Saud University (Science), 24, 211-214. http://dx.doi.org/10.1016/j.jksus.2010.08.015
- Shadia, M. (1992) Numerical Solution of Delay Differential and Neutral Differential Equations Using Spline Methods. Ph.D. Thesis, Assuit University, Assuit.
- El-Hawary, H.M. and Mahmoud, S.M. (2003) Spline Collocation Methods for Solving Delay-Differential Equations. Applied Mathematics and Computation, 146, 359-372. http://dx.doi.org/10.1016/S0096-3003(02)00586-6
- Evans, D.J. and Raslan, K.R. (2005) The Adomian Decomposition Method for Solving Delay Differential Equation. International Journal of Computer Mathematics, 82, 49-54. http://dx.doi.org/10.1080/00207160412331286815
- Vanani, S.K. and Aminataei, A. (2008) On the Numerical Solution of Neutral Delay Differential Equations Using Multiquadric Approximation Scheme. Bulletin of the Korean Mathematical Society, 45, 663-670. http://dx.doi.org/10.4134/BKMS.2008.45.4.663
- Karako, F. and Bereketoğlu, H. (2009) Solutions of Delay Differential Equations by Using Differential Transform Method. International Journal of Computer Mathematics, 86, 914-923. http://dx.doi.org/10.1080/00207160701750575
- Lainscsek, C. and Sejnowski, T.J. (2015) Delay Differential Analysis of Time Series. Neural Computation, 27, 594- 614. http://dx.doi.org/10.1162/NECO_a_00706
- Abazari, R. and Kilicman, A. (2014) Application of Differential Transform Method on Nonlinear Integro-Differential Equations with Proportional Delay. Neural Computing and Applications, 24, 391-397. http://dx.doi.org/10.1007/s00521-012-1235-4
- Driver, R.D. (1977) Ordinary and Delay Differential Equations. Applied Mathematical Sciences, 20. http://dx.doi.org/10.1007/978-1-4684-9467-9
- Odibat, Z.M., et al. (2010) A Multi-Step Differential Transform Method and Application to Non-Chaotic or Chaotic Systems. Computers & Mathematics with Applications, 59, 1462-1472. http://dx.doi.org/10.1016/j.camwa.2009.11.005
- Hairer, E., et al. (1993) Solving Ordinary Differential Equations I: Nonstiff Problems. 2nd Edition, Springer-Verlag, Berlin.