Research Article  Open Access
Şuayip Yüzbaşı, Mehmet Sezer, "Exponential Collocation Method for Solutions of Singularly Perturbed Delay Differential Equations", Abstract and Applied Analysis, vol. 2013, Article ID 493204, 9 pages, 2013. https://doi.org/10.1155/2013/493204
Exponential Collocation Method for Solutions of Singularly Perturbed Delay Differential Equations
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 boundaryvalue problems of singularly perturbed delay differential equations, for example, the study of bistable devices [1], description of the human pupillight 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 piecewiseuniform mesh, Kadalbajoo and Kumar [13] have applied the fitted mesh Bspline 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 boundaryvalue problems for differential difference equations and boundaryvalue 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 meanperiodic 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 singularlyperturbed 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 .

(a)
(b)
(c)
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 .

(a)
(b)
(c)
(d)
Example 3 (see [9]). Let us consider the boundaryvalue 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).
(a)
(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).
(a)
(b)
(c)
(d)
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 Site  Google Scholar
 A. Longtin and J. G. Milton, “Complex oscillations in the human pupil light reflex with “mixed” and delayed feedback,” Mathematical Biosciences, vol. 90, no. 12, pp. 183–199, 1988. View at: Publisher Site  Google Scholar  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
 M. WażewskaCzyż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  MathSciNet
 V. Y. Glizer, “Asymptotic solution of a singularly perturbed set of functionaldifferential 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 Site  Google Scholar  MathSciNet
 V. Y. Glizer, “Asymptotic solution of a boundaryvalue problem for linear singularlyperturbed 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 Site  Google Scholar  Zentralblatt MATH  MathSciNet
 K. C. Patidar and K. K. Sharma, “$\epsilon $uniformly convergent nonstandard 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 Site  Google Scholar  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 Site  Google Scholar  Zentralblatt MATH  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 Site  Google Scholar  Zentralblatt MATH  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 Site  Google Scholar  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 Site  Google Scholar  Zentralblatt MATH  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 Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. K. Kadalbajoo and D. Kumar, “Fitted mesh $B$spline collocation method for singularly perturbed differentialdifference equations with small delay,” Applied Mathematics and Computation, vol. 204, no. 1, pp. 90–98, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  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 Site  Google Scholar  Zentralblatt MATH  MathSciNet
 P. Rai and K. K. Sharma, “Numerical study of singularly perturbed differentialdifference equation arising in the modeling of neuronal variability,” Computers & Mathematics with Applications, vol. 63, no. 1, pp. 118–132, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 C. G. Lange and R. M. Miura, “Singular perturbation analysis of boundary value problems for differentialdifference equations. V. Small shifts with layer behavior,” SIAM Journal on Applied Mathematics, vol. 54, no. 1, pp. 249–272, 1994. View at: Publisher Site  Google Scholar  Zentralblatt MATH  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 Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. K. Kadalbajoo and K. K. Sharma, “An $\epsilon $uniform fitted operator method for solving boundaryvalue 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 Site  Google Scholar  MathSciNet
 F. Alharbi, “Predefined exponential basis set for halfbounded 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 expfunction 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 Site  Google Scholar  Zentralblatt MATH  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 Site  Google Scholar  Zentralblatt MATH  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 Site  Google Scholar  Zentralblatt MATH  MathSciNet
 H. Ouerdiane and M. Ounaies, “Expansion in series of exponential polynomials of meanperiodic functions,” Complex Variables and Eliptic Equations, vol. 57, no. 5, pp. 469–487, 2012. View at: Publisher Site  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  Zentralblatt MATH  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 Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. Shahmorad, “Numerical solution of the general form linear FredholmVolterra integrodifferential equations by the Tau method with an error estimation,” Applied Mathematics and Computation, vol. 167, no. 2, pp. 1418–1429, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Ş. Yüzbaşı, “An efficient algorithm for solving multipantograph equation systems,” Computers & Mathematics with Applications, vol. 64, no. 4, pp. 589–603, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Ş. Yüzbaşı and M. Sezer, “An improved Bessel collocation method with a residual error function to solve a class of LaneEmden differential equations,” Mathematical and Computer Modelling, vol. 57, no. 56, 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 Site  Google Scholar  Zentralblatt MATH  MathSciNet
 F. A. Oliveira, “Collocation and residual correction,” Numerische Mathematik, vol. 36, no. 1, pp. 27–31, 1980. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
Copyright
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.