Abstract and Applied Analysis

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

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

Research Article

Exponential Collocation Method for Solutions of Singularly Perturbed Delay Differential Equations

Şuayip Yüzbaşı¹ and Mehmet Sezer²

¹Department of Mathematics, Faculty of Science, Akdeniz University, Antalya 07058, Turkey

²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 .

Table 1: Estimated maximum absolute errors for various values of , , and of (38).

Estimated absolute errors

(a)

(b)

(c)

Figure 1: For and various values of and of (38),    (a) graph of the approximate solutions for , (b) graph of the approximate solutions for , and (c) graph of the approximate solutions for .

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 .

Table 2: Estimated maximum absolute errors for various values of , , and of (40).

Estimated absolute errors

(a)

(b)

(c)

(d)

Figure 2: Graphs of the approximate solutions of (40), (a) for , , and various values of and , (b) for , , and various values of and , (c) for , , and various values of and , and (d) for , and , , , and .

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).

(a)

(b)

Figure 3: For (42), (a) comparison of the exact solution and the approximate solutions for , , and various values of and and (b) comparison of the exact solution for , and , and the approximate solutions for , and , , , and .

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)

Figure 4: For (45), (a) graph of the approximate solutions for , , and various values of and , (b) comparison of the estimated error functions for , , and various values of and , (c) graph of the approximate solutions for ,  , , , , , and , and (d) graph of the estimated error functions for , , , and , , , and .

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.

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.

M. C. Mackey and L. Glass, “Oscillation and chaos in physiological control systems,” Science, vol. 197, no. 4300, pp. 287–289, 1977.

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.

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.

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.

K. C. Patidar and K. K. Sharma, “-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.

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.

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.

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.

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.

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.

M. K. Kadalbajoo and D. Kumar, “Fitted mesh -spline collocation method for singularly perturbed differential-difference equations with small delay,” Applied Mathematics and Computation, vol. 204, no. 1, pp. 90–98, 2008.

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.

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.

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.

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.

M. K. Kadalbajoo and K. K. Sharma, “An -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.

F. Alharbi, “Predefined exponential basis set for half-bounded multi domain spectral method,” Applied Mathematics, vol. 1, pp. 146–152, 2010.

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.

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.

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.

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.

L. Székelyhidi, “On the extension of exponential polynomials,” Mathematica Bohemica, vol. 125, no. 3, pp. 365–370, 2000.

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.

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.

Ş. Yüzbaşı, “An efficient algorithm for solving multi-pantograph equation systems,” Computers & Mathematics with Applications, vol. 64, no. 4, pp. 589–603, 2012.

Ş. 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.

İ. Çelik, “Collocation method and residual correction using Chebyshev series,” Applied Mathematics and Computation, vol. 174, no. 2, pp. 910–920, 2006.

F. A. Oliveira, “Collocation and residual correction,” Numerische Mathematik, vol. 36, no. 1, pp. 27–31, 1980.