Abstract and Applied Analysis

Volume 2012, Article ID 925134, 13 pages

http://dx.doi.org/10.1155/2012/925134

## Numerical Algorithms for Computing Eigenvalues of Discontinuous Dirac System Using Sinc-Gaussian Method

^{1}Department of Mathematics, Faculty of Science, King Abdulaziz University, Jeddah 21589, Saudi Arabia^{2}Department of Mathematics, Faculty of Science, Beni-Suef University, Beni-Suef 62511, Egypt

Received 8 March 2012; Accepted 26 May 2012

Academic Editor: Zhenya Yan

Copyright © 2012 A. H. Bhrawy et al. 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

The eigenvalues of a discontinuous regular Dirac systems with transmission conditions at the point of discontinuity are computed using the sinc-Gaussian method. The error analysis of this method for solving discontinuous regular Dirac system is discussed. It shows that the error decays exponentially in terms of the number of involved samples. Therefore, the accuracy of the new method is higher than the classical sinc-method. Numerical results indicating the high accuracy and effectiveness of these algorithms are presented. Comparisons with the classical sinc-method are given.

#### 1. Introduction

The mathematical modeling of many practical problems in mechanics and other areas of mathematical physics requires solutions of boundary value problems (see, e.g., [1–7]). Boundary value problems with discontinuity conditions arise in different branches of mathematics, mechanics, radio, electronics, geophysics, and other fields of natural science and technology. For example, discontinuous conditions inside an interval are connected with discontinuous or nonsmooth properties of media (see [8–11] and references there in).

Sampling theory is one of the most important mathematical tools used in communication engineering since it enables engineers to reconstruct signals from some of their sampled data. A fundamental result in information theory is the Whittaker-Kotel'nikov-Shannon (WKS) sampling theorem [12–14]. It states that any , where is the space of all entire functions of exponential type which lie in when restricted to , can be reconstructed from its sampled values by the formula where Series (1.1) converges absolutely and uniformly on compact subsets of , and uniform on , cf. [15]. Expansion (1.1) is used in several approximation problems which are known as sinc methods, see, for example, [16–19]. In particular the sinc-method is used to approximate eigenvalues of boundary value problems, see, for examples, [20–24]. The sinc-method has a slow rate of decay at infinity, which is as slow as . There are several attempts to improve the rate of decay. One of the interesting ways is to multiply the sinc-function in (1.1) by a kernel function, see, for example, [25–27]. Let and . Assume that such that , then for we have the expansion [28], The speed of convergence of the series in (1.3) is determined by the decay of . But the decay of an entire function of exponential type cannot be as fast as as , for some positive , [28]. In [29], Qian has introduced the following regularized sampling formula. For , and , where , which is called the Gaussian function, , and denotes the integer part of , see also [30, 31]. Qian also derived the following error bound. If , and , then In [28], Schmeisser and Stenger extended the operator (1.4) to the complex domain . For , and , they defined the operator [28], where and . Note that the summation limits in (1.6) depend on the real part of . Schmeisser and Stenger, [28], proved that if is an entire function of exponential type , then for , , , , we have where

The amplitude error arises when the exact values of (1.6) are replaced by the approximations . We assume that are close to , that is, there is , sufficiently small such that Let , and be fixed numbers. The authors in [9] proved that if (1.9) is held, then for , we have where

We are concerned with the computation of eigenvalues the Dirac system with boundary conditions and transmission conditions where ; the real valued function and are continuous in and , and have finite limits , ; , and .

In this paper we will use the sinc-Gaussian sampling formula (1.6) to compute eigenvalues of the Dirac system (1.12)–(1.16). As expected, the new method reduced the error bounds remarkably, see examples at the end of this paper. Special attention is given to the comparison of the numerical results obtained by the new method with those found by classical sinc-method. We would like to mention that works in direction of computing eigenvalues with the sinc-Gaussian, are few, see for example, [9]. Also papers in computing of eigenvalues with discontinuous are few, see [10, 32]. However, the computing of eigenvalues by sinc-Gaussian technique which has discontinuity conditions, do not exist as for as we know.

The paper is organized as follows: Section 2 contains some preliminary results and the approximated values of the eigenvalues of the Dirac system with discontinuous. The method with error estimates are contained in Section 3. The last section involves some illustrative examples for showing the high accuracy of the proposed technique.

#### 2. The Approximated Eigenvalues of Dirac System

In this section we derive approximate values of the eigenvalues of problem (1.12)–(1.16). Recall that problem (1.12)–(1.16) has denumerable set of real and simple eigenvalues, compare with [33, 34]. Let be the solution of (1.12) satisfying the following initial conditions: Since satisfies (1.13), then the eigenvalues of the problem (1.12)–(1.16) are the zeros of the function Notice that both and are entire functions of and satisfies the system of integral equations, see [32, 35], where , , and , , are the Volterra integral operators defined by

For convenience, we define the constants Define and , , to be

Lemma 2.1. *The functions and are entire in for any fixed and satisfy the growth condition
*

*Proof. *Since , then from (2.4) and (2.5) we obtain
Using the inequalities and for , leads for to
The above inequality can be reduced to
Similarly, we can prove that
Then from (2.12) and (2.15) and Lemma 3.1 of [34, pp. 204], we obtain (2.12).

In a similar manner, we will prove the following lemma for and .

Lemma 2.2. *The functions and are entire in for any fixed and satisfy the growth condition
*

*Proof. *Since , then from (2.6) and (2.7) we obtain
Then from (2.4) and (2.5) and Lemma 2.1, we get
Similarly, we can prove that

#### 3. The Method and Error Analysis

In this section we derive the method of computing the eigenvalues of problem (1.12)–(1.16) numerically. We aim to approximate and hence its zeros, that is, the eigenvalues. The idea is to split into two parts, one is known and the other is unknown, but is an entire function of exponential type. Then we approximate the unknown part using (1.6) to get the approximate and then compute the approximate zeros. Now, let us split into where is the unknown part involving integral operators and is the known part Then, from Lemmas 2.1 and 2.2, we have the following result.

Lemma 3.1. *The function is entire in and the following estimate holds:
**
where
*

*Proof. *From (3.2), we have
Using the inequalities and for , and Lemmas 2.1 and 2.2 imply (3.4).

Then is an entire function of exponential type . In the following we let since all eigenvalues are real. Now we approximate the function using the operator (1.6) where and and then we obtain where The samples , cannot be computed explicitly in the general case. We approximate these samples numerically by solving the initial-value problems defined by (1.12) and (2.2) to obtain the approximate values , , that is, . Accordingly we have the explicit expansion Therefore we get, compare with (1.10), Now let . From (3.7) and (3.10) we obtain Let be an eigenvalue and be its desired approximation, that is, and . From (3.11) we have . Now we define an enclosure interval for . Define the curves The curves , trap the curve of for suitably large . Hence the closure interval is determined by solving , gives an interval . Next we estimate the error for the eigenvalue .

Theorem 3.2. *Let be an eigenvalue of (1.12)–(1.16) and be its approximation. Then, for , one has the following estimate:
**
where the interval is defined above.*

*Proof. *Replacing by in (3.11) we obtain
where we have used . Using the mean value theorem yields that for some ,
Since is simple and is sufficiently large, then and we get (3.13).

#### 4. Numerical Examples and Comparisons

This section includes two examples illustrating the sinc-Gaussian method. In the following examples, we consider being the *k*th root of . Also, it is observed that the approximation and the exact solution are all inside the interval . We indicate in these two examples the effect of the amplitude error in the proposed method by determining enclosure intervals for different values of . All examples are computed in [32] by using the classical sinc method. We see that the sinc-Gaussian method gives remarkably better results.

*Example 4.1. *Consider the system
subject to
Here
and .

Direct calculations give
therefore the eigenvalues are, .

Let and denote the absolute errors associated with the results of classical sinc method and sinc-Gaussian method, respectively. In Table 1, we give comparison between the absolute error of sinc-Gaussian and the classical sinc-method.

In Table 2, we observe that the approximation and the exact solution are all inside the interval for different values of .

*Example 4.2. *In this example we consider the system
subject to
where
, and .

Direct calculations give
therefore the eigenvalues are , .

Table 3 gives comparison between the absolute error of sinc-Gaussian and the classical sinc-method. Moreover, in Table 4 the approximation and the exact solution are all inside the interval for different values of .

#### Acknowledgment

This paper was funded by the Deanship of Scientific Research (DSR), King Abdulaziz University, Jeddah. The authors, therefore, acknowledge with thanks DSR technical and financial support.

#### References

- R. P. Agarwal,
*Boundary Value Problems for Higher Order Differential Equations*, World Scientific, Singapore, Singapore, 1986. View at Zentralblatt MATH - A. H. Bhrawy and A. S. Alofi, “A Jacobi-Gauss collocation method for solving nonlinear Lane-Emden type equations,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 17, no. 1, pp. 62–70, 2012. View at Publisher · View at Google Scholar - A. H. Bhrawy, A. S. Alofi, and S. I. El-Soubhy, “Spectral shifted Jacobi tau and collocation methods for solving fifth-order boundary value problems,”
*Abstract and Applied Analysis*, Article ID 823273, 14 pages, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - A. H. Bhrawy and W. M. Abd-Elhameed, “New algorithm for the numerical solutions of nonlinear third-order differential equations using Jacobi-Gauss collocation method,”
*Mathematical Problems in Engineering*, Article ID 837218, 14 pages, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - E. H. Doha and W. M. Abd-Elhameed, “Efficient solutions of multidimensional sixth-order
boundary value problems using symmetric generalized Jacobi-Galerkin method,”
*Abstract and Applied Analysis*, vol. 2012, Article ID Article ID 749370, 19 pages, 2012. View at Publisher · View at Google Scholar - E. H. Doha, A. H. Bhrawy, and R. M. Hafez, “A Jacobi dual-Petrov-Galerkin method for solving some odd-order ordinary differential equations,”
*Abstract and Applied Analysis*, Article ID 947230, 21 pages, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - A. Imani, A. Aminataei, and A. Imani, “Collocation method via Jacobi polynomials for solving nonlinear ordinary differential equations,”
*International Journal of Mathematics and Mathematical Sciences*, Article ID 673085, 11 pages, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - O. N. Litvinenko and V. I. Soshnikov,
*The Theory of Heterogeneous Lines and their Applications in Radio Engineering*, Radio, Moscow, Russia, 1964. - M. H. Annaby and R. M. Asharabi, “Computing eigenvalues of boundary-value problems using sinc-Gaussian method,”
*Sampling Theory in Signal and Image Processing*, vol. 7, no. 3, pp. 293–311, 2008. View at Google Scholar · View at Zentralblatt MATH - M. H. Annaby and R. M. Asharabi, “Approximating eigenvalues of discontinuous problems by sampling theorems,”
*Journal of Numerical Mathematics*, vol. 16, no. 3, pp. 163–183, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - R. Kh. Amirov, “On a system of Dirac differential equations with discontinuity conditions inside an interval,”
*Ukrainian Mathematical Journal*, vol. 57, no. 5, pp. 601–613, 2005. View at Publisher · View at Google Scholar - V. Kotelnikov,
*On the Carrying Capacity of the Either and Wire in Telecommunications, Material for the First All-Union Conference on Questions of Communications*, Izd. Red. Upr. Svyazi RKKA, Moscow, Russia, 1933. - C. E. Shannon, “Communication in the presence of noise,” vol. 37, pp. 10–21, 1949. View at Google Scholar
- E. T. Whittaker, “On the functions which are represented by the expansion of the interpolation
theory,”
*Proceedings of the Royal Society of Edinburgh A*, vol. 35, pp. 181–194, 1915. View at Google Scholar - P. L. Butzer, G. Schmeisser, and R. L. Stens, “An introduction to sampling analysis,” in
*Non Uniform Sampling: Theory and Practices*, F. Marvasti, Ed., pp. 17–121, Kluwer, New York, NY, USA, 2001. View at Google Scholar - M. A. Kowalski, K. A. Sikorski, and F. Stenger,
*Selected Topics in Approximation and Computation*, SIAM, Philadelphia, Pa, USA; Oxford University Press, New York, NY, USA, 1995. - J. Lund and K. L. Bowers,
*Sinc Methods for Quadrature and Differential Equations*, Society for Industrial and Applied Mathematics (SIAM), Philadelphia, Pa, USA, 1992. View at Publisher · View at Google Scholar - F. Stenger, “Numerical methods based on Whittaker cardinal, or sinc functions,”
*SIAM Review*, vol. 23, no. 2, pp. 165–224, 1981. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - F. Stenger,
*Numerical Methods Based On Sinc and Analytic Functions*, Springer-Verlag, New York, NY, USA, 1993. View at Publisher · View at Google Scholar - M. H. Annaby and M. M. Tharwat, “On computing eigenvalues of second-order linear pencils,”
*IMA Journal of Numerical Analysis*, vol. 27, no. 2, pp. 366–380, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - M. H. Annaby and M. M. Tharwat, “Sinc-based computations of eigenvalues of Dirac systems,”
*Behaviour & Information Technology*, vol. 47, no. 4, pp. 699–713, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - M. H. Annaby and M. M. Tharwat, “On the computation of the eigenvalues of Dirac systems,”
*Calcolo*. In press. View at Publisher · View at Google Scholar - A. Boumenir, “Higher approximation of eigenvalues by the sampling method,”
*Behaviour & Information Technology*, vol. 40, no. 2, pp. 215–225, 2000. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - A. Boumenir, “Sampling and eigenvalues of non-self-adjoint Sturm-Liouville problems,”
*SIAM Journal on Scientific Computing*, vol. 23, no. 1, pp. 219–229, 2001. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - P. L. Butzer and R. L. Stens, “A modification of the Whittaker-Kotelnikov-Shannon sampling series,”
*Aequationes Mathematicae*, vol. 28, no. 3, pp. 305–311, 1985. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - R. Gervais, Q. I. Rahman, and G. Schmeisser, “A bandlimited function simulating a duration-limited one,” in
*Anniversary Volume on Approximation Theory and Functional Analysis*, pp. 355–362, Birkhäuser, Basel, Switzerland, 1984. View at Google Scholar · View at Zentralblatt MATH - R. L. Stens, “Sampling by generalized kernels,” in
*Sampling Theory in Fourier and Signal Analysis: Advanced Topics*, J. R. Higgins and R. L. Stens, Eds., pp. 130–157, Oxford University Press, Oxford, UK, 1999. View at Google Scholar - G. Schmeisser and F. Stenger, “Sinc approximation with a Gaussian multiplier,”
*Sampling Theory in Signal and Image Processing*, vol. 6, no. 2, pp. 199–221, 2007. View at Google Scholar · View at Zentralblatt MATH - L. Qian, “On the regularized Whittaker-Kotel'nikov-Shannon sampling formula,”
*Proceedings of the American Mathematical Society*, vol. 131, no. 4, pp. 1169–1176, 2003. View at Publisher · View at Google Scholar - L. Qian and D. B. Creamer, “A modification of the sampling series with a Gaussian multiplier,”
*Sampling Theory in Signal and Image Processing*, vol. 5, no. 1, pp. 1–19, 2006. View at Google Scholar · View at Zentralblatt MATH - L. Qian and D. B. Creamer, “Localized sampling in the presence of noise,”
*Applied Mathematics Letters*, vol. 19, no. 4, pp. 351–355, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - M. M. Tharwat, A. H. Bhrawy, and A. Yildirim, “Computing eigenvalues of discontinuous Dirac system by using sinc method,”
*International Journal of Computer Mathematics*. In press. View at Publisher · View at Google Scholar - B. M. Levitan and I. S. Sargsjan, “Introduction to spectral theory: self adjoint ordinary differential operators,” in
*Translation of Mthematical Monographs*, vol. 39, American Mathematical Society, Providence, RI, USA, 1975. View at Google Scholar - B. M. Levitan and I. S. Sargsjan,
*Sturm-Liouville and Dirac Operators*, Kluwer Academic, Dordrecht, The Netherlands, 1991. - M. M. Tharwat, “Discontinuous Sturm-Liouville problems and associated sampling theories,”
*Abstract and Applied Analysis*, Article ID 610232, 30 pages, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH