SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

Open Access Research

Approximation of eigenvalues of discontinuous Sturm-Liouville problems with eigenparameter in all boundary conditions

Mohammed M Tharwat1*, Ali H Bhrawy12 and Abdulaziz S Alofi1

Author Affiliations

1 Department of Mathematics, Faculty of Science, King Abdulaziz University, Jeddah, Saudi Arabia

2 Permanent address: Department of Mathematics, Faculty of Science, Beni-Suef University, Beni-Suef, Egypt

For all author emails, please log on.

Boundary Value Problems 2013, 2013:132  doi:10.1186/1687-2770-2013-132

The electronic version of this article is the complete one and can be found online at: http://www.boundaryvalueproblems.com/content/2013/1/132


Received:12 March 2013
Accepted:4 May 2013
Published:20 May 2013

© 2013 Tharwat et al.; licensee Springer

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

In this paper, we apply a sinc-Gaussian technique to compute approximate values of the eigenvalues of Sturm-Liouville 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; Sturm-Liouville problems; transmission conditions; sinc-Gaussian; sinc method; truncation and amplitude errors

1 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 band-limited, the sampling process can be done via the celebrated Whittaker-Kotel’nikov-Shannon (WKS) sampling theorem [1-3]. By a band-limited signal with band width τ, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M1">View MathML</a>, we mean a function in the Paley-Wiener space

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M2">View MathML</a>

(1.1)

The WKS sampling theorem is a fundamental result in information theory. It states that any <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M3">View MathML</a> can be reconstructed from its sampled values <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M4">View MathML</a>, where <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M5">View MathML</a> and <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M6">View MathML</a>, by the formula

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M7">View MathML</a>

(1.2)

where

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M8">View MathML</a>

(1.3)

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., [4-7]. In particular the sinc-method is used to approximate eigenvalues of boundary value problems; see, for example, [8-14]. The sinc-method has a slow rate of decay at infinity, which is as slow as <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M9">View MathML</a>. There are several attempts to improve the rate of decay. One of the interesting ways is to multiply the sinc-function in (1.2) by a kernel function; see, e.g., [15-17]. Let <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M10">View MathML</a> and <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M11">View MathML</a>. Assume that <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M12">View MathML</a> such that <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M13">View MathML</a>, then for <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M3">View MathML</a> we have the expansion, [18]

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M15">View MathML</a>

(1.4)

The speed of convergence of the series in (1.4) is determined by the decay of <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M16">View MathML</a>. But the decay of an entire function of exponential type cannot be as fast as <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M17">View MathML</a> as <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M18">View MathML</a>, for some positive c[18]. In [19], Qian has introduced the following regularized sampling formula. For <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M10">View MathML</a>, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M20">View MathML</a> and <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M21">View MathML</a>, Qian defined the operator [19]

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M22">View MathML</a>

(1.5)

where <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M23">View MathML</a>, which is called the Gaussian function, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M24">View MathML</a>, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M25','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M25">View MathML</a> and <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M26','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M26">View MathML</a> denotes the integer part of <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M27">View MathML</a>; see also [20,21]. Qian also derived the following error bound. If <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M28','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M28">View MathML</a>, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M10">View MathML</a> and <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M30','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M30">View MathML</a>, then [19,20]

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M31','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M31">View MathML</a>

(1.6)

In [18] Schmeisser and Stenger extended the operator (1.5) to the complex domain ℂ. For <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M32','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M32">View MathML</a>, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M10">View MathML</a> and <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M34','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M34">View MathML</a>, they defined the operator [18]

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M35','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M35">View MathML</a>

(1.7)

where <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M36">View MathML</a> and <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M20">View MathML</a>. 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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M38','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M38">View MathML</a>

(1.8)

where ϕ is a non-decreasing, non-negative function on <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M39','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M39">View MathML</a> and <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M40','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M40">View MathML</a>, then for <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M41">View MathML</a>, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M42">View MathML</a>, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M20">View MathML</a>, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M44">View MathML</a>, we have

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M45','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M45">View MathML</a>

(1.9)

where

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M46','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M46">View MathML</a>

(1.10)

The amplitude error arises when the exact values <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M47','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M47">View MathML</a> of (1.7) are replaced by the approximations <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M48','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M48">View MathML</a>. We assume that <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M48','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M48">View MathML</a> are close to <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M47','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M47">View MathML</a>, i.e., there is <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M51','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M51">View MathML</a> sufficiently small such that

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M52','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M52">View MathML</a>

(1.11)

Let <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M53','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M53">View MathML</a>, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M42">View MathML</a> and <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M20">View MathML</a> be fixed numbers. The authors in [22] proved that if (1.11) holds, then for <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M56','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M56">View MathML</a>, we have

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M57','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M57">View MathML</a>

(1.12)

where

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M58','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M58">View MathML</a>

(1.13)

It is well known that many topics in mathematical physics require the investigation of the eigenvalues and eigenfunctions of Sturm-Liouville 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 [23-27] 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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M59','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M59">View MathML</a>

(1.14)

with boundary conditions

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M60','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M60">View MathML</a>

(1.15)

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M61','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M61">View MathML</a>

(1.16)

and transmission conditions

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M62','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M62">View MathML</a>

(1.17)

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M63','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M63">View MathML</a>

(1.18)

where μ is a complex spectral parameter; <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M64','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M64">View MathML</a> is a given real-valued function, which is continuous in <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M65','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M65">View MathML</a> and <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M66','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M66">View MathML</a> and has a finite limit <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M67','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M67">View MathML</a>; <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M68','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M68">View MathML</a>, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M69','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M69">View MathML</a>, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M70','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M70">View MathML</a>, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M71','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M71">View MathML</a>, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M72','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M72">View MathML</a>, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M73','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M73">View MathML</a> (<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M74','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M74">View MathML</a>) are real numbers; <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M75','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M75">View MathML</a>, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M76','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M76">View MathML</a> (<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M74','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M74">View MathML</a>); <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M78','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M78">View MathML</a> and

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M79','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M79">View MathML</a>

(1.19)

The eigenvalue problem (1.14)-(1.18) when <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M80','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M80">View MathML</a> is a Sturm-Liouville 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 <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M81','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M81">View MathML</a> which can be written in the form <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M82','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M82">View MathML</a>, where <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M83','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M83">View MathML</a> (known part) is the function for the case <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M84','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M84">View MathML</a>. The ingenuity of the approach is in trying to choose the function <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M81','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M81">View MathML</a> so that <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M86','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M86">View MathML</a> (unknown part) and can be approximated by the WKS sampling theorem if its values at some equally spaced points are known; see [8-14].

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 sinc-Gaussian 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 <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M87','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M87">View MathML</a> is an entire function of exponential type and satisfies (1.8), that is, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M87','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M87">View MathML</a> is not necessary <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M89','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M89">View MathML</a>-function. Then we approximate the <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M87','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M87">View MathML</a> using (1.7) and obtain better results. We would like to mention that the papers in computing eigenvalues by the sinc-Gaussian method are few; see [22,32,33]. In Section 2 we derive the sinc-Gaussian 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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M91','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M91">View MathML</a>

denote the solution of (1.14) satisfying the following initial conditions:

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M92','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M92">View MathML</a>

(2.1)

Since <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M93','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M93">View MathML</a> 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],

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M94','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M94">View MathML</a>

(2.2)

According to [30], see also [34-40], the function <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M95','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M95">View MathML</a> is an entire function of μ where zeros are real and simple. We aim to approximate <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M96','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M96">View MathML</a> and hence its zeros, i.e., the eigenvalues, by using (1.7). The idea is to split <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M96','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M96">View MathML</a> 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 <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M96','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M96">View MathML</a> and then compute the approximate zeros. By using the method of variation of constants, we can see that the solution <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M93','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M93">View MathML</a> satisfies the Volterra integral equations, cf.[30],

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M100','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M100">View MathML</a>

(2.3)

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M101','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M101">View MathML</a>

(2.4)

where <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M102','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M102">View MathML</a> and <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M103','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M103">View MathML</a> are the Volterra operators

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M104','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M104">View MathML</a>

(2.5)

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M105','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M105">View MathML</a>

(2.6)

Differentiating (2.3) and (2.4), we obtain

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M106','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M106">View MathML</a>

(2.7)

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M107','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M107">View MathML</a>

(2.8)

where <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M108','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M108">View MathML</a> and <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M109','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M109">View MathML</a> are the Volterra-type integral operators

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M110','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M110">View MathML</a>

(2.9)

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M111','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M111">View MathML</a>

(2.10)

Define <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M112','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M112">View MathML</a> and <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M113','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M113">View MathML</a>, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M74','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M74">View MathML</a>, to be

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M115','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M115">View MathML</a>

(2.11)

In the following, we make use of the known estimates [41]

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M116','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M116">View MathML</a>

(2.12)

where <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M117','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M117">View MathML</a> is some constant (we may take <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M118','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M118">View MathML</a>cf.[41]). For convenience, we define the constants

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M119','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M119">View MathML</a>

(2.13)

As in [14], we split <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M96','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M96">View MathML</a> into two parts via

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M121','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M121">View MathML</a>

(2.14)

where <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M122','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M122">View MathML</a> is the known part

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M123','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M123">View MathML</a>

(2.15)

and <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M124','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M124">View MathML</a> is the unknown one

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M125','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M125">View MathML</a>

(2.16)

Then the function <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M124','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M124">View MathML</a> is entire in μ for each <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M127','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M127">View MathML</a> for which, cf.[14],

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M128','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M128">View MathML</a>

(2.17)

where

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M129','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M129">View MathML</a>

(2.18)

and

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M130','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M130">View MathML</a>

(2.19)

Then <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M124','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M124">View MathML</a> is an entire function of exponential type <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M132','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M132">View MathML</a>. In the following, we let <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M133','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M133">View MathML</a> since all eigenvalues are real. Now we approximate the function <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M124','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M124">View MathML</a> using the operator (1.7) where <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M135','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M135">View MathML</a> and <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M136','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M136">View MathML</a> and then, from (1.9), we obtain

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M137','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M137">View MathML</a>

(2.20)

where

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M138','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M138">View MathML</a>

(2.21)

The samples <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M139','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M139">View MathML</a>, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M140','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M140">View MathML</a> cannot be computed explicitly in the general case. We approximate these samples numerically by solving the initial-value problems defined by (1.14) and (2.1) to obtain the approximate values <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M141','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M141">View MathML</a>, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M140','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M140">View MathML</a>, i.e., <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M143','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M143">View MathML</a>. 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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M144','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M144">View MathML</a>

(2.22)

Therefore we get, cf. (1.12),

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M145','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M145">View MathML</a>

(2.23)

Now let <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M146','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M146">View MathML</a>. From (2.20) and (2.23) we obtain

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M147','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M147">View MathML</a>

(2.24)

Let <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M148','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M148">View MathML</a> be an eigenvalue and <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M149','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M149">View MathML</a> be its desired approximation, i.e., <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M150','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M150">View MathML</a> and <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M151','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M151">View MathML</a>. From (2.24) we have <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M152','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M152">View MathML</a>. Define the curves

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M153','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M153">View MathML</a>

(2.25)

The curves <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M154','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M154">View MathML</a>, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M155','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M155">View MathML</a> trap the curve of <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M96','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M96">View MathML</a> for suitably large N. Hence the closure interval is determined by solving <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M157','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M157">View MathML</a>, which gives an interval

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M158','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M158">View MathML</a>

It is worthwhile to mention that the simplicity of the eigenvalues guarantees the existence of approximate eigenvalues, i.e., the <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M149','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M149">View MathML</a>’s for which <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M151','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M151">View MathML</a>. Next we estimate the error <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M161','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M161">View MathML</a> for the eigenvalue <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M162','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M162">View MathML</a>.

Theorem 2.1Let<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M148','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M148">View MathML</a>be an eigenvalue of (1.14)-(1.18) and<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M149','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M149">View MathML</a>be its approximation. Then, for<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M133','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M133">View MathML</a>, we have the following estimate:

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M166','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M166">View MathML</a>

(2.26)

where the interval<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M167','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M167">View MathML</a>is defined above.

Proof Replacing μ by <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M149','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M149">View MathML</a> in (2.24), we obtain

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M169','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M169">View MathML</a>

(2.27)

where we have used <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M170','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M170">View MathML</a>. Using the mean value theorem yields that for some <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M171','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M171">View MathML</a>,

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M172','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M172">View MathML</a>

(2.28)

Since <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M162','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M162">View MathML</a> is simple and N is sufficiently large, then <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M174','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M174">View MathML</a>, and we get (2.26). □

3 Examples

This section includes two examples illustrating the sinc-Gaussian method. All examples are computed in [14] with the classical sinc-method. It is clearly seen that the sinc-Gaussian 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

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M175','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M175">View MathML</a>

(3.1)

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M176','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M176">View MathML</a>

(3.2)

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M177','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M177">View MathML</a>

(3.3)

Here <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M178','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M178">View MathML</a>, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M179','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M179">View MathML</a>, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M180','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M180">View MathML</a>, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M181','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M181">View MathML</a> and

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M182','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M182">View MathML</a>

(3.4)

The characteristic function is

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M183','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M183">View MathML</a>

(3.5)

The function <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M122','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M122">View MathML</a> will be

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M185','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M185">View MathML</a>

(3.6)

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 <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M96','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M96">View MathML</a> computed by MATHEMATICA.

Table 1. The approximation<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M187','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M187">View MathML</a>and the exact solution<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M188','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M188">View MathML</a>for different choices ofhandN

Table 2. Absolute error<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M196','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M196">View MathML</a>

Table 3. For<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M202','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M202">View MathML</a>and<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M203','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M203">View MathML</a>, the exact solutions<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M204','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M204">View MathML</a>are all inside the interval<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M205','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M205">View MathML</a>for different values ofε

Figures 1 and 2 illustrate the enclosure intervals dominating <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M214','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M214">View MathML</a> for <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M202','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M202">View MathML</a>, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M203','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M203">View MathML</a> and <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M217','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M217">View MathML</a>, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M218','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M218">View MathML</a>, respectively. The middle curve represents <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M96','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M96">View MathML</a>, while the upper and lower curves represent the curves of <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M154','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M154">View MathML</a>, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M155','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M155">View MathML</a>, respectively. We notice that when <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M218','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M218">View MathML</a>, all two curves are almost identical. Similarly, Figures 3 and 4 illustrate the enclosure intervals dominating <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M223','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M223">View MathML</a> for <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M203','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M203">View MathML</a>, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M202','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M202">View MathML</a> and <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M217','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M217">View MathML</a>, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M218','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M218">View MathML</a>, respectively.

thumbnailFigure 1. The enclosure interval dominating<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M228','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M228">View MathML</a>for<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M229','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M229">View MathML</a>,<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M230','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M230">View MathML</a>and<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M231','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M231">View MathML</a>.

thumbnailFigure 2. The enclosure interval dominating<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M228','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M228">View MathML</a>for<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M229','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M229">View MathML</a>,<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M230','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M230">View MathML</a>and<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M235','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M235">View MathML</a>.

thumbnailFigure 3. The enclosure interval dominating<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M236','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M236">View MathML</a>for<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M229','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M229">View MathML</a>,<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M230','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M230">View MathML</a>and<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M231','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M231">View MathML</a>.

thumbnailFigure 4. The enclosure interval dominating<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M236','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M236">View MathML</a>for<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M229','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M229">View MathML</a>,<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M230','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M230">View MathML</a>and<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M235','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M235">View MathML</a>.

Example 2 Consider the boundary value problem

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M244','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M244">View MathML</a>

(3.7)

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M245','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M245">View MathML</a>

(3.8)

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M246','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M246">View MathML</a>

(3.9)

where <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M247','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M247">View MathML</a>, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M248','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M248">View MathML</a>, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M249','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M249">View MathML</a>, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M250','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M250">View MathML</a>, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M251','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M251">View MathML</a> and

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M252','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M252">View MathML</a>

(3.10)

The function <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M122','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M122">View MathML</a> will be

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M254','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M254">View MathML</a>

(3.11)

The characteristic determinant of the problem is

<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M255','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M255">View MathML</a>

(3.12)

where <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M256','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M256">View MathML</a> and <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M257','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M257">View MathML</a> are Airy functions, and <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M258','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M258">View MathML</a> and <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M259','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M259">View MathML</a> 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 approximation<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M187','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M187">View MathML</a>and the exact solution<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M188','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M188">View MathML</a>for different choices ofhandN

Table 5. Absolute error<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M196','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M196">View MathML</a>

Table 6. For<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M275','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M275">View MathML</a>and<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M276','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M276">View MathML</a>, the exact solutions<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M204','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M204">View MathML</a>are all inside the interval<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M205','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M205">View MathML</a>for different values ofε

Here Figures {5, 6}, {7, 8} illustrate the enclosure intervals dominating <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M287','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M287">View MathML</a> and <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M288','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M288">View MathML</a> for <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M276','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M276">View MathML</a>, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M275','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M275">View MathML</a> and <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M217','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M217">View MathML</a>, <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M218','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M218">View MathML</a>, respectively.

thumbnailFigure 5. The enclosure interval dominating<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M293','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M293">View MathML</a>for<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M294','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M294">View MathML</a>,<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M295','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M295">View MathML</a>and<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M231','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M231">View MathML</a>.

thumbnailFigure 6. The enclosure interval dominating<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M293','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M293">View MathML</a>for<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M294','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M294">View MathML</a>,<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M295','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M295">View MathML</a>and<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M235','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M235">View MathML</a>.

thumbnailFigure 7. The enclosure interval dominating<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M301','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M301">View MathML</a>for<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M294','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M294">View MathML</a>,<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M295','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M295">View MathML</a>and<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M231','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M231">View MathML</a>.

thumbnailFigure 8. The enclosure interval dominating<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M301','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M301">View MathML</a>for<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M294','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M294">View MathML</a>,<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M295','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M295">View MathML</a>and<a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M235','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M235">View MathML</a>.

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 Sturm-Liouville 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 miss-distance is defined at the right end of the interval of integration whose roots are the eigenvalues to be computed. The unknown part <a onClick="popup('http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M124','MathML',630,470);return false;" target="_blank" href="http://www.boundaryvalueproblems.com/content/2013/1/132/mathml/M124">View MathML</a> of the miss-distance 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. (130-065-D1433). The authors, therefore, acknowledge with thanks DSR technical and financial support.

References

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

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

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

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

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

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

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

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

  9. Boumenir, A: Sampling and eigenvalues of non-self-adjoint Sturm-Liouville problems. SIAM J. Sci. Comput.. 23, 219–229 (2001). Publisher Full Text OpenURL

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

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

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

  13. 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 OpenURL

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

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

  16. Butzer, PL, Stens, RL: A modification of the Whittaker-Kotel’nikov-Shannon sampling series. Aequ. Math.. 28, 305–311 (1985). Publisher Full Text OpenURL

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

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

  19. Qian, L: On the regularized Whittaker-Kotel’nikov-Shannon sampling formula. Proc. Am. Math. Soc.. 131, 1169–1176 (2002)

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

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

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

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

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

  25. 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 OpenURL

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

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

  28. Likov, AV, Mikhailov, YA: The Theory of Heat and Mass Transfer, Qosenergaizdat, Moscow-Leningrad (1963) (in Russian)

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

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

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

  32. Bhrawy, AH, Tharwat, MM, Al-Fhaid, A: Numerical algorithms for computing eigenvalues of discontinuous Dirac system using sinc-Gaussian method. Abstr. Appl. Anal. doi:10.1155/2012/925134 (2012)

  33. Annaby, MH, Tharwat, MM: A sinc-Gaussian technique for computing eigenvalues of second-order linear pencils. Appl. Numer. Math.. 63, 129–137 (2013)

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

  35. 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 OpenURL

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

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

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

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

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

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