\documentclass[reqno]{amsart}
\usepackage{hyperref}

\AtBeginDocument{{\noindent\small
\emph{Electronic Journal of Differential Equations},
Vol. 2017 (2017), No. 11, pp. 1--12.\newline
ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu}
\thanks{\copyright 2017 Texas State University.}
\vspace{8mm}}

\begin{document}
\title[\hfilneg EJDE-2017/11\hfil Nonlocal Sturm-Liouville problems]
{Nonlocal Sturm-Liouville problems with integral terms in
the boundary conditions}

\author[M. Kandemir, O. Sh. Mukhtarov \hfil EJDE-2017/11\hfilneg]
{Mustafa Kandemir, Oktay Sh. Mukhtarov}

\address{Mustafa Kandemir \newline
Department of Mathematics,
Education Faculty,
 Amasya University,
Amasya, Turkey}
\email{mkandemir5@yahoo.com}

\address{Oktay Sh. Mukhtarov \newline
Department of Mathematics,
Faculty of Science and Arts,
Gaziosmanpasa University, 60100 Tokat, Turkey. \newline
Institute of Mathematics and Mechanics,
Azerbaijan National Academy of Sciences,
Baku, Azerbaijan}
\email{omukhtarov@yahoo.com}

\dedicatory{Communicated by Ludmila S. Pulkina}

\thanks{Submitted October 11, 2016. Published January 12, 2017.}
\subjclass[2010]{34A36, 34B08, 34B24}
\keywords{ Sturm-Liouville problem; nonlocal boundary conditions;
\hfill\break\indent coercive; solvability; Fredholmness}

\begin{abstract}
 We consider a new type Sturm-Liouville problems whose main 
 feature is the nature of boundary conditions.
 Namely, we study the nonhomogeneous Sturm-Liouville equation
 $$
 p(x)u''(x)+(q(x)-\lambda )u=f(x)
 $$
 on two disjoint intervals $[-1,0)$  and  $(0,1]$, subject to the
 nonlocal boundary-transmission conditions
\begin{align*}
 &\alpha _ku^{(m_k)}(-1)+\beta _ku^{(m_k)}(-0)+\eta
 _ku^{(m_k)}(+0)+\gamma _ku^{(m_k)}(1)   \\
 & +\sum_{j=1}^{n_k}\delta _{kj}u^{(m_k)}(x_{kj})+\sum_{\upsilon
 =1}^{2}\sum_{j=0}^{m_k}\int_{\Omega _{\upsilon }}\mathcal{K}
 _{k\upsilon j}(t)u^{(j)}(t)dt=f_k,\quad k=1,2,3,4.
 \end{align*}
 where $\Omega _1:=[-1,0)$, $\Omega _2:=(0,1]$ and
 $x_{kj}\in (-1,0)\cup (0,1)$ are internal points. By using our own
 approaches we establish such important properties as Fredholmness,
 coercive solvability and isomorphism with respect to the spectral
 parameter $\lambda$.
\end{abstract}

\maketitle
\numberwithin{equation}{section}
\newtheorem{theorem}{Theorem}[section]
\allowdisplaybreaks

\section{Introduction}

Various generalizations of classical Sturm-Liouville problems for
ordinary linear differential equations have attracted a lot of
attention because of the appearance of new important
applications in physical sciences and applied mathematics. For
instance, theoretical investigations have  become interested in
the discontinuous Sturm-Liouville problems for its application in
physics. The discontinuity of the coefficients of the equations in
the Sturm-Liouville problems corresponds to the fact that the
heterogeneous media consists of two different materials. On the
other hand, transmission problems appear frequently  in various
fields of physics  such as in electrostatics, magnetostatics and in
solid mechanic for discontinuous problems (in these regard see,
\cite{bor,mar}). Solvability and some spectral properties of
nonlocal Sturm-Liouville problems have been investigated by many
authors; see for example, \cite{all,all1,ima,ism,kad, mam,sha,ugu}).
An important special case of the nonlocal Sturm-Liouville problems
are so-called multipoint Sturm-Liouville problems. Such problems
have been extensively studied by many authors; see for
example,\cite{gar,gen,huy} and references therein.

 In general, the mathematical problems encountered in the
study of boundary value transmission problems or nonclassical
problems cannot be treated with the usual techniques within the
standard framework of Sturm-Liouville problems. In classical theory
of boundary-value problems for ordinary differential equations is
usually considered for equations with continuous coefficients and
for boundary conditions which contain only endpoints of the
considered interval. This article deals with one nonclassical
boundary-value problem for a second-order ordinary differential
equation with discontinuous coefficients and boundary conditions
containing not only endpoints of the considered interval, but also
a finite number of internal points and integral terms. Namely, we
consider the differential equation
\begin{equation}\label{1}
L(\lambda )u:=p(x)u''(x)+(q(x)-\lambda )u(x)=f(x),\quad x\in
[ -1,0)\cup (0,1]
\end{equation}
together with new type boundary conditions
\begin{equation}\label{2}
\begin{aligned}
L_ku&:=\alpha _ku^{(m_k)}(-1)+\beta _ku^{(m_k)}(-0)+\eta
_ku^{(m_k)}(+0)+\gamma _ku^{(m_k)}(1)   \\
&\quad +\sum_{j=1}^{n_k}\delta _{kj}u^{(m_k)}(x_{kj})+\sum_{\upsilon
=1}^{2}\sum_{j=0}^{m_k}\int_{\Omega _{\upsilon }}\mathcal{K}
_{k\upsilon j}(t)u^{(j)}(t)dt=f_k,
\end{aligned}
\end{equation}
for $k=1,2,3,4$, where $p(x)$ is piecewise constant function, $p(x)=p_1$ for
$x\in[-1,0)$, $p(x)=p_2$ for $x\in (0,1]$; $\lambda $-complex parameter;
$p_i$ $(i=1,2)$, $\alpha _k$, $\beta _k$, $\eta _k$,
$\gamma_k,\delta _{ki} $ ($i=1,2$, $k=1,2,3,4$) are complex coefficients;
$m_k$ $(k=1,2,3,4)$ are integers;
$\Omega _1:=(-1,0)$, $\Omega _2:=(0,1)$;
$\mathcal{K}_{k\upsilon j}\in W_q^{m_k}(-1,0)\dot{+}W_q^{m_k}(0,1)$;
$x_{kj}\in (-1,0)\cup (0,1)$ are internal points and $q(x)$ is measurable
function on $[-1,0)\cup (0,1]$. Naturally, we shall assume that, $p_1\neq 0$,
$p_2\neq 0$ and $| \alpha _k| +|\beta _k| +| \eta _k| +|\gamma _k| \neq 0$
$(k=1,2,3,4)$. Some special cases of
the considered Sturm-Liouville problem \eqref{1}--\eqref{2} arise
after an application of the method of separation of variables to the
varied assortment of physical problems, namely, in heat and mass
transfer problems (see, for example, \cite{lik}), in diffraction
problems (for example, \cite{agr}), in vibrating string problems,
when the string loaded additionally with point masses (see,
\cite{voi}) and etc. Some problems with transmission conditions
which arise in mechanics were studied in \cite{mar,voi}.
Investigation of various spectral properties of some nonlocal
boundary-value problems can be found in some works of  Imanbaev
\cite{ima}, Sadybekov \cite{sad},
 Shakhmurov \cite{sha}, Aliyev \cite{ali} and  Rasulov  \cite{ras}.
 Note that some new type Sturm-Liouville problems with nonlocal
 boundary conditions
were investigated by authors of this paper and some others
 \cite{ayd,ayd1,ayd2,kan,kan1,muk1,muk2,osh5}.

\section{Homogeneous equation with nonhomogeneous transmission conditions}

For convenience we denote
\[
S_ku :=\sum_{j=1}^{n_k}\delta _{kj}u^{(m_k)}(x_{kj}), \quad
\mathcal{F}_ku :=\sum_{\upsilon =1}^{2}
 \sum_{j=0}^{m_k}\int_{\Omega _{\upsilon }}\mathcal{K}_{k\upsilon j}(t)u^{(j)}(t)dt,
\quad k=1,2,3,4\,.
\]
We consider the homogeneous differential equation
\begin{equation}\label{3}
L_0(\lambda )u:=p(x)u''(x)-\lambda u(x)=0
\end{equation}
 with the nonlocal and nonhomogeneous  boundary conditions
\begin{equation}\label{4}
\begin{aligned}
L_{k0}u &:=\alpha _ku^{(m_k)}(-1)+\beta
_ku^{(m_k)}(-0)+\eta_ku^{(m_k)}(+0)   \\
&\quad +\gamma _ku^{(m_k)}(1)+S_ku=f_k,\quad k=1,2,3,4.
\end{aligned}
\end{equation}
For convenience we shall use the notation
\begin{gather*}
\omega_1:=-(p_1^{-1}\lambda )^{1/2},\quad
\omega_2:=(p_1^{-1}\lambda )^{1/2}, \quad
\omega _3:=-(p_2^{-1}\lambda )^{1/2},\quad
\omega _4=(p_2^{-1}\lambda )^{1/2}, \\
\underline{\omega}:=\min \{ \arg p_1,\,\arg p_2\}, \quad
\bar{\omega}:=\max \{ \arg p_1,\,\arg p_2\}, \\
\theta :=
\begin{vmatrix}
{\alpha _1\omega _1^{m_1}} & {\beta _1\omega _2^{m_1}} & {\eta
_1\omega _3^{m_1}} & {\gamma _1\omega _4^{m_1}} \\
{\alpha _2\omega _1^{m_2}} & {\beta _2\omega _2^{m_2}} & {\eta
_2\omega _3^{m_2}} & {\gamma _2\omega _4^{m_2}} \\
{\alpha _3\omega _1^{m_3}} & {\beta _3\omega _2^{m_3}} & {\eta
_3\omega _3^{m_3}} & {\gamma _3\omega _4^{m_3}} \\
{\alpha _4\omega _1^{m_4}} & {\beta _4\omega _2^{m_4}} & {\eta
_4\omega _3^{m_4}} & {\gamma _4\omega _4^{m_4}}
\end{vmatrix},
\\
B_{\varepsilon }(\underline{\omega},\bar{\omega}):=\{ \lambda
\in \mathbb{C}:  \pi +\bar{\omega}+\varepsilon <\arg \lambda
<3\pi +\underline{\omega}-\varepsilon \}
\end{gather*}
for real $\varepsilon >0$ small enough.

The direct sum of Sobolev spaces $W_q^k(-1,0)\dot{+}W_q^k(0,1)$
(for an integer $k\geq 0$ and real $q>1$) is defined as Banach space of
complex-valued functions $u=u(x)$ defined on $[-1,0)\cup (0,1]$ which belong
to $W_q^k(-1,0)$ and $W_q^k(0,1)$ on intervals $(-1,0)$ and $(0,1)$
respectively, with the norm
\begin{equation*}
\| u\| _{q,k}=\| u\|_{W_q^k(-1,0)}+\| u\| _{W_q^k(0,1)}.
\end{equation*}
Here, as usual, $W_q^k(a,b)$ is the Sobolev space, i.e. the
Banach space consisting of all measurable functions $u(x)$ that have
generalized derivatives on the interval $(a,b)$  up to $k$-th order
inclusive with the finite norm
\begin{equation*}
\| u\| _{W_q^k(a,b)}=\sum_{i=0}^k\Big(
\int_{a}^{b}| u^{( i) }(x)| ^{q}dx\Big) ^{1/q}.
\end{equation*}

\begin{theorem} \label{thm1}
If $\theta \neq 0$ then for any $\varepsilon >0$ there exist
$\rho _{\varepsilon }>0$ such that for all
$\lambda \in B_{\varepsilon }(\underline{\omega},\bar{\omega})$ for which
$| \lambda| >\rho _{\varepsilon }$, the problem
\eqref{3}-\eqref{4} has a unique solution $u(x,\lambda )$ that belongs to
$W_q^{l}(-1,0)\dot{+}W_q^{l}(0,1)$ for arbitrary
$l\geq \max \{2,\max
\left\{ m_1,m_2,m_3,m_4\right\} +1\}$ and for these $\lambda $ the
coercive estimate
\begin{equation}\label{5}
\sum_{k=0}^{l}| \lambda | ^{{l-k}}\| u\|_{q,k}
\leq C(\varepsilon )\sum_{j=0}^{4}| \lambda | ^{{l-m_{j}-}
\frac{{1}}{q}}| f_{\upsilon }|
\end{equation}
is valid.
\end{theorem}

\begin{proof}
Let $\lambda =\mu ^{2}$. Let us define four
basic solutions $u_i=u_i(x,\mu )$ ($i=1,2,3,4$) of  \eqref{3} as
\begin{equation*}%\label{6}
u_i(x,\mu ):=
\begin{cases}
\exp(\omega _i\mu (x-\xi _i))& \text{for } x\in I_i \\
0 & \text{for } x\notin I_i,
\end{cases}
\end{equation*}
where, $\xi _1=-1$, $\xi _2=\xi _3=0$, $\xi _4=1$; $j=1$ for
$i=1,2$ and $j=2$ for $i=3,4;$ $I_1=I_2=[-1,0)$,
$I_3=I_4=(0,1]$. Then the general
solution of  \eqref{3} can be written in the form
\begin{equation} \label{6}
u(x,\mu )=\sum_{k=1}^{4}C_ku_k(x,\mu ).
\end{equation}
Substituting this expression into  \eqref{4}
yields the following  system of linear homogeneous equations with respect
to variables  $C_1,C_2,C_3,C_4$:
\begin{equation}\label{7}
\begin{aligned}
&(\omega _1\mu )^{m_k}(\alpha _k+\beta _ke^{\omega _1\mu
})C_1+(\omega _2\mu )^{m_k}(\alpha _ke^{-\omega _2\mu
}+\beta
_k)C_2   \\
&+(\omega _3\mu )^{m_k}(\eta _k+\gamma _ke^{\omega _3\mu
})C_3+(\omega _4\mu )^{m_k}(\eta _ke^{-\omega _4\mu }+\gamma
_k)C_4
=f_k,\quad k=1,2,3,4.
\end{aligned}
\end{equation}
From $\lambda \in B_{\varepsilon }(\underline{\omega},\bar{\omega})$ it
follows that
\begin{gather*}
\frac{\pi +\varepsilon }{2}<\arg (\omega _i\mu )<\frac{3\pi -\varepsilon }{2}
\quad\text{for } i=1,3;\\
-\frac{\pi -\varepsilon }{2}<\arg (\omega _i\mu )<\frac{\pi -\varepsilon }{
2} \quad\text{for } i=2,4.
\end{gather*}
Consequently, for these $\lambda $ and for  $\varepsilon >0$ (small enough), we have
\begin{equation*}
(-1)^{k+1}\operatorname{Re}(\omega _k\mu )\leq -| \lambda | |
\omega _k| \sin \frac{\varepsilon }{2},\quad k=1,2,3,4.
\end{equation*}
Hence, the determinant of the system \eqref{7} has the form
\begin{align*}
\Delta (\lambda )
&=\lambda^{\frac{1}{2}\sum_{i=1}^{4}m_i}
\Bigg( \begin{vmatrix}
{\alpha _1\omega _1^{m_1}} & {\beta _1\omega _2^{m_1}} & {\eta
_1\omega _3^{m_1}} & {\gamma _1\omega _4^{m_1}} \\
{\alpha _2\omega _1^{m_2}} & {\beta _2\omega _2^{m_2}} & {\eta
_2\omega _3^{m_2}} & {\gamma _2\omega _4^{m_2}} \\
{\alpha _3\omega _1^{m_3}} & {\beta _3\omega _2^{m_3}} & {\eta
_3\omega _3^{m_3}} & {\gamma _3\omega _4^{m_3}} \\
{\alpha _4\omega _1^{m_4}} & {\beta _4\omega _2^{m_4}} & {\eta
_4\omega _3^{m_4}} & {\gamma _4\omega _4^{m_4}}
\end{vmatrix} \\
&\quad + e^{\lambda^{1/2}\sum_{i=1}^{4}(-1)^{i+1}\omega _i}
\begin{vmatrix}
{\beta _1\omega _1^{m_1}} & {\alpha _1\omega _2^{m_1}} & {\gamma
_1\omega _3^{m_1}} & {\eta _1\omega _4^{m_1}} \\
{\beta _2\omega _1^{m_2}} & {\alpha _2\omega _2^{m_2}} & {\gamma
_2\omega _3^{m_2}} & {\eta _2\omega _4^{m_2}} \\
{\beta _3\omega _1^{m_3}} & {\alpha _3\omega _2^{m_3}} & {\gamma
_3\omega _3^{m_3}} & {\eta _3\omega _4^{m_3}} \\
{\beta _4\omega _1^{m_4}} & {\alpha _4\omega _2^{m_4}} & {\gamma
_4\omega _3^{m_4}} & {\eta _4\omega _4^{m_4}}
\end{vmatrix}
\Bigg) \\
&=\lambda ^{{m}}( \theta +r(\lambda ))
\end{align*}
where $m=m_1+m_2+m_3+m_4$ and
$ r(\lambda )\to 0$ as  $| \lambda |\to \infty $
in the angle $ B_{\varepsilon }(\underline{\omega},\bar{\omega})$.
Since $\theta \neq 0$, there
exist $\rho _{\varepsilon }>0$ such that for all complex numbers
$\lambda $ satisfying
$\lambda \in B_{\varepsilon }(\underline{\omega},\bar{\omega})$ and
$| \lambda| >\rho _{\varepsilon }$ we have
$\Delta (\lambda )\neq 0$. So,
for these $\lambda $, system \eqref{7} has a unique solution
\begin{equation*}
C_i(\lambda )=\frac{1}{\Delta (\lambda )}\sum_{k=1}^{4}\Delta
_{ik}(\lambda )f_k,\quad i=1,2,3,4
\end{equation*}
where $\Delta _{ik}(\lambda )$ is an algebraic complement of
 $(i,k) $-th element of the determinant $\Delta (\lambda )$. It
is easy to see that each of the
determinant $\Delta _{ik}(\lambda )$  has the representation
\begin{equation*}
\Delta _{ik}(\lambda )=( \theta _{ik}+r_{ik}(\lambda )) \lambda ^{{m-m_k}}
\end{equation*}
where $\theta _{ik}$ are complex numbers and $r_{ik}\to 0$ as
$| \lambda | \to \infty $ in the angle
$sB_{\varepsilon }(\underline{\omega},\bar{\omega})$. Then we have
\begin{equation*}
C_i(\lambda )=\sum_{k=1}^{4}\lambda ^{-{m_k}}\frac{\theta
_{ik}+r_{ik}(\lambda )}{\theta +r(\lambda )}f_k,\quad i=1,2,3,4.
\end{equation*}
Therefore, the solution of problem \eqref{3}-\eqref{4} has the form
\begin{equation*}
u(x,\lambda )=\sum_{i=1}^{4}\sum_{k=1}^{4}\lambda ^{-{m_k}}\frac{\theta
_{ik}+r_{ik}(\lambda )}{\theta +r(\lambda )}f_ku_i(x,\lambda ).
\end{equation*}
From this it follows that for each integer $l\geq 0$
\begin{equation}\label{8}
\| u^{(l)}\| _{L_q(-1,1)}\leq
C\sum_{k=1}^{4}\Big( | \lambda |
^{l-m_k}| f_k| \sum_{i=1}^{4}\|
u_i(.,\lambda )\| _{L_q(I _i)}\Big) .
\end{equation}
Further, by \eqref{6} we have the inequality
\begin{align*}
\| u_1(.,\lambda )\|_{L_q(-1,0)}^{q}
&=\int_{-1}^{0}e^{q\operatorname{Re}(\omega _1\lambda
)(x+1)}dx\leq \int_{-1}^{0}e^{-q| \lambda |
| \omega _1| \sin \frac{{\varepsilon
}}{2}(x+1)}dx  \\
&=\big( -q| \lambda |
| \omega _1| \sin \frac{{\varepsilon}}{2}\big) ^{-1}
\big( e^{-q| \lambda
| | \omega _1| \sin \frac{{\varepsilon }}{2}}-1\big) \\
&\leq C(\varepsilon )| \lambda |^{-1}
\end{align*}
as $| \lambda | \to \infty $ in the angle
 $ \lambda \in B_{\varepsilon }(\underline{\omega},\bar{\omega})$. In a
similar way we have
\begin{equation*}
\| u_1(\cdot,\lambda )\| _{L_q(I_i)}^{q}
\leq C(\varepsilon )| \lambda | ^{{-1}},\quad i=2,3,4
\end{equation*}
as $| \lambda | \to \infty $ in the angle $
\lambda \in B_{\varepsilon }(\underline{\omega},\bar{\omega})$.
Substituting these inequalities in \eqref{8} we have
\begin{equation*}
\| u^{(l)}\| _{L_q(-1,1)}\leq C(\varepsilon)\sum_{k=1}^{4}|
\lambda | ^{{l-m_k-\frac{{1}}{q}}}| f_k|
\end{equation*}
which, in turn, gives us the needed estimation \eqref{5}. The proof is
complete.
\end{proof}

\section{Fredholm property of problem with multipoint and functional
conditions}

 Let us consider problem \eqref{1}-\eqref{2} and the operator
$\mathcal{L}$ corresponding to this problem. Suppose that
$l\geq \max \{ 2,\max \{ m_1,m_2,m_3,m_4\} +1\} $ and
define a linear operator $\mathcal{L}$ from
$W_q^{l}(-1,0)\dot{+}W_q^{l}(0,1)$ into
$W_q^{l-2}(-1,0)\dot{+}W_q^{l-2}(0,1)+\mathbb{C}^{4}$ by action low
\begin{equation*}
\mathcal{L}u=(L(\lambda )u,L_1u,L_2u,L_3u,L_4u).
\end{equation*}

\begin{theorem} \label{thm2}
Let the following conditions be satisfied:
\begin{itemize}
\item[(1)] $p_1\neq 0$, $p_2\neq 0$;

\item[(2)] the functionals $\mathcal{F}_k,k=1,2,3,4,$ in $W_q^{m_k}(-1,0)\dot{+
}W_q^{m_k}(0,1)$ are continuous;

\item[(3)] $q(x)$ is measurable function on $[-1,0)\cup (0,1]$.
\end{itemize}
Then the linear operator $\mathcal{L}$ is bounded and Fredholm.
\end{theorem}

\begin{proof}
The operator $\mathcal{L}$ can be rewritten in the form
\begin{gather*}
\mathcal{L}_0u=\big(L_0(\lambda)u,L_{10}u,L_{20}u,L_{30}u,L_{40}u\big), \\
\mathcal{L}_1u=\big( q(x)u+\lambda _0u,\mathcal{F}_1u,
\mathcal{F}_2u,\mathcal{F}_3u,\mathcal{F}_4u\big)
\end{gather*}
where $\lambda _0\in B_{\varepsilon }(\underline{\omega},\bar{\omega})$
is some complex number sufficiently large in modulus.
By Theorem \ref{thm1} the operator $\mathcal{L}_0$ is an isomorphism from
$W_q^{l}(-1,0)\dot{+}W_q^{l}(0,1)$ onto
$W_q^{l-2}(-1,0)\dot{+}W_q^{l-2}(0,1)\dot{+}\mathbb{C}^{4}$.
Further, it is easy to see that the
linear operator $\mathcal{L}_1$ acts compactly from
$W_q^{l}(-1,0)\dot{+}W_q^{l}(0,1)$ onto
$W_q^{l-2}(-1,0)\dot{+}W_q^{l-2}(0,1)\dot{+}\mathbb{C}^{4}$.

Consequently, we can apply the theorem of Fredholm operator perturbation
\cite[p. 238]{muk2} to the operator
$\mathcal{L}=\mathcal{L} _0+\mathcal{L}_1$, which follows that
$\mathcal{L}$ is Fredholm. Moreover, it is obvious that the operator
$\mathrm{{\mathcal{L}}}$ is bounded. So, the proof of the theorem is complete.
\end{proof}

\section{Isomorphism and coerciveness of the principal part of the problem}

 Consider problem \eqref{1}-\eqref{2} without internal points, namely,
\begin{gather}\label{9}
L_0(\lambda )u:=p(x)u''(x)-\lambda u(x)=f(x), \\
\label{10}
L_{k0}u :=\alpha _ku^{(m_k)}(-1)+\beta_ku^{(m_k)}(-0)+\eta_ku^{(m_k)}(+0)
+\gamma _ku^{(m_k)}(1)=f_k,
\end{gather}
for $k=1,2,3,4$.
The corresponding operator is
\begin{equation*}
\widetilde{\mathcal{L}}_0u=(L_0(\lambda
)u,L_{10}u,L_{20}u,L_{30}u,L_{40}u).
\end{equation*}

\begin{theorem} \label{thm3}
Let the following conditions be satisfied:
\begin{itemize}
\item[(1)] $\theta \neq 0$;
\item[(2)] $l\geq \max \{ 2,\max \{m_1,m_2,m_3,m_4\} +1\} $.
\end{itemize}
Then for each $\varepsilon >0$ there exist $\rho _{\varepsilon }>0$ such that
for all complex numbers $\lambda $ satisfying
$\lambda \in B_{\varepsilon }(\underline{\omega},\bar{\omega})$,
$| \lambda | >\rho _{\varepsilon }$ the operator
$\widetilde{\mathcal{L}}_0(\lambda )$ from
$W_q^{l}(-1,0)\dot{+}W_q^{l}(0,1)$ onto
$W_q^{l-2}(-1,0)\dot{+}W_q^{l-2}(0,1)\dot{+}\mathbb{C}^{4}$ is an
isomorphism and for these $\lambda $ the following inequality holds for the
solution of \eqref{9}--\eqref{10},
\begin{equation}\label{11}
\begin{aligned}
&\sum_{k=0}^{l}| \lambda | ^{\frac{{l-k}}{2}}\| u\| _{W_{q,k}} \\
&\leq C(\varepsilon )\Big( \| f\|_{W_{q,l-2}}+| \lambda |
^{\frac{{l-2}}{2}}\|f\| _{L_{q,0}}
 + \sum_{\upsilon =1}^{4}| \lambda |
^{( l-m_{\upsilon }-\frac{{1}}{q})/2}| f_{\upsilon }| \Big) .
\end{aligned}
\end{equation}
\end{theorem}

\begin{proof}
 It is obvious that the linear operator $\widetilde{\mathcal{L}}_0(\lambda )$
 is continuous from the space $W_q^{l}(-1,0)\dot{+}W_q^{l}(0,1)$ to
$W_q^{l-2}(-1,0)\dot{+}W_q^{l-2}(0,1)\dot{+}\mathbb{C}^{4}$.
Let $( f(x),f_1,f_2,f_3,f_4) \in W_q^{l-2}(-1,0)\dot{+}W_q^{l-2}(0,1)
\dot{+}\mathbb{C}^{4}$  be any element. We shall seek the solution
$u(x,\lambda )$ of problem \eqref{9}-\eqref{10} in the form of the sum
$u(x,\lambda)=u_1(x,\lambda )+u_2(x,\lambda )$ as follows. By
$f_{\upsilon}(x)$ $(\upsilon =1,2)$ we shall denote the
restriction of $f(x)$ on the interval $\Omega_{\upsilon }$. Let
$\widetilde{f}_{\upsilon }(\cdot)\in W_q^{l-2}(\mathbb{R})$
be an extension of $f_{\upsilon }(\cdot)\in W_q^{l-2}(I_{\upsilon })$ such
that the extension operator $S_{\upsilon }f_{\upsilon }:=\tilde{f}_{\upsilon
}$ from $W_q^{l-2}(I_{\upsilon })$ to
$W_q^{l-2}(\mathbb{R})$ is bounded for $\upsilon =1,2$.
\cite[Lemma 1.7.6]{yak}, where as usual $\mathbb{R}=(-\infty ,\infty )$.
First consider the equations
\begin{equation*}
-p_{\upsilon }(x)u''(x)+\lambda u(x)=\widetilde{f}_{\upsilon
}(x),x\in \mathbb{R}
\end{equation*}
for $\upsilon =1,2$. By applying the \cite[Theorem 3.2.1]{yak} we
see that this equation has a unique solution
 $\tilde{u}_{1\upsilon}=\tilde{u}_{1\upsilon}(\cdot,\lambda )\in W_q^{l}
(\mathbb{R})$ and for $u_{1\upsilon }(x,\lambda )$ (i.e. the restriction of
$\tilde{u}_{1\upsilon }(x,\lambda $ on interval)
$\Omega_{\upsilon }$) the estimate
\begin{equation}\label{12}
\sum_{k=0}^{l}| \lambda |
^{\frac{{l-k}}{2}}\| u_{1\upsilon }\|
_{W_q^k(I_{\Omega_{\upsilon }})}\leq C(\varepsilon )(
\| f\| _{W_q^{l-2}(I_{\upsilon })}+|
\lambda | ^{\frac{{l-2}}{2}}\| f\|
_{L_q(\Omega_{\upsilon })}) ,
\end{equation}
for $\upsilon =1,2$,
is valid for all complex numbers $\lambda $ satisfying
$\lambda \in B_{\varepsilon }(\underline{\omega},\bar{\omega})$.
Consequently, the function
\begin{equation*}
u_1(x,\lambda )=\begin{cases}
u_{11}(x,\lambda ),&\text{for }x\in (-1,0) \\
u_{12}(x,\lambda ),&\text{for }x\in (0,1)
\end{cases}
\end{equation*}
satisfies  equation \eqref{9}. In terms of this solution, we construct the
 boundary-value problem
\begin{gather*}
p(x)u''(x)-\lambda u(x)=0,\quad x\in (-1,0)\cup (0,1), \\
L_{k0}u=f_k-L_{k0}u_1(.,\lambda ),k=1,2,3,4.
\end{gather*}
By Theorem \ref{thm1}, this problem has a unique solution
$u_2=u_2(x,\lambda )$ that belongs to
$W_q^{l}(-1,0)\dot{+}W_q^{l}(0,1)$  for all complex numbers $\lambda $
satisfying $\lambda \in B_{\varepsilon }(\underline{\omega},\bar{\omega})$,
sufficiently large in modulus, and for these $\lambda $ the estimate
\begin{equation}\label{13}
\sum_{k=0}^{l}| \lambda |
^{\frac{l-k}{2}}\| u_2\| _{q,k}\leq C(\varepsilon
)\sum_{\upsilon =1}^{4}|
\lambda | ^{{( l-m_{\upsilon }-\,\frac{{1}}{q}) }\frac{{
1}}{2}}( | f_{\upsilon }| +| L_{\upsilon 0}u_1| )
\end{equation}
holds. By applying the of Theorem \ref{thm1} and taking into account
\cite[Theorem 1.7.7/2]{sha}, we have that for all
$\lambda \in B_{\varepsilon }(\underline{\omega},\bar{\omega})$ and
 $l\geq \max \{ 2,\max \{m_1,m_2,m_3,m_4\} +1\} $ that the following estimates hold.
\begin{equation}\label{14}
\begin{aligned}
| \lambda | ^{( l-m_{\upsilon }-\frac{{1}}{q})/2}
| L_{\upsilon 0}u_1|
&\leq C| \lambda | ^{( l-m_{\upsilon }-\frac{{1}}{q})/2}
\| u_1\| _{C^{m_{\upsilon }}[-1,0]+C^{m_{\upsilon }}[0,1]} \\
&  \leq  C( | \lambda | ^{\frac{{l}}{2}}\|
u_1\| _{q,0}+\| u_1\| _{q,l}) \\
&  \leq  C(\varepsilon )( \| f\|
_{q,l-2}+| \lambda | ^{{\frac{l-2}{2}}}\|
f\| _{q,0}) .
\end{aligned}
\end{equation}
From \eqref{13} and \eqref{14} we have the inequality
\begin{equation}\label{15}
\begin{aligned}
&\sum_{k=0}^{l}| \lambda |^{{\frac{l-k}{2}}}\| u_2\| _{q,k}\\
&\leq C(\varepsilon )\Big( \| f\| _{q,l-2}+|\lambda | ^{{\frac{l-2}{2}}}\|
f\| _{q,0}  
 + \sum_{\upsilon =1}^{4}| \lambda |
^{( l-m_{\upsilon }-\frac{{1}}{2})/2}| f_{\upsilon }| \Big) .
\end{aligned}
\end{equation}
It is easy to see that the function $u(x,\lambda )$ defined as
$u(x,\lambda)=u_1(x,\lambda )+u_2(x,\lambda )$ is the solution of the considered
problem \eqref{9}-\eqref{10}. Taking into account the estimates \eqref{12} and
\eqref{15}, we see that for this solution the needed estimation \eqref{11}
 is valid. Moreover,
from  estimate \eqref{11} it follows the uniqueness of the solution. On the
other hand by Theorem \ref{thm2} the operator
$\widetilde{{\mathcal{L}}}$ is Fredholm from $W_q^{l}(-1,0)\dot{+}W_q^{l}(0,1)$
to $W_q^{l-2}(-1,0)\dot{+}W_q^{l-2}(0,1)\dot{+}\mathbb{C}^{4}$.
Now, isomorphism of this operator follows from the fact that it is a
Fredholm and one-to-one operator. So,the proof of the theorem is
complete.
\end{proof}

\section{Solvability and coerciveness of the main problem with nonlocal
boundary conditions}

Now, we can study the main problem \eqref{1}-\eqref{2}

\begin{theorem} \label{thm4}
 Let the following conditions be satisfied:
\begin{itemize}
\item[(1)] $\theta \neq 0$;

\item[(2)] $l\geq \max \{ 2,\max \{m_1,m_2,m_3,m_4\} +1\} $,

\item[(3)] The functionals $\mathcal{F}_{\upsilon }$ are continuous in
$W_q^{m_{\upsilon }}(-1,0)\dot{+}W_q^{m_{\upsilon }}(0,1)$.
\end{itemize}
Then for each $\varepsilon >0$ there exist $\rho _{\varepsilon }>0$ such that
for all complex numbers $\lambda \in B_{\varepsilon }(\underline{\omega},
\bar{\omega})$ for which $| \lambda | >\rho _{\varepsilon}$ the operator
\begin{equation*}
\widetilde{\mathcal{L}}(\lambda )u:=(L(\lambda
)u,L_1u,L_2u,L_3u,L_4u)
\end{equation*}
is an isomorphism from $W_q^{l}(-1,0)\dot{+}W_q^{l}(0,1)$ onto $
W_q^{l-2}(-1,0)\dot{+}W_q^{l-2}(0,1)\dot{+}\mathbb{C}^{4}$ and for these
$\lambda $ the following coercive estimate holds for the solution of problem
\eqref{1}-\eqref{2}
\begin{equation}\label{16}
\sum_{k=0}^{l}| \lambda |
^{{\frac{l-k}{2}}}\| u\| _{q,k}
\leq C(\varepsilon )\Big( \| f\| _{q,l-2}+|
\lambda | ^{{\frac{l-2}{2}}}\|f\| _{q,0}
  + \sum_{\upsilon =1}^{4}| \lambda |
^{( l-m_{\upsilon }-\frac{{1}}{q})/2} | f_{\upsilon }| \Big)
\end{equation}
where $C(\varepsilon )$ is a constant which depends only on
$\varepsilon $.
\end{theorem}

\begin{proof}
Let $( f(x),f_1,f_2,f_3,f_4) $ be any element of
$W_q^{l-2}(-1,0)\dot{+}W_q^{l-2}(0,1)\dot{+}\mathbb{C}^{4}$. Assume that
there exists a solution $u=u(x,\lambda )$ of problem \eqref{1}-\eqref{2}
corresponding to this element. Then this solution satisfies the equalities
\begin{gather}\label{17}
L_0(\lambda )u=L(\lambda )u-q(x)u, \\
\label{18}
L_{k0}u=L_ku-S_ku-\mathcal{F}_ku,~k=1,2,3,4.
\end{gather}
By applying  Theorem \ref{thm3} to the problem \eqref{17}-\eqref{18} we have that for this
solution the following a priory estimate hold
\begin{equation}\label{19}
\begin{aligned}
\sum_{k=0}^{l}| \lambda |^{{\frac{l-k}{2}}}\| u\| _{q,k}
&\leq C(\varepsilon )\Big( \| L(\lambda)u-q(x)u\| _{q,l-2}
 +| \lambda | ^{{\frac{l-2}{2}} }\| L(\lambda )u-q(x)u\| _{q,0}  \\
&\quad +\sum_{\upsilon =1}^{4}| \lambda |
^{( l-m_{\upsilon }-\frac{{1}}{q})/2}| L_{\upsilon
}u-S_ku-\mathcal{F}_{\upsilon }u| \Big)  \\
&\quad + C(\varepsilon )\Big( \| f\|
_{q,l-2}+| \lambda | ^{{\frac{l-2}{2}}}\|f\|_{q,0}+\| q(x)u\| _{q,l-2}   \\
&\quad  +| \lambda | ^{{\frac{l-2}{2}}}\|
q(x)u\| _{q,0}+\sum_{\upsilon =1}^{4}|
\lambda| ^{( l-m_{\upsilon }-\frac{{1}}{q})/2} |f_{\upsilon }|   \\
&\quad + \sum_{\upsilon =1}^{4}| \lambda |
^{( l-m_{\upsilon }-\frac{{1}}{q})/2}
 (| S_ku| +|\mathcal{F}_{\upsilon }u| )\Big)
\end{aligned}
\end{equation}
Let $\delta $ be any real number satisfying
\begin{equation*}
0<\delta <\min \big\{  \frac{1}{2},1+x_{ki},|
x_{ki}| ,1-x_{ki} :
k=1,2,3,4,\;i=1,2,\dots,n_k\big\} .
\end{equation*}
By applying the same approach as in  \cite[sec. 2.8.3]{muk1} it is easy to
construct a function $\psi _{\delta }(x)\in C_0^{\infty }[-1,1]$ such that
\begin{gather*}
\psi _{\delta }(x)=1\quad \text{for }
x\in [ -1+\delta ,-\delta ]\cup[ \delta ,1-\delta ]\,,\\
\psi _{\delta }(x)=0\quad \text{for }x\in [ -1,-1+\frac{\delta }{2}]\cup
[ -\frac{\delta }{2},\frac{\delta }{2}]\cup [ 1-\frac{\delta }{2},1]
\end{gather*}
and $0\leq \psi _{\delta }(x)\leq 1$ for all $x\in [ -1,1]$. It is
obvious that
\begin{equation}\label{20}
| S_ku| \leq C\| ( \psi _{\delta}u) ^{(m_k)}\| _{C[-1,1]}.
\end{equation}
By \cite[Theorem 3.10.4]{ras}, for $u\in W_q^{l}(-1,0)\dot{+}
W_q^{l}(0,1)$ the following estimate holds,
\begin{equation}\label{21}
| \lambda | ^{( l-m_{\upsilon }-\frac{{1}}{q})/2}\| u^{(m_{\upsilon })}\|
_{C[-1,1]}
\leq C( \| u\| _{q,l}+| \lambda | ^{\frac{{l}}{2}}\| u\| _{q,0})\,.
\end{equation}
By  Theorem \ref{thm4}, from \eqref{20} and \eqref{21} it follows that for all
$\lambda \in B_{\varepsilon }(\underline{\omega},\bar{\omega})$
sufficiently large in modulus the following estimate holds,
\begin{equation} \label{22}
\begin{aligned}
| \lambda | ^{( l-m_{\upsilon }-\frac{{1}}{q})/2}| S_{\upsilon }u|
&\leq C| \lambda | ^{( l-m_{\upsilon }-\frac{{1}}{q})/2}\|
  (\psi _{\delta }u)^{(m_{\upsilon })}\| _{C[-1,1]} \\
&\leq  C\big( \| \psi _{\delta }u\| _{q,l}+|
\lambda | ^{\frac{{l}}{2}}\| \psi _{\delta }u\|_{q,0}\big) \\
&  \leq C(\varepsilon )\big( \| L_0(\lambda )(\psi_{\delta
}u)\| _{q,l-2}+| \lambda | ^{{\frac{l-2}{2}}
}\| L_0(\lambda )(\psi _{\delta }u)\| _{q,0}\big) \\
&  \leq C(\varepsilon )\Big( \| L_0(\lambda)u\| _{q,l-2}+| \lambda |^{{\frac{l-2}{2}}}
 \|L_0(\lambda )u\| _{q,0}   \\
&\quad +  \| q(x)u\| _{q,l-2}+| \lambda| ^{{\frac{l-2}{2}}}\| q(x)u\|_{q,0}
 +\sum_{k=0}^{l-1}| \lambda | ^{{\frac{l-1-k}{2}}}\| u\| _{q,k}\Big)\\
&\leq C(\varepsilon )\Big( \| f\|_{q,l-2}+| \lambda | ^{{\frac{l-2}{2}}}\|
f\| _{q,0}\\
&\quad  +  \| q(x)u\| _{q,l-2}+| \lambda| ^{{\frac{l-2}{2}}}
 \| q(x)u\| _{q,0}+\sum_{k=0}^{l-1}| \lambda | ^{{\frac{l-1-k}{2}}}
 \| u\| _{q,k}\Big)
\end{aligned}
\end{equation}
By \cite[Theorem 1.3.3]{ayd}  there is a positive constant $C$ such that for
all $u$ in the set $W_q^{l}(-1,0)\dot{+}W_q^{l}(0,1)$ and for each $k=0,1,\dots ,l-1$
the following inequality is valid
\begin{equation}\label{23}
\| u\| _{q,k}\leq C\| u\| _{q,k+1}^{\frac{k}{k+1}}\| u\| _{q,0}^{\frac{1}{k+1}}.
\end{equation}
Applying the well-known Young inequality
\begin{equation*}
ab\leq \frac{1}{p}( \alpha a) ^{p}+\frac{1}{q}( \frac{b}{\alpha }) ^{q}
\end{equation*}
where $a>0$, $b>0$, $\alpha >0$, $1<p,q<\infty $, $\frac{1}{p}+\frac{1}{q}=1$
to the right-hand of \eqref{22} for
\begin{equation*}
a=\| u\| _{q,k+1}^{\frac{k}{k+1}},\quad
b=\|u\| _{q,0}^{\frac{1}{k+1}}, \quad
p=\frac{k+1}{k},
\end{equation*}
we have
\begin{equation*}
\| u\| _{q,k}\leq C\Big( \frac{k}{k+1}\alpha ^{{\frac{k+1}{
k}}}\| u\| _{q,k+1}+\frac{1}{k+1}\alpha ^{-(k+1)}\|
u\| _{q,0}\Big)
\end{equation*}
for $k=0,1,\dots ,l-1$.
We denote
\begin{gather*}
A(\alpha )=\max \big\{ C\frac{k}{k+1}\alpha ^{{\frac{k+1}{k}}}:
 k=0,1,\dots ,l-1\big\},\\
B(\alpha )=\max \big\{C\frac{1}{k+1}\alpha ^{-(k+1)}:
k=0,1,\dots ,l-1\big\} .
\end{gather*}
Then from  inequality \eqref{21}, we have
\begin{equation}\label{24}
\begin{aligned}
| \lambda | ^{( l-m_{\upsilon }-\frac{{1}}{q})/2}| S_{\upsilon }u| 
&\leq C(\varepsilon )( \| f\| _{q,l-2}+|\lambda | ^{{\frac{l-2}{2}}}
 \| f\|_{q,0}) \\
&\quad + C(\varepsilon )\sum_{k=0}^{l-1}| \lambda | ^{{\frac{l-1-k
}{2}}}( A(\alpha )\| u\| _{q,k+1}+B(\alpha )\|u\| _{q,0})\\
&  \leq \big( C(\varepsilon )A(\alpha )+D(\varepsilon ,\alpha
)| \lambda | ^{-1/2}\big)
\sum_{k=0}^{l}| \lambda |
^{{\frac{l-k}{2}}}\| u\| _{q,k}
\end{aligned}
\end{equation}
where $D(\varepsilon ,\alpha )$ is a constant which depends only on
$\varepsilon $ and $\alpha $. In view of \cite[Theorem 1.7.7/2]{yak}, for any
$\zeta >0$ we obtain
\begin{equation*}
\| u\| _{q,k}\leq \zeta \| u\|_{q,k+1}+C(\zeta )\| u\| _{q,0}.
\end{equation*}
On the other hand, from \cite[Lemma 1.8]{ayd}
 and \cite[Theorem 8.19]{ras} we have
\begin{equation}\label{25}
\begin{aligned}
| \mathcal{F}_ku|
&\leq \sum_{j=0}^{m_k}\Big(
| \int_{\Omega _1}\mathcal{K}_{k1j}(t)u^{(j)}(t)dt|
 +| \int_{\Omega _1}\mathcal{K}_{k2j}(t)u^{(j)}(t)dt
| \Big)   \\
&\leq \sup_k \Big( \sum_{j=0}^{m_k}\int_{\Omega
_1}| \mathcal{K}_{k1j}(t)u^{(j)}(t)|
dt+\sum_{j=0}^{m_k}\int_{\Omega _1}| \mathcal{K}
_{k2j}(t)u^{(j)}(t)| dt\Big)   \\
&\leq \sup_k \Big( \sum_{j=0}^{m_k}\int_{\Omega
_1}| \mathcal{K}_{k1j}(t)u(t)|
dt+\sum_{j=0}^{m_k}\int_{\Omega _1}| \mathcal{K}
_{k2j}(t)u(t)| dt\Big)   \\
&\leq C_1\| u\| _{q,k}+C_2\| u\|_{q,k}   \\
&\leq C\| u\| _{q,k}.
\end{aligned}
\end{equation}
From  \eqref{23} and  \eqref{24} we have
\begin{equation}\label{26}
\begin{aligned}
&\| q(x)u\| _{q,l-2}+| \lambda | ^{\frac{l-2}{2}}
\| q(x)u\| _{q,0}+\sum_{\upsilon=1}^{4}|
\lambda | ^{( l-m_{\upsilon }-\frac{{1}}{q})/2}
(| S_ku| +| \mathcal{F}_{\upsilon }u| )  \\
&\leq C(\varepsilon )(\| f\| _{q,l-2}+|
\lambda | ^{{\frac{l-2}{2}}}\| f\|_{q,0})
+\zeta (\| u\| _{q,l}+| \lambda | ^{\frac{l-2}{2}}\| u\| _{q,0})   \\
&\quad +C(\zeta )| \lambda | ^{\frac{l-2}{2}}\|
u\| _{q,0}+( C(\varepsilon )A(\alpha )+D(\varepsilon ,\alpha
)| \lambda | ^{-1/2})
\sum_{k=0}^{l}| \lambda | ^{{\frac{l-k}{2}}}\|u\| _{q,k}   \\
&\quad +C\sum_{\upsilon =1}^{4}| \lambda |
^{(l-m_{\upsilon }-\frac{{1}}{q})/2} \| u\|_{q,k}   \\
&\leq C(\varepsilon )(\| f\| _{q,l-2}+|\lambda| ^{{\frac{l-2}{2}}}
\| f\| _{q,0})   \\
&\quad +\big( C(\varepsilon )A(\alpha )+D(\varepsilon ,\alpha )|
\lambda | ^{{-}\frac{{1}}{2q}}\big) \sum_{k=0}^{l}|
\lambda | ^{{\frac{l-k}{2}}}\| u\| _{q,k}
\end{aligned}
\end{equation}
Substituting \eqref{25} into \eqref{19} we obtain
\begin{align*}
\sum_{k=0}^{l}| \lambda |^{{\frac{l-k}{2}}}\| u\| _{q,k}
&\leq C(\varepsilon )\Big( \| f\| _{q,l-2}+|
\lambda | ^{{\frac{l-2}{2}}}\| f\|_{q,0}
+\sum_{\upsilon =1}^{4}| \lambda | ^{( l-m_{\upsilon }-\frac{{1}}{q})/2}
| f_{\upsilon }| \Big) \\
&\quad +( C(\varepsilon )A(\alpha )+D(\varepsilon ,\alpha
)| \lambda | ^{{-}\frac{{1}}{2q}})
\sum_{k=0}^{l}| \lambda | ^{{\frac{l-k}{2}}}\| u\| _{q,k}.
\end{align*}
For a fixed $\varepsilon >0$ we can choose $\alpha >0$ so small, and
$| \lambda | $ so large that
\begin{equation*}
C(\varepsilon )A(\alpha )+D(\varepsilon ,\alpha )| \lambda| ^{-1/2q} <1.
\end{equation*}
Thus, for $\lambda \in B_{\varepsilon }(\underline{\omega},\bar{\omega})$
sufficiently large in modulus we obtain a priori estimate (5.1). From this
estimate it follows the uniqueness property of the solution of problem
\eqref{1}-\eqref{2}, i.e. the operator $\widetilde{\mathcal{L}}(\lambda )$
is one-to-one operator. Moreover, by  Theorem \ref{thm2} the operator
$\widetilde{\mathcal{L}}(\lambda )$ from
$W_q^{l}(-1,0)\dot{+} W_q^{l}(0,1)$ to
$W_q^{l-2}(-1,0)\dot{+}W_q^{l-2}(0,1)\dot{+}\mathbb{C}^{4}$ is Fredholm.
 Consequently, the existence of a solution results in
its uniqueness. So, the proof of the theorem is complete.
\end{proof}

\begin{thebibliography}{99}

\bibitem{agr} M. S. Agranovich;
\emph{Spectral properties of diffraction
problems, in: The Generalized method of eigenoscillations in the
theory of diffraction}, N. N. Voitovich, B. Z. Kazenelenbaum, A. N.
Sivov, Nauka, Moscow, 1977 (in Russian: translated into English
Wiley-VCH, Berlin), 1999.

\bibitem{ali} Z. S. Aliyev;
\emph{Basis properties of a fourth order differential
operator with spectral parameter in the boundary condition}, Open
Mathematics, 8(2) (2010), 378-388.

\bibitem{all}  B. P. Allahverdiev,   E. Bairamov,   E., Ugurlu;
\emph{Eigenparameter dependent Sturm-Liouville problems in
boundary conditions with transmission conditions}, Journal of
Mathematical Analysis and Applications, 401, No. 1 (2013), 388-396.

\bibitem{all1} B. P. Allahverdiev,  E. Ugurlu;
\emph{On dilation, scattering and spectral theory for
two-interval singular differential operators}, Bulletin Mathematique
 de la Societe des Sciences Mathematiques de Roumanie, 58(106), No. 4 (2015), 3
83-392.

\bibitem{ayd} K. Aydemir;
\emph{Boundary value problems with eigenvalue depending
boundary and transmission conditions}, Boundary Value Problems,
2014/1/131 (2014).

\bibitem{ayd1}  K. Aydemir,  O. Sh. Mukhtarov;
\emph{Second-order differential operators with interior singularity},
Advances in Difference Equations  (2015)
2015:26 DOI 10.1186/s13662-015-0360-7.

\bibitem{ayd2}  K. Aydemir,  O. Sh. Mukhtarov;
\emph{Spectrum and Green's Function of a Many-Interval Sturm-Liouville
Problem}, Z. Naturforsch. 70(5)a(2015), 301-308

\bibitem{bor} M. Borsuk;
\emph{Transmission problems for elliptic
second-order equations in non-smooth domains}, Springer Basel AG
2010.

\bibitem{gar} M. Garcia-Huidobra, P. Gupta Chaitan, R. Manasevich;
\emph{Some multipoint boundary value problems of Neumann-Dirichlet type
involving a multipoint p-Laplace like operator}, J. Math. Anal. Appl.
333 (2007), 247-264.

\bibitem{gen} F. Genoud, B. P. Rynne;
\emph{Some recent results on the spectrum
of multi-point eigenvalue problems for the p-Laplacian}, Commun Appl.
Anl., 15 (2011), 413-434.

\bibitem{huy} H. Huy Pham,  Palencia E. Sanchez;
\emph{Phenomenes des transports a travers descouches minces de conductivite eleve},
J. Math. Anal. Appl.,  47 (1974), 284-309.

\bibitem{ima} N. S. Imanbaev, M. A. Sadybekov;
\emph{Characteristic determinant of the spectral problem for the ordinary
differential operator with the boundary load, International conference on
analysis and applied mathematics (ICAAM 2014)} - AIP Conference Proceedings,
1611 (2014), 261-265.

\bibitem{ism} Z. Ismailov, P. Ipek;
\emph{Selfadjoint singular differential operators
of first order and their spectrum}, Electronic Journal of
Differential Equations, 21 (2016), 1-9

\bibitem{kad} M. Kadakal, O. Sh. Mukhtarov;
\emph{Discontinuous Sturm-Liouville
problems containing eigenparameter in the boundary conditions}, Acta
Mathematica Sinica, 22(5) (2006), 1519-1528.

\bibitem{kan} M. Kandemir;
\emph{Irregular boundary value problems for elliptic
differential-operator equations with discontinuous coefficients and
transmission conditions}, Kuwait Journal of Science and Engineering,
39(1A)(2012), 71-97.

\bibitem{kan1} M. Kandemir, Ya.Yakubov;
\emph{Regular boundary value problems with a discontinuous coefficient,
 functional-multipoint conditions, and a
linear spectral parameter}, Israel Journal of Mathematics, 180(2010),
255-270.

\bibitem{kat} T. Kato;
\emph{Perturbation theory for linear operators},
Sipringer-Verlag, New York Inc., 1966.

\bibitem{ker} N. B. Kerimov, Kh. R.  Mamedov;
\emph{On a boundary value problem with Spectral parameter in the
boundary conditions}, (Russian) Sibirsk. Mat. Zh.,  40(2) (1999),
325-335.
\bibitem{lik} A. V. Likov, Yu. A. Mikhailov;
\emph{The theory of heat and mass transfer}, Qosenergaizdat, (Russian), 1963.

\bibitem{mam}  K. R. Mamedov, F. A. Cetinkaya;
\emph{Eigenparameter dependent inverse
boundary value problem for a class of Sturm-Liouville operator},
Boundary Value Problems, 2014: 194 (2014).

\bibitem{mar} S. A. Margareth, J. E. Rivera, S. Mauricio, O. V. Villagran;
\emph{Transmission problem in thermoelasticity}, Hindawi Pub.
Cor. B. V. P., ID 190548(2011), 33.

\bibitem{muk2} O. Sh. Mukhtarov, K. Aydemir;
\emph{Eigenfunction expansion for Sturm-Liouville problems with
transmission conditions at one interior point},
 Acta Mathematica Scienta, 35(3) (2015), 639-649.

\bibitem{osh5} O. Sh. Mukhtarov, K. Aydemir;
\emph{New type Sturm-Liouville problems in associated Hilbert spaces}, 
J. Funct. Spaces Appl., vol. 2014 (2014), Article ID 606815, 7 pages,
 doi:10.1155/2014/606815. 

\bibitem{muk1} O. Sh. Mukhtarov, S. Yakubov;
\emph{Problems for ordinary differential equations with transmission conditions},
Applicable Analysis, 81(2002), 1033-1064.

\bibitem{ras} M. L. Rasulov;
\emph{Application of contour integral method},
(in Russian) Navka, Moskow, 1997.

\bibitem{sad} M. A. Sadybekov, B. Kh. Turmetov, B. T. Torebek;
\emph{Solvability of nonlocal boundary-value problems for the Laplace equation
in the ball}, Electronic Journal of Differential Equations, 157 (2014), 1-14.

\bibitem{sha} V. B. Shakhmurov;
\emph{Linear and nonlinear abstract elliptic
equations with VMO coefficients and applications}, Fixed point theory
and applications, 6(2013), 1-21.


\bibitem{ugu} E. U\v{g}urlu, E. Bairamov;
\emph{Dissipative operators with impulsive conditions}, Math Chem,
51 1670-1680(2013).

\bibitem{voi} N. N. Voitovich, B. Z. Katsenelbaum, A. Sivov;
\emph{Generalized method of eigen-vibration in the theory of diffraction},
(Russian) Nauka, Moscow, 1997.

\bibitem{yak} S. Yakubov, Ya. Yakubov;
\emph{Differential-operator equation ordinary and partial differential equation},
Chapman and Hall/CRC, Boca Raton, London State New York Washington, State D. C.
1999.


\end{thebibliography}

\end{document}
