\documentclass[reqno]{amsart}
\usepackage{hyperref}

\AtBeginDocument{{\noindent\small
\emph{Electronic Journal of Differential Equations},
Vol. 2016 (2016), No. 161, pp. 1--12.\newline
ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu}
\thanks{\copyright 2016 Texas State University.}
\vspace{8mm}}

\begin{document}
\title[\hfilneg EJDE-2016/161\hfil Solvability of  boundary-value problems]
{Solvability of boundary-value problems for Poisson equations
 with Hadamard type \\ boundary operator}

\author[B. Kh. Turmetov, M. D. Koshanova, K. I. Usmanov \hfil EJDE-2016/161\hfilneg]
{Batirkhan Turmetov, Maira Koshanova, Kayrat Usmanov}

\address{Batirkhan Kh. Turmetov \newline
Department of Mathematics,
Akhmet Yasawi International Kazakh-Turkish University,
161200 Turkistan, Kazakhistan. \newline
Institute of Mathematics and Mathematical Modeling,
Ministry of Education and Science Republic of Kazakhstan,
050010 Almaty, Kazakhistan}
\email{batirkhan.turmetov@ayu.edu.kz}

\address{Maira D. Koshanova \newline
Department of Mathematics,
Akhmet Yasawi International Kazakh-Turkish University,
161200 Turkistan, Kazakhistan}
\email{maira\_koshanova@mail.ru}

\address{Kayrat I. Usmanov \newline
Department of Mathematics,
Akhmet Yasawi International Kazakh-Turkish University,
161200 Turkistan, Kazakhistan}
\email{y\_kairat@mail.ru}

\thanks{Submitted December 21, 2015. Published June 27, 2016.}
\subjclass[2010]{35J05, 35J25, 34B10, 26A33}
\keywords{Fractional derivative; Hadamard operator; Poisson equation;
\hfill\break\indent  Neumann problem}

\begin{abstract}
 In this article we study properties of some integro-differential operators
 of fractional order. As an application of the properties of these operators
 for Poisson equation we examine questions on solvability of a fractional
 analogue of Neumann problem and analogues of periodic boundary-value
 problems for circular domains. The exact conditions for solvability of
 these problems are found.
\end{abstract}

\maketitle
\numberwithin{equation}{section}
\newtheorem{theorem}{Theorem}[section]
\newtheorem{lemma}[theorem]{Lemma}
\newtheorem{remark}[theorem]{Remark}
\newtheorem{problem}[theorem]{Problem}
\allowdisplaybreaks


\section{Introduction}\label{intr}

Let $Q$ be a bounded domain from $\mathbb{R}^n$  with a smooth boundary $S$.
 It is known that classical problems for the Poisson equation.
\begin{equation}\label{e1-1}
\Delta u(x) = f(x),x \in Q,
\end{equation}
are Dirichlet and Neumann problems. Let $\nu $ be a normal vector to $S$,
and ${D_\nu } = \frac{d}{{d\nu }}$ be an operator of differentiation
along the normal, $D_\nu ^0 = I$  be a unit vector.
Then Dirichlet and Neumann boundary conditions can be given in the  form
\begin{equation}\label{e1-2}
D_\nu ^\alpha u(x) = {g_\alpha }(x),x \in S,
\end{equation}
where $\alpha  = 0$ or $\alpha  = 1$, $D_\nu ^0 u(x) = u(x)$.
It is known that the Dirichlet problem is unconditionally solvable,
and for solvability of the Neumann condition the following condition
is necessary \cite{Bitsazde}:	
\begin{equation}\label{e1-3}
\int_Q {f(x)dx = } \int_S {{g_1}(x)dx}.
\end{equation}

In this article, we introduce fractional analogues of the boundary operators
$D_\nu ^\alpha $, and for the equation \eqref{e1-1} we study the
boundary-value problem with the boundary condition \eqref{e1-2} for all
values of the parameter $\alpha  \in (0,\infty )$.
Moreover, we investigate solvability of some analogues of periodic
boundary-value problems for circular domains.

The structure of this paper is as follows.
In Introduction we provide an overview of some papers published on the subject.
Further, we give concepts of Hadamard type integral-differential operators
of fractional order. In the second section we study properties of
integral-differential operators of fractional order in the class of
smooth functions. Properties of these operators in H\"{o}lder class are studied.
Propositions about reversibility of the operators are proved.
In the third section we present some auxiliary statements related to the
properties of solutions of the Dirichlet problem for the Poisson equation.
In the fourth and fifth sections we consider applications of these
integral - differential operators of fractional order to examine questions on
 solvability of some boundary-value problems with boundary operators of
fractional order. In the fourth section we study questions about solvability
of a fractional analogue of the Neumann problem. The problem is solved by
reducing it to an equivalent Dirichlet problem with the additional condition
at the point $x = 0$. In the fifth section we also study analogues of
periodic problems for circular domains. The problem is reduced to
two auxiliary problems: Dirichlet problem and an analogue of Neumann problem.

Note that the local and nonlocal boundary-value problems with boundary
operators of fractional order for the second order elliptic equations were
studied in \cite{berdyshev-smj,kirane-AMS,karachik-sam,kirane-smj1,kirane-smj2,
kirane-MMAS,Krasnoschok,Muratbekova,torebek-bvp,turmetov-de1,
turmetov-sam,turmetov-de2,umarov-dm,umarov-fcaa}
and for higher-order equations in
\cite{Bekaeva,berdyshev-ams,berdyshev-MM,turmetov-ejde}.
As the boundary operators in
\cite{kirane-AMS,karachik-sam,kirane-smj1,kirane-smj2,kirane-MMAS,Krasnoschok,
torebek-bvp,turmetov-de1,turmetov-sam,turmetov-de2,umarov-dm,umarov-fcaa}
operators with Riemann-Liouville and Caputo type derivatives, and in
\cite{Bekaeva,berdyshev-smj,Muratbekova} the Hadamard - Marchaud type operators
were considered.	
We also note that applications of boundary-value problems for elliptic equations
with boundary operators of fractional order have been considered in
 \cite{AkhVI1,AkhVI2,Veliev}.
	Now let us turn to the definitions of integration and differentiation
operator of fractional order.	

Let $\Omega  = \{ {x \in {\mathbb{R}^n}:| x | < 1}\}$ be a unit ball,
$n \ge 2$, $\partial \Omega  = \{ {x \in {\mathbb{R}^n}:| x | = 1} \}$
 - unit sphere. Suppose further that, $u(x)$ is a smooth function in the
 domain $\Omega $, $r = |x|,\theta  = x/r$, $\delta  = r\frac{d}{{dr}}$-
Dirac operator, where
\[
r\frac{d}{{dr}} = \sum_{j = 1}^n {{x_j}\frac{\partial }{{\partial {x_j}}}} ,
\quad \alpha  > 0.
\]
Further, let $0 < \alpha  < \infty $. The expression
$$
{J^\alpha }[u](x) = \frac{1}{{\Gamma (\alpha )}}
\int_0^r {{{\big( {\ln \frac{r}{s}} \big)}^{\alpha  - 1}}u(s\theta )\frac{{ds}}{s}}
 $$
is called integration operator of the $\alpha $ order in the Hadamard sense
(see e.g. \cite{kilbas-book}). Furthermore, we assume that ${J^0}[u](x) = u(x)$.

Note that, if $u(0) \ne 0$, then in the class of continuous functions the
operator ${J^\alpha }$ is not defined, since the integral
$\int_0^1 {{{( {\ln \frac{1}{s}} )}^{\alpha  - 1}}{s^{ - 1}}ds} $
 diverges. Therefore, as the differentiation operator we consider the
 Hadamard - Caputo type operator. Namely, differentiation operator
of the $\alpha  > 0$ order is the expression:
$$
D^\alpha[u](x) = \frac{1}{\Gamma (\ell  - \alpha )}
\int_0^r {{{( {\ln \frac{r}{s}} )}^{\ell  - 1 - \alpha }}{{( {s\frac{d}{{ds}}} )}
^\ell }u(s\theta )\frac{{ds}}{s}} ,\ell  - 1 < \alpha  \le \ell ,\;\ell\geq 1.
$$

\section{Properties of $J^\alpha $ and $D ^\alpha $ operators}\label{s2}

In this section we study properties of ${J^\alpha }$ and ${D^\alpha }$ operators.
Further, by the symbol $C$ we denote the constant whose value can be different.

\begin{lemma}\label{lem2-1}
Let $\alpha  > 0,\quad 0 < \lambda  < 1$ and
$u(x) \in {C^{\lambda  + p}}( {\overline \Omega  } ), p \geq 0$.
If the condition $u(0) = 0$ holds, then
${J^\alpha }[ u ](x) \in {C^{\lambda  + p}}( {\overline \Omega  } )$
and ${J^\alpha }[ u ](0) = 0$.
\end{lemma}

\begin{proof}
If $u( 0 ) = 0$, then
$$
| {{J^\alpha }[ u ](x)} |
\le \frac{1}{{\Gamma ( \alpha  )}}{\int_0^1 {( {\ln \frac{1}{s}} )}
^{\alpha  - 1}}\frac{{|u( {s\theta } )|}}{s}
\le \frac{C}{{\Gamma ( \alpha  )}}
{\int_0^1 {( {\ln \frac{1}{s}} )} ^{\alpha  - 1}}{s^{\lambda  - 1}}ds.
$$
Since the last integral converges, the function
 ${J^\alpha }[ u ](x)$ is defined in the domain $\overline \Omega  $.
 Let $x^{(1)},\;x^{(2)}$ be arbitrary points of the domain $\overline \Omega$.
Denote $h(x) = J^\alpha[u](x)$.
Then
\begin{align*}
| {h({x^{( 1 )}}) - h({x^{( 2 )}})} |
& \le \frac{1}{{\Gamma (\alpha )}}
\int_0^1 {( {\ln \frac{1}{s}} ){s^{\mu  - 1}}|u(s{x^{( 1 )}})
- u(s{x^{( 2 )}})|ds}  \\
&\le \frac{{C|{x^{( 1 )}} - {x^{( 2 )}}{|^\lambda }}}{{\Gamma (\alpha )}}
\int_0^1 {{{( {\ln \frac{1}{s}} )}^{\alpha  - 1}}{s^{\mu  + \lambda  - 1}}ds}  \\
&\le C|{x^{( 1 )}} - {x^{( 2 )}}{|^\lambda },
\end{align*}
i.e. $h(x) \in {C^\lambda }( {\overline \Omega  } )$.
	Further, if $\beta  = ( {{\beta _1},{\beta _2}, \ldots {\beta _n}} )$
is a multi-index and
$\partial _x^\beta  = \frac{{{\partial ^{| \beta  |}}}}{{\partial x_1^{{\beta _1}}
\cdots \partial x_n^{{\beta _n}}}}$,
then for all  $\beta $ with length $|\beta | \le p$ and
${x^{( 1 )}},\;{x^{( 2 )}} \in \overline \Omega  $ we have
\begin{align*}
| {\partial _x^\beta h({x^{( 1 )}}) - \partial _x^\beta h({x^{( 2 )}})} |
&\le \frac{1}{{\Gamma (\alpha )}}
 \int_0^1 {{{( {\ln \frac{1}{s}} )}^{\alpha  - 1}}{s^{\mu  + | \beta  | - 1}}|
 {\partial _x^\beta u({y_1}) - \partial _x^\beta h({y_2})} |ds} \\
&\le C|{x^{(1)}} - {x^{(2)}}{|^\lambda },
\end{align*}
where  $y = sx = ( {s{x_1},s{x_2}, \ldots s{x_n}} )$, and, consequently
 ${J^\alpha }[ u ](x) \in {C^{\lambda  + p}}( {\overline \Omega  } )$.
Lastly,
$$\lim _{x \to 0} \;|J_0^\alpha [u](x)|
\leq C \lim _{x \to 0} |x|^{\lambda} = 0.
$$
Then $J_0^\alpha [u](0)=0$.
\end{proof}

Similarly we can prove the following statement.

\begin{lemma}\label{lem2-2}
Let $\ell  - 1 < \alpha  \le \ell$, $\ell  = 1,2, \ldots $,
$0 < \lambda  < 1$ and $u(x) \in {C^{\lambda  + p}}( {\bar \Omega } ),p \ge \ell $.
Then ${D^\alpha }[u](x) \in {C^{\lambda  + p - \ell }}( {\bar \Omega } )$
and the equality ${D^\alpha }[u](0) = 0$. holds.
\end{lemma}

\begin{lemma}\label{lem2-3}
Let  $\ell  - 1 < \alpha  \le \ell$, $\ell  = 1,2, \ldots $,
$0 < \lambda  < 1$ and $u(x) \in {C^{\lambda  + p}}( {\bar \Omega } )$,
$p \ge \ell $, $p = 1,2,\dots $. Then for any $x \in \bar \Omega$:
\begin{equation}\label{e2-1}
{J^\alpha }[ {{D^\alpha }[u]} ]( x ) = u(x) - u(0),
\end{equation}
and if $u(0) = 0$,  then we obtain
\begin{equation}\label{e2-2}
{D^\alpha }[ {{J^\alpha }[u]} ](x) = u(x).
\end{equation}
\end{lemma}

\begin{proof}
If $u(x) \in {C^{\lambda  + p}}( {\bar \Omega } ),\;\;p \ge \ell $,
then by  Lemma \ref{lem2-2} we obtain
${D^\alpha }[u](x) \in {C^{\lambda  + p - \ell }}( {\bar \Omega } )$
and ${D^\alpha }[u](0) = 0$. Let us prove  equality \eqref{e2-1} for
the case $\alpha  = \ell $ - integer.
Since ${D^\alpha }[u](0) = 0$, then in the class of these functions
the operator ${J^\alpha }$ is defined, and in this case:
$$
{J^\ell }[ {{D^\ell }[u]} ](x) = \frac{1}{{( {\ell  - 1} )!}}
\int_0^r {{s^{ - 1}}{{( {\ln \frac{r}{s}} )}^{\ell  - 1}}{{( {s\frac{d}{{ds}}} )}
^\ell }u(s\theta )ds}.
$$
Integrating by parts the last integral  $\ell  - 1$ times,  we obtain
$$
{J^\ell }[ {{D^\ell }[u]} ](x)
= \int_0^r \frac{d}{{ds}}[ {u(s\theta )} ]ds =
 {u(s\theta )} |_{s = 0}^{s = r} = u(x) - u(0).
$$
Let now $\ell  - 1 < \alpha  < \ell $, $\ell  = 1,2, \ldots $.
Then
$$
{J^\alpha }[ {{D^\alpha }[u]} ](x)
= {J^\alpha }[ {{J^{\ell  - \alpha }}[{\delta ^\ell }[u]]} ](x).
$$
Further, since ${J^\alpha } \cdot {J^{\ell  - \alpha }} = {J^\ell }$
(see e.g. \cite{kilbas-book}, page 114),  it follows that
$$
{J^\alpha }[ {{D^\alpha }[u]} ](x)
= {J^\ell }[ {{\delta ^\ell }[u]} ](x) = u(x) - u(0).
$$
The equality \eqref{e2-1} is proved. Let us turn to the proof of
the equality \eqref{e2-2}.
If $\alpha  = \ell $, then
\begin{align*}
{D^\ell }[ {{J^\ell }[u]} ](x)
& = {\delta ^\ell }\Big\{ {\frac{1}{{( {\ell  - 1} )!}}
\int_0^r {{s^{ - 1}}{{( {\ln \frac{r}{s}} )}^{\ell  - 1}}u(s\theta )ds} } \Big\} \\
&= {\delta ^{\ell  - 1}}\Big\{ {\frac{1}{{( {\ell  - 2} )!}}
\int_0^r {{s^{ - 1}}{{( {\ln \frac{r}{s}} )}^{\ell  - 2}}u(s\theta )ds} } \Big\} \\
&= r\frac{d}{{dr}}\{ {\int_0^r {{s^{ - 1}}u(s\theta )ds} } \} = u(x).
\end{align*}
Further, in the case  $\ell  - 1 < \alpha  < \ell$, $\ell  = 1,2, \ldots$
by the definitions of ${D^\alpha }$  and  ${J^\alpha }$ operators, we obtain
\begin{align*}
D^\alpha [ {{J^\alpha }[u]} ](x)
&= \frac{1}{{\Gamma ( {\ell  - \alpha } )}}
 \int_0^r {\frac{1}{{\ell  - \alpha }}{{( {\ln \frac{r}{s}} )}^{\ell
 - \alpha  - 1}}{\delta ^\ell }[ {{J^\alpha }[ u ]} ](sx)\frac{{ds}}{s}} \\
& = \frac{1}{{\Gamma ( {\ell  - \alpha } )}}r\frac{d}{{dr}}
 \int_0^r {\frac{1}{{\ell  - \alpha }}{{( {\ln \frac{r}{s}} )}^{\ell
 - \alpha }}\frac{d}{{ds}}[ {{\delta ^{\ell  - 1}}[ {{J^\alpha }[ u ]} ]} ]ds } \\
 &= \frac{1}{{\Gamma ( {\ell  - \alpha } )}}r\frac{d}{{dr}}
 \int_0^r {{( {\ln \frac{r}{s}} )}^{\ell  - 1 - \alpha }}
 {\delta ^{\ell  - 1}}[ {{J^\alpha }[ u ]} ]\frac{{ds}}{s}.
\end{align*}
Performing this operation again $( {\ell  - 1} )$ times, we have
$$
{D^\alpha }[ {{J^\alpha }[u]} ](x)
= {\delta ^\ell }[ {{J^{\ell  - \alpha }}{J^\alpha }[ u ]} ](x)
= {\delta ^\ell }[ {{J^\ell }[ u ]} ](x) = u(x).
$$
\end{proof}

\begin{lemma}\label{lem2-4}
Let $\ell  - 1 < \alpha  \le \ell$, $\ell  = 1,2, \ldots $,
$0 < \lambda  < 1$, $f(x)$ be a smooth function in the domain $\bar \Omega $
and $\Delta u(x) = f(x)$, $x \in \Omega $. Then
\begin{equation}\label{e2-3}
\Delta {D^\alpha }[u](x) = |x{|^{ - 2}}{D^\alpha }[|x{|^2}f](x),\;x \in \Omega.
\end{equation}
\end{lemma}

\begin{proof}
We represent the function ${D^\alpha }[ u ](x)$ in the following form:
$$
{D^\alpha }[ u ](x) = \frac{1}{{\Gamma ( {\ell  - \alpha } )}}
\int_0^1 {{{( {\ln \frac{1}{\xi }} )}^{\ell  - 1 - \alpha }}
{{( {\xi \frac{d}{{d\xi }}} )}^\ell }[ {u(\xi x)} ]\frac{{d\xi }}{\xi }} .
$$
Further, since
$$
\Delta {( {\xi \frac{d}{{d\xi }}} )^\ell }[u](\xi x)
= {\xi ^2}{( {\xi \frac{d}{{d\xi }} + 2} )^\ell }f(\xi x)
 = {\xi ^2}{( {\delta  + 2} )^\ell }[f](\xi x),
$$
it follows that
\begin{align*}
\Delta {D^\alpha }[u](x)
&= \frac{1}{{\Gamma (\ell  - \alpha )}}\int_0^1 {{{( {\ln \frac{1}{\xi }} )}^{\ell
- 1 - \alpha }}{\xi ^2}{{( {\delta  + 2} )}^\ell }[ f ](\xi x)\frac{{d\xi }}{\xi } }\\
&={r^{ - 2}}{J^{\ell  - \alpha }}[ {{\delta ^\ell }[ {{r^2}f} ]} ](x) \\
&= {r^{ - 2}}{D^\alpha }[ {{r^2}f} ](x).
\end{align*}
\end{proof}

\begin{remark} \label{rmk2.1} \rm
It is easy to prove that for the function
$F(x) = |x{|^{ - 2}}{D^\alpha }[|x{|^2}f](x)$ the following representation holds:
\begin{equation}\label{e2-4}
F(x) = \big( {r\frac{d}{{dr}} + 2} \big){f_{\ell  - \alpha }}(x),
\end{equation}
where
\begin{equation}\label{e2-5}
{f_{\ell  - \alpha }}(x) = {r^{ - 2}}{J^{\ell  - \alpha }}[ {{r^2}{\delta ^{\ell  - 1}}[f]} ](x).
\end{equation}
\end{remark}

\section{A property of the Dirichlet problem solution}

In the domain $\Omega $ we consider the  Dirichlet problem
\begin{equation}\label{e3-1}
\begin{gathered}
\Delta v( x ) = F( x ), \quad x \in \Omega, \\
v( x ) = g( x ),\quad  x \in \partial \Omega.
\end{gathered}
\end{equation}
It is known \cite{Gilbarg} that if $0 < \lambda  < 1$,
$F( x ) \in {C^{\lambda  + p}}(\overline \Omega  )$,
$g( x ) \in {C^{\lambda  + p + 2}}(\partial \Omega ),\,p \geq 1$,
then a solution of the problem exists, is unique, belongs to the class
 ${C^{\lambda  + 2}}$ and  can be represented in the form:
\begin{equation}\label{e3-2}
v( x ) =  - \frac{1}{{{\omega _n}}}\int_\Omega  {G( {x,y} )F( y )dy}
 + \frac{1}{{{\omega _n}}}\int_{\partial \Omega }
 {\frac{{1 - {{| x |}^2}}}{{{{| {x - y} |}^n}}}g( y )d{s_y}},
\end{equation}
where ${\omega _n}$ is a square of the unit sphere, $G( {x,y} )$ is the
 Green function of the problem \eqref{e3-1}. Moreover,  $G( {x,y} )$
is represented in the form \cite{Bitsazde}:
$$
G( {x,y} ) = \begin{cases}
\frac{1}{{n - 2}}[ {{{| {x - y} |}^{2 - n}} - {{| {x|y|
- \frac{y}{{| y |}}} |}^{2 - n}}} ], & n \ge 3\\
\ln \frac{1}{{| {x - y} |}}, &n = 2.
\end{cases}
$$

Let $\rho  = | y |$.

\begin{lemma}\label{lem3-1}
Let  $F(y), g(y)$ be smooth enough functions and $F( y )$ be represented
in the form $F( y ) = ( {\rho \frac{\partial }{{\partial \rho }}
+ 2} ){f_1}( y ),v(x)$ be a solution of the problem \eqref{e3-1}.
Then the condition $v(0) = 0$ holds if and only if
\begin{equation}\label{e3-3}
\int_\Omega  {{f_1}( y )dy}  = \int_{\partial \Omega } {g( y )d{s_y}}.
\end{equation}
\end{lemma}

\begin{proof}
Since $F( y )$ and $g( y )$ are smooth enough functions, then solution of the
problem \eqref{e3-1} exists and can be represented as \eqref{e3-2}.
Then in the case $n \ge 3$ we have
\begin{equation}\label{e3-4}
v(0) =  - \frac{1}{{{\omega _n}}}\int_\Omega  {\frac{1}{{n - 2}}[ {{{| y |}^{2 - n}}
 - 1} ]F( y )dy}  + \frac{1}{{{\omega _n}}}\int_{\partial \Omega } {g( y )dy} .
\end{equation}
We consider the following two integrals:
\begin{gather*}
{I_1}(\rho ,\xi ) = \int_0^1 {{\rho ^{n - 1}}[ {{\rho ^{2 - n}} - 1} ]
\rho \frac{\partial }{{\partial \rho }}{f_1}( {\rho ,\xi } )d\rho } , \\
{I_2}(\rho ,\xi ) = 2\int_0^1 {{\rho ^{n - 1}}[ {{\rho ^{2 - n}} - 1} ]{f_1}
( {\rho ,\xi } )d\rho } .
\end{gather*}
	Integrating ${I_1}$ by parts, we obtain
\begin{align*}
{I_1}(\rho ,\xi )
&= \int_0^1 {[ {{\rho ^2} - {\rho ^n}} ]\frac{\partial }{{\partial \rho }}{f_1}
( {\rho ,\xi } )d\rho } \\
& =  - \int_0^1 {[ {2\rho  - n{\rho ^{n - 1}}} ]{f_1}( {\rho ,\xi } )d\rho } \\
&= \int_0^1 {{\rho ^{n - 1}}[ {n - 2{\rho ^{2 - n}}} ]{f_1}( {\rho ,\xi } )d\rho }.
\end{align*}
Since  $F( y )$ has the form
$( {\rho \frac{\partial }{{\partial \rho }} + 2} ){f_1}( y )$,
 moving to spherical coordinates for the first integral in the right-hand
side of \eqref{e3-4}, we have
\begin{align*}
& - \frac{1}{{{\omega _n}}}\int_\Omega  {\frac{1}{{n - 2}}
 [ {{{| y |}^{2 - n}} - 1} ]( {\rho \frac{\partial }
 {{\partial \rho }} + 2} ){f_1}( y )dy}  \\
& =  - \frac{1}{{( {n - 2} ){\omega _n}}}\int_{| \xi  | = 1}
 {\int_0^1 {{\rho ^{n - 1}}[ {{\rho ^{2 - n}} - 1} ]
 ( {\rho \frac{\partial }{{\partial \rho }} + 2} ){f_1}
 ( {\rho ,\xi } )d\rho d\xi  } } \\
&=  - \frac{1}{{( {n - 2} ){\omega _n}}}\int_{| \xi  | = 1} 
 {[ {{I_1}(\rho ,\xi ) + {I_2}(\rho ,\xi )} ]d\xi } \\
&=  - \frac{1}{{{\omega _n}}}\int_{| \xi  | = 1} 
\int_0^1 {\rho ^{n - 1}}{f_1}( {\rho ,\xi } )d\rho d\xi  
=    - \frac{1}{{{\omega _n}}}\int_{ \Omega  } {{f_1}( y )} dy.
\end{align*}
Consequently, if $v(0) = 0$, then the equality \eqref{e3-3} holds. 
Hence, necessity of the condition \eqref{e3-3} is proved. Sufficiency 
is proved in reverse order.
\end{proof}

\section{Neumann type problem}

In this section we consider a fractional analogue of the Neumann problem with 
the boundary operator ${D^\alpha }$.

\begin{problem}\label{probl4.1} \rm
Let $0 < \alpha $. Find a function 
$u(x) \in {C^2}(\Omega ) \cap C(\overline \Omega  )$ such that 
${D^\alpha }[ u ](x) \in C(\overline \Omega  )$, and satisfying the
equation
\begin{equation}\label{e4-1}
  \Delta u(x) = f(x), x \in \Omega ,
\end{equation}
and the boundary value condition
\begin{equation}\label{e4-2}
 {D^\alpha }[ u ](x) = g(x), x \in \partial \Omega .
\end{equation}
\end{problem}

Since ${J^0} = I$ ,  when $\alpha  = 1$ we have
$$
 {{D^1}u(x)} \big|_{\partial \Omega }
 = { {{J^0}[\delta [u]](x)} \big|_{\partial \Omega }} 
= { {r\frac{{du(x)}}{{dr}}} \big|_{\partial \Omega }} 
= { {\frac{{\partial u(x)}}{{\partial \nu }}} \big|_{\partial \Omega }}.
$$
Therefore, when $\alpha  = 1$ the problem  \eqref{e4-1} - \eqref{e4-2}
 coincides with the classical Neumann problem.
	
\begin{theorem} \label{thm4.2}
Let $\ell - 1 < \alpha  \le \ell$, $\ell  = 1,2,\dots $,  $0 < \lambda  < 1$, 
$f(x) \in {C^{\lambda  + 2\ell  - 1}}(\overline \Omega  )$, 
$g(x) \in {C^{\lambda  + \ell  + 1}}(\partial \Omega )$. 
Then for solvability of the problem \ref{probl4.1} it is necessary and sufficient 
the condition
\begin{equation}\label{e4-3}
\int_\Omega  {{f_{\ell  - \alpha }}(y)dy}  = \int_{\partial \Omega } g(y)dy.
\end{equation}
where the function ${f_{\ell  - \alpha }}(x)$ is defined by the equality 
\eqref{e2-5}.

	If a solution of the problem exists, then it is unique up to a constant term, 
belongs to the class ${C^{\lambda  + \ell  + 1}}(\overline \Omega  )$ 
and can be represented in the  form
\begin{equation}\label{e4-4}
  u(x) = C + {J^\alpha }[v](x),
\end{equation}
where $v(x)$ is a solution of problem \eqref{e3-1} with the function
 $F(x) =r^{-2}D^\alpha [{r^2}f](x)$ and satisfies the condition $v(0) = 0$.
\end{theorem}

\begin{proof}
Let $u(x)$ be a solution of  problem \ref{probl4.1}. Apply the operator ${D^\alpha }$ 
to the function $u(x)$, and denote $v(x) = {D^\alpha }[u](x)$. 
Find conditions, which the function $v(x)$ satisfies.  It is obvious that 
${ {v(x)} |_{\partial \Omega }} = { {{D^\alpha }[u](x)} |_{\partial \Omega }} 
= g(x)$. Applying the operator $\Delta $ to the equality
 $v(x) = {D^\alpha }[u](x)$, due to \eqref{e2-3},  we obtain
$$
\Delta v(x) = {r^{ - 2}}{D^\alpha }[{r^2}f](x). 
$$
Therefore, if $u(x)$ is a solution of the problem \ref{probl4.1}, then 
 $v(x) = {D^\alpha }[u](x)$ will be a solution of \eqref{e3-1} 
with the function $F(x)=r^{-2}D^\alpha[{r^2}f]$. Moreover, according to 
Lemma \ref{lem2-2}, the function $v(x)$ satisfies the condition $v(0) = 0$.
By \eqref{e2-4}, the function $F( x )$ can be represented in the form 
$$
F(x) = \Big( r\frac{d}{dr} + 2\Big) f_{\ell  - \alpha }(x),
$$ 
where ${f_{\ell  - \alpha }}(x)$ is defined by the equality \eqref{e2-5}. 
Then, by Lemma \ref{lem3-1} for the equality $v(0) = 0$ the following condition is necessary:
$$
\int_\Omega  {{f_{\ell  - \alpha }}(y)dy}  = \int_{\partial \Omega } {g(y)d{S_y}}.
$$
Therefore, necessity of the condition \eqref{e4-3} is proved. 
Applying the operator ${J^\alpha }$ to the equality 
$v(x) = {D^\alpha }[u](x)$, because of \eqref{e2-1}, we obtain
$$
u(x) - u(0) = {J^\alpha }[v](x).
$$
Hence, if a solution of the problem \ref{probl4.1}
 exists, then it can be represented 
as \eqref{e4-4}.
We show that the condition \eqref{e4-3} is sufficient for the existence 
of any solution of the problem \ref{probl4.1}.

Indeed, let $v(x)$ be a solution of the problem \eqref{e3-1} with 
$F( x ) = {r^{ - 2}}{D^\alpha }[{r^2}f](x)$. 
If  $f(x) \in {C^{\lambda  + 2\ell  - 1}}(\overline \Omega  )$, 
then $F( x ) \in {C^{\lambda  + \ell  - 1}}(\overline \Omega  )$, 
and since $g(x) \in {C^{\lambda  + \ell  + 1}}(\partial \Omega )$, 
 a solution of  \eqref{e3-1} exists, is unique and belongs to the class 
${C^{\lambda  + p + 1}}( {\bar \Omega } )$ (see e.g. \cite{Gilbarg}).
 We represent the function $F( x ) = |x{|^{ - 2}}{D^\alpha }[|x{|^2}f](x)$ as 
$F( x ) = ( {r\frac{d}{{dr}} + 2} ){f_{\ell  - \alpha }}(x)$.
If for the function ${f_{\ell  - \alpha }}(x)$ the condition \eqref{e4-3} holds, 
then corresponding solution of the problem \eqref{e3-1} satisfies the 
condition $v(0) = 0$. Then we should to consider the function
 $u(x) = C + {J^\alpha }[v](x)$, which satisfies all conditions of
 problem \ref{probl4.1}. By Lemma \ref{lem2-1} this function belongs to the class 
${C^{\lambda  + p + 1}}( {\bar \Omega } )$.
	Further, using \eqref{e2-2}, we obtain
\[
{ {{D^\alpha }[u](x)} |_{\partial \Omega }} 
= { {{D^\alpha }[C] + {D^\alpha }[{J^\alpha }[v]](x)} |_{\partial \Omega }} 
= {\ {v(x)} |_{\partial \Omega }} = g(x).
\]
Moreover,
\begin{align*}
\Delta u(x) 
&= \Delta \Big[ {\frac{1}{{\Gamma ( {1 - \alpha } )}}\int_0^r {{{
( {\ln \frac{r}{s}} )}^{1 - \alpha }}v(s\theta )\frac{{ds}}{s}} } \Big]\\
&= \Delta \Big[ {\frac{1}{{\Gamma ( {1 - \alpha } )}}
\int_0^1 {{{( {\ln \frac{1}{\xi }} )}^{1 - \alpha }}v(\xi x)
\frac{{d\xi }}{\xi }} } \Big] \\
&= \frac{1}{{\Gamma ( {1 - \alpha } )}}
\int_0^1 {{{( {\ln \frac{1}{\xi }} )}^{1 - \alpha }}{\xi ^2}F(\xi x)
\frac{{d\xi }}{\xi }}  \\
&= \frac{1}{{\Gamma ( {1 - \alpha } )}}\int_0^1 {{{( {\ln \frac{1}{\xi }} 
 )}^{1 - \alpha }}{\xi ^2}|\xi x{|^{ - 2}}{D^\alpha }[|\xi x{|^2}f(\xi x)]
 \frac{{d\xi }}{\xi }}  \\
&=\frac{{|x{|^{ - 2}}}}{{\Gamma ( {1 - \alpha } )}}
 \int_0^r {{{( {\ln \frac{r}{s}} )}^{1 - \alpha }}{D^\alpha }[{s^2}f(s\theta )]
\frac{{ds}}{s}}  \\
&= {r^{ - 2}}{J^\alpha }[ {{D^\alpha }[{r^2}f]} ](x) 
= {r^{ - 2}} \cdot {r^2}f(x) = f(x).
\end{align*}
Thus, the function $u(x) = C + {J^\alpha }[v](x)$ satisfies all conditions 
of  problem \ref{probl4.1}.
\end{proof}

\begin{remark} \label{rmk4.1} \rm
If  $ \alpha  = 1$, then ${f_1}(x) = {r^{ - 2}}{J^0}[{r^2}f](x) = f(x)$ and 
\eqref{e4-3} coincides with the condition of solvability of the Neumann 
problem \eqref{e1-3}.
\end{remark}

\section{Boundary-value problems with periodic conditions}

In this section we study some analogues of periodic problems in $\Omega $.
Let $x = ({x_1},\tilde x) \in \Omega , \tilde x = ({x_2},\dots ,{x_n})$
 For any $x = ({x_1},\tilde x) \in \Omega $ we put ``opposite''
 point ${x^*} = (-x_1,a\tilde x) \in \Omega $, where 
$a = ({a_2},{a_3},\dots ,{a_n})$ and ${a_j},j = 2,\dots ,n$ 
take one of the values $ \pm 1$.
Denote 
$$
\partial {\Omega _ + } = \{ {x \in \partial \Omega :{x_1} \ge 0} \},\quad
\partial {\Omega _ - } = \{ {x \in \partial \Omega :{x_1} \le 0} \},\quad
I = \{ {x \in \partial \Omega :{x_1} = 0} \}.
$$
Let $0 < \alpha  \le 1$. Consider in $\Omega $ the following problem:

\begin{problem} \label{problem5.1}\rm
Find a function $u(x) \in {C^2}(\Omega ) \cap C(\bar \Omega )$, such that 
${D^\alpha }[u](x) \in C(\bar \Omega )$ and
\begin{gather}\label{e5-1}
   \Delta u(x) = f(x), \quad x \in \Omega, \\
\label{e5-2}
   u(x) - {( - 1)^k}u({x^*}) = {g_0}(x), \quad x \in \partial {\Omega _ + }, \\
\label{e5-3}
   {D^\alpha }[u](x) + {( - 1)^k}{D^\alpha }[u](x^{*}) 
= {g_1}(x), \quad x \in \partial {\Omega _ + },
\end{gather}
where $k = 1,2$.
\end{problem}

The problem \eqref{e5-1}--\eqref{e5-3} in the case $\alpha  = 1$ have been 
studied in \cite{sadybekov-de},\cite{sadybekov-emj} and in the case 
$0 < \alpha  < 1$ for the Riemann - Liouville and Caputo operators 
in \cite{sadybekov-ejd}.
	
If $x = (0,\tilde x) \in I$, then $x^{*} = (0,\alpha \tilde x) \in I$, 
therefore, a necessary condition for  existence of a solution from the 
class $u(x) \in {C^2}(\Omega ) \cap C(\bar \Omega )$,
 ${D^\alpha }[u](x) \in C(\bar \Omega )$ is the fulfillment of the 
 conditions
\begin{gather}\label{e5-4}
  {g_0}(0,\tilde x) =  - {( - 1)^k}{g_0}(0,a\tilde x),\\
  \frac{{\partial {g_0}(0,\tilde x)}}{{\partial {x_j}}} 
=  - {( - 1)^k}\frac{{\partial {g_0}(0,a\tilde x)}}{{\partial {x_j}}},\quad
j = 1,\dots ,n,(0,\tilde x) \in I, \\
\label{e5-5}
  {g_1}(0,\tilde x) = {( - 1)^k}{g_1}(0,a\tilde x),\quad (0,\tilde x) \in I.
   \end{gather}

\begin{theorem} \label{thm5.2}
Let $0 < \lambda  < 1$, $f(x) \in {C^{\lambda  + 1}}(\bar \Omega )$,
${g_0}(x) \in {C^{\lambda  + 2}}(\partial {\Omega _ + })$, 
${g_1}(x) \in {C^{\lambda  + 2}}(\partial {\Omega _ + })$ and the matching 
conditions \eqref{e5-4},\eqref{e5-5} hold. Then
\begin{itemize}
\item[(1)] if $k = 1$, then a solution of the problem 
\eqref{e5-1}--\eqref{e5-3} exists, is unique and belongs to the 
class ${C^{\lambda  + 2}}(\bar \Omega )$;

\item[(2)] if $k = 2$ , then for solvability of the problem
 \eqref{e5-1} - \eqref{e5-3} the following condition is necessary and sufficient:
\begin{equation}\label{e5-6}
  \int_\Omega  {{f_{1 - \alpha }}(y)dy}  = \int_{\partial {\Omega _ + }} {{g_1}(y)d{S_y}}.
  \end{equation}
\end{itemize}
  If a solution exists, then it is unique up to a constant term, and belongs 
to the class ${C^{\lambda  + 2}}(\bar \Omega )$.
\end{theorem}

\begin{proof}
First we prove uniqueness. Let  $u(x)$ be a solution of the homogenous 
problem \eqref{e5-1} - \eqref{e5-3}. Putting the function  $u(x)$ 
into the boundary value conditions of the problem \eqref{e5-1}--\eqref{e5-3}, 
we have
\begin{gather}\label{e5-7}
  u(x) = {( - 1)^k}u(x^{*}),\quad x \in \partial {\Omega _ + }, \\
\label{e5-8}
  {D^\alpha }[u](x) =  - {( - 1)^k}{D^\alpha }[u](x^{*}),\quad
x \in \partial {\Omega _ + }.
\end{gather}

If $x \in \partial {\Omega _ - }$, then $x^{*} \in \partial {\Omega _ + }$. 
Then the condition \eqref{e5-7} implies $u(x^{*}) = (-1)^k u(x)$, 
$ x \in \partial {\Omega _ - }$, and \eqref{e5-8} yields 
${D^\alpha }[u](x^{*}) =  - {( - 1)^k}{D^\alpha }[u](x),x \in \partial {\Omega _ - }$ Consequently,  the equalities  \eqref{e5-7} and \eqref{e5-8} hold for all points $x \in \partial \Omega $, i.e.
$$
u(x) = {( - 1)^k}u(x^{*}), {D^\alpha }[u](x)
 =  - {( - 1)^k}{D^\alpha }[u](x^{*}), \quad x \in \partial \Omega. 
$$
Since ${D^\alpha }[u](x) \in C(\bar \Omega )$, then from the equality 
$u(x) = {( - 1)^k}u(x^{*}), x \in \partial \Omega $ it follows:
 ${D^\alpha }u(x) = {( - 1)^k}{D^\alpha }u(x^{*}), x \in \partial \Omega $.
Consequently, ${D^\alpha }u(x) = 0, x \in \partial \Omega $, i.e. solution 
of the homogeneous problem \eqref{e5-1} - \eqref{e5-3} is also solution 
of the homogeneous Problem \ref{probl4.1}. 
Then by Theorem \ref{thm4.2}:  $u(x) \equiv C,x \in \bar \Omega $.
Hence, putting $u(x) \equiv C$ into \eqref{e5-7}, when $k = 1$ we have 
$u(x) \equiv 0$. Therefore, solution of the problem \eqref{e5-1}--\eqref{e5-3} 
when  $k = 1$ is unique, and when $k = 2$ it is unique up to a constant term. 
Uniqueness is proved.
Now let us turn to study existence of a solution. 
Consider the  auxiliary functions
$$
v(x) = \frac{1}{2}( {u(x) + u(x^{*})} ), \quad
w(x) = \frac{1}{2}( {u(x) - u(x^{*})} ).
$$
It is obvious that $u(x) = v(x) + w(x)$. Moreover, 
$v(x) = v(x^{*})$, $w(x) =  - w(x^{*})$.

We find problems, which these functions satisfy. Let $k = 1$.
Applying the operator  $\Delta $ to the functions $v(x)$ and $w(x)$, we have
\begin{gather*}
\Delta v(x) = \frac{1}{2}[ {\Delta u(x) + \Delta u(x^{*})} ] 
= \frac{1}{2}[ {f(x) + f(x^{*})} ] \equiv {f^ + }(x),\quad  x \in \Omega , \\
\Delta w(x) = \frac{1}{2}[ {\Delta u(x) - \Delta u(x^{*})} ] 
= \frac{1}{2}[ {f(x) - f(x^{*})} ] \equiv {f^ - }(x), x \in \Omega .
\end{gather*}
Further, from the boundary value conditions \eqref{e5-2} and \eqref{e5-3} 
we obtain
\begin{gather*}
{ {v(x)} \big|_{\partial {\Omega _ + }}} 
= { {\frac{1}{2}[ {u(x) + u(x^{*})} ]} \big|_{\partial {\Omega _ + }}} 
= \frac{{{g_0}(x)}}{2},\\
{ {{D^\alpha }w(x)} \big|_{\partial \Omega }} 
= { {\frac{1}{2}[ {{D^\alpha }u(x) - {D^\alpha }u(x^{*})} ]} 
\big|_{\partial {\Omega _ + }}} = \frac{{{g_1}(x)}}{2}.
\end{gather*}
If $x \in \partial {\Omega _ - }$, then $x^{*} \in \partial {\Omega _ + }$, 
so the following equalities hold:
\begin{gather*}
{ {v(x)} \big|_{\partial {\Omega _ - }}} = { {\frac{1}{2}[ {u(x^{*}) + u(x)} ]} 
\big|_{\partial {\Omega _ + }}} = \frac{{{g_0}(x^{*})}}{2},\\
{ {{D^\alpha }w(x)} \big|_{\partial \Omega }} 
= { { - \frac{1}{2}[ {{D^\alpha }u(x^{*}) - {D^\alpha }u(x)} ]}
 \big|_{\partial {\Omega _ + }}} =  - \frac{{{g_1}(x^{*})}}{2}.
\end{gather*}
We introduce the functions
\[
2{\tilde g_0}(x) =  \begin{cases}
{g_0}(x),x \in \partial {\Omega _ + }\\
{g_0}(x^{*}), x \in \partial {\Omega _ - }
\end{cases} ,\quad
2{\tilde g_1}(x) = \begin{cases}
{g_1}(x),x \in \partial {\Omega _ + }\\
 - {g_1}(x^{*}), x \in \partial {\Omega _ - }
\end{cases}. 
\]
Therefore, functions $v(x)$ and  $w(x)$ are solutions of the two problems:
\begin{gather}\label{e5-9}
  \Delta v(x) = {f^ + }(x), \quad x \in \Omega ;\quad 
{{v(x)} |_{\partial \Omega }} = {\tilde g_0}(x), \\
\label{e5-10}
  \Delta w(x) = {f^ - }(x),x \in \Omega ;\quad 
{{{D^\alpha }w(x)} |_{\partial \Omega }} = {\tilde g_1}(x).
\end{gather}

If for the functions $f(x),{g_0}(x)$ and ${g_1}(x)$ the conditions of the theorem 
hold, then ${f^ \pm }(x) \in {C^{\lambda  + 1}}(\bar \Omega )$,
${\tilde g_0}(x) \in {C^{\lambda  + 2}}(\partial \Omega )$,
${\tilde g_1}(x) \in {C^{\lambda  + 2}}(\partial \Omega )$. 
Then a solution of the Dirichlet problem \eqref{e5-9} exists, is unique, 
belongs to the class ${C^{\lambda  + 2}}(\bar \Omega )$. 
By Theorem \ref{thm4.2}, for solvability of the problem \eqref{e5-10} 
it is necessary and sufficient the following condition:
\begin{equation}\label{e5-11}
  \int_\Omega  {f_{1 - \alpha }^ - (y)dy}  
= \int_{\partial \Omega } {{{\tilde g}_1}(y)dy} ,
\end{equation}
where $f_{1 - \alpha }^ - (y) = {r^{ - 2}}{J^{1 - \alpha }}[{r^2}{f^ - }](x)$.

Since
\begin{gather*}
\int_\Omega  {f_{1 - \alpha }^ - (y)dy}  
= \frac{1}{2}\int_\Omega  {f_{1 - \alpha }^ - (y)dy}  
- \frac{1}{2}\int_\Omega  {f_{1 - \alpha }^ - (y*)dy = 0,} \\
\int_{\partial \Omega } {{{\tilde g}_1}(y)d{S_y} = } 
\frac{1}{2}\int_{\partial \Omega } {{g_1}(y)d{S_y} 
- \frac{1}{2}\int_{\partial \Omega } {{g_1}(y*)d{S_y} = 0} } ,
\end{gather*}
it follows that the condition for solvability of \eqref{e5-11}  always holds, 
and therefore, in this case 
${f^ - }(x) \in {C^{\lambda  + 1}}(\bar \Omega )$,
${\tilde g_1}(x) \in {C^{\lambda  + 2}}(\partial \Omega )$ a solution of
 problem \eqref{e5-10} exists and belongs to the class
 ${C^{\lambda  + 1}}(\bar \Omega )$. Note that a solution of the problem 
\eqref{e5-10} is unique up to a constant term $C$. 
Since the function $w(x)$ should have the property $w(x) =  - w(x^{*})$, 
we obtain $C \equiv 0$. Therefore, the existence of a solution of 
 problem \eqref{e5-1}--\eqref{e5-3} for the case $k = 1$ is proved.

Let $k = 2$. In this case for auxiliary functions $v(x)$ and $w(x)$
 we obtain the following problems:
\begin{gather}\label{e5-12}
  \Delta w(x) = {f^ - }(x),\quad x \in \Omega ;\quad 
{ {w(x)} |_{\partial \Omega }} = {\tilde g_0}(x), \\
\label{e5-13}
  \Delta v(x) = {f^ + }(x),\quad x \in \Omega ;\quad 
{ {{D^\alpha }v(x)} |_{\partial \Omega }} = {\tilde g_1}(x).
\end{gather}
Here 
\[
2{\tilde g_0}(x) =  \begin{cases}
{g_0}(x), & x \in \partial {\Omega _ + }\\
 - {g_0}(x^{*}), & x \in \partial {\Omega _ - }
\end{cases} ,\quad
2{\tilde g_1}(x) =  \begin{cases}
{g_1}(x), & x \in \partial {\Omega _ + }\\
{g_1}(x^{*}),& x \in \partial {\Omega _ - }.
\end{cases} 
\]
When the conditions of the theorem hold, a solution of 
 problem \eqref{e5-12} exists, it is unique and belongs to the 
class ${C^{\lambda  + 1}}(\bar \Omega )$. 
And for solvability of the problem \eqref{e5-13} it is necessary 
and sufficient the condition
\begin{equation}\label{e5-14}
  \int_\Omega  {f_{1 - \alpha }^ + (y)dy}  
= \int_{\partial \Omega } {{{\tilde g}_1}(y)d{S_y}}.
\end{equation}
Since
\begin{gather*}
\int_\Omega  {f_{1 - \alpha }^ + (y)dy}  
= \frac{1}{2}\int_\Omega  {{f_{1 - \alpha }}(y)dy}  
+ \frac{1}{2}\int_\Omega  {{f_{1 - \alpha }}(y*)dy 
= } \int_\Omega  {{f_{1 - \alpha }}(y)dy} ,\\
\int_{\partial \Omega } {{{\tilde g}_1}(y)d{S_y}}  
= \frac{1}{2}\int_{\partial {\Omega _ + }} {{g_1}(y)d{S_y}}  
+ \frac{1}{2}\int_{\partial {\Omega _ - }} {{g_1}(y*)d{S_y}} 
 = \int_{\partial {\Omega _ + }} {{g_1}(y)d{S_y}} ,
\end{gather*}
it follows that \eqref{e5-14} can be rewritten as \eqref{e5-6}.
 Under this condition, a solution of  problem \eqref{e5-13} exists, 
is unique up to a constant term, and belongs to the class
${C^{\lambda  + 1}}(\bar \Omega )$.
\end{proof}

\begin{remark} \label{rmk5.1} \rm
When $\alpha  = 1$  the propositions  in Theorems  \ref{thm4.2} and 
\ref{thm5.2} coincide  with the results in \cite{sadybekov-de,sadybekov-emj}.
\end{remark}

\subsection*{Acknowledgements}
This research is supported by a grant from the Ministry
of Science and Education of the Republic of Kazakhstan
(Grant No.0819 /GF4).

\begin{thebibliography}{99}

\bibitem{AkhVI1} T. M. Akhmedov, E. Veliev, M. Ivakhnychenko; 
\emph{Fractional operators approach in electromagnetic wave reflection problems},
Journal of electromagnetic waves and applications. 21, No.13 (2007), 1787--1802.

\bibitem{AkhVI2} T. M. Akhmedov, E. Veliev, M. Ivakhnychenko; 
\emph{Description of the boundaries in scattering problems using fractional
operators}, Radio Physics and Electronics. 14,(2009), 133--141.

\bibitem{Bekaeva} A. E. Bekaeva, V. V. Karachik, B. Kh. Turmetov; 
\emph{On solvability of some boundary value problems for polyharmonic 
equation with Hadamard--Marchaud boundary operator}, Russian Math., 58, No.
7 (2014), 11-–24.

\bibitem{berdyshev-smj} A. S. Berdyshev, B. Kh. Turmetov, B. J. Kadirkulov;
\emph{Some properties and applications of the integrodifferential
operators of Hadamard-Marchaud type in the class of harmonic
functions}, Siberian Mathematical Journal. 53, No.4 (2012), 600--610.

\bibitem{berdyshev-ams} A. S. Berdyshev, A. Cabada, B. Kh. Turmetov;
\emph{On solvability of a boundary value problem for a nonhomogeneous 
biharmonic equation with a boundary operator of a fractional order}, 
Acta Mathematica Scientia. 34B, No.6 (2014),  1695--1706.

\bibitem{berdyshev-MM} A. S. Berdyshev, A. Cabada, B. Kh. Turmetov;
\emph{On solvability of some boundary value problem for polyharmonic equation
 with boundary operator of a fractional order}, 
Applied Mathematical Modelling.  39 (2015), 4548--4569.

\bibitem{Bitsazde} A. V. Bitsazde;
\emph{Equations of Mathematical Physics}, Mir Publishers,
Moscow, 1980. (Russian).

\bibitem{Gilbarg} D. Gilbarg, N. S. Trudinger;
\emph{Elliptic Partial Differential Equations of Second Order},
 Grundlehren, Vol. 224, Springer-Verlag, Berlin, 1983.

\bibitem{kirane-AMS} B. J. Kadirkulov, M. Kirane;
\emph{On solvability of a boundary value problem for the Poisson equation 
with a nonlocal boundary operator},Acta Mathematica Scientia.  
35, No.5 (2015),  970-–980.

\bibitem{karachik-sam} V. V. Karachik, B. Kh. Turmetov, B. T. Torebek;
\emph{On some integro-differential operators in the class of
harmonic functions and their applications}, Siberian Advances in Mathematics.
 22,  No.2 (2012),  115--134.

\bibitem{kilbas-book} A. A. Kilbas, H. M. Srivastava, J. J. Trujillo;
\emph{Theory and Applications of Fractional Differential
Equations}, Elsevier. North-Holland. Mathematics Studies. (2006) 539p.

\bibitem{kirane-smj1} M. Kirane, N. -e. Tatar;
\emph{Absence of local and global solutions to an elliptic system
with time-fractional dynamical boundary conditions}, Siberian
Mathematical Journal. 48, No.3 (2007), 593--605.

\bibitem{kirane-smj2} M. Kirane, N. -e. Tatar;
\emph{Nonexistence for the Laplace equation with a dynamical
boundary condition of fractional type}, Siberian Mathematical
Journal 48 No.5 (2007), 1056--1064.

\bibitem{kirane-MMAS} M. Kirane, B. T. Torebek; 
\emph{On a nonlocal problem for the Laplace equation in the unit ball 
with fractional boundary conditions}, Mathematical Methods in the Applied Sciences. 
 2015, doi: 10.1002/mma.3554.

\bibitem{Krasnoschok} M. Krasnoschok,  N. Vasylyeva;
\emph{On a nonclassical fractional boundary-value problem for the Laplace operator},
Journal of Differential Equations. 257 (2014),   1814--1839.

\bibitem{Muratbekova} M. A. Muratbekova, K. M. Shinaliyev, B. Kh. Turmetov; 
\emph{On solvability of a nonlocal problem for the Laplace equation with the 
fractional-order boundary operator},Boundary Value Problems. 
2014, No.29 (2014), doi:10.1186/1687-2770-2014-29.

\bibitem{sadybekov-de} M. A. Sadybekov, B. Kh. Turmetov; 
\emph{On an analog of periodic boundary value problems for the Poisson 
equation in the disk},Differential Equations.  50, No.2(2014),  268--273.

\bibitem{sadybekov-emj}     M. A. Sadybekov, B. Kh. Turmetov; 
\emph{On analogues of periodic boundary value problems for the Laplace 
operator in a ball}, Eurasian Math. J., 3, No.1 (2012), 143-–146.

\bibitem{sadybekov-ejd} 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.
 2014, No.157 (2014), 1--14.

\bibitem{torebek-bvp} B. T. Torebek, B. Kh. Turmetov;
\emph{On solvability of a boundary value problem for the Poisson
equation with the boundary operator of a fractional order},
Boundary Value Problems. 2013, No.93 (2013),
doi:10.1186/1687-2770-2013-93.

\bibitem{turmetov-de1} B. Kh. Turmetov;
\emph{On a boundary value problem for a harmonic equation},
Differential Equations. 32, No.8 (1996),  1093--1096.

\bibitem{turmetov-sam} B. Kh. Turmetov;
\emph{On Smoothness of a Solution to a Boundary Value Problem with
Fractional Order Boundary Operator}, Siberian Advances in Mathematics. 15,
No.2 (2005), 115--125.

\bibitem{turmetov-ejde} B. Kh. Turmetov;
\emph{Solvability of fractional analogues of the neumann problem for 
a nonhomogeneous biharmonic equation}, Electronic Journal of 
Differential Equations. 2015, No. 82 (2015),  1--21.

\bibitem{turmetov-de2} B. Kh. Turmetov, B. T. Torebek;
\emph{Modified Bavrin operators
and their applications}, Differential Equations.  51, No.2 (2015),  243--254.

\bibitem{umarov-dm} S. R. Umarov; 
\emph{ On some boundary value problems for elliptic equations with
a boundary operator of fractional order}, Doklady Russian Academy
of Science. 333,  No.6 (1993), 708--710.

\bibitem{umarov-fcaa} S. R. Umarov, Yu. F. Luchko, R. Gorenflo;
\emph{On boundary value problems for elliptic equations with
boundary operators of fractional order}, Fractional Calculus and
Applied Analysis. 2,  No.4 (2000), 454--468.

\bibitem{Veliev} E. I. Veliev, M. V. Ivakhnychenko ;
\emph{ Fractional boundary conditions in plane wave diffraction on a strip}, Progress In
Electromagnetics Research. 79 (2008),   443--462.


\end{thebibliography}

\end{document}
