Abstract and Applied Analysis

Volume 2014 (2014), Article ID 794716, 6 pages

http://dx.doi.org/10.1155/2014/794716

Research Article

Solving Singularly Perturbed Multipantograph Delay Equations Based on the Reproducing Kernel Method

F. Z. Geng and S. P. Qian

Department of Mathematics, Changshu Institute of Technology, Changshu, Jiangsu 215500, China

Received 13 November 2013; Accepted 20 January 2014; Published 27 February 2014

Academic Editor: Valery Y. Glizer

Copyright © 2014 F. Z. Geng and S. P. Qian. 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.

A numerical method is presented for solving the singularly perturbed multipantograph delay equations with a boundary layer at one end point. The original problem is reduced to boundary layer and regular domain problems. The regular domain problem is solved by combining the asymptotic expansion and the reproducing kernel method (RKM). The boundary layer problem is treated by the method of scaling and the RKM. Two numerical examples are provided to illustrate the effectiveness of the present method. The results from the numerical example show that the present method can provide very accurate analytical approximate solutions.

1. Introduction

In this paper, we consider the following singularly perturbed multipantograph delay equation:

where , , is a positive integer and and are assumed to be sufficiently smooth, such that (1) has a unique solution with a boundary layer at .

Singularly perturbed problems arise frequently in applications including geophysical fluid dynamics, oceanic and atmospheric circulation, chemical reactions, and optimal control. These problems are characterized by the presence of a small parameter that multiplies the highest order derivative, and they are stiff and there exist boundary layers where the solutions change rapidly.

Functional differential equations with proportional delays are usually referred to as pantograph equations. These equations arise in a variety of applications, such as number theory, electrodynamics, astrophysics, nonlinear dynamical systems, probability theory on algebraic structure, quantum mechanics, and cell growth.

Recently, singularly perturbed delayed differential equations have attracted significant attention. The numerical treatment of such problems presents some major computational difficulties, and therefore discussion on numerical solutions of singularly perturbed delayed differential equation is rare. Amiraliyev et al. [1, 2] proposed a uniform numerical method for dealing with singularly perturbed delay initial value problems. Kadalbajoo et al. [3–6] presented some effective methods for solving singularly perturbed delay boundary value problems. In [7], Amiraliyev and Cimen also introduced a numerical method for singularly perturbed delay boundary value problems. In [8], Rai and Sharma described a numerical method based on fitted operator finite difference scheme for the boundary value problems for singularly perturbed delay differential equations with turning point and mixed shifts.

Reproducing kernel theory has important applications in numerical analysis, differential equations, probability, and statistics, amongst other fields [9–20]. Recently, using the reproducing kernel method (RKM), the authors have discussed various differential equations [11–20]. However, it is very difficult to expand the application of the RKM to singularly perturbed delayed differential equations. Geng [13] developed a method for solving a class of singularly perturbed boundary value problems based on the RKM and a proper transformation. Nevertheless, this method fails to solve singularly perturbed delayed differential equations.

In this paper, based on the RKM presented in [9, 11], an effective numerical method will be presented for solving singularly perturbed delayed boundary value problem (1).

The rest of the paper is organized as follows. In the next section, the numerical technique for (1) is introduced. Error analysis is introduced in Section 2. The numerical example is given in Section 3. Section  4 ends this paper with a brief conclusion.

2. Numerical Method

We divide the domain into two subdomains, namely, and , where is a positive real number. The asymptotic approximation technique and the RKM are combined to solve (1) in the regular domain . And then the value of asymptotic approximation in the regular domain is used as the boundary condition at the so-called transition point . In the boundary layer domain , (1) is solved by combining the method of scaling and the RKM. After solving both the regular and boundary layer domain problems their solutions are combined to obtain an approximate solution to the original problem over the entire domain .

2.1. Solution of the Regular Domain

We seek the regular region solution as an asymptotic expansion of the form

where are unknown functions to be determined.

Substituting into (1) and equating the coefficients of like powers of , we obtain

By the RKM, the solutions of the above equations can be approximated by

where , , .

Therefore, the solution of regular region can be approximated by

In the following, we will show how to solve (3) using the RKM in detail.

Consider the following operator equation:

Under the assumption that (3) has a unique solution, we will give the approximate solution of (3) in the reproducing kernel space .

The reproducing kernel space is defined as are absolutely continuous real valued functions, . The inner product and norm in are given, respectively, by

Its reproducing kernel is

For detailed method of obtaining reproducing kernel , please refer to [9, 15].

In (6), it is clear that is a bounded linear operator (please see [9] for the definition of ). Put and , where is the reproducing kernel of , is the adjoint operator of . The orthonormal system of can be derived from the Gram-Schmidt orthogonalization process applied to ,

Theorem 1. For (6), if is dense on , then is the complete system of and .

Proof. Note here that

Clearly, .

For each fixed , let , which means that

Since is dense on , . It follows that from the existence of . So the proof of the Theorem 1 is complete.

Theorem 2. If is dense on and the solution of (6) is unique, then the solution of (6) is

Proof. Applying Theorem 1, it is easy to see that is the complete orthonormal basis of .

Note that for each ; hence we have

and the proof of the theorem is complete.

The approximate solution can be obtained by taking finitely many terms in the series representation of and

2.2. Solution of the Boundary Layer Domain

Consider

For the boundary layer domain, we scale with ; then (15) becomes

Using the similar method for solving (3), we can obtain the approximate solution of (16) in the reproducing kernel space is an absolutely continuous real valued function, , in which the inner product and norm are given, respectively, by

The approximate solution of (16) can be represented by

where . Therefore, the approximation of solution of boundary layer domain problem (15) can be obtained by

From (5) and (19), the approximate solution of (1) on the entire region is immediately obtained as follows:

3. Numerical Examples

Example 1. Consider the following singular perturbation problem with pantograph delay

where is given such that its exact solution is . Using the present method and taking , , (regular domain), (boundary layer domain), and , , , respectively, the numerical results are given in Figures 1, 2, and 3. All computations are performed by using Mathematica 5.1.

(a)

(b)

Figure 1: Absolute errors of Example 1 for ((a): regular domain; (b): boundary later domain).

(a)

(b)

Figure 2: Absolute errors of Example 1 for ((a): regular domain; (b): boundary later domain).

(a)

(b)

Figure 3: Absolute errors of Example 1 for ((a): regular domain; (b): boundary later domain).

Example 2. For comparison, we consider the following singular perturbation problem without delay [21–25]

Its exact solution is . Using the present method and taking , , (regular domain), (boundary layer domain), , the numerical results compared with other methods are given in Tables 1, 2, and 3. It is shown that the present method can yield better results than existing methods.

Table 1: Comparison of maximum absolute error with other methods for .

[21]
[22]
Present method

Table 2: Comparison of maximum absolute error with other methods for .

[21]
[22]
Present method

Table 3: Comparison of computed solution with [22–25] for Example  4.2 with and .

Exact solution
Present method
[22]
[23]
[24, 25]

0.01
−0.987875
−0.987921
−0.989854
−1.398575
−2.607798

0.03
−0.967160
−0.967161
−0.969099
−1.603239
−1.940047

0.05
−0.945600
−0.945602
−0.947500
−1.579587
−1.698810

0.07
−0.923240
−0.923242
−0.925100
−1.550175
−1.598463

0.09
−0.900080
−0.900082
−0.901900
−1.519601
−1.544179

0.10
−0.888200
−0.888202
−0.890000
−1.503912
−1.509999

0.30
−0.608600
−0.608603
−0.610000
−1.134145
−1.144100

0.50
−0.249000
−0.249005
−0.250000
−0.657717
−0.664829

0.70
0.190600
0.190594
0.190000
−0.074626
−0.078895

0.90
0.710200
0.710193
0.710000
0.615126
0.613702

1.00
1.000000
1.000000
1.000000
1.000000
1.000000

4. Conclusion

In this paper, a new method is proposed for solving singularly perturbed multipantograph delay equations. The present method is based on the RKM, the asymptotic expansion technique, and the method of scaling. The major advantage of the method is that it can produce good globally continuous approximate solutions. The results from the numerical example show that the present method is an accurate and reliable analytical technique for treating singularly perturbed multipantograph delay equations.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors would like to express their thanks to the unknown referees for their careful reading and helpful comments. The work was supported by the NSFC (Grant nos. 11201041 and 11026200), the Special Funds of the National Natural Science Foundation of China (Grant no. 11141003), and Qing Lan Project of Jiangsu Province.

References

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.

I. G. Amiraliyeva, F. Erdogan, and G. M. Amiraliyev, “A uniform numerical method for dealing with a singularly perturbed delay initial value problem,” Applied Mathematics Letters, vol. 23, no. 10, pp. 1221–1225, 2010.

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.

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, “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.

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.

G. M. Amiraliyev and E. Cimen, “Numerical method for a singularly perturbed convection-diffusion problem with delay,” Applied Mathematics and Computation, vol. 216, no. 8, pp. 2351–2359, 2010.

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. Cui and Y. Lin, Nonlinear Numerical Analysis in the Reproducing Kernel Space, Nova Science, New York, NY, USA, 2009.

A. Berlinet and C. Thomas-Agnan, Reproducing Kernel Hilbert Spaces in Probability and Statistics, Kluwer Academic, Boston, Mass, USA, 2004.

F. Z. Geng and M. Cui, “Solving a nonlinear system of second order boundary value problems,” Journal of Mathematical Analysis and Applications, vol. 327, no. 2, pp. 1167–1181, 2007.

F. Z. Geng and M. Cui, “New method based on the HPM and RKHSM for solving forced Duffing equations with integral boundary conditions,” Journal of Computational and Applied Mathematics, vol. 233, no. 2, pp. 165–172, 2009.

F. Z. Geng, “A novel method for solving a class of singularly perturbed boundary value problems based on reproducing kernel method,” Applied Mathematics and Computation, vol. 218, no. 8, pp. 4211–4215, 2011.

F. Z. Geng, S. P. Qian, and S. Li, “A numerical method for singularly perturbed turning point problems with an interior layer,” Journal of Computational and Applied Mathematics, vol. 255, pp. 97–105, 2014.

B. Wu and X. Li, “Iterative reproducing kernel method for nonlinear oscillator with discontinuity,” Applied Mathematics Letters, vol. 23, no. 10, pp. 1301–1304, 2010.

W. Wang, B. Han, and M. Yamamoto, “Inverse heat problem of determining time-dependent source parameter in reproducing kernel space,” Nonlinear Analysis: Real World Applications, vol. 14, no. 1, pp. 875–887, 2013.

Y. Wang, M. Du, F. Tan, Z. Li, and T. Nie, “Using reproducing kernel for solving a class of fractional partial differential equation with non-classical conditions,” Applied Mathematics and Computation, vol. 219, no. 11, pp. 5918–5925, 2013.

M. Inc, A. Akgül, and A. Kiliçman, “A novel method for solving KdV equation based on reproducing kernel Hilbert space method,” Abstract and Applied Analysis, vol. 2013, Article ID 578942, 11 pages, 2013.

M. Mohammadi and R. Mokhtari, “Solving the generalized regularized long wave equation on the basis of a reproducing kernel space,” Journal of Computational and Applied Mathematics, vol. 235, no. 14, pp. 4003–4014, 2011.

G. Akram and H. U. Rehman, “Numerical solution of eighth order boundary value problems in reproducing kernel space,” Numerical Algorithms, vol. 62, no. 3, pp. 527–540, 2013.

M. K. Kadalbajoo and P. Arora, “B-spline collocation method for the singular-perturbation problem using artificial viscosity,” Computers & Mathematics with Applications, vol. 57, no. 4, pp. 650–663, 2009.

A. Kaushik, V. Kumar, and A. K. Vashishth, “An efficient mixed asymptotic-numerical scheme for singularly perturbed convection diffusion problems,” Applied Mathematics and Computation, vol. 218, no. 17, pp. 8645–8658, 2012.

A. Andargie and Y. N. Reddy, “Fitted fourth-order tridiagonal finite difference method for singular perturbation problems,” Applied Mathematics and Computation, vol. 192, no. 1, pp. 90–100, 2007.

M. M. Chawla, “A fourth-order tridiagonal finite difference method for general non-linear two-point boundary value problems with mixed boundary conditions,” IMA Journal of Applied Mathematics, vol. 21, no. 1, pp. 83–93, 1978.

M. M. Chawla and P. N. Shivakumar, “An efficient finite difference method for two-point boundary value problems,” Neural, Parallel & Scientific Computations, vol. 4, no. 3, pp. 387–395, 1996.