Abstract
In this paper, we apply a sincGaussian technique to compute approximate values of the eigenvalues of SturmLiouville problems which contain an eigenparameter appearing linearly in two boundary conditions, in addition to an internal point of discontinuity. The error of this method decays exponentially in terms of the number of involved samples. Therefore the accuracy of the new technique is higher than that of the classical sinc method. Numerical worked examples with tables and illustrative figures are given at the end of the paper.
MSC: 34L16, 94A20, 65L15.
Keywords:
sampling theory; SturmLiouville problems; transmission conditions; sincGaussian; sinc method; truncation and amplitude errors1 Introduction
By a sampling theorem we mean a representation of a certain function in terms of its values at a discrete set of points. In communication theory, it means a reconstruction of a signal (information) in terms of a discrete set of data. This has several applications, especially in the transmission of information. If the signal is bandlimited, the sampling process can be done via the celebrated WhittakerKotel’nikovShannon (WKS) sampling theorem [13]. By a bandlimited signal with band width τ, , we mean a function in the PaleyWiener space
The WKS sampling theorem is a fundamental result in information theory. It states that any can be reconstructed from its sampled values , where and , by the formula
where
and the series converges absolutely and uniformly on any finite interval of ℝ. Expansion (1.2) is used in several approximation problems which are known as sinc methods; see, e.g., [47]. In particular the sincmethod is used to approximate eigenvalues of boundary value problems; see, for example, [814]. The sincmethod 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 sincfunction in (1.2) by a kernel function; see, e.g., [1517]. Let and . Assume that such that , then for we have the expansion, [18]
The speed of convergence of the series in (1.4) 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 c[18]. In [19], Qian has introduced the following regularized sampling formula. For , and , Qian defined the operator [19]
where , which is called the Gaussian function, , and denotes the integer part of ; see also [20,21]. Qian also derived the following error bound. If , and , then [19,20]
In [18] Schmeisser and Stenger extended the operator (1.5) to the complex domain ℂ. For , and , they defined the operator [18]
where and . Note that the summation limits in (1.7) depend on the real part of z. Schmeisser and Stenger [18] proved that if f is an entire function such that
where ϕ is a nondecreasing, nonnegative function on and , then for , , , , we have
where
The amplitude error arises when the exact values of (1.7) are replaced by the approximations . We assume that are close to , i.e., there is sufficiently small such that
Let , and be fixed numbers. The authors in [22] proved that if (1.11) holds, then for , we have
where
It is well known that many topics in mathematical physics require the investigation of the eigenvalues and eigenfunctions of SturmLiouville type boundary value problems. Therefore, the Sturmian theory is one of the most actual and extensively developing fields of theoretical and applied mathematics. Particularly, in recent years, highly important results in this field have been obtained for the case when the eigenparameter appears not only in the differential equation but also in the boundary conditions. The literature on such results is voluminous, and we refer to [2327] and corresponding bibliography cited therein. In particular, [24,26,28,29] contain many references to problems in physics and mechanics. Our task is to use formula (1.7) to compute the eigenvalues numerically of the differential equation
with boundary conditions
and transmission conditions
where μ is a complex spectral parameter; is a given realvalued function, which is continuous in and and has a finite limit ; , , , , , () are real numbers; , (); and
The eigenvalue problem (1.14)(1.18) when is a SturmLiouville problem which contains an eigenparameter μ in two boundary conditions, in addition to an internal point of discontinuity. In [30], Tharwat proved that the eigenvalue problem (1.14)(1.18) has a denumerable set of real and simple eigenvalues using techniques similar to those established in [23,24,31], where also sampling theorems have been established. Tharwat et al., in [14], computed the eigenvalues of the problem (1.14)(1.18) by using the sinc method. In the sinc method, the basic idea is as follows: The eigenvalues are characterized as the zeros of an analytic function which can be written in the form , where (known part) is the function for the case . The ingenuity of the approach is in trying to choose the function so that (unknown part) and can be approximated by the WKS sampling theorem if its values at some equally spaced points are known; see [814].
Our goal in this paper is to improve the results presented in Tharwat et al.[14] with the least conditions. In this paper we use the sincGaussian sampling formula (1.7) to compute eigenvalues of (1.14)(1.18) numerically. As is expected, the new method reduced the error bounds remarkably; see the examples at the end of this paper. Also here, we use the same idea but the unknown part is an entire function of exponential type and satisfies (1.8), that is, is not necessary function. Then we approximate the using (1.7) and obtain better results. We would like to mention that the papers in computing eigenvalues by the sincGaussian method are few; see [22,32,33]. In Section 2 we derive the sincGaussian technique to compute the eigenvalues of (1.14)(1.18) with error estimates. The last section involves some illustrative examples.
2 Treatment of the eigenvalue problem (1.14)(1.18)
In this section we derive approximate values of the eigenvalues of the eigenvalue problem (1.14)(1.18). Recall that the problem (1.14)(1.18) has a denumerable set of real and simple eigenvalues, cf.[30]. Let
denote the solution of (1.14) satisfying the following initial conditions:
Since satisfies (1.15), (1.17) and (1.18), then the eigenvalues of problem (1.14)(1.18) are the zeros of the characteristic determinant, cf.[30],
According to [30], see also [3440], the function is an entire function of μ where zeros are real and simple. We aim to approximate and hence its zeros, i.e., the eigenvalues, by using (1.7). The idea is to split into two parts, one is known and the other is unknown, but is an entire function of exponential type and satisfies (1.8). Then we approximate the unknown part using (1.7) to get the approximate and then compute the approximate zeros. By using the method of variation of constants, we can see that the solution satisfies the Volterra integral equations, cf.[30],
where and are the Volterra operators
Differentiating (2.3) and (2.4), we obtain
where and are the Volterratype integral operators
In the following, we make use of the known estimates [41]
where is some constant (we may take cf.[41]). For convenience, we define the constants
As in [14], we split into two parts via
Then the function is entire in μ for each for which, cf.[14],
where
and
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.7) where and and then, from (1.9), we obtain
where
The samples , cannot be computed explicitly in the general case. We approximate these samples numerically by solving the initialvalue problems defined by (1.14) and (2.1) to obtain the approximate values , , i.e., . Here we use a computer algebra system, MATHEMATICA, to obtain the approximate solutions with the required accuracy. However, a separate study for the effect of different numerical schemes and the computational costs would be interesting. Accordingly, we have the explicit expansion
Therefore we get, cf. (1.12),
Now let . From (2.20) and (2.23) we obtain
Let be an eigenvalue and be its desired approximation, i.e., and . From (2.24) we have . Define the curves
The curves , trap the curve of for suitably large N. Hence the closure interval is determined by solving , which gives an interval
It is worthwhile to mention that the simplicity of the eigenvalues guarantees the existence of approximate eigenvalues, i.e., the ’s for which . Next we estimate the error for the eigenvalue .
Theorem 2.1Letbe an eigenvalue of (1.14)(1.18) andbe its approximation. Then, for, we have the following estimate:
where the intervalis defined above.
Proof Replacing μ by in (2.24), we obtain
where we have used . Using the mean value theorem yields that for some ,
Since is simple and N is sufficiently large, then , and we get (2.26). □
3 Examples
This section includes two examples illustrating the sincGaussian method. All examples are computed in [14] with the classical sincmethod. It is clearly seen that the sincGaussian method gives remarkably better results. We indicate in these two examples the effect of the amplitude error in the method by determining enclosure intervals for different values of ε. We also indicate the effect of N and h by several choices. We would like to mention that MATHEMATICA has been used to obtain the exact values for these examples where eigenvalues cannot computed concretely. MATHEMATICA is also used in rounding the exact eigenvalues, which are square roots. Each example is exhibited via figures that accurately illustrate the procedure near to some of the approximated eigenvalues. More explanations are given below.
Example 1 Consider the boundary value problem
The characteristic function is
As is clearly seen, eigenvalues cannot be computed explicitly. Tables 1, 2, 3 indicate the application of our technique to this problem and the effect of ε. By exact we mean the zeros of computed by MATHEMATICA.
Table 1. The approximationand the exact solutionfor different choices ofhandN
Table 2. Absolute error
Table 3. Forand, the exact solutionsare all inside the intervalfor different values ofε
Figures 1 and 2 illustrate the enclosure intervals dominating for , and , , respectively. The middle curve represents , while the upper and lower curves represent the curves of , , respectively. We notice that when , all two curves are almost identical. Similarly, Figures 3 and 4 illustrate the enclosure intervals dominating for , and , , respectively.
Figure 1. The enclosure interval dominatingfor,and.
Figure 2. The enclosure interval dominatingfor,and.
Figure 3. The enclosure interval dominatingfor,and.
Figure 4. The enclosure interval dominatingfor,and.
Example 2 Consider the boundary value problem
The characteristic determinant of the problem is
where and are Airy functions, and and are derivatives of Airy functions. As in the above example, the three tables (Tables 4, 5, 6) indicate the application of our technique to this problem and the effect of ε.
Table 4. The approximationand the exact solutionfor different choices ofhandN
Table 5. Absolute error
Table 6. Forand, the exact solutionsare all inside the intervalfor different values ofε
Here Figures {5, 6}, {7, 8} illustrate the enclosure intervals dominating and for , and , , respectively.
4 Conclusion
With a simple analysis, and with values of solutions of initial value problems computed at a few values of the eigenparameter, we have computed the eigenvalues of discontinuous SturmLiouville problems which contain an eigenparameter appearing linearly in two boundary conditions, with a certain estimated error. The method proposed is a shooting procedure, i.e., the problem is reformulated as two initial value ones, due to the interior discontinuity, of size two and a missdistance is defined at the right end of the interval of integration whose roots are the eigenvalues to be computed. The unknown part of the missdistance can be written in terms of a function which is an entire function of exponential type. Therefore, we propose to approximate such term by means of a truncated cardinal series with sampling values approximated by solving numerically corresponding suitable initial value problems. Finally, in Section 3 we introduced two instructive examples. The computations show that, as compared to the classical sampling expansion in [14], the variant with the Gaussian multiplier provides a strikingly high improvement of the accuracy.
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
The authors have equal contributions to each part of this article. All the authors read and approved the final manuscript.
Acknowledgements
This work was funded by the Deanship of Scientific Research (DSR), King Abdulaziz University, Jeddah, under grant No. (130065D1433). The authors, therefore, acknowledge with thanks DSR technical and financial support.
References

Kotel’nikov, V: On the carrying capacity of the ‘ether’ and wire in telecommunications. Material for the First AllUnion Conference on Questions of Communications, pp. 55–64. Izd. Red. Upr. Svyazi RKKA, Moscow (1933)

Shannon, CE: Communications in the presence of noise. Proc. IRE. 37, 10–21 (1949)

Whittaker, ET: On the functions which are represented by the expansion of the interpolation theory. Proc. R. Soc. Edinb., Sect. A. 35, 181–194 (1915)

Stenger, F: Numerical methods based on Whittaker cardinal, or sinc functions. SIAM Rev.. 23, 156–224 (1981)

Lund, J, Bowers, K: Sinc Methods for Quadrature and Differential Equations, SIAM, Philadelphia (1992)

Stenger, F: Numerical Methods Based on Sinc and Analytic Functions, Springer, New York (1993)

Kowalski, M, Sikorski, K, Stenger, F: Selected Topics in Approximation and Computation, Oxford University Press, New York (1995)

Boumenir, A: Higher approximation of eigenvalues by sampling. BIT Numer. Math.. 40, 215–225 (2000). Publisher Full Text

Boumenir, A: Sampling and eigenvalues of nonselfadjoint SturmLiouville problems. SIAM J. Sci. Comput.. 23, 219–229 (2001). Publisher Full Text

Annaby, MH, Tharwat, MM: On computing eigenvalues of secondorder linear pencils. IMA J. Numer. Anal.. 27, 366–380 (2007)

Annaby, MH, Tharwat, MM: Sincbased computations of eigenvalues of Dirac systems. BIT Numer. Math.. 47, 699–713 (2007). Publisher Full Text

Annaby, MH, Tharwat, MM: On the computation of the eigenvalues of Dirac systems. Calcolo. 49, 221–240 (2012). Publisher Full Text

Tharwat, MM, Bhrawy, AH, Yildirim, A: Numerical computation of eigenvalues of discontinuous Dirac system using Sinc method with error analysis. Int. J. Comput. Math.. 89, 2061–2080 (2012). Publisher Full Text

Tharwat, MM, Bhrawy, AH, Yildirim, A: Numerical computation of eigenvalues of discontinuous SturmLiouville problems with parameter dependent boundary conditions using Sinc method. Numer. Algorithms. 63, 27–48 (2013). Publisher Full Text

Gervais, R, Rahman, QI, Schmeisser, G: A bandlimited function simulating a durationlimited one. In: Butzer PL, Stens RL (eds.) Approximation Theory and Functional Analysis, pp. 355–362. Birkhäuser, Basel (1984)

Butzer, PL, Stens, RL: A modification of the WhittakerKotel’nikovShannon sampling series. Aequ. Math.. 28, 305–311 (1985). Publisher Full Text

Stens, RL: Sampling by generalized kernels. In: Higgins JR, Stens RL (eds.) Sampling Theory in Fourier and Signal Analysis: Advanced Topics, pp. 130–157. Oxford University Press, Oxford (1999)

Schmeisser, G, Stenger, F: Sinc approximation with a Gaussian multiplier. Sampl. Theory Signal. Image Process, Int. J.. 6, 199–221 (2007)

Qian, L: On the regularized WhittakerKotel’nikovShannon sampling formula. Proc. Am. Math. Soc.. 131, 1169–1176 (2002)

Qian, L, Creamer, DB: A modification of the sampling series with a Gaussian multiplie. Sampl. Theory Signal. Image Process, Int. J.. 5, 1–20 (2006)

Qian, L, Creamer, DB: Localized sampling in the presence of noise. Appl. Math. Lett.. 19, 351–355 (2006). Publisher Full Text

Annaby, MH, Asharabi, RM: Computing eigenvalues of boundary value problems using sincGaussian method. Sampl. Theory Signal. Image Process, Int. J.. 7, 293–312 (2008)

Walter, J: Regular eigenvalue problems with eigenvalue parameter in the boundary condition. Math. Z.. 133, 301–312 (1973). Publisher Full Text

Fulton, CT: Twopoint boundary value problems with eigenvalue parameter contained in the boundary conditions. Proc. R. Soc. Edinb., Sect. A. 77, 293–308 (1977). Publisher Full Text

Hinton, DB: An expansion theorem for an eigenvalue problem with eigenvalue parameter in the boundary condition. Q. J. Math.. 30, 33–42 (1979). Publisher Full Text

Shkalikov, AA: Boundary value problems for ordinary differential equations with a parameter in boundary conditions. Tr. Semin. Im. I.G. Petrovskogo. 9, 190–229 (in Russian) (1983)

Binding, PA, Browne, PJ, Watson, BA: StrumLiouville problems with boundary conditions rationally dependent on the eigenparameter II. J. Comput. Appl. Math.. 148, 147–169 (2002). Publisher Full Text

Likov, AV, Mikhailov, YA: The Theory of Heat and Mass Transfer, Qosenergaizdat, MoscowLeningrad (1963) (in Russian)

Tikhonov, AN, Samarskii, AA: Equations of Mathematical Physics, Macmillan Co., New York (1963)

Tharwat, MM: Discontinuous SturmLiouville problems and associated sampling theories. Abstr. Appl. Anal. doi:10.1155/2011/610232 (2011)

Titchmarsh, EC: Eigenfunction Expansions Associated with Second Order Differential Equations. Part I, Clarendon, Oxford (1962)

Bhrawy, AH, Tharwat, MM, AlFhaid, A: Numerical algorithms for computing eigenvalues of discontinuous Dirac system using sincGaussian method. Abstr. Appl. Anal. doi:10.1155/2012/925134 (2012)

Annaby, MH, Tharwat, MM: A sincGaussian technique for computing eigenvalues of secondorder linear pencils. Appl. Numer. Math.. 63, 129–137 (2013)

Annaby, MH, Tharwat, MM: On sampling theory and eigenvalue problems with an eigenparameter in the boundary conditions. SUT J. Math.. 42, 157–176 (2006)

Annaby, MH, Tharwat, MM: On sampling and Dirac systems with eigenparameter in the boundary conditions. J. Appl. Math. Comput.. 36, 291–317 (2011). Publisher Full Text

Kandemir, M, Mukhtarov, OS: Discontinuous Sturm Liouville problems containing eigenparameter in the boundary conditions. Acta Math. Sin.. 34, 1519–1528 (2006)

Mukhtarov, OS, Kadakal, M, Altinisik, N: Eigenvalues and eigenfunctions of discontinuous SturmLiouville problems with eigenparameter in the boundary conditions. Indian J. Pure Appl. Math.. 34, 501–516 (2003)

Tharwat, MM, Bhrawy, AH: Computation of eigenvalues of discontinuous Dirac system using Hermite interpolation technique. Adv. Differ. Equ. doi:10.1186/16871847201259 (2012)

Tharwat, MM, Yildirim, A, Bhrawy, AH: Sampling of discontinuous Dirac systems. Numer. Funct. Anal. Optim.. 34, 323–348 (2013). Publisher Full Text

Tharwat, MM: On sampling theories and discontinuous Dirac systems with eigenparameter in the boundary conditions. Bound. Value Probl. doi:10.1186/16872770201365 (2013)

Chadan, K, Sabatier, PC: Inverse Problems in Quantum Scattering Theory, Springer, Berlin (1989)