Abstract and Applied Analysis

Volume 2013 (2013), Article ID 493204, 9 pages

http://dx.doi.org/10.1155/2013/493204

## Exponential Collocation Method for Solutions of Singularly Perturbed Delay Differential Equations

^{1}Department of Mathematics, Faculty of Science, Akdeniz University, Antalya 07058, Turkey^{2}Department of Mathematics, Faculty of Science, Celal Bayar University, Manisa 45040, Turkey

Received 16 June 2013; Accepted 22 July 2013

Academic Editor: Valery Y. Glizer

Copyright © 2013 Şuayip Yüzbaşı and Mehmet Sezer. 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

This paper deals with the singularly perturbed delay differential equations under boundary conditions. A numerical approximation based on the exponential functions is proposed to solve the singularly perturbed delay differential equations. By aid of the collocation points and the matrix operations, the suggested scheme converts singularly perturbed problem into a matrix equation, and this matrix equation corresponds to a system of linear algebraic equations. Also, an error analysis technique based on the residual function is introduced for the method. Four examples are considered to demonstrate the performance of the proposed scheme, and the results are discussed.

#### 1. Introduction

The mathematical models of many practical phenomena in many areas of sciences often result in boundary-value problems of singularly perturbed delay differential equations, for example, the study of bistable devices [1], description of the human pupil-light reflex [2], a variety of models for physiological processes or diseases [3, 4], evolutionary biology [4], variational problems in control theory [5, 6], and so forth. These problems mainly depend on a small positive parameter and a delay parameter in such a way that the solution varies rapidly in some parts of the domain and varies slowly in some other parts of the domain. Also, this class of problems possesses boundary layers, that is, regions of rapid change in the solution near one of the boundary points.

In the recent years, many researchers have a great interest in singularly perturbed delay differential equation problems. For example, for these problems, Patidar and Sharma [7] have studied nonstandard finite difference methods, Kadalbajoo and Ramesh [8] have applied the hybrid method, Kadalbajoo and Sharma [9] have presented a numerical study, Rai and Sharma [10] have worked on a numerical scheme based fitted operator methods, Kadalbajoo and Ramesh [11] have studied the finite difference scheme, the hybrid method, and the fitted mesh methods, Amiraliyev and Erdogan [12] have presented a numerical study based on finite difference scheme and piecewise-uniform mesh, Kadalbajoo and Kumar [13] have applied the fitted mesh B-spline collocation method, and Kadalbajoo and Sharma [14] have presented a numerical study involved which finite difference scheme. In addition, Rai and Sharma [15] have solved the singularly perturbed differential difference equation arising in the modeling of neuronal variability by using fitted operator scheme, Lange and Miura [16] have given the singular perturbation analysis of boundary-value problems for differential difference equations and boundary-value problems, and Kadalbajoo and Sharma [17, 18] have studied the numerical solutions of singularly perturbed delay differential equations by various methods.

On the other hand, exponential polynomials or exponential functions have interesting applications in many optical and quantum electronics [19], some nonlinear phenomena modeled by partial differential equations [20], many statistical discussions (especially in data analysis) [21], the safety analysis of control synthesis [22], the problem of expressing mean-periodic functions [23], and the study of spectral synthesis [24, 25]. These polynomials are based on the exponential base set .

Recently, Yüzbaşı and Sezer have studied the exponential polynomial solutions of the systems of linear differential equations in [26].

In this study, we consider the singularly perturbed delay differential equation with the boundary conditions where is a small positive parameter (), is a small shifting parameter , and are given constants, is an unknown function, and and are the known functions defined on interval .

The aim of this paper is to give an approximate solution of the problems (1)-(2) in the form where the exponential basis set is defined by and , are unknown coefficients.

To find a solution in the form (3) of (1) under the conditions (2), we will use the equally spaced collocation points

#### 2. Matrix Relations for Exponential Functions

In this section, we construct the matrix relations related to the exponential solution (3). Note that these relations will be used in Section 3.

Firstly, the approximate solution defined by (3) of (1) can be written in the matrix form where Also, the relation between and its first derivative is given by and the relation between and its second derivative is in the form where By placing the relation (7) into first derivative of (5), we have Similarly, from the relations (5) and (8), we obtain the matrix form By writing in (10), we get the relation The relation between and is as follows: where By substituting (13) into (12), we have the matrix form

#### 3. Exponential Collocation Method

In this section, to compute the unknown coefficients in the approximate solution (3), we use the following procedure by using the matrix relations in Section 2.

Firstly, let us substitute the matrix relations (5), (11), (12), and (15) into (1) as follows: The collocation points defined by (4) are placed into (16), and we have the system of the matrix equations as The system can be written in the matrix form where We note that the matrix equation (18) corresponds to a system of algebraic equations with the unknown coefficients .

Briefly, (18) can be expressed in the form where From the relation (5), the matrix forms of the conditions (2) are written as Briefly, we write the above matrix forms of the conditions as follows: where To obtain the solution of (1) under conditions (2), we replace the row matrices (23) with any two rows of the matrix (20), and thus we have the augmented matrix For simplicity, if the last two rows of the matrix (20) are replaced, the augmented matrix (25) becomes However, we do not have to replace the last rows. For example, if the matrix is singular, then the rows that have the same factor or all zeros are replaced.

If , then the coefficients are uniquely determined by where By substituting the determined coefficients into (3), we obtain the approximate solution On the other hand, when , if , then we may find a particular solution. Otherwise if , then there is not a solution.

#### 4. Error Estimation Based on Residual Function and Improvement of Solution

In this section, we apply the error estimation technique in [27–29] and the residual correction method in [30, 31] for our method and the problems (1)-(2). For our purpose, let us define the residual function for the present method as where denotes the approximate solution (29) of the problems (1)-(2). Therefore, satisfies and the conditions If is the exact solution of the problems (1)-(2), then becomes the error function. By substituting (33) into the problem (1)-(2) and by using (30), we obtain the error differential equation By using (33), the inhomogeneous conditions (2) and (32) are reduced to the homogeneous conditions From (34) and (35), we can clearly write the error problem By solving the problem (36) in the same way as Section 3, the approximation is obtained for .

Consequently, by summing the exponential polynomial solution and the estimated error function , we obtain the corrected exponential solution We note that the errors , () can be estimated by the error function when the exact solution of (1) is unknown.

#### 5. Numerical Examples

In this section, we apply the presented method to some examples. In examples, the terms , , , and , respectively, represent the exact solution, the approximate solution, the corrected approximate solution, and the estimated absolute error function. Also, denotes the estimated maximum error for the values , , , and .

*Example 1 (see [13]). *Firstly, let us consider the singularly-perturbed delay differential equation
with the boundary conditions
Firstly, we obtain the approximate solutions for various values of by the presented method in Section 3. Secondly, the approximate solutions are corrected by the residual correction technique for various values of . Hence, the corrected approximate solutions are obtained. In Table 1, we give the estimated maximum absolute errors for various values of , , and . Figures 1(a), 1(b), and 1(c) display the corrected approximate solutions for some values of , , and .

*Example 2 (see [7]). *Now, we consider the singularly perturbed delay differential equation
with the boundary conditions
In Table 2, the estimated maximum absolute errors for various values of , , and are presented. For various values of , , and , Figures 2(a), 2(b), 2(c), and 2(d) show the graphs of the corrected approximate solutions .

*Example 3 (see [9]). *Let us consider the boundary-value problem
The exact solution of this problem is given by
where
For some values of , the corrected approximate solutions are compared with the exact solution in Figures 3(a) and 3(b).

*Example 4 (see [14]). *Finally, we consider the problem
Figure 4(a) displays the corrected approximate solutions for , , and various values of and . The estimated error functions for these approximate solutions are shown in Figure 4(b). It is seen from Figure 4(b) that the absolute errors decrease while value of increases. In Figure 4(c), we show the corrected approximate solutions for , , and different values of . Figure 4(d) shows the estimated error functions for the approximate solutions in Figure 4(c).

#### 6. Conclusions

In this paper, a numerical scheme based on the exponential functions and the collocation points is presented for the singularly perturbed delay differential equations. Numerical examples are given to demonstrate the applicability and the efficiency of the method. Since the exact solutions of the problems in Examples 1, 2, and 4 are not available, we have computed the estimated maximum absolute errors. For Examples 1 and 2, the maximum absolute errors for some values of , , , and are tabulated in Tables 1 and 2. Also, the approximate solutions for different values of , , , and are compared with the exact solution. It is seen from Figure 3(a) that the accuracies of the approximate solutions increase while and increase. However, when values of and are increased, it is observed from Tables 1 and 2 and Figure 4(c) that the errors usually increase. In addition, the estimated absolute error functions are shown in Figure 4(b) for Example 3. It is seen from Figure 4(b) that the errors decrease as and increase. Moreover, the approximate solutions can be very easily obtained using the software programs.

#### Acknowledgments

The authors would like to thank the referees for their constructive comments and suggestions to improve the paper. The first author is supported by *the Scientific Research Project Administration* of Akdeniz University.

#### References

- M. W. Derstine, H. M. Gibbs, F. A. Hopf, and D. L. Kaplan, “Bifurcation gap in a hybrid optically bistable system,”
*Physical Review A*, vol. 26, no. 6, pp. 3720–3722, 1982. View at Publisher · View at Google Scholar · View at Scopus - A. Longtin and J. G. Milton, “Complex oscillations in the human pupil light reflex with “mixed” and delayed feedback,”
*Mathematical Biosciences*, vol. 90, no. 1-2, pp. 183–199, 1988. View at Publisher · View at Google Scholar · View at MathSciNet - M. C. Mackey and L. Glass, “Oscillation and chaos in physiological control systems,”
*Science*, vol. 197, no. 4300, pp. 287–289, 1977. View at Google Scholar · View at Scopus - M. Ważewska-Czyżewska and A. Lasota, “Mathematical problems of the dynamics of a system of red blood cells,”
*Matematyka Stosowana*, vol. 6, pp. 23–40, 1976. View at Google Scholar · View at MathSciNet - V. Y. Glizer, “Asymptotic solution of a singularly perturbed set of functional-differential equations of Riccati type encountered in the optimal control theory,”
*Nonlinear Differential Equations and Applications*, vol. 5, no. 4, pp. 491–515, 1998. View at Publisher · View at Google Scholar · View at MathSciNet - V. Y. Glizer, “Asymptotic solution of a boundary-value problem for linear singularly-perturbed functional differential equations arising in optimal control theory,”
*Journal of Optimization Theory and Applications*, vol. 106, no. 2, pp. 309–335, 2000. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - K. C. Patidar and K. K. Sharma, “$\epsilon $-uniformly convergent non-standard finite difference methods for singularly perturbed differential difference equations with small delay,”
*Applied Mathematics and Computation*, vol. 175, no. 1, pp. 864–890, 2006. View at Publisher · View at Google Scholar · View at MathSciNet - M. K. Kadalbajoo and V. P. Ramesh, “Hybrid method for numerical solution of singularly perturbed delay differential equations,”
*Applied Mathematics and Computation*, vol. 187, no. 2, pp. 797–814, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. K. Kadalbajoo and K. K. Sharma, “Numerical analysis of singularly perturbed delay differential equations with layer behavior,”
*Applied Mathematics and Computation*, vol. 157, no. 1, pp. 11–28, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - P. Rai and K. K. Sharma, “Numerical analysis of singularly perturbed delay differential turning point problem,”
*Applied Mathematics and Computation*, vol. 218, no. 7, pp. 3483–3498, 2011. View at Publisher · View at Google Scholar · View at MathSciNet - M. K. Kadalbajoo and V. P. Ramesh, “Numerical methods on Shishkin mesh for singularly perturbed delay differential equations with a grid adaptation strategy,”
*Applied Mathematics and Computation*, vol. 188, no. 2, pp. 1816–1831, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - G. M. Amiraliyev and F. Erdogan, “Uniform numerical method for singularly perturbed delay differential equations,”
*Computers & Mathematics with Applications*, vol. 53, no. 8, pp. 1251–1259, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. K. Kadalbajoo and D. Kumar, “Fitted mesh $B$-spline collocation method for singularly perturbed differential-difference equations with small delay,”
*Applied Mathematics and Computation*, vol. 204, no. 1, pp. 90–98, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. K. Kadalbajoo and K. K. Sharma, “A numerical method based on finite difference for boundary value problems for singularly perturbed delay differential equations,”
*Applied Mathematics and Computation*, vol. 197, no. 2, pp. 692–707, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - P. Rai and K. K. Sharma, “Numerical study of singularly perturbed differential-difference equation arising in the modeling of neuronal variability,”
*Computers & Mathematics with Applications*, vol. 63, no. 1, pp. 118–132, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - C. G. Lange and R. M. Miura, “Singular perturbation analysis of boundary value problems for differential-difference equations. V. Small shifts with layer behavior,”
*SIAM Journal on Applied Mathematics*, vol. 54, no. 1, pp. 249–272, 1994. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. K. Kadalbajoo and K. K. Sharma, “Numerical treatment of boundary value problems for second order singularly perturbed delay differential equations,”
*Computational & Applied Mathematics*, vol. 24, no. 2, pp. 151–172, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. K. Kadalbajoo and K. K. Sharma, “An $\epsilon $-uniform fitted operator method for solving boundary-value problems for singularly perturbed delay differential equations: layer behavior,”
*International Journal of Computer Mathematics*, vol. 80, no. 10, pp. 1261–1276, 2003. View at Publisher · View at Google Scholar · View at MathSciNet - F. Alharbi, “Predefined exponential basis set for half-bounded multi domain spectral method,”
*Applied Mathematics*, vol. 1, pp. 146–152, 2010. View at Google Scholar - M. M. Alipour, G. Domairry, and A. G. Davodi, “An application of exp-function method to approximate general and explicit solutions for nonlinear Schrödinger equations,”
*Numerical Methods for Partial Differential Equations*, vol. 27, no. 5, pp. 1016–1025, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - R. Shanmugam, “Generalized exponential and logarithmic polynomials with statistical applications,”
*International Journal of Mathematical Education in Science and Technology*, vol. 19, no. 5, pp. 659–669, 1988. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Xu, L. Chen, Z. Zeng, and Z.-B. Li, “Reachability analysis of rational eigenvalue linear systems,”
*International Journal of Systems Science*, vol. 41, no. 12, pp. 1411–1419, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - H. Ouerdiane and M. Ounaies, “Expansion in series of exponential polynomials of mean-periodic functions,”
*Complex Variables and Eliptic Equations*, vol. 57, no. 5, pp. 469–487, 2012. View at Publisher · View at Google Scholar - L. Székelyhidi, “On the extension of exponential polynomials,”
*Mathematica Bohemica*, vol. 125, no. 3, pp. 365–370, 2000. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - K. Ross,
*Abstract Harmonic Analysis I, II*, Springer, Berlin, Germany, 1963. - Ş. Yüzbaşı and M. Sezer, “An exponential matrix method for solving systems of linear differential equations,”
*Mathematical Methods in the Applied Sciences*, vol. 36, no. 3, pp. 336–348, 2013. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Shahmorad, “Numerical solution of the general form linear Fredholm-Volterra integro-differential equations by the Tau method with an error estimation,”
*Applied Mathematics and Computation*, vol. 167, no. 2, pp. 1418–1429, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Ş. Yüzbaşı, “An efficient algorithm for solving multi-pantograph equation systems,”
*Computers & Mathematics with Applications*, vol. 64, no. 4, pp. 589–603, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Ş. Yüzbaşı and M. Sezer, “An improved Bessel collocation method with a residual error function to solve a class of Lane-Emden differential equations,”
*Mathematical and Computer Modelling*, vol. 57, no. 5-6, pp. 1298–1311, 2013. View at Google Scholar - İ. Çelik, “Collocation method and residual correction using Chebyshev series,”
*Applied Mathematics and Computation*, vol. 174, no. 2, pp. 910–920, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - F. A. Oliveira, “Collocation and residual correction,”
*Numerische Mathematik*, vol. 36, no. 1, pp. 27–31, 1980. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet