Skip to main content

Sturm-Liouville problem with moving discontinuity points

Abstract

In this paper, we present a new discontinuous Sturm-Liouville problem with symmetrically located discontinuities which are defined depending on a parameter in the neighborhood of an interior point in the interval. Also the problem contains an eigenparameter in a boundary condition. We investigate some spectral properties of the eigenvalues, obtain asymptotic formulae for the eigenvalues and the corresponding eigenfunctions and construct Green’s function for the problem. We give an illustrative example with tables and figures at the end of the paper.

1 Introduction

The Sturm-Liouville theory plays an important role in solving many mathematical physics problems [1]. Such research is motivated by the theory of heat and mass transfer, or vibrating string problems when the string is loaded additionally with point masses (see [2, 3]; see also the references therein). When using Fourier’s method for heat transfer problems, it becomes necessary to solve a related Sturm-Liouville problem in which the separation constant plays the role of a spectral parameter (or eigenparameter). In the case of inhomogeneous materials and composite walls, the Sturm-Liouville equation has variable and not necessarily continuous coefficients, so that transmission conditions across the interfaces should be added in the problem. Heat conduction problems in composite walls have been analysed by several researchers [46]. In all of these works, the temperature distribution of composite walls involving two or more layers are investigated spectrally and the results are formulated by the eigenvalues and eigenfunctions of the auxiliary spectral problem. In [7], the presence of an infinite number of eigenvalues has been shown and an asymptotic formula has been obtained for the eigenvalues of the spectral problem for the temperature distribution in composite walls. The problem in that work included a two-layer composite wall consisting of different materials, having a common contact surface. Physically, we can say that the problem in the present paper is derived from a heat conduction problem of a three-layer composite wall consisting of two interfaces that are located symmetrically. Thus we add transmission conditions at these interfaces i.e. at the points of discontinuity to the problem.

Eigenvalue problems have been studied by some authors in the continuous case; see [810]. In recent years, there were also some studies on the discontinuous eigenvalue problems [1119] and construction of Green’s functions for discontinuous Sturm-Liouville problems which contain an eigenparameter under one or two boundary conditions [14, 16].

We consider the following Sturm-Liouville problem:

$$ \tau ( y ): = - y\prime\prime + q ( x )y = \lambda y,\quad x \in I, $$
(1)

with an eigenparameter dependent on a boundary condition:

$$\begin{aligned}& B_{a} ( y ): = \beta_{1}y ( a ) + \beta_{2}y' ( a ) = 0, \end{aligned}$$
(2)
$$\begin{aligned}& B_{b} ( y ): = \lambda \bigl( \alpha_{1}'y ( b ) - \alpha_{2}'y' ( b ) \bigr) + \bigl( \alpha_{1}y ( b ) - \alpha_{2}y' ( b ) \bigr) = 0, \end{aligned}$$
(3)

and coupled transmission conditions at the points of discontinuity, \(\theta_{ - \varepsilon}\) and \(\theta_{ + \varepsilon}\):

$$\begin{aligned}& T_{ - \varepsilon} ( y ): = y ( \theta_{ - \varepsilon} + ) - \bigl( \mu_{1}y ( \theta_{ - \varepsilon} - ) + \mu_{2}y' ( \theta_{ - \varepsilon} - ) \bigr) = 0, \end{aligned}$$
(4)
$$\begin{aligned}& T_{ - \varepsilon} ' ( y ): = y' ( \theta_{ - \varepsilon} + ) - \bigl( \mu_{1}'y ( \theta_{ - \varepsilon} - ) + \mu_{2}'y' ( \theta_{ - \varepsilon} - ) \bigr) = 0, \end{aligned}$$
(5)
$$\begin{aligned}& T_{ + \varepsilon} ( y ): = y ( \theta_{ + \varepsilon} + ) - \bigl( \eta_{1}y ( \theta_{ + \varepsilon} - ) + \eta_{2}y' ( \theta_{ + \varepsilon} - ) \bigr) = 0, \end{aligned}$$
(6)
$$\begin{aligned}& T_{ + \varepsilon} ' ( y ): = y' ( \theta_{ + \varepsilon} + ) - \bigl( \eta_{1}'y ( \theta_{ + \varepsilon} - ) + \eta_{2}'y' ( \theta_{ + \varepsilon} - ) \bigr) = 0, \end{aligned}$$
(7)

where \(I: = [ a,\theta_{ - \varepsilon} ) \cup ( \theta_{ - \varepsilon},\theta_{ + \varepsilon} ) \cup ( \theta_{ + \varepsilon},b ]\); \(\theta \in ( a,b )\) and

$$ 0 < \varepsilon < \min \{ \theta - a,b - \theta \}, $$
(8)

λ is a spectral parameter; \(q ( x )\) is a given real valued function which is continuous in \([ a,\theta_{ - \varepsilon} )\), \(( \theta_{ - \varepsilon},\theta_{ + \varepsilon} )\) and \(( \theta_{ + \varepsilon},b ]\) and has finite limits \(q ( \theta_{ - \varepsilon} \pm )\), \(q ( \theta_{ + \varepsilon} \pm )\); \(\beta_{i}\), \(\alpha_{i}\), \(\alpha_{i}'\), \(\mu_{i}\), \(\mu_{i}'\), \(\eta_{i}\), \(\eta_{i}'\) (\(i = 1,2 \)) are real numbers such that \(\vert \beta_{1} \vert + \vert \beta_{2} \vert \ne 0\) and

$$ \begin{aligned} &\rho: = \det \begin{pmatrix} \alpha_{1}' & \alpha_{2}' \\ \alpha_{1} & \alpha_{2} \end{pmatrix} > 0, \\ &D_{1}: = \det \begin{pmatrix} \mu_{1} & \mu_{2} \\ \mu_{1}' & \mu_{2}' \end{pmatrix} > 0,\quad\quad D_{2}: = \det \begin{pmatrix} \eta_{1} & \eta_{2} \\ \eta_{1}' & \eta_{2}' \end{pmatrix} > 0. \end{aligned} $$
(9)

Also for convenience we will use the notations \(\theta_{ \pm \varepsilon} \pm: = ( \theta \pm \varepsilon ) \pm 0\).

We introduce a new Sturm-Liouville problem with discontinuities which are defined depending on a parameter in the neighborhood of an interior point of θ. ε is a parameter controlling the variation of the neighborhood process and by using the variation of this parameter, it is possible to determine points of discontinuity. Thus, the points of discontinuity can be placed anywhere. Accordingly, they can be called ‘moving discontinuity points.’ In the special case, when our problem is not with an eigenparameter in the boundary condition and discontinuities are \(d_{1} = d\) and \(d_{2} = \pi - d\) in the interval \([ 0,\pi ]\), was derived in [11, 12].

The main aim of the present work is to present moving discontinuity points, to emphasise some of the distinguishing properties of such problems from the ones of discontinuous Sturm-Liouville problems, and to study the spectral properties of the problem. To achieve our aim we extend some classic results of Sturm-Liouville theory to the new moving discontinuous case. First, we define a linear operator A in a suitable Hilbert space H such that the eigenvalues of the problem (1)-(7) coincide with those of A and construct a special fundamental system of solutions. Then we obtain the asymptotic formulae for the eigenvalues and the corresponding eigenfunctions depending on the parameter ε and construct Green’s function for the problem (1)-(7). Finally, an illustrative example, which shows how to determine discontinuities for different values ofε, is given.

2 An operator formulation

In this section we will introduce the special inner product in the Hilbert space \(L_{2} ( a,b ) \oplus \mathbb{C}\) and a symmetric linear operator A is defined in this Hilbert space in such a way that the problem (1)-(7) can be considered as the eigenvalue problem of this operator.

Definition 1

We define a Hilbert space H of two component vectors by \(H: = L_{2} ( a,b ) \oplus \mathbb{C}\) with the inner product:

$$\begin{aligned} \langle U,V \rangle_{H} : =& \int_{a}^{\theta_{ - \varepsilon}} u ( x )\overline{v} ( x )\,dx + \frac{1}{D_{1}} \int_{\theta_{ - \varepsilon}}^{\theta_{ + \varepsilon}} u ( x ) \overline{v} ( x )\,dx + \frac{1}{D_{1}D_{2}} \int_{\theta_{ + \varepsilon}}^{b} u ( x )\overline{v} ( x )\,dx \\ &{}+ \frac{1}{\rho D_{1}D_{2}}h\overline{k}, \end{aligned}$$
(10)

where \(U = \binom{u ( x ) }{ h} \), \(V = \binom{v ( x )}{ k } \in H\), \(u ( x ), v ( x ) \in L_{2} ( a,b )\) and \(h, k \in \mathbb{C}\).

For a function \(u ( x )\), which is defined on I and has finite limits \(u ( \theta_{ - \varepsilon} \pm ): = \lim_{x \to \theta_{ - \varepsilon} \pm} u ( x )\), \(u ( \theta_{ + \varepsilon} \pm ): = \lim_{x \to \theta_{ + \varepsilon} \pm} u ( x )\), by \(u_{ ( i )} ( x ) \) (\(i = \overline{1,4} \)) we denote the function

$$\begin{aligned}& u_{ ( 1 )} ( x ): = \textstyle\begin{cases} u ( x ),& x \in [ a,\theta_{ - \varepsilon} ), \\ u ( \theta_{ - \varepsilon} - ),& x = \theta_{ - \varepsilon}, \end{cases}\displaystyle \quad\quad u_{ ( 2 )} ( x ): = \textstyle\begin{cases} u ( \theta_{ - \varepsilon} + ),& x = \theta_{ - \varepsilon}, \\ u ( x ),& x \in ( \theta_{ - \varepsilon},\theta_{ + \varepsilon} ), \end{cases}\displaystyle \\& u_{ ( 3 )} ( x ): = \textstyle\begin{cases} u ( x ),& x \in ( \theta_{ - \varepsilon},\theta_{ + \varepsilon} ), \\ u ( \theta_{ + \varepsilon} - ),& x = \theta_{ + \varepsilon}, \end{cases}\displaystyle \quad\quad u_{ ( 4 )} ( x ): = \textstyle\begin{cases} u ( \theta_{ + \varepsilon} + ),& x = \theta_{ + \varepsilon}, \\ u ( x ),& x \in ( \theta_{ + \varepsilon},b ], \end{cases}\displaystyle \end{aligned}$$

which are defined on \(I_{1} = [ a,\theta_{ - \varepsilon} ]\), \(I_{2} = [ \theta_{ - \varepsilon},\theta_{ + \varepsilon} ] \) and \(I_{3} = [ \theta_{ + \varepsilon},b ]\), respectively.

For convenience we will use the notations

$$ R ( u ): = \alpha_{1}u ( b ) - \alpha_{2}u' ( b ),\quad\quad R' ( u ): = \alpha_{1}'u ( b ) - \alpha_{2}'u' ( b ). $$
(11)

Definition 2

We define a linear operator \(A: D ( A ) \to H\) by

$$ A \begin{pmatrix} u ( x ) \\ R' ( u ) \end{pmatrix}: = \begin{pmatrix} \tau ( u ) \\ - R ( u ) \end{pmatrix}, $$
(12)

where the domain \(D ( A )\) of the linear operator A is defined as the set of all \(U = \binom{ u ( x ) }{ R' ( u ) }\) which satisfies the conditions (i) \(\tau ( u ) \in L_{2} ( a,b )\), (ii) \(u_{ ( i )} (\cdot )\), \(u_{ ( i )}' (\cdot )\) are absolutely continuous functions in \(I_{i} \) (\(i = 1,2,3 \)), (iii) \(B_{a} ( u ) = 0\), (iv) \(T_{ \pm \varepsilon} ( u ) = T_{ \pm \varepsilon} ' ( u ) = 0\).

Now we can rewrite the problem (1)-(7) in the operator form as \(AU = \lambda U\) where \(U = \binom{ u ( x ) }{ R' ( u ) }\in D ( A )\).

The eigenvalues and eigenfunctions of the problem (1)-(7) are defined as the eigenvalues and the first components of the corresponding eigenelements of the operator A, respectively.

Lemma 1

The operator A in H is symmetric.

Proof

For \(U,V \in D ( A )\),

$$\begin{aligned} \langle AU,V \rangle_{H} =& \int_{a}^{\theta_{ - \varepsilon}} \tau ( u )\overline{v} ( x )\,dx + \frac{1}{D_{1}} \int_{\theta_{ - \varepsilon}}^{\theta_{ + \varepsilon}} \tau ( u )\overline{v} ( x )\,dx + \frac{1}{D_{1}D_{2}} \int_{\theta_{ + \varepsilon}}^{b} \tau ( u )\overline{v} ( x )\,dx \\ &{}- \frac{1}{\rho D_{1}D_{2}}R ( u )R' ( \overline{v} ). \end{aligned}$$
(13)

By two partial integrations, we obtain

$$\begin{aligned} \langle AU,V \rangle_{H} =& \langle U,AV \rangle_{H} + W ( u,\overline{v};\theta_{ - \varepsilon} - ) - W ( u,\overline{v};a ) + \frac{1}{D_{1}}W ( u,\overline{v};\theta_{ + \varepsilon} - ) \\ &{}- \frac{1}{D_{1}}W ( u,\overline{v};\theta_{ - \varepsilon} + ) + \frac{1}{D_{1}D_{2}}W ( u,\overline{v};b ) - \frac{1}{D_{1}D_{2}}W ( u,\overline{v}; \theta_{ + \varepsilon} + ) \\ &{}- \frac{1}{\rho D_{1}D_{2}} \bigl( R ( u )R' ( \overline{v} ) - R' ( u )R ( \overline{v} ) \bigr), \end{aligned}$$
(14)

where, as usual, by \(W ( u,v;x )\) we denote the Wronskian of the functions \(u ( x )\) and \(v ( x )\):

$$W ( u,v;x ) = u ( x )v' ( x ) - u' ( x )v ( x ). $$

Since u and satisfy (2), it follows that

$$ W ( u,\overline{v};a ) = 0, $$
(15)

from (4)-(7), we get

$$\begin{aligned}& W ( u,\overline{v};\theta_{ - \varepsilon} + ) = D_{1}W ( u, \overline{v};\theta_{ - \varepsilon} - ), \end{aligned}$$
(16)
$$\begin{aligned}& W ( u,\overline{v};\theta_{ + \varepsilon} + ) = D_{2}W ( u, \overline{v};\theta_{ + \varepsilon} - ). \end{aligned}$$
(17)

Further, from (11) it is easy to verify that

$$ R(u)R'(\overline{v}) - R'(u)R(\overline{v}) = \rho W ( u,\overline{v};b ). $$
(18)

Finally, substituting (15)-(18) into (14), we have

$$\langle AU,V \rangle_{H} = \langle U,AV \rangle_{H}, $$

so A is symmetric. □

Corollary 1

All eigenvalues of the problem (1)-(7) are real.

3 Construction of fundamental solutions

We will define the two solutions

$$ \phi_{\lambda} ( x ) = \textstyle\begin{cases} \phi_{ - \varepsilon,\lambda} ( x ),& x \in [ a,\theta_{ - \varepsilon} ), \\ \phi_{\varepsilon,\lambda} ( x ),& x \in ( \theta_{ - \varepsilon},\theta_{ + \varepsilon} ), \\ \phi_{ + \varepsilon,\lambda} ( x ),& x \in ( \theta_{ + \varepsilon},b ], \end{cases}\displaystyle \quad\quad\chi_{\lambda} ( x ) = \textstyle\begin{cases} \chi_{ - \varepsilon,\lambda} ( x ),& x \in [ a,\theta_{ - \varepsilon} ), \\ \chi_{\varepsilon,\lambda} ( x ),& x \in ( \theta_{ - \varepsilon},\theta_{ + \varepsilon} ), \\ \chi_{ + \varepsilon,\lambda} ( x ),& x \in ( \theta_{ + \varepsilon},b ], \end{cases} $$
(19)

of (1) as follows: Let \(\phi_{ - \varepsilon,\lambda} ( x ) = \phi_{ - \varepsilon} ( x,\lambda )\) be the solution of (1) on \([ a,\theta_{ - \varepsilon} ]\), which satisfies the initial conditions

$$ y ( a ) = \beta_{2},\quad\quad y' ( a ) = - \beta_{1}. $$
(20)

By virtue of Theorem 1.5 in [20], after defining this solution, we define the solution \(\phi_{\varepsilon,\lambda} ( x ) = \phi_{\varepsilon} ( x,\lambda )\) of (1) on \([ \theta_{ - \varepsilon},\theta_{ + \varepsilon} ]\) by means of the solution \(\phi_{ - \varepsilon,\lambda} ( x )\) by the nonstandard initial conditions

$$ \begin{aligned} &y ( \theta_{ - \varepsilon} + ) = \mu_{1}\phi_{ - \varepsilon,\lambda} ( \theta_{ - \varepsilon} - ) + \mu_{2}\phi_{ - \varepsilon,\lambda} ' ( \theta_{ - \varepsilon} - ), \\ &y' ( \theta_{ - \varepsilon} + ) = \mu_{1}' \phi_{ - \varepsilon,\lambda} ( \theta_{ - \varepsilon} - ) + \mu_{2}' \phi_{ - \varepsilon,\lambda} ' ( \theta_{ - \varepsilon} - ). \end{aligned} $$
(21)

After defining this solution, we may define the solution \(\phi_{ + \varepsilon,\lambda} ( x ) = \phi_{ + \varepsilon} ( x,\lambda )\) of (1) on \([ \theta_{ + \varepsilon},b ]\) by means of the solution \(\phi_{\varepsilon,\lambda} ( x )\) by the nonstandard initial conditions

$$ \begin{aligned} &y ( \theta_{ + \varepsilon} + ) = \eta_{1}\phi_{\varepsilon,\lambda} ( \theta_{ + \varepsilon} - ) + \eta_{2}\phi_{\varepsilon,\lambda} ' ( \theta_{ + \varepsilon} - ), \\ &y' ( \theta_{ + \varepsilon} + ) = \eta_{1}' \phi_{\varepsilon,\lambda} ( \theta_{ + \varepsilon} - ) + \eta_{2}' \phi_{\varepsilon,\lambda} ' ( \theta_{ + \varepsilon} - ). \end{aligned} $$
(22)

Hence, \(\phi_{\lambda} ( x ) = \phi ( x,\lambda )\) satisfies of (1), (2) and (4)-(7) on I.

Analogously, first we define the solution \(\chi_{ + \varepsilon,\lambda} ( x ) = \chi_{ + \varepsilon} ( x,\lambda )\) on \([ \theta_{ + \varepsilon},b ]\) by the initial conditions

$$ y ( b ) = \alpha_{2}'\lambda + \alpha_{2},\quad\quad y' ( b ) = \alpha_{1}'\lambda + \alpha_{1}. $$
(23)

Again, after defining this solution, we define the solution \(\chi_{\varepsilon,\lambda} ( x ) = \chi_{\varepsilon} ( x,\lambda )\) of (1) on \([ \theta_{ - \varepsilon},\theta_{ + \varepsilon} ]\) by the initial conditions

$$ \begin{aligned} &y ( \theta_{ + \varepsilon} - ) = \frac{1}{D_{2}} \bigl( \eta_{2}' \chi_{ + \varepsilon,\lambda} ( \theta_{ + \varepsilon} + ) - \eta_{2} \chi_{ + \varepsilon,\lambda} ' ( \theta_{ + \varepsilon} + ) \bigr), \\ &y' ( \theta_{ + \varepsilon} - ) = \frac{1}{D_{2}} \bigl( - \eta_{1}'\chi_{ + \varepsilon,\lambda} ( \theta_{ + \varepsilon} + ) + \eta_{1}\chi_{ + \varepsilon,\lambda} ' ( \theta_{ + \varepsilon} + ) \bigr). \end{aligned} $$
(24)

After defining this solution, we define the solution \(\chi_{ - \varepsilon,\lambda} ( x ) = \chi_{ - \varepsilon} ( x,\lambda )\) of (1) on \([ a,\theta_{ - \varepsilon} ]\) by the initial conditions

$$ \begin{aligned} &y ( \theta_{ - \varepsilon} - ) = \frac{1}{D_{1}} \bigl( \mu_{2}' \chi_{\varepsilon,\lambda} ( \theta_{ - \varepsilon} + ) - \mu_{2} \chi_{\varepsilon,\lambda} ' ( \theta_{ - \varepsilon} + ) \bigr), \\ &y' ( \theta_{ - \varepsilon} - ) = \frac{1}{D_{1}} \bigl( - \mu_{1}'\chi_{\varepsilon,\lambda} ( \theta_{ - \varepsilon} + ) + \mu_{1}\chi_{\varepsilon,\lambda} ' ( \theta_{ - \varepsilon} + ) \bigr). \end{aligned} $$
(25)

Hence, \(\chi_{\lambda} ( x ) = \chi ( x,\lambda )\) satisfies of (1), (3) and (4)-(7) on I.

Let us consider the Wronskians

$$\begin{aligned}& W ( \phi_{ \pm \varepsilon,\lambda},\chi_{ \pm \varepsilon,\lambda};x ) = \phi_{ \pm \varepsilon,\lambda} ( x ) \chi_{ \pm \varepsilon,\lambda} ' ( x ) - \phi_{ \pm \varepsilon,\lambda} ' ( x )\chi_{ \pm \varepsilon,\lambda} ( x ), \\& W ( \phi_{\varepsilon,\lambda},\chi_{\varepsilon,\lambda};x ) = \phi_{\varepsilon,\lambda} ( x ) \chi_{\varepsilon,\lambda} ' ( x ) - \phi_{\varepsilon,\lambda} ' ( x )\chi_{\varepsilon,\lambda} ( x ), \end{aligned}$$

which are independent of \(x \in I_{i} \) (\(i = 1,2,3 \)) and entire functions where \(I_{i} \) (\(i = 1,2,3 \)). After a short calculation we see that \(D_{1}D_{2}\omega_{ - \varepsilon} ( \lambda ) = D_{2}\omega_{\varepsilon} ( \lambda ) = \omega_{ + \varepsilon} ( \lambda )\). Now we may introduce the characteristic function \(\omega ( \lambda )\) as

$$ \omega ( \lambda ): = \omega_{ - \varepsilon} ( \lambda ) = \frac{1}{D_{1}} \omega_{\varepsilon} ( \lambda ) = \frac{1}{D_{1}D_{2}}\omega_{ + \varepsilon} ( \lambda ). $$
(26)

Corollary 2

The zeros of the functions \(\omega_{ - \varepsilon} ( \lambda )\), \(\omega_{\varepsilon} ( \lambda )\) and \(\omega_{ + \varepsilon} ( \lambda )\) coincide.

Theorem 1

The eigenvalues of the problem (1)-(7) are the zeros of the function \(\omega ( \lambda )\).

Proof

Let \(\omega ( \lambda_{0} ) = 0\). Then \(W ( \phi_{ - \varepsilon,\lambda_{0}},\chi_{ - \varepsilon,\lambda_{0}};x ) = 0\) and therefore the functions \(\phi_{ - \varepsilon,\lambda_{0}} ( x )\) and \(\chi_{ - \varepsilon,\lambda_{0}} ( x )\) are linearly dependent, i.e. \(\chi_{ - \varepsilon,\lambda_{0}} ( x ) = k_{1}\phi_{ - \varepsilon,\lambda_{0}} ( x )\), \(x \in [ a,\theta_{ - \varepsilon} ]\) for some \(k_{1} \ne 0\). From this, it follows that \(\chi_{\lambda_{0}} ( x )\) satisfies also the first boundary condition (2), so \(\chi_{\lambda_{0}} ( x )\) is an eigenfunction for the eigenvalue \(\lambda_{0}\).

Now let \(y ( x,\lambda_{0} )\) be any eigenfunction correspond to eigenvalue \(\lambda_{0}\), but \(\omega ( \lambda_{0} ) \ne 0\). Then the pair of the functions \(( \phi_{ - \varepsilon} (\cdot ),\chi_{ - \varepsilon} (\cdot ) )\), \(( \phi_{\varepsilon} (\cdot ),\chi_{\varepsilon} (\cdot ) )\) and \(( \phi_{ + \varepsilon} (\cdot ),\chi_{ + \varepsilon} (\cdot ) )\) would be linearly independent on \(I_{i} \) (\(i = 1,2,3 \)), respectively. Therefore \(y ( x,\lambda_{0} )\) may be represented as

$$ y ( x,\lambda_{0} ) = \textstyle\begin{cases} c_{1}\phi_{ - \varepsilon,\lambda_{0}} ( x ) + c_{2}\chi_{ - \varepsilon,\lambda_{0}} ( x ),& x \in [ a,\theta_{ - \varepsilon} ), \\ c_{3}\phi_{\varepsilon,\lambda_{0}} ( x ) + c_{4}\chi_{\varepsilon,\lambda_{0}} ( x ),& x \in ( \theta_{ - \varepsilon},\theta_{ + \varepsilon} ), \\ c_{5}\phi_{ + \varepsilon,\lambda_{0}} ( x ) + c_{6}\chi_{ + \varepsilon,\lambda_{0}} ( x ),& x \in ( \theta_{ + \varepsilon},b ], \end{cases} $$
(27)

where at least one of the constants \(c_{i} \) (\(i = \overline{1,6} \)) is not zero. Considering the equations

$$ \begin{aligned} &B_{a} \bigl( y ( x,\lambda_{0} ) \bigr) = 0,\quad\quad B_{b} \bigl( y ( x,\lambda_{0} ) \bigr) = 0, \\ &T_{ \pm \varepsilon} \bigl( y ( x,\lambda_{0} ) \bigr) = 0,\quad\quad T_{ \pm \varepsilon} ' \bigl( y ( x,\lambda_{0} ) \bigr) = 0, \end{aligned} $$
(28)

as a system of linear equations of the variables \(c_{i}\) (\(i = \overline{1,6} \)) and taking (21), (22), (24) and (25) into account, it follows that the determinant of this system is

$$\begin{aligned}& \left| \begin{matrix} 0 & \omega_{ - \varepsilon} ( \lambda_{0} ) & 0 & 0 & 0 & 0 \\ \phi_{\varepsilon,\lambda_{0}} ( \theta_{ - \varepsilon} + ) & \chi_{\varepsilon,\lambda_{0}} ( \theta_{ - \varepsilon} + ) & - \phi_{\varepsilon,\lambda_{0}} ( \theta_{ - \varepsilon} + ) & - \chi_{\varepsilon,\lambda_{0}} ( \theta_{ - \varepsilon} + ) & 0 & 0 \\ \phi_{\varepsilon,\lambda_{0}}' ( \theta_{ - \varepsilon} + ) & \chi_{\varepsilon,\lambda_{0}}' ( \theta_{ - \varepsilon} + ) & - \phi_{\varepsilon,\lambda_{0}}' ( \theta_{ - \varepsilon} + ) & - \chi_{\varepsilon,\lambda_{0}}' ( \theta_{ - \varepsilon} + ) & 0 & 0 \\ 0 & 0 & \phi_{ + \varepsilon,\lambda_{0}} ( \theta_{ + \varepsilon} + ) & \chi_{ + \varepsilon,\lambda_{0}} ( \theta_{ + \varepsilon} + ) & - \phi_{ + \varepsilon,\lambda_{0}} ( \theta_{ + \varepsilon} + ) & - \chi_{ + \varepsilon,\lambda_{0}} ( \theta_{ + \varepsilon} + ) \\ 0 & 0 & \phi_{ + \varepsilon,\lambda_{0}}' ( \theta_{ + \varepsilon} + ) & \chi_{ + \varepsilon,\lambda_{0}}' ( \theta_{ + \varepsilon} + ) & - \phi_{ + \varepsilon,\lambda_{0}}' ( \theta_{ + \varepsilon} + ) & - \chi_{ + \varepsilon,\lambda_{0}}' ( \theta_{ + \varepsilon} + ) \\ 0 & 0 & 0 & 0 & \omega_{ + \varepsilon} ( \lambda_{0} ) & 0 \end{matrix}\right| \\& \quad= - \omega_{ - \varepsilon} ( \lambda_{0} )\omega_{\varepsilon} ( \lambda_{0} )\omega_{ + \varepsilon}^{2} ( \lambda_{0} ) \ne 0. \end{aligned}$$

Therefore, the system has only the trivial solution \(c_{i} = 0\) (\(i = \overline{1,6} \)). Thus we get a contradiction, which completes the proof. □

4 Asymptotic approximate formulae

Now we will derive asymptotic formulae of the eigenvalues and eigenfunctions in a way similar to the techniques of [9, 20] and [13, 15].

Lemma 2

Let \(\phi_{\lambda} ( x )\) be the solution of equation (1) defined in Section  3, and let \(\lambda = s^{2}\). Then the following integral equations hold for \(k = 0 \) and \(k = 1\):

$$\begin{aligned}& \phi_{ - \varepsilon,\lambda}^{ ( k )} ( x ) = \beta_{2} \bigl( \cos \bigl( s ( x - a ) \bigr) \bigr)^{ ( k )} - \frac{\beta_{1}}{s} \bigl( \sin \bigl( s ( x - a ) \bigr) \bigr)^{ ( k )} \\& \hphantom{\phi_{ - \varepsilon,\lambda}^{ ( k )} ( x ) =}{}+ \frac{1}{s} \int_{a}^{x} \sin \bigl( s ( x - t ) \bigr)^{ ( k )} q ( t )\phi_{ - \varepsilon,\lambda} ( t )\,dt, \end{aligned}$$
(29)
$$\begin{aligned}& \phi_{\varepsilon,\lambda}^{ ( k )} ( x ) = \bigl( \mu_{1} \phi_{ - \varepsilon,\lambda} ( \theta_{ - \varepsilon} - ) + \mu_{2} \phi_{ - \varepsilon,\lambda} ' ( \theta_{ - \varepsilon} - ) \bigr) \bigl( \cos \bigl( s ( x - \theta_{ - \varepsilon} ) \bigr) \bigr)^{ ( k )} \\& \hphantom{\phi_{\varepsilon,\lambda}^{ ( k )} ( x ) =}{}+ \frac{1}{s} \bigl( \mu_{1}'\phi_{ - \varepsilon,\lambda} ( \theta_{ - \varepsilon} - ) + \mu_{2}'\phi_{ - \varepsilon,\lambda} ' ( \theta_{ - \varepsilon} - ) \bigr) \bigl( \sin \bigl( s ( x - \theta_{ - \varepsilon} ) \bigr) \bigr)^{ ( k )} \\& \hphantom{\phi_{\varepsilon,\lambda}^{ ( k )} ( x ) =}{}+ \frac{1}{s} \int_{\theta_{ - \varepsilon}}^{x} \bigl( \sin \bigl( s ( x - t ) \bigr) \bigr)^{ ( k )}q ( t )\phi_{\varepsilon,\lambda} ( t )\,dt, \end{aligned}$$
(30)
$$\begin{aligned}& \phi_{ + \varepsilon,\lambda}^{ ( k )} ( x ) = \bigl( \eta_{1} \phi_{\varepsilon,\lambda} ( \theta_{ + \varepsilon} - ) + \eta_{2} \phi_{\varepsilon,\lambda} ' ( \theta_{ + \varepsilon} - ) \bigr) \bigl( \cos \bigl( s ( x - \theta_{ + \varepsilon} ) \bigr) \bigr)^{ ( k )} \\& \hphantom{\phi_{ + \varepsilon,\lambda}^{ ( k )} ( x ) =}{}+ \frac{1}{s} \bigl( \eta_{1}'\phi_{\varepsilon,\lambda} ( \theta_{ + \varepsilon} - ) + \eta_{2}'\phi_{\varepsilon,\lambda} ' ( \theta_{ + \varepsilon} - ) \bigr) \bigl( \sin \bigl( s ( x - \theta_{ + \varepsilon} ) \bigr) \bigr)^{ ( k )} \\& \hphantom{\phi_{ + \varepsilon,\lambda}^{ ( k )} ( x ) =}{}+ \frac{1}{s} \int_{\theta_{ + \varepsilon}}^{x} \bigl( \sin \bigl( s ( x - t ) \bigr) \bigr)^{ ( k )}q ( t )\phi_{ + \varepsilon,\lambda} ( t )\,dt. \end{aligned}$$
(31)

Proof

For proving it is enough to substitute \(s^{2}\phi_{ - \varepsilon,\lambda} ( t ) + \phi_{ - \varepsilon,\lambda}'' ( t )\), \(s^{2}\phi_{\varepsilon,\lambda} ( t ) + \phi_{\varepsilon,\lambda}'' ( t )\) and \(s^{2}\phi_{ + \varepsilon,\lambda} ( t ) + \phi_{ + \varepsilon,\lambda}'' ( t )\) instead of \(q ( t )\phi_{ - \varepsilon,\lambda} ( t )\), \(q ( t )\phi_{\varepsilon,\lambda} ( t )\) and \(q ( t )\phi_{ + \varepsilon,\lambda} ( t )\) in the integral terms of (29), (30) and (31), respectively, and integrate by parts twice. □

Lemma 3

Let \(\lambda = s^{2}\). \(\operatorname{Im} s = \ell\). Then the function \(\phi_{\lambda} ( x )\) has the following asymptotic representations for \(\vert \lambda \vert \to \infty\), which hold uniformly for \(x \in I_{ i } \) (\(i = 1,2,3 \)):

$$\begin{aligned}& \phi_{ - \varepsilon,\lambda}^{ ( k )} ( x ) = \beta_{2} \bigl( \cos \bigl( s ( x - a ) \bigr) \bigr)^{ ( k )} + O \bigl( \vert s \vert ^{k - 1}e^{\vert \ell \vert ( x - a )} \bigr), \end{aligned}$$
(32)
$$\begin{aligned}& \phi_{\varepsilon,\lambda}^{ ( k )} ( x ) = - s\mu_{2}\beta_{2}\sin \bigl( s ( \theta_{ - \varepsilon} - a ) \bigr) \bigl( \cos \bigl( s ( x - \theta_{ - \varepsilon} ) \bigr) \bigr)^{ ( k )} + O \bigl( \vert s \vert ^{k}e^{\vert \ell \vert ( x - a )} \bigr), \end{aligned}$$
(33)
$$\begin{aligned}& \phi_{ + \varepsilon,\lambda}^{ ( k )} ( x ) = s^{2}\mu_{2} \eta_{2}\beta_{2}\sin \bigl( s ( \theta_{ - \varepsilon} - a ) \bigr)\sin \bigl( s ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} ) \bigr) \bigl( \cos \bigl( s ( x - \theta_{ + \varepsilon} ) \bigr) \bigr)^{ ( k )} \\& \hphantom{\phi_{ + \varepsilon,\lambda}^{ ( k )} ( x ) =}{}+ O \bigl( \vert s \vert ^{k + 1}e^{\vert \ell \vert ( x - a )} \bigr), \end{aligned}$$
(34)

for \(\beta_{2} \ne 0\),

$$\begin{aligned}& \phi_{ - \varepsilon,\lambda}^{ ( k )} ( x ) = - \frac{\beta_{1}}{s} \bigl( \sin \bigl( s ( x - a ) \bigr) \bigr)^{ ( k )} + O \bigl( \vert s \vert ^{k - 2}e^{\vert \ell \vert ( x - a )} \bigr), \end{aligned}$$
(35)
$$\begin{aligned}& \phi_{\varepsilon,\lambda}^{ ( k )} ( x ) = - \mu_{2} \beta_{1}\cos \bigl( s ( \theta_{ - \varepsilon} - a ) \bigr) \bigl( \cos \bigl( s ( x - \theta_{ - \varepsilon} ) \bigr) \bigr)^{ ( k )} + O \bigl( \vert s \vert ^{k - 1}e^{\vert \ell \vert ( x - a )} \bigr), \end{aligned}$$
(36)
$$\begin{aligned}& \phi_{ + \varepsilon,\lambda}^{ ( k )} ( x ) = s\mu_{2}\eta_{2} \beta_{1}\cos \bigl( s ( \theta_{ - \varepsilon} - a ) \bigr)\sin \bigl( s ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} ) \bigr) \bigl( \cos \bigl( s ( x - \theta_{ + \varepsilon} ) \bigr) \bigr)^{ ( k )} \\& \hphantom{\phi_{ + \varepsilon,\lambda}^{ ( k )} ( x ) =}{}+ O \bigl( \vert s \vert ^{k}e^{\vert \ell \vert ( x - a )} \bigr), \end{aligned}$$
(37)

for \(\beta_{2} = 0\).

Proof

These formulae can be proven similar to Titchmarsh’s proof [20] and also the techniques in [13, 15]. □

Lemma 4

Let \(\lambda = s^{2}\). \(\operatorname{Im} s = \ell\). Then the characteristic function \(\omega ( \lambda )\) has the following asymptotic representations:

Case 1. If \(\beta_{2} \ne 0\), \(\alpha_{2}' \ne 0\), then

$$\begin{aligned} \omega ( \lambda ) =& \frac{1}{D_{1}D_{2}}s^{5}\beta_{2} \alpha_{2}'\mu_{2}\eta_{2}\sin \bigl( s ( \theta_{ - \varepsilon} - a ) \bigr)\sin \bigl( s ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} ) \bigr)\sin \bigl( s ( b - \theta_{ + \varepsilon} ) \bigr) \\ &{}+ O \bigl( \vert s \vert ^{4}e^{\vert \ell \vert ( b - a )} \bigr). \end{aligned}$$
(38)

Case 2. If \(\beta_{2} \ne 0\), \(\alpha_{2}' = 0\), then

$$\begin{aligned} \omega ( \lambda ) =& \frac{1}{D_{1}D_{2}}s^{4}\beta_{2} \alpha_{1}'\mu_{2}\eta_{2}\sin \bigl( s ( \theta_{ - \varepsilon} - a ) \bigr)\sin \bigl( s ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} ) \bigr)\cos \bigl( s ( b - \theta_{ + \varepsilon} ) \bigr) \\ &{}+ O \bigl( \vert s \vert ^{3}e^{\vert \ell \vert ( b - a )} \bigr). \end{aligned}$$
(39)

Case 3. If \(\beta_{2} = 0\), \(\alpha_{2}' \ne 0\), then

$$\begin{aligned} \omega ( \lambda ) =& \frac{1}{D_{1}D_{2}}s^{4}\beta_{1} \alpha_{2}'\mu_{2}\eta_{2}\cos \bigl( s ( \theta_{ - \varepsilon} - a ) \bigr)\sin \bigl( s ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} ) \bigr)\sin \bigl( s ( b - \theta_{ + \varepsilon} ) \bigr) \\ &{}+ O \bigl( \vert s \vert ^{3}e^{\vert \ell \vert ( b - a )} \bigr). \end{aligned}$$
(40)

Case 4. If \(\beta_{2} = 0\), \(\alpha_{2}' = 0\), then

$$\begin{aligned} \omega ( \lambda ) = &\frac{1}{D_{1}D_{2}}s^{3}\beta_{1} \alpha_{1}'\mu_{2}\eta_{2}\cos \bigl( s ( \theta_{ - \varepsilon} - a ) \bigr)\sin \bigl( s ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} ) \bigr)\cos \bigl( s ( b - \theta_{ + \varepsilon} ) \bigr) \\ &{}+ O \bigl( \vert s \vert ^{2}e^{\vert \ell \vert ( b - a )} \bigr). \end{aligned}$$
(41)

Proof

The proof is immediate by substituting (34) and (37) into the representation

$$\omega ( \lambda ) = \frac{1}{D_{1}D_{2}} \bigl\{ \bigl( \lambda \alpha_{1}' + \alpha_{1} \bigr) \phi_{ + \varepsilon,\lambda} ( b ) - \bigl( \lambda \alpha_{2}' + \alpha_{2} \bigr)\phi_{ + \varepsilon,\lambda} ' ( b ) \bigr\} . $$

 □

Corollary 3

The eigenvalues of the problem (1)-(7) is bounded from below.

We are now ready to find the asymptotic approximation formulae for the eigenvalues of the problem (1)-(7). Since the eigenvalues coincide with the zeros of the entire functions \(\omega ( \lambda )\), it follows that they have no finite accumulation point. Moreover, all eigenvalues are real and bounded below by Corollaries 1 and 3. Therefore, we may renumber them as \(\lambda_{0} \le \lambda_{1} \le \lambda_{2} \le\cdots\), which are counted according to their multiplicity. Below we shall denote \(s_{n}^{2} = \lambda_{n}\).

Theorem 2

The problem (1)-(7) has a precisely denumerable number of real eigenvalues, whose behaviour may be expressed by the three sequences \(\{ \lambda_{n}' \}\), \(\{ \lambda_{n}'' \} \) and \(\{ \lambda_{n} ''' \}\) with the following asymptotics representations for \(n \to \infty\):

Case 1. If \(\beta_{2} \ne 0\), \(\alpha_{2}' \ne 0\), then

$$ \begin{aligned} &s_{n}' = \frac{ ( n - 1 )\pi}{ ( \theta_{ - \varepsilon} - a )} + O \biggl( \frac{1}{n} \biggr),\quad\quad s_{n}^{\prime\prime} = \frac{ ( n - 1 )\pi}{ ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} )} + O \biggl( \frac{1}{n} \biggr), \\ &s_{n} ''' = \frac{ ( n - 2 )\pi}{ ( b - \theta_{ + \varepsilon} )} + O \biggl( \frac{1}{n} \biggr). \end{aligned} $$
(42)

Case 2. If \(\beta_{2} \ne 0\), \(\alpha_{2}' = 0\), then

$$ \begin{aligned} &s_{n}' = \frac{ ( n - 1 )\pi}{ ( \theta_{ - \varepsilon} - a )} + O \biggl( \frac{1}{n} \biggr),\quad\quad s_{n}^{\prime\prime} = \frac{ ( n - 1 )\pi}{ ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} )} + O \biggl( \frac{1}{n} \biggr),\\ & s_{n} ''' = \frac{ ( n - 1/2 )\pi}{ ( b - \theta_{ + \varepsilon} )} + O \biggl( \frac{1}{n} \biggr). \end{aligned} $$
(43)

Case 3. If \(\beta_{2} = 0\), \(\alpha_{2}' \ne 0\), then

$$ \begin{aligned} &s_{n}' = \frac{ ( n - 1/2 )\pi}{ ( \theta_{ - \varepsilon} - a )} + O \biggl( \frac{1}{n} \biggr), \quad\quad s_{n}^{\prime\prime} = \frac{ ( n - 1 )\pi}{ ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} )} + O \biggl( \frac{1}{n} \biggr),\\ & s_{n} ''' = \frac{ ( n - 1 )\pi}{ ( b - \theta_{ + \varepsilon} )} + O \biggl( \frac{1}{n} \biggr). \end{aligned} $$
(44)

Case 4. If \(\beta_{2} = 0\), \(\alpha_{2}' = 0\), then

$$ \begin{aligned} &s_{n}' = \frac{ ( n - 1/2 )\pi}{ ( \theta_{ - \varepsilon} - a )} + O \biggl( \frac{1}{n} \biggr), \quad\quad s_{n}^{\prime\prime} = \frac{ ( n - 1 )\pi}{ ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} )} + O \biggl( \frac{1}{n} \biggr), \\ &s_{n} ''' = \frac{ ( n - 1/2 )\pi}{ ( b - \theta_{ + \varepsilon} )} + O \biggl( \frac{1}{n} \biggr). \end{aligned} $$
(45)

Proof

We will only consider the first case. By applying the well-known Rouche theorem on a sufficiently large contour, it follows that \(\omega ( \lambda )\) has the same number of zeros inside the contour as the leading term in (38). Hence, if \(\lambda_{0}' \le \lambda_{1}' \le \lambda_{2}' \le\cdots\) are the zeros of \(\omega ( \lambda )\) and \(s_{n}^{2'} = \lambda_{n}'\) we have

$$ s_{n}' = \frac{ ( n - 1 )\pi}{ ( \theta_{ - \varepsilon} - a )} + \delta_{n}', $$
(46)

for sufficiently large n, where \(\vert \delta_{n}' \vert \le \frac{\pi}{2 ( \theta_{ - \varepsilon} - a )}\). By using (38) we have \(\delta_{n}' = O ( \frac{1}{n} )\), which completes the proof for the first formula of Case 1. The proof for the other cases are similar. □

Then from (32)-(37) (for \(k = 0\)) and the above theorem, the asymptotic behaviour of the eigenfunctions of the problem (1)-(7) is given by:

Case 1. If \(\beta_{2} \ne 0\), \(\alpha_{2}' \ne 0\), then

$$\begin{aligned}& \phi_{\lambda_{n}'} ( x ) = \textstyle\begin{cases} \beta_{2}\cos ( \frac{ ( n - 1 )\pi}{ ( \theta_{ - \varepsilon} - a )} ( x - a ) ) + O ( \frac{1}{n} ),& x \in [ a,\theta_{ - \varepsilon} ), \\ O ( \frac{1}{n} ),& x \in ( \theta_{ - \varepsilon},\theta_{ + \varepsilon} ), \\ O ( \frac{1}{n} ),& x \in ( \theta_{ + \varepsilon},b ], \end{cases}\displaystyle \\& \phi_{\lambda_{n}''} ( x ) = \textstyle\begin{cases} \beta_{2}\cos ( \frac{ ( n - 1 )\pi}{ ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} )} ( x - a ) ) + O ( \frac{1}{n} ),& x \in [ a,\theta_{ - \varepsilon} ), \\ - \frac{ ( n - 1 )\pi}{ ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} )}\mu_{2}\beta_{2}\sin ( \frac{ ( n - 1 )\pi}{ ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} )} ( \theta_{ - \varepsilon} - a ) )\cos ( \frac{ ( n - 1 )\pi}{ ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} )} ( x - \theta_{ - \varepsilon} ) ) \\ \quad{}+ O ( 1 ),& x \in ( \theta_{ - \varepsilon},\theta_{ + \varepsilon} ), \\ O ( \frac{1}{n} ),& x \in ( \theta_{ + \varepsilon},b ], \end{cases}\displaystyle \\& \phi_{\lambda_{n} ''' } ( x ) = \textstyle\begin{cases} \beta_{2}\cos ( \frac{ ( n - 2 )\pi}{ ( b - \theta_{ + \varepsilon} )} ( x - a ) ) + O ( \frac{1}{n} ),& x \in [ a,\theta_{ - \varepsilon} ), \\ - \frac{ ( n - 2 )\pi}{ ( b - \theta_{ + \varepsilon} )}\mu_{2}\beta_{2}\sin ( \frac{ ( n - 2 )\pi}{ ( b - \theta_{ + \varepsilon} )} ( \theta_{ - \varepsilon} - a ) )\cos ( \frac{ ( n - 2 )\pi}{ ( b - \theta_{ + \varepsilon} )} ( x - \theta_{ - \varepsilon} ) ) \\ \quad{}+ O ( 1 ),& x \in ( \theta_{ - \varepsilon},\theta_{ + \varepsilon} ), \\ ( \frac{ ( n - 2 )\pi}{ ( b - \theta_{ + \varepsilon} )} )^{2}\mu_{2}\eta_{2}\beta_{2}\sin ( \frac{ ( n - 2 )\pi}{ ( b - \theta_{ + \varepsilon} )} ( \theta_{ - \varepsilon} - a ) )\sin ( \frac{ ( n - 2 )\pi}{ ( b - \theta_{ + \varepsilon} )} ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} ) ) \\ \quad{}\times \cos ( \frac{ ( n - 2 )\pi}{ ( b - \theta_{ + \varepsilon} )} ( x - \theta_{ + \varepsilon} ) ) + O ( n ),& x \in ( \theta_{ + \varepsilon},b ]. \end{cases}\displaystyle \end{aligned}$$

Case 2. If \(\beta_{2} \ne 0\), \(\alpha_{2}' = 0\), then

$$\begin{aligned}& \phi_{\lambda_{n}'} ( x ) = \textstyle\begin{cases} \beta_{2}\cos ( \frac{ ( n - 1 )\pi}{ ( \theta_{ - \varepsilon} - a )} ( x - a ) ) + O ( \frac{1}{n} ),& x \in [ a,\theta_{ - \varepsilon} ), \\ O ( \frac{1}{n} ),& x \in ( \theta_{ - \varepsilon},\theta_{ + \varepsilon} ), \\ O ( \frac{1}{n} ),& x \in ( \theta_{ + \varepsilon},b ], \end{cases}\displaystyle \\& \phi_{\lambda_{n}''} ( x ) = \textstyle\begin{cases} \beta_{2}\cos ( \frac{ ( n - 1 )\pi}{ ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} )} ( x - a ) ) + O ( \frac{1}{n} ),& x \in [ a,\theta_{ - \varepsilon} ), \\ - \frac{ ( n - 1 )\pi}{ ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} )}\mu_{2}\beta_{2}\sin ( \frac{ ( n - 1 )\pi}{ ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} )} ( \theta_{ - \varepsilon} - a ) )\cos ( \frac{ ( n - 1 )\pi}{ ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} )} ( x - \theta_{ - \varepsilon} ) ) \\ \quad{}+ O ( 1 ),& x \in ( \theta_{ - \varepsilon},\theta_{ + \varepsilon} ), \\ O ( \frac{1}{n} ),& x \in ( \theta_{ + \varepsilon},b ], \end{cases}\displaystyle \\& \phi_{\lambda_{n} ''' } ( x ) = \textstyle\begin{cases} \beta_{2}\cos ( \frac{ ( n - 1/2 )\pi}{ ( b - \theta_{ + \varepsilon} )} ( x - a ) ) + O ( \frac{1}{n} ),&\hspace{-3pt} x \in [ a,\theta_{ - \varepsilon} ), \\ - \frac{ ( n - 1/2 )\pi}{ ( b - \theta_{ + \varepsilon} )} \mu_{2}\beta_{2}\sin ( \frac{ ( n - 1/2 )\pi}{ ( b - \theta_{ + \varepsilon} )} ( \theta_{ - \varepsilon} - a ) )\cos ( \frac{ ( n - 1/2 )\pi}{ ( b - \theta_{ + \varepsilon} )} ( x - \theta_{ - \varepsilon} ) ) \\ \quad{}+ O ( 1 ),& \hspace{-3pt}x \in ( \theta_{ - \varepsilon},\theta_{ + \varepsilon} ), \\ ( \frac{ ( n - 1/2 )\pi}{ ( b - \theta_{ + \varepsilon} )} )^{2}\mu_{2}\eta_{2}\beta_{2}\sin ( \frac{ ( n - 1/2 )\pi}{ ( b - \theta_{ + \varepsilon} )} ( \theta_{ - \varepsilon} - a ) )\sin ( \frac{ ( n - 1/2 )\pi}{ ( b - \theta_{ + \varepsilon} )} ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} ) )\\ \quad{}\times \cos ( \frac{ ( n - 1/2 )\pi}{ ( b - \theta_{ + \varepsilon} )} ( x - \theta_{ + \varepsilon} ) ) + O ( n ), &\hspace{-3pt} x \in ( \theta_{ + \varepsilon},b ]. \end{cases}\displaystyle \end{aligned}$$

Case 3. If \(\beta_{2} = 0\), \(\alpha_{2}' \ne 0\), then

$$\begin{aligned}& \phi_{\lambda_{n}'} ( x ) = \textstyle\begin{cases} - \frac{\beta_{1} ( \theta_{ - \varepsilon} - a )}{ ( n - 1/2 )\pi} \sin ( \frac{ ( n - 1/2 )\pi}{ ( \theta_{ - \varepsilon} - a )} ( x - a ) ) + O ( \frac{1}{n^{2}} ),& x \in [ a,\theta_{ - \varepsilon} ), \\ O ( \frac{1}{n} ), &x \in ( \theta_{ - \varepsilon},\theta_{ + \varepsilon} ), \\ O ( 1 ),& x \in ( b,\theta_{ + \varepsilon} ], \end{cases}\displaystyle \\& \phi_{\lambda_{n}''} ( x ) = \textstyle\begin{cases} - \frac{\beta_{1} ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} )}{ ( n - 1 )\pi} \sin ( \frac{ ( n - 1 )\pi}{ ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} )} ( x - a ) ) + O ( \frac{1}{n^{2}} ),& x \in [ a,\theta_{ - \varepsilon} ), \\ - \mu_{2}\beta_{1}\cos ( \frac{ ( n - 1 )\pi}{ ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} )} ( \theta_{ - \varepsilon} - a ) )\cos ( \frac{ ( n - 1 )\pi}{ ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} )} ( x - \theta_{ - \varepsilon} ) ) \\ \quad{}+ O ( \frac{1}{n} ), &x \in ( \theta_{ - \varepsilon},\theta_{ + \varepsilon} ), \\ O ( 1 ),& x \in ( b,\theta_{ + \varepsilon} ], \end{cases}\displaystyle \\& \phi_{\lambda_{n} ''' } ( x ) = \textstyle\begin{cases} - \frac{\beta_{1} ( b - \theta_{ + \varepsilon} )}{ ( n - 1 )\pi} \sin ( \frac{ ( n - 1 )\pi}{ ( b - \theta_{ + \varepsilon} )} ( x - a ) ) + O ( \frac{1}{n^{2}} ),& x \in [ a,\theta_{ - \varepsilon} ), \\ - \mu_{2}\beta_{1}\cos ( \frac{ ( n - 1 )\pi}{ ( b - \theta_{ + \varepsilon} )} ( \theta_{ - \varepsilon} - a ) )\cos ( \frac{ ( n - 1 )\pi}{ ( b - \theta_{ + \varepsilon} )} ( x - \theta_{ - \varepsilon} ) ) \\ \quad{}+ O ( \frac{1}{n} ),& x \in ( \theta_{ - \varepsilon},\theta_{ + \varepsilon} ), \\ \frac{ ( n - 1 )\pi}{ ( b - \theta_{ + \varepsilon} )}\mu_{2}\beta_{1}\eta_{2}\cos ( \frac{ ( n - 1 )\pi}{ ( b - \theta_{ + \varepsilon} )} ( \theta_{ - \varepsilon} - a ) )\sin ( \frac{ ( n - 1 )\pi}{ ( b - \theta_{ + \varepsilon} )} ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} ) ) \\ \quad{}\times \cos ( \frac{ ( n - 1 )\pi}{ ( b - \theta_{ + \varepsilon} )} ( x - \theta_{ - \varepsilon} ) ) + O ( 1 ),& x \in ( b,\theta_{ + \varepsilon} ]. \end{cases}\displaystyle \end{aligned}$$

Case 4. If \(\beta_{2} = 0\), \(\alpha_{2}' = 0\), then

$$\begin{aligned}& \phi_{\lambda_{n}'} ( x ) = \textstyle\begin{cases} - \frac{\beta_{1} ( \theta_{ - \varepsilon} - a )}{ ( n - 1/2 )\pi} \sin ( \frac{ ( n - 1/2 )\pi}{ ( \theta_{ - \varepsilon} - a )} ( x - a ) ) + O ( \frac{1}{n^{2}} ), &x \in [ a,\theta_{ - \varepsilon} ), \\ O ( \frac{1}{n} ),& x \in ( \theta_{ - \varepsilon},\theta_{ + \varepsilon} ), \\ O ( 1 ),& x \in ( b,\theta_{ + \varepsilon} ], \end{cases}\displaystyle \\& \phi_{\lambda_{n}''} ( x ) = \textstyle\begin{cases} - \frac{\beta_{1} ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} )}{ ( n - 1 )\pi} \sin ( \frac{ ( n - 1 )\pi}{ ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} )} ( x - a ) ) + O ( \frac{1}{n^{2}} ),& x \in [ a,\theta_{ - \varepsilon} ), \\ - \mu_{2}\beta_{1}\cos ( \frac{ ( n - 1 )\pi}{ ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} )} ( \theta_{ - \varepsilon} - a ) )\cos ( \frac{ ( n - 1 )\pi}{ ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} )} ( x - \theta_{ - \varepsilon} ) ) \\ \quad{}+ O ( \frac{1}{n} ),& x \in ( \theta_{ - \varepsilon},\theta_{ + \varepsilon} ), \\ O ( 1 ),& x \in ( b,\theta_{ + \varepsilon} ], \end{cases}\displaystyle \\& \phi_{\lambda_{n} ''' } ( x ) = \textstyle\begin{cases} - \frac{\beta_{1} ( b - \theta_{ + \varepsilon} )}{ ( n - 1/2 )\pi} \sin ( \frac{ ( n - 1/2 )\pi}{ ( b - \theta_{ + \varepsilon} )} ( x - a ) ) + O ( \frac{1}{n^{2}} ),& x \in [ a,\theta_{ - \varepsilon} ), \\ - \mu_{2}\beta_{1}\cos ( \frac{ ( n - 1/2 )\pi}{ ( b - \theta_{ + \varepsilon} )} ( \theta_{ - \varepsilon} - a ) )\cos ( \frac{ ( n - 1/2 )\pi}{ ( b - \theta_{ + \varepsilon} )} ( x - \theta_{ - \varepsilon} ) ) \\ \quad{}+ O ( \frac{1}{n} ),& x \in ( \theta_{ - \varepsilon},\theta_{ + \varepsilon} ), \\ \frac{ ( n - 1/2 )\pi}{ ( b - \theta_{ + \varepsilon} )}\mu_{2}\beta_{1}\eta_{2}\cos ( \frac{ ( n - 1/2 )\pi}{ ( b - \theta_{ + \varepsilon} )} ( \theta_{ - \varepsilon} - a ) )\sin ( \frac{ ( n - 1/2 )\pi}{ ( b - \theta_{ + \varepsilon} )} ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} ) ) \\ \quad{}\times \cos ( \frac{ ( n - 1/2 )\pi}{ ( b - \theta_{ + \varepsilon} )} ( x - \theta_{ - \varepsilon} ) ) + O ( 1 ),& x \in ( b,\theta_{ + \varepsilon} ]. \end{cases}\displaystyle \end{aligned}$$

All these asymptotic formulae hold uniformly for \(x \in I\).

5 Green’s function

Let \(F = \binom{ f ( x ) }{ h }\) be a continuous function. To study the completeness of the eigenelements of A, and hence the completeness of the eigenfunctions (1)-(7), we derive Green’s function of the problem (1)-(7) as well as the resolvent of A. Indeed let \(\lambda \in \mathbb{C}\) not be an eigenvalue of A and consider the inhomogeneous problem

$$ ( \lambda \mathrm{I} - A )Y = F,\quad x \in I, $$
(47)

where I is the identity operator and \(Y = \binom{ y ( x ) }{ R' ( y ) } \in D ( A )\). Since

$$ ( \lambda \mathrm{I} - A )Y = \lambda \begin{pmatrix} y ( x ) \\ R' ( y ) \end{pmatrix} - \begin{pmatrix} \tau ( y ) \\ - R ( y ) \end{pmatrix} = \begin{pmatrix} f ( x ) \\ h \end{pmatrix}, $$
(48)

we have

$$\begin{aligned}& ( \lambda - \tau )y ( x ) = f ( x ),\quad x \in I, \end{aligned}$$
(49)
$$\begin{aligned}& \lambda R' ( y ) + R ( y ) = h. \end{aligned}$$
(50)

Now we can represent the general solution of homogeneous differential equation (1), appropriate to (49) in the following form:

$$y ( x,\lambda ) =\textstyle\begin{cases} c_{1}\phi_{ - \varepsilon,\lambda} ( x ) + c_{2}\chi_{ - \varepsilon,\lambda} ( x ),& x \in [ a,\theta_{ - \varepsilon} ), \\ c_{3}\phi_{\varepsilon,\lambda} ( x ) + c_{4}\chi_{\varepsilon,\lambda} ( x ),& x \in ( \theta_{ - \varepsilon},\theta_{ + \varepsilon} ), \\ c_{5}\phi_{ + \varepsilon,\lambda} ( x ) + c_{6}\chi_{ + \varepsilon,\lambda} ( x ),& x \in ( \theta_{ + \varepsilon},b ], \end{cases} $$

in which \(c_{i}\) (\(i = \overline{1,6}\)) are arbitrary constants. By applying the method of variation of the constants, we shall search the general solution of the non-homogeneous linear differential equation (49) in the following form:

$$ y ( x,\lambda ) = \textstyle\begin{cases} c_{1} ( x,\lambda )\phi_{ - \varepsilon,\lambda} ( x ) + c_{2} ( x,\lambda )\chi_{ - \varepsilon,\lambda} ( x ), &x \in [ a,\theta_{ - \varepsilon} ), \\ c_{3} ( x,\lambda )\phi_{\varepsilon,\lambda} ( x ) + c_{4} ( x,\lambda )\chi_{\varepsilon,\lambda} ( x ), &x \in ( \theta_{ - \varepsilon},\theta_{ + \varepsilon} ), \\ c_{5} ( x,\lambda )\phi_{ + \varepsilon,\lambda} ( x ) + c_{6} ( x,\lambda )\chi_{ + \varepsilon,\lambda} ( x ),& x \in ( \theta_{ + \varepsilon},b ], \end{cases} $$
(51)

where the functions \(c_{i} ( x,\lambda )\) (\(i = \overline{1,6} \)) satisfy the linear system of equations

$$\begin{aligned}& \textstyle\begin{cases} c_{1}' ( x,\lambda )\phi_{ - \varepsilon,\lambda} ( x ) + c_{2}' ( x,\lambda )\chi_{ - \varepsilon,\lambda} ( x ) = 0, \\ c_{1}' ( x,\lambda )\phi '_{ - \varepsilon,\lambda} ( x ) + c_{2}' ( x,\lambda )\chi '_{ - \varepsilon,\lambda} ( x ) = f ( x ), \end{cases}\displaystyle \mbox{for } x \in \bigl[ a,\theta_{ - \varepsilon} ), \end{aligned}$$
(52)
$$\begin{aligned}& \textstyle\begin{cases} c_{3}' ( x,\lambda )\phi_{\varepsilon,\lambda} ( x ) + c_{4}' ( x,\lambda )\chi_{\varepsilon,\lambda} ( x ) = 0, \\ c_{3}' ( x,\lambda )\phi '_{\varepsilon,\lambda} ( x ) + c_{4}' ( x,\lambda )\chi '_{\varepsilon,\lambda} ( x ) = f ( x ), \end{cases}\displaystyle \mbox{for } x \in ( \theta_{ - \varepsilon},\theta_{ + \varepsilon} ), \end{aligned}$$
(53)
$$\begin{aligned}& \textstyle\begin{cases} c_{5}' ( x,\lambda )\phi_{ + \varepsilon,\lambda} ( x ) + c_{6}' ( x,\lambda )\chi_{ + \varepsilon,\lambda} ( x ) = 0, \\ c_{5}' ( x,\lambda )\phi '_{ + \varepsilon,\lambda} ( x ) + c_{6}' ( x,\lambda )\chi '_{ + \varepsilon,\lambda} ( x ) = f ( x ), \end{cases}\displaystyle \mbox{for }x \in ( \theta_{ - \varepsilon},b \bigr]. \end{aligned}$$
(54)

Since λ is not an eigenvalue and \(\omega_{ - \varepsilon} ( \lambda ) \ne 0\), \(\omega_{\varepsilon} ( \lambda ) \ne 0\) and \(\omega_{ + \varepsilon} ( \lambda ) \ne 0\), each of the linear systems in (52)-(54) has a unique solution, which leads to

$$ \begin{aligned} &c_{1} ( x,\lambda ) = \frac{1}{\omega_{ - \varepsilon} ( \lambda )} \int_{x}^{\theta_{ - \varepsilon}} \chi_{ - \varepsilon,\lambda} ( t )f ( t )\,dt + c_{1} ( \lambda ), \\ &c_{2} ( x,\lambda ) = \frac{1}{\omega_{ - \varepsilon} ( \lambda )} \int_{a}^{x} \phi_{ - \varepsilon,\lambda} ( t )f ( t )\,dt + c_{2} ( \lambda ), \\ &c_{3} ( x,\lambda ) = \frac{1}{\omega_{\varepsilon} ( \lambda )} \int_{x}^{\theta_{ + \varepsilon}} \chi_{\varepsilon,\lambda} ( t )f ( t )\,dt + c_{3} ( \lambda ), \\ &c_{4} ( x,\lambda ) = \frac{1}{\omega_{\varepsilon} ( \lambda )} \int_{\theta_{ - \varepsilon}}^{x} \phi_{\varepsilon,\lambda} ( t )f ( t )\,dt + c_{4} ( \lambda ), \\ &c_{5} ( x,\lambda ) = \frac{1}{\omega_{ + \varepsilon} ( \lambda )} \int_{x}^{b} \chi_{ + \varepsilon,\lambda} ( t )f ( t )\,dt + c_{5} ( \lambda ), \\ &c_{6} ( x,\lambda ) = \frac{1}{\omega_{ + \varepsilon} ( \lambda )} \int_{\theta_{ + \varepsilon}}^{x} \phi_{ + \varepsilon,\lambda} ( t )f ( t )\,dt + c_{6} ( \lambda ), \end{aligned} $$
(55)

where \(c_{i} ( \lambda ) \) (\(i = \overline{1,6} \)) are arbitrary constants. Substituting (55) into (51), we obtain the solution of (49),

$$ y ( x,\lambda ) = \textstyle\begin{cases} \frac{\phi_{ - \varepsilon,\lambda} ( x )}{\omega_{ - \varepsilon} ( \lambda )}\int_{x}^{\theta_{ - \varepsilon}} \chi_{ - \varepsilon,\lambda} ( t )f ( t )\,dt + \frac{\chi_{ - \varepsilon,\lambda} ( x )}{\omega_{ - \varepsilon} ( \lambda )}\int_{a}^{x} \phi_{ - \varepsilon,\lambda} ( t )f ( t )\,dt \\ \quad{}+ c_{1} ( \lambda )\phi_{ - \varepsilon,\lambda} ( x ) + c_{2} ( \lambda )\chi_{ - \varepsilon,\lambda} ( x ),& x \in [ a,\theta_{ - \varepsilon} ), \\ \frac{\phi_{\varepsilon,\lambda} ( x )}{\omega_{\varepsilon} ( \lambda )}\int_{x}^{\theta_{ + \varepsilon}} \chi_{\varepsilon,\lambda} ( t )f ( t )\,dt + \frac{\chi_{\varepsilon,\lambda} ( x )}{\omega_{\varepsilon} ( \lambda )}\int_{\theta_{ - \varepsilon}}^{x} \phi_{\varepsilon,\lambda} ( t )f ( t )\,dt \\ \quad{}+ c_{3} ( \lambda )\phi_{\varepsilon,\lambda} ( x ) + c_{4} ( \lambda )\chi_{\varepsilon,\lambda} ( x ),& x \in ( \theta_{ - \varepsilon},\theta_{ + \varepsilon} ), \\ \frac{\phi_{ + \varepsilon,\lambda} ( x )}{\omega_{ + \varepsilon} ( \lambda )}\int_{x}^{b} \chi_{ + \varepsilon,\lambda} ( t )f ( t )\,dt + \frac{\chi_{ + \varepsilon,\lambda} ( x )}{\omega_{ + \varepsilon} ( \lambda )}\int_{\theta_{ + \varepsilon}}^{x} \phi_{ + \varepsilon,\lambda} ( t )f ( t )\,dt \\ \quad{}+ c_{5} ( \lambda )\phi_{ + \varepsilon,\lambda} ( x ) + c_{6} ( \lambda )\chi_{ + \varepsilon,\lambda} ( x ), &x \in ( \theta_{ + \varepsilon},b ]. \end{cases} $$
(56)

Then from (2), (50) and (4)-(7), we get

$$ \begin{aligned} &c_{1} ( \lambda ) = \frac{1}{\omega_{\varepsilon} ( \lambda )} \int_{\theta_{ - \varepsilon}}^{\theta_{ + \varepsilon}} \chi_{\varepsilon,\lambda} ( t )f ( t )\,dt + \frac{1}{\omega_{ + \varepsilon} ( \lambda )} \int_{\theta_{ + \varepsilon}}^{b} \chi_{ + \varepsilon,\lambda} ( t )f ( t )\,dt + \frac{h}{\omega_{ + \varepsilon} ( \lambda )}, \\ &c_{2} ( \lambda ) = 0,\\ &c_{3} ( \lambda ) = \frac{1}{\omega_{ + \varepsilon} ( \lambda )} \int_{\theta_{ + \varepsilon}}^{b} \chi_{ + \varepsilon,\lambda} ( t )f ( t )\,dt + \frac{h}{\omega_{ + \varepsilon} ( \lambda )},\\ &c_{4} ( \lambda ) = \frac{1}{\omega_{ - \varepsilon} ( \lambda )} \int_{a}^{\theta_{ - \varepsilon}} \phi_{ - \varepsilon,\lambda} ( t )f ( t )\,dt, \\ &c_{5} ( \lambda ) = \frac{h}{\omega_{ + \varepsilon} ( \lambda )},\\ &c_{6} ( \lambda ) = \frac{1}{\omega_{ - \varepsilon} ( \lambda )} \int_{a}^{\theta_{ - \varepsilon}} \phi_{ - \varepsilon,\lambda} ( t )f ( t )\,dt + \frac{1}{\omega_{\varepsilon} ( \lambda )} \int_{\theta_{ - \varepsilon}}^{\theta_{ + \varepsilon}} \phi_{\varepsilon,\lambda} ( t )f ( t )\,dt. \end{aligned} $$
(57)

Substituting (57) and (26) into (56), then (56) can be written as

$$ y ( x,\lambda ) = \textstyle\begin{cases}\frac{\phi_{ - \varepsilon,\lambda} ( x )}{\omega ( \lambda )}\int_{x}^{\theta_{ - \varepsilon}} \chi_{ - \varepsilon,\lambda} ( t )f ( t )\,dt + \frac{\chi_{ - \varepsilon,\lambda} ( x )}{\omega ( \lambda )}\int_{a}^{x} \phi_{ - \varepsilon,\lambda} ( t )f ( t )\,dt \\ \quad{}+ \frac{\phi_{ - \varepsilon,\lambda} ( x )}{D_{1}\omega ( \lambda )}\int_{\theta_{ - \varepsilon}}^{\theta_{ + \varepsilon}} \chi_{\varepsilon,\lambda} ( t )f ( t )\,dt \\ \quad{}+ \frac{\phi_{ - \varepsilon,\lambda} ( x )}{D_{1}D_{2}\omega ( \lambda )}\int_{\theta_{ + \varepsilon}}^{b} \chi_{ + \varepsilon,\lambda} ( t )f ( t )\,dt + \frac{h}{D_{1}D_{2}\omega ( \lambda )}\phi_{ - \varepsilon,\lambda} ( x ),& x \in [ a,\theta_{ - \varepsilon} ), \\ \frac{\phi_{\varepsilon,\lambda} ( x )}{D_{1}\omega ( \lambda )}\int_{x}^{\theta_{ + \varepsilon}} \chi_{\varepsilon,\lambda} ( t )f ( t )\,dt + \frac{\chi_{\varepsilon,\lambda} ( x )}{D_{1}\omega ( \lambda )}\int_{\theta_{ - \varepsilon}}^{x} \phi_{\varepsilon,\lambda} ( t )f ( t )\,dt \\ \quad{}+ \frac{\chi_{\varepsilon,\lambda} ( x )}{\omega ( \lambda )}\int_{a}^{\theta_{ - \varepsilon}} \phi_{ - \varepsilon,\lambda} ( t )f ( t )\,dt \\ \quad{}+ \frac{\phi_{\varepsilon,\lambda} ( x )}{D_{1}D_{2}\omega ( \lambda )}\int_{\theta_{ + \varepsilon}}^{b} \chi_{ + \varepsilon,\lambda} ( t )f ( t )\,dt + \frac{h}{D_{1}D_{2}\omega ( \lambda )}\phi_{\varepsilon,\lambda} ( x ),& x \in ( \theta_{ - \varepsilon},\theta_{ + \varepsilon} ), \\ \frac{\phi_{ + \varepsilon,\lambda} ( x )}{D_{1}D_{2}\omega ( \lambda )}\int_{x}^{b} \chi_{ + \varepsilon,\lambda} ( t )f ( t )\,dt + \frac{\chi_{ + \varepsilon,\lambda} ( x )}{D_{1}D_{2}\omega ( \lambda )}\int_{\theta_{ + \varepsilon}}^{x} \phi_{ + \varepsilon,\lambda} ( t )f ( t )\,dt \\ \quad{}+ \frac{\chi_{ + \varepsilon,\lambda} ( x )}{\omega ( \lambda )}\int_{a}^{\theta_{ - \varepsilon}} \phi_{ - \varepsilon,\lambda} ( t )f ( t )\,dt \\ \quad{}+ \frac{\chi_{ + \varepsilon,\lambda} ( x )}{D_{1}\omega ( \lambda )}\int_{\theta_{ - \varepsilon}}^{\theta_{ + \varepsilon}} \phi_{\varepsilon,\lambda} ( t )f ( t )\,dt+ \frac{h}{D_{1}D_{2}\omega ( \lambda )}\phi_{ + \varepsilon,\lambda} ( x ),& x \in ( \theta_{ + \varepsilon},b ]. \end{cases} $$
(58)

Then (58) can be rewritten in the form

$$\begin{aligned} y ( x,\lambda ) =& \int_{a}^{\theta_{ - \varepsilon}} G ( x,t;\lambda )f ( t )\,dt + \frac{1}{D_{1}} \int_{\theta_{ - \varepsilon}}^{\theta_{ + \varepsilon}} G ( x,t;\lambda )f ( t )\,dt \\ &{}+ \frac{1}{D_{1}D_{2}} \int_{\theta_{ + \varepsilon}}^{b} G ( x,t;\lambda )f ( t )\,dt + \frac{h\phi_{\lambda} ( x )}{D_{1}D_{2}\omega ( \lambda )}, \end{aligned}$$
(59)

where

$$ G ( x,t;\lambda ) = \textstyle\begin{cases} \frac{\phi_{\lambda} ( t )\chi_{\lambda} ( x )}{\omega ( \lambda )},& a \le t \le x \le b, x \ne \theta_{ - \varepsilon},\theta_{ + \varepsilon}; t \ne \theta_{ - \varepsilon},\theta_{ + \varepsilon}, \\ \frac{\phi_{\lambda} ( x )\chi_{\lambda} ( t )}{\omega ( \lambda )},& a \le x \le t \le b, x \ne \theta_{ - \varepsilon},\theta_{ + \varepsilon}; t \ne \theta_{ - \varepsilon},\theta_{ + \varepsilon}, \end{cases} $$
(60)

is the Green’s function of the problem (1)-(7).

Hence, we have

$$\begin{aligned} Y = &( \lambda \mathrm{I} - A )^{ - 1}F \\ =& \begin{pmatrix} \int_{a}^{\theta_{ - \varepsilon}} G ( x,t;\lambda )f ( t )\,dt + \frac{1}{D_{1}}\int_{\theta_{ - \varepsilon}}^{\theta_{ + \varepsilon}} G ( x,t;\lambda )f ( t )\,dt + \frac{1}{D_{1}D_{2}}\int_{\theta_{ + \varepsilon}}^{b} G ( x,t;\lambda )f ( t )\,dt + \frac{h\phi_{\lambda} ( x )}{D_{1}D_{2}\omega ( \lambda )} \\ R' ( y ) \end{pmatrix} \end{aligned}$$

the resolvent operator of the problem (1)-(7).

6 Example

We indicate in this example the effect on determining the discontinuities of different values of θ and ε. For each value of θ and ε, we shall display the characteristic function and give the first five eigenvalues of the problem.

Example

Consider the boundary value problem

$$\begin{aligned}& - y'' = \lambda y, \end{aligned}$$
(61)
$$\begin{aligned}& y' ( - 2 ) = 0, \end{aligned}$$
(62)
$$\begin{aligned}& \lambda y ( 4 ) - y' ( 4 ) = 0, \end{aligned}$$
(63)
$$\begin{aligned}& y ( \theta_{ - \varepsilon} + ) = 2y ( \theta_{ - \varepsilon} - ), \end{aligned}$$
(64)
$$\begin{aligned}& y' ( \theta_{ - \varepsilon} + ) = y ( \theta_{ - \varepsilon} - ) + 2^{ - 1}y' ( \theta_{ - \varepsilon} - ), \end{aligned}$$
(65)
$$\begin{aligned}& y ( \theta_{ + \varepsilon} + ) = 2^{ - 1}y ( \theta_{ + \varepsilon} - ), \end{aligned}$$
(66)
$$\begin{aligned}& y' ( \theta_{ + \varepsilon} + ) = y ( \theta_{ + \varepsilon} - ) + 2y' ( \theta_{ + \varepsilon} - ), \end{aligned}$$
(67)

where \(I = [ - 2,4 ]\).

Let \(\lambda = s^{2}\). The eigenvalues of the problem (61)-(67) are the squares of the zeros of the characteristic function of \(\omega ( \lambda )\), given by

$$\begin{aligned} \omega ( \lambda ) =& s^{2}\biggl\{ \cos \bigl( s ( \theta_{ - \varepsilon} + 2 ) \bigr)\cos \bigl( s ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} ) \bigr)\cos \bigl( s ( 4 - \theta_{ + \varepsilon} ) \bigr) \\ &{} - \frac{1}{4} \sin \bigl( s ( \theta_{ - \varepsilon} + 2 ) \bigr)\sin \bigl( s ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} ) \bigr) \cos \bigl( s ( 4 - \theta_{ + \varepsilon} ) \bigr) \\ &{}- 4\cos \bigl( s ( \theta_{ - \varepsilon} + 2 ) \bigr)\sin \bigl( s ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} ) \bigr) \sin \bigl( s ( 4 - \theta_{ + \varepsilon} ) \bigr) \\ &{}- \sin \bigl( s ( \theta_{ - \varepsilon} + 2 ) \bigr)\cos \bigl( s ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} ) \bigr)\sin \bigl( s ( 4 - \theta_{ + \varepsilon} ) \bigr) \biggr\} \\ &{}+ s\biggl\{ 5\cos \bigl( s ( \theta_{ - \varepsilon} + 2 ) \bigr)\cos \bigl( s ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} ) \bigr)\sin \bigl( s ( 4 - \theta_{ + \varepsilon} ) \bigr) \\ &{}- \frac{3}{4}\sin \bigl( s ( \theta_{ - \varepsilon} + 2 ) \bigr) \sin \bigl( s ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} ) \bigr) \sin \bigl( s ( 4 - \theta_{ + \varepsilon} ) \bigr) \\ &{}+ \frac{9}{2}\cos \bigl( s ( \theta_{ - \varepsilon} + 2 ) \bigr)\sin \bigl( s ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} ) \bigr) \cos \bigl( s ( 4 - \theta_{ + \varepsilon} ) \bigr) \\ &{} + \sin \bigl( s ( \theta_{ - \varepsilon} + 2 ) \bigr)\cos \bigl( s ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} ) \bigr)\cos \bigl( s ( 4 - \theta_{ + \varepsilon} ) \bigr) \biggr\} \\ &{}+ \frac{3}{2}\cos \bigl( s ( \theta_{ - \varepsilon} + 2 ) \bigr)\sin \bigl( s ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} ) \bigr)\sin \bigl( s ( 4 - \theta_{ + \varepsilon} ) \bigr) \\ &{}- 4\cos \bigl( s ( \theta_{ - \varepsilon} + 2 ) \bigr) \cos \bigl( s ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} ) \bigr) \cos \bigl( s ( 4 - \theta_{ + \varepsilon} ) \bigr) \\ &{} + \frac{1}{2}\sin \bigl( s ( \theta_{ - \varepsilon} + 2 ) \bigr)\sin \bigl( s ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} ) \bigr) \cos \bigl( s ( 4 - \theta_{ + \varepsilon} ) \bigr) \\ &{}- \frac{1}{s} \cos \bigl( s ( \theta_{ - \varepsilon} + 2 ) \bigr)\sin \bigl( s ( \theta_{ + \varepsilon} - \theta_{ - \varepsilon} ) \bigr)\cos \bigl( s ( 4 - \theta_{ + \varepsilon} ) \bigr) \\ =& 0. \end{aligned}$$
(68)

(1) Let \(\theta = - 1\) be an interior point in \(I = [ - 2,4 ]\). From (8), we get \(0 < \varepsilon < 1\).

(i) If \(\varepsilon = 1/4\), then the points of discontinuity are \(\theta_{ - \varepsilon} = - 5/4\) and \(\theta_{ + \varepsilon} = - 3/4\). Equation (68) is then reduced to

$$\begin{aligned} \omega ( \lambda ) =& s^{2}\biggl\{ \cos \biggl( \frac{3}{4}s \biggr)\cos \biggl( \frac{1}{2}s \biggr)\cos \biggl( \frac{19}{4}s \biggr) - \frac{1}{4} \sin \biggl( \frac{3}{4}s \biggr)\sin \biggl( \frac{1}{2}s \biggr)\cos \biggl( \frac{19}{4}s \biggr) \\ &{}- 4\cos \biggl( \frac{3}{4}s \biggr)\sin \biggl( \frac{1}{2}s \biggr)\sin \biggl( \frac{19}{4}s \biggr) - \sin \biggl( \frac{3}{4}s \biggr)\cos \biggl( \frac{1}{2}s \biggr)\sin \biggl( \frac{19}{4}s \biggr) \biggr\} \\ &{}+ s\biggl\{ 5\cos \biggl( \frac{3}{4}s \biggr)\cos \biggl( \frac{1}{2}s \biggr)\sin \biggl( \frac{19}{4}s \biggr) - \frac{3}{4}\sin \biggl( \frac{3}{4}s \biggr)\sin \biggl( \frac{1}{2}s \biggr)\sin \biggl( \frac{19}{4}s \biggr) \\ &{}+ \frac{9}{2}\cos \biggl( \frac{3}{4}s \biggr)\sin \biggl( \frac{1}{2}s \biggr) \cos \biggl( \frac{19}{4}s \biggr) + \sin \biggl( \frac{3}{4}s \biggr)\cos \biggl( \frac{1}{2}s \biggr)\cos \biggl( \frac{19}{4}s \biggr) \biggr\} \\ &{}+ \frac{3}{2}\cos \biggl( \frac{3}{4}s \biggr)\sin \biggl( \frac{1}{2}s \biggr)\sin \biggl( \frac{19}{4}s \biggr) - 4\cos \biggl( \frac{3}{4}s \biggr)\cos \biggl( \frac{1}{2}s \biggr)\cos \biggl( \frac{19}{4}s \biggr) \\ &{}+ \frac{1}{2}\sin \biggl( \frac{3}{4}s \biggr)\sin \biggl( \frac{1}{2}s \biggr)\cos \biggl( \frac{19}{4}s \biggr) - \frac{1}{s}\cos \biggl( \frac{3}{4}s \biggr)\sin \biggl( \frac{1}{2}s \biggr)\cos \biggl( \frac{19}{4}s \biggr) \\ =& 0. \end{aligned}$$
(69)

The graph of the characteristic function for \(\varepsilon = 1/4\) is displayed in Figure 1.

Figure 1
figure 1

The graph of the characteristic function for \(\pmb{\theta = - 1 }\) and \(\pmb{\varepsilon = 1/4}\) .

(ii) If \(\varepsilon = 1/2\), then the points of discontinuity are \(\theta_{ - \varepsilon} = - 3/2\) and \(\theta_{ + \varepsilon} = - 1/2\). Equation (68) is then reduced to

$$\begin{aligned} \omega ( \lambda ) =& s^{2}\biggl\{ \cos \biggl( \frac{1}{2}s \biggr)\cos ( s )\cos \biggl( \frac{9}{2}s \biggr) - \frac{1}{4} \sin \biggl( \frac{1}{2}s \biggr)\sin ( s )\cos \biggl( \frac{9}{2}s \biggr) \\ &{}- 4\cos \biggl( \frac{1}{2}s \biggr)\sin ( s )\sin \biggl( \frac{9}{2}s \biggr) - \sin \biggl( \frac{1}{2}s \biggr)\cos ( s ) \sin \biggl( \frac{9}{2}s \biggr) \biggr\} \\ &{}+ s\biggl\{ 5\cos \biggl( \frac{1}{2}s \biggr)\cos ( s )\sin \biggl( \frac{9}{2}s \biggr) - \frac{3}{4}\sin \biggl( \frac{1}{2}s \biggr)\sin ( s )\sin \biggl( \frac{9}{2}s \biggr) \\ &{}+ \frac{9}{2}\cos \biggl( \frac{1}{2}s \biggr)\sin ( s ) \cos \biggl( \frac{9}{2}s \biggr) + \sin \biggl( \frac{1}{2}s \biggr)\cos ( s )\cos \biggl( \frac{9}{2}s \biggr) \biggr\} \\ &{}+ \frac{3}{2}\cos \biggl( \frac{1}{2}s \biggr)\sin ( s )\sin \biggl( \frac{9}{2}s \biggr) - 4\cos \biggl( \frac{1}{2}s \biggr) \cos ( s )\cos \biggl( \frac{9}{2}s \biggr) \\ &{}+ \frac{1}{2}\sin \biggl( \frac{1}{2}s \biggr)\sin ( s )\cos \biggl( \frac{9}{2}s \biggr) - \frac{1}{s}\cos \biggl( \frac{1}{2}s \biggr)\sin ( s )\cos \biggl( \frac{9}{2}s \biggr) \\ =& 0. \end{aligned}$$
(70)

The graph of the characteristic function for \(\varepsilon = 1/2\) is displayed in Figure 2.

Figure 2
figure 2

The graph of the characteristic function for \(\pmb{\theta = - 1 }\) and \(\pmb{\varepsilon = 1/2}\) .

(2) Let \(\theta = 1\) be midpoint of the interval \(I = [ - 2,4 ]\). From (8), we get \(0 < \varepsilon < 3\).

(i) If \(\varepsilon = 1\), then the points of discontinuity are \(\theta_{ - \varepsilon} = 0\) and \(\theta_{ + \varepsilon} = 2\). Equation (68) is then reduced to

$$\begin{aligned} \omega ( \lambda ) =& s^{2}\biggl\{ \cos^{3} ( 2s ) - \frac{21}{4}\cos ( 2s )\sin^{2} ( 2s ) \biggr\} \\ &{}+ s\biggl\{ \frac{21}{2}\sin ( 2s )\cos^{2} ( 2s )- \frac{3}{4}\sin^{3} ( 2s ) \biggr\} \\ &{}+ 2\cos ( 2s ) \sin^{2} ( 2s ) - 4\cos^{3} ( 2s ) - \frac{1}{s}\sin ( 2s )\cos^{2} ( 2s ) \\ =& 0. \end{aligned}$$
(71)

The graph of the characteristic function for \(\varepsilon = 1\) is displayed in Figure 3.

Figure 3
figure 3

The graph of the characteristic function for \(\pmb{\theta = 1 }\) and \(\pmb{\varepsilon = 1}\) .

(ii) If \(\varepsilon = 3/2\), then the points of discontinuity are \(\theta_{ - \varepsilon} = - 1/2\) and \(\theta_{ + \varepsilon} = 5/2\). The equation (68) is then reduced to

$$\begin{aligned} \omega ( \lambda ) =& s^{2}\biggl\{ \cos^{2} ( 3s ) - \frac{17}{8}\sin^{2} ( 3s ) \biggr\} \\ &{}+ s\biggl\{ \frac{3}{2}\sin ( 6s ) - \frac{3}{4}\sin ( 3s )\sin^{2} \biggl( \frac{3}{2}s \biggr)+ \frac{9}{2}\sin ( 3s )\cos^{2} \biggl( \frac{3}{2}s \biggr) \biggr\} \\ &{}+ \sin^{2} ( 3s )- 4\cos ( 3s )\cos^{2} \biggl( \frac{3}{2}s \biggr)- \frac{1}{s}\sin ( 3s )\cos^{2} \biggl( \frac{3}{2}s \biggr) \\ =& 0. \end{aligned}$$
(72)

The graph of the characteristic function for \(\varepsilon = 3/2\) is displayed in Figure 4.

Figure 4
figure 4

The graph of the characteristic function for \(\pmb{\theta = 1 }\) and \(\pmb{\varepsilon = 3/2}\) .

(3) Let \(\theta = 3\) be an interior point in \(I = [ - 2,4 ]\). From (8), we get \(0 < \varepsilon < 1\).

(i) If \(\varepsilon = 1/4\), then the points of discontinuity are \(\theta_{ - \varepsilon} = 11/4\) and \(\theta_{ + \varepsilon} = 13/4\). Equation (68) is then reduced to

$$\begin{aligned} \omega ( \lambda ) =& s^{2}\biggl\{ \cos \biggl( \frac{19}{4}s \biggr)\cos \biggl( \frac{1}{2}s \biggr)\cos \biggl( \frac{3}{4}s \biggr) - \frac{1}{4} \sin \biggl( \frac{19}{4}s \biggr)\sin \biggl( \frac{1}{2}s \biggr)\cos \biggl( \frac{3}{4}s \biggr) \\ &{}- 4\cos \biggl( \frac{19}{4}s \biggr)\sin \biggl( \frac{1}{2}s \biggr)\sin \biggl( \frac{3}{4}s \biggr) - \sin \biggl( \frac{19}{4}s \biggr)\cos \biggl( \frac{1}{2}s \biggr)\sin \biggl( \frac{3}{4}s \biggr) \biggr\} \\ &{}+ s\biggl\{ 5\cos \biggl( \frac{19}{4}s \biggr)\cos \biggl( \frac{1}{2}s \biggr)\sin \biggl( \frac{3}{4}s \biggr) - \frac{3}{4}\sin \biggl( \frac{19}{4}s \biggr)\sin \biggl( \frac{1}{2}s \biggr)\sin \biggl( \frac{3}{4}s \biggr) \\ &{}+ \frac{9}{2}\cos \biggl( \frac{19}{4}s \biggr)\sin \biggl( \frac{1}{2}s \biggr) \cos \biggl( \frac{3}{4}s \biggr) + \sin \biggl( \frac{19}{4}s \biggr)\cos \biggl( \frac{1}{2}s \biggr)\cos \biggl( \frac{3}{4}s \biggr) \biggr\} \\ &{}+ \frac{3}{2}\cos \biggl( \frac{19}{4}s \biggr)\sin \biggl( \frac{1}{2}s \biggr)\sin \biggl( \frac{3}{4}s \biggr) - 4\cos \biggl( \frac{19}{4}s \biggr)\cos \biggl( \frac{1}{2}s \biggr)\cos \biggl( \frac{3}{4}s \biggr) \\ &{}+ \frac{1}{2}\sin \biggl( \frac{19}{4}s \biggr)\sin \biggl( \frac{1}{2}s \biggr)\cos \biggl( \frac{3}{4}s \biggr) - \frac{1}{s}\cos \biggl( \frac{19}{4}s \biggr)\sin \biggl( \frac{1}{2}s \biggr)\cos \biggl( \frac{3}{4}s \biggr) \\ =& 0. \end{aligned}$$
(73)

The graph of the characteristic function for \(\varepsilon = 1/4\) is displayed in Figure 5.

Figure 5
figure 5

The graph of the characteristic function for \(\pmb{\theta = 3 }\) and \(\pmb{\varepsilon = 1/4}\) .

(ii) If \(\varepsilon = 1/2\), then the points of discontinuity are \(\theta_{ - \varepsilon} = 5/2\) and \(\theta_{ + \varepsilon} = 7/2\). Equation (68) is then reduced to

$$\begin{aligned} \omega ( \lambda ) =& s^{2}\biggl\{ \cos \biggl( \frac{9}{2}s \biggr)\cos ( s )\cos \biggl( \frac{1}{2}s \biggr) - \frac{1}{4} \sin \biggl( \frac{9}{2}s \biggr)\sin ( s )\cos \biggl( \frac{1}{2}s \biggr) \\ &{}- 4\cos \biggl( \frac{9}{2}s \biggr)\sin ( s )\sin \biggl( \frac{1}{2}s \biggr) - \sin \biggl( \frac{9}{2}s \biggr)\cos ( s ) \sin \biggl( \frac{1}{2}s \biggr) \biggr\} \\ &{}+ s\biggl\{ 5\cos \biggl( \frac{9}{2}s \biggr)\cos ( s )\sin \biggl( \frac{1}{2}s \biggr) - \frac{3}{4}\sin \biggl( \frac{9}{2}s \biggr)\sin ( s )\sin \biggl( \frac{1}{2}s \biggr) \\ &{}+ \frac{9}{2}\cos \biggl( \frac{9}{2}s \biggr)\sin ( s ) \cos \biggl( \frac{1}{2}s \biggr) + \sin \biggl( \frac{9}{2}s \biggr)\cos ( s )\cos \biggl( \frac{1}{2}s \biggr) \biggr\} \\ &{}+ \frac{3}{2}\cos \biggl( \frac{9}{2}s \biggr)\sin ( s )\sin \biggl( \frac{1}{2}s \biggr) - 4\cos \biggl( \frac{9}{2}s \biggr) \cos ( s )\cos \biggl( \frac{1}{2}s \biggr) \\ &{}+ \frac{1}{2}\sin \biggl( \frac{9}{2}s \biggr)\sin ( s )\cos \biggl( \frac{1}{2}s \biggr) - \frac{1}{s}\cos \biggl( \frac{9}{2}s \biggr)\sin ( s )\cos \biggl( \frac{1}{2}s \biggr) \\ =& 0. \end{aligned}$$
(74)

The graph of the characteristic function for \(\varepsilon = 1/2\) is displayed in Figure 6.

Figure 6
figure 6

The graph of the characteristic function for \(\pmb{\theta = 3 }\) and \(\pmb{\varepsilon = 1/2}\) .

References

  1. Buschmann, D, Stolz, G, Weidmann, J: One-dimensional Schrödinger operators with local point interactions. J. Reine Angew. Math. 467, 169-186 (1995)

    MathSciNet  MATH  Google Scholar 

  2. Titeux, I, Yakubov, Y: Completeness of root functions for thermal conduction in a strip with piecewise continuous coefficients. Math. Models Methods Appl. Sci. 7(7), 1035-1050 (1997)

    Article  MathSciNet  MATH  Google Scholar 

  3. Khmelnytskaya, KV, Serroukh, I: The heat transfer problem for inhomogeneous material in photoacoustic applications and spectral parameter power series. Math. Methods Appl. Sci. (2013). doi:10.1002/mma.2732

    MathSciNet  Google Scholar 

  4. Tittle, CW: Boundary value problems in composite media. J. Appl. Phys. 36, 1486-1488 (1965)

    Article  MATH  Google Scholar 

  5. Mulholland, GP, Cobble, MH: Diffusion through composite media. Int. J. Heat Mass Transf. 15(1), 147-160 (1972)

    Article  Google Scholar 

  6. Özisik, MN: Heat Conduction. Wiley, New York (1980)

    Google Scholar 

  7. Oturanç, G, Sahin, AZ: Eigenvalue analysis of temperature distribution in composite walls. Int. J. Energy Res. 25, 1189-1196 (2001)

    Article  Google Scholar 

  8. Hinton, BD: An expansion theorem for an eigenvalue problem with eigenvalue parameter in the boundary conditions. Q. J. Math. 30, 33-42 (1979)

    Article  MathSciNet  MATH  Google Scholar 

  9. Fulton, CT: Two-point boundary value problems with eigenvalues parameter contained in the boundary conditions. Proc. R. Soc. Edinb. A 77A, 293-308 (1977)

    Article  MathSciNet  Google Scholar 

  10. Walter, J: Regular eigenvalue problems with eigenvalue parameter in the boundary conditions. Math. Z. 133, 301-312 (1973)

    Article  MathSciNet  MATH  Google Scholar 

  11. Willis, C: Inverse Sturm Liouville problems with two discontinuities. Inverse Probl. 1, 263-289 (1985)

    Article  MathSciNet  MATH  Google Scholar 

  12. Kobayashi, M: Eigenfunction expansions: a discontinuous version. SIAM J. Appl. Math. 50(3), 910-917 (1990)

    Article  MathSciNet  MATH  Google Scholar 

  13. Altınışık, N, Kadakal, M, Mukhtarov, OS: Eigenvalues and eigenfunctions of discontinuous Sturm Liouville problems with eigenparameter dependent boundary conditions. Acta Math. Hung. 102, 159-175 (2004)

    Article  MATH  Google Scholar 

  14. Mukhtarov, OS, Kadakal, M: Some spectral properties of one Sturm-Liouville type problem with discontinuous weight. Sib. Math. Zh. 46(4), 860-875 (2004) (in Russian). Translation in Sib. Math. J. 46(4) 681-694 (2005)

    MathSciNet  Google Scholar 

  15. Kadakal, M, Mukhtarov, OS: Sturm-Liouville problems with discontinuities at two points. Comput. Math. Appl. 54, 1367-1379 (2007)

    Article  MathSciNet  MATH  Google Scholar 

  16. Akdoğan, Z, Demirci, M, Mukhtarov, OS: Green function of discontinuous boundary value problem with transmission conditions. Math. Models Methods Appl. Sci. 30, 1719-1738 (2007)

    Article  MATH  Google Scholar 

  17. Wang, A, Sun, J, Hao, X, Yao, S: Completeness of eigenfunctions of Sturm Liouville problems with transmission conditions. Methods Appl. Anal. 16(3), 299-312 (2009)

    Article  MathSciNet  MATH  Google Scholar 

  18. Hıra, F, Altınışık, N: Sturm Liouville problem with discontinuity conditions at several points. AIP Conf. Proc. 1470, 126-129 (2012)

    Article  Google Scholar 

  19. Tharwat, MM, Bhrawy, AH, Alofi, AS: Approximation of eigenvalues of discontinuous Sturm-Liouville problems with eigenparameter in all boundary conditions. Bound. Value Probl. (2013). doi:10.1186/1687-2770-2013-132

    MathSciNet  Google Scholar 

  20. Titchmarsh, EC: Eigenfunctions Expansion Associated with Second Order Differential Equations I, 2nd edn. Oxford University Press, London (1962)

    Google Scholar 

Download references

Acknowledgements

All authors are very grateful to the anonymous referees for their valuable suggestions. This paper was supported as a scientific research project (PYO.FEN.1904.11.010) by Ondokuz Mayıs University.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Fatma Hıra.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

All authors contributed equally to the writing of this paper. All authors read and approved the final manuscript.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hıra, F., Altınışık, N. Sturm-Liouville problem with moving discontinuity points. Bound Value Probl 2015, 237 (2015). https://doi.org/10.1186/s13661-015-0502-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13661-015-0502-6

MSC

Keywords