\documentclass[reqno]{amsart}
\usepackage{hyperref}

\AtBeginDocument{{\noindent\small
\emph{Electronic Journal of Differential Equations},
Vol. 2017 (2017), No. 287, pp. 1--18.\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/287\hfil Direct and inverse Robin-Regge problems]
{Direct and inverse Robin-Regge problems}

\author[M. M\"oller, V. Pivovarchik \hfil EJDE-2017/287\hfilneg]
{Manfred M\"oller, Vyacheslav Pivovarchik}

\address{Manfred  M\"oller \newline
The John Knopfmacher Centre for Applicable
Analysis and Number Theory,
School of Mathematics,
University of the Witwatersrand,
Johannesburg, WITS 2050, South Africa}
\email{manfred.moller@wits.ac.za}

\address{ Vyacheslav Pivovarchik \newline
Department of Algebra and Geometry,
South-Ukrainian National Pedagogocal University,
Staroportofrankovskaya str. 26,
Odessa, 65020, Ukraine}
\email{vpivovarchik@gmail.com}

\dedicatory{Communicated by Tuncay Aktosun}

\thanks{Submitted August 8, 2017. Published November 17, 2017.}
\subjclass[2010]{34B09, 34B07, 34B24, 34L20, 34A55}
\keywords{String problem; viscous damping; Robin boundary condition;
\hfill\break\indent  generalized Regge boundary condition; spectral asymptotics;
 inverse problems}

\begin{abstract}
 The Sturm-Liouville equation is considered on a bounded interval
 with a Robin boundary condition at the left end and a generalized
 Regge condition at the right end. Properties of the spectrum of
 this problem are derived, including the location of the spectrum
 in the lower half-plane  and the asymptotic distribution of the
 eigenvalues. The inverse problem is solved, given one spectrum or
 parts of two spectra.
\end{abstract}

\maketitle
\numberwithin{equation}{section}
\newtheorem{theorem}{Theorem}[section]
\newtheorem{lemma}[theorem]{Lemma}
\newtheorem{proposition}[theorem]{Proposition}
\newtheorem{remark}[theorem]{Remark}
\newtheorem{definition}[theorem]{Definition}
\newtheorem{corollary}[theorem]{Corollary}
\allowdisplaybreaks



\section{Introduction}

It is well known that the spectra of the Neumann-Dirichlet
(or the Dirichlet-Neumann) and the Dirichlet-Dirichlet boundary-value problems
 for the Sturm-Liouville equation
\begin{equation}\label{eq:1}
-y''+q(x)y=\lambda ^2y, \quad x\in[0,a],
\end{equation}
 with spectral parameter $\lambda $,  generated by the same potential $q$, 
uniquely determine this potential in $L_2(0,a)$. It is also known that the 
spectra of two boundary-value problems with the same Robin boundary condition
at one of the ends and different Robin conditions at the other end of the
interval uniquely determine the potential and the constants
in the boundary conditions. These results are due to Borg \cite{B} 
(see also \cite{L,LG,M}). If the boundary conditions are given data, 
then the problem of recovering the potential appears to be overdetermined 
in case of Robin conditions. One does not need to know all the
eigenvalues of the two spectra in this case.  In a further development 
of this theory one spectrum together with knowledge of a part of the potential
 is used, see \cite{GRS,Ha,HL,MP,Pi4,Su,WX}.

Another direction of generalization of the above results is to use eigenvalues 
of more than two spectra to determine the potential.
In \cite{GS} it was shown that two thirds of the union of  three spectra of 
boundary-value problems with the same boundary condition at one of the
ends uniquely determines the potential. In \cite{H1} a similar but more 
general sufficient condition of unique solvability  was given for the case 
when the known eigenvalues were taken from $n$ different spectra, 
see \cite{H2} for a topical review.

It is also known that if we impose the so called generalized Regge condition
\begin{equation}\label{eq:2}
y'(a)+(i\lambda\alpha+\beta) y(a)=0
\end{equation}
with $\alpha>0$ at the right end of the interval, then we need roughly 
speaking only one spectrum instead of two spectra. 
This was shown first in \cite{KN1,KN2} for the boundary-value problem generated 
by the string equation. For the Sturm-Liouville problem this was
shown in \cite{PvM} where the condition at the left was $y(0)=0$. 
In the case of the Dirichlet-Regge problem with $\alpha =1$ and $\beta =0$, 
Regge \cite{Reg} described how to recover the potential from the spectrum 
using the Jost function and the Gelfand-Levitan equation. A detailed account of
the spectral asymptotics for the Dirichlet-Regge problem, 
including the case $\alpha =1$, is given in \cite[Section 7.2]{MPb}.

In this paper we consider the Sturm-Liouville problem \eqref{eq:1} 
with (generalized) Regge boundary condition \eqref{eq:2} 
at the right end and Robin boundary condition
\begin{equation}\label{eq:3}
y'(0)-\delta y(0)=0
\end{equation}
 at the left end.
Our aim in this paper is to describe the spectrum of this problem and to solve 
the inverse problems by one spectrum and by parts of two spectra with two 
different $\delta $ in the Robin condition \eqref{eq:3}.

In Section 2  we show that the problem of small transverse vibrations of a 
smooth inhomogeneous string damped at the right end can be reduced to 
\eqref{eq:1}--\eqref{eq:3} by means of the Liouville transformation. 
Also in that section we show that all eigenvalues of this string problem 
lie in the upper half-plane.

In Section 3 we investigate the spectrum of \eqref{eq:1}--\eqref{eq:3} 
admitting eigenvalues in the lower half-plane, and in particular the 
asymptotic distribution of eigenvalues is given.

 In Section 4 we solve the corresponding inverse problem, namely we show 
that if a sequence of complex numbers $(\lambda_k)$ satisfies certain conditions 
(which appear to be necessary and sufficient), then there exist a unique  real
 $L_2(0,a)$ potential and a unique set of real parameters $\delta$, $\alpha>0$  
and $\beta$ which generate a Robin-Regge problem with the spectrum   $(\lambda_k)$.
 A corresponding existence and uniqueness result is obtained for the string problem.

 In Section 5 we show that the characteristic functions of two Robin-Regge 
problems with the same Regge condition at the right end and different Robin 
conditions at the left satisfy a functional equation which was investigated 
in \cite{MP3}. Using the results of \cite{MP3} and the results of Section 4 
we solve the following inverse problem: given two sets of complex numbers, 
given $\delta_1-\delta_2 $, find the   potential and a set of real 
parameters $\delta_1$, $\alpha>0$  and $\beta$ which generate a 
Robin-Regge problems  with the same Regge condition at the right end and 
with different coefficients $\delta_2$ and $\delta_1$ in the  Robin condition 
at the left end for which the given sets of complex numbers are parts of 
spectra with $\delta _1$ and $\delta _2$, respectively.


\section{String problem}

Small transversal vibrations of a smooth inhomogeneous string subject to 
viscous damping are described by the boundary-value problem
\begin{gather}
\label{2.1}
\frac{\partial^2 u}{\partial s^2}-\rho(s)\frac{\partial^2 u}{\partial t^2}=0,\\
\label{2.2}
\Big(\frac{\partial u}{\partial s}-pu\Big)\Big|_{s=0}=0,\\
\label{2.3}
\Big(\frac{\partial u}{\partial s}+\nu\frac{\partial u}{\partial t}\Big)\Big|_{s=l}=0.
\end{gather}
Here $u=u(s,t)$ is the transversal displacement of a point of the string 
which is as far as $s$ from the left end of the string at time $t$, $l$ 
is the length of the string, $\rho\geq \varepsilon>0$ is its density, and 
the constant $p>0$ is the Robin coefficient. We will assume 
$\rho \in  L_\infty (0,l)$. The left end of the string is constrained 
elastically while the right end is free to move in the direction
orthogonal to the equilibrium position of the string subject to damping 
with damping coefficient $\nu>0$. Substituting $u( s,t)=v(s)e^{i\lambda t}$ 
we arrive at
\begin{gather} \label{2.4}
\frac{\partial^2 v}{\partial s^2}+\lambda^2\rho(s) v= 0, \\
\label{2.5}
\Big(\frac{\partial v}{\partial s}-pv\Big)\Big|_{s=0}=0,\\
\label{2.6}
\Big(\frac{\partial v}{\partial s}+i\lambda\nu v\Big)\Big|_{s=l}=0.
\end{gather}
The eigenvalue problem  \eqref{2.4}--\eqref{2.6} is described by
the operator pencil
\begin{equation}\label{eq:2.7a}
L(\lambda )=\lambda ^2M-i\lambda K-A,
\end{equation}
where
 the operators $A$, $K$ and $M$ act in the Hilbert space 
$H=L_2(0,l)\oplus \mathbb{C}$ according to
\begin{gather} \label{2.7}
A\begin{pmatrix} v \\ c \end{pmatrix}=\begin{pmatrix} - v'' \\ v'(l) \end{pmatrix},\\
\label{2.8}
D(L):=D(A)=\Big\{\begin{pmatrix} v \\ c \end{pmatrix}: 
 v\in W_2^2(0,l), \ v'(0)-p v(0)=0, \ c=v(l)\Big\},\\
\label{2.9}
M=\begin{pmatrix} \rho  & 0 \\ 0 & 0 \end{pmatrix}, \quad   
 K=\begin{pmatrix}0  & 0 \\ 0 & \nu \end{pmatrix},
\end{gather}
where $W_2^2(0,l)$ is the Sobolev space of order $2$ over the interval $(0,l)$.

Recall that an operator $A$ is said to be strictly positive, written $A\gg 0$, 
if there is $\varepsilon>0$ such that $(AY,Y)\ge \varepsilon (Y,Y)$ for 
all $Y$ in the domain of $A$.

\begin{proposition} \label{prop2.1}
The operators $A$, $K$ and $M$ are self-adjoint, $M\ge 0$ and $K\ge 0$ are
 bounded, $M+K\gg 0$, and $A\gg0$  has a compact resolvent.
\end{proposition}

\begin{proof}
The statements about $M$ and $K$ are obvious. The proof of the statements 
about $A$ is similar to that of \cite[Proposition 2.1.1]{MPb}.
To prove that $A$ is self-adjoint, we will use \cite[Theorem 10.3.5]{MPb}, 
so that we have to find $U_3$, $U_1$ and $U$ defined in
\cite[(10.3.12) and (10.3.13)]{MPb}. Since the differential expression 
in $A$ is $-v''=(-v')'$, the first two quasi-derivatives of $v$ are 
$v^{[0]}=v$ and  $v^{[1]}=-v'$, see  \cite[Definition 10.2.1]{MPb}. 
It is easy to see that, in the notation of  \cite[Section 10.3]{MPb},
\[
 U_1=\begin{pmatrix} -p&-1&0&0 \end{pmatrix} , \quad
 U_2=\begin{pmatrix} 0&0&1&0 \end{pmatrix} , \quad
 V=\begin{pmatrix} 0&0&\tilde\beta &-1 \end{pmatrix} ,
\]
with $\tilde \beta =0$, so that
\[
U_3 =\begin{pmatrix} 0&-1&0&0\\1&0&0&0\\0&0&0&1\\0&0&-1&0\\0&0&\tilde\beta &-1\\
0&0&-1&0 \end{pmatrix},\quad
U=\begin{pmatrix} -p&-1&0&0&0&0\\0&0&1&0 &-1&0\\ 0&0&\tilde\beta &-1&0&-1\end{pmatrix}.
\]
We have introduced the parameter $\tilde \beta $ since we will encounter the
 matrix $V$ with $\tilde \beta \not=0$ in the next section.
Denoting the $j$-th standard basis eigenvector in $\mathbb{C}^n$ by $e_j$, 
it follows that
\begin{gather*}
N(U_1)=\operatorname{span}\{e_1-pe_2,e_3,e_4\},\\
U_3(N(U_1))=\operatorname{span}\{pe_1+e_2,e_4-\tilde\beta e_5 +e_6,e_3-e_5\},\\
R(U^*)=\operatorname{span}\{pe_1+e_2,e_3-e_5,e_4-\tilde\beta e_3+e_6\},
\end{gather*}
so that   $U_3(N(U_1))=R(U^*)$. Hence $A$ is self-adjoint by 
Theorem  \cite[Theorem 10.3.5]{MPb}.
Now it follows from \cite[Theorem 10.3.8]{MPb} that $A$ has a compact 
resolvent and that $A$ is bounded below.

Finally, for $Y=(v,c)^{\rm T}\in D(A)$ we conclude in view of
\cite[(10.2.5)]{MPb} that
\begin{equation}
\begin{aligned}
(AY,Y)&=\int_0^a|v'(x)|^2dx-v'(x)\overline{v(x)}|_0^l+v'(l)\overline{c}\\
&=\int_0^a|v'(x)|^2dx-v'(l)\overline{v(l)}+v'(0)\overline{v(0)}+v'(l)\overline{c}\\
&=\int_0^a|v'(x)|^2dx+p|v(0)|^2\ge 0.
\end{aligned} \label{eq:Apos}
\end{equation}
Furthermore, if $(AY,Y)=0$, then $v'=0$ and $v(0)=0$, which gives $v=0$
and then $c=v(0)=0$.
Therefore $A>0$ has been shown. Since $A$ has a compact resolvent and
thus a discrete spectrum consisting only of eigenvalues, it follows
 that $A\gg0$.
\end{proof}

\begin{proposition}\label{prop2.2}
All eigenvalues of \eqref{2.4}--\eqref{2.6}  lie in the open upper half-plane.
\end{proposition}

\begin{proof}
Let  $Y=(v,c)^{\rm T}\in D(A)\cap N(K)$ and $\lambda \in \mathbb{R}$
such that $\lambda ^2MY-AY=0$. It follows that $v''+\lambda ^2\rho y=0$, 
$v(l)=c=0$ and $v'(l)=0$, and the
unique solution of this initial value problem is $v=0$.
 An application of   \cite[Lemma 1.2.4, 3]{MPb} completes the proof.
\end{proof}

For the spectral asymptotics and for the inverse problem it is more convenient 
to write \eqref{2.4} as a Sturm-Liouville equation. Hence
we give an alternate approach under the additional
 assumption that  $\rho\in W_2^2(0,l)$. 
Then the Liouville transform \cite[p.\ 292]{CH}
\begin{gather}
\label{2.10}
x(s)=\int_0^s\rho^{1/2}(r)\,dr,\quad 0\le s\le l,\\
\label{2.11}
y(\lambda,x)=\rho^{1/4}(s(x))v(\lambda,s(x)),\quad
0\le x\le a,\; \lambda \in \mathbb{C},
\end{gather}
leads to the equivalent boundary-value problem
\begin{gather}
\label{2.12}
y'' -q(x)y+\lambda^2y=0,\\
\label{2.13}
y'(\lambda, 0)-\delta y(\lambda,0)=0,\\
\label{2.14}
y'(\lambda,a)+(i\lambda\alpha  +\beta)y(\lambda,a)=0,
\end{gather}
where
\begin{gather}
\label{2.15}
q(x)=\rho^{-1/4}(s(x))\frac{d^2}{dx^2}\rho^{1/4}(s(x)),\\
\label{2.16}
a=\int_0^l\rho^{1/2}(r)\,dr,\\
\label{2.17}
\beta=-\rho^{-1/4}(s(a))\frac{d\rho^{1/4}(s(x))}{dx}\Big|_{x=a},\\
\label{2.18}
\delta=p\rho^{-\frac{1}{2}}(0)+\rho^{-1/4}(0)
\frac{d\rho^{1/4}(s(x))}{d x}\Big|_{x=0},\\
\label{2.19}
\alpha=\rho^{-\frac{1}{2}}(s(a))\nu>0.
\end{gather}

The eigenvalue problem \eqref{2.12}--\eqref{2.14} is described by
the operator pencil
\begin{equation}\label{eq:Ltilde}
\tilde{L}(\lambda )=\lambda^2\tilde{M}-i\lambda\tilde{K}-\tilde{A},
\end{equation}
where
 the operators $\tilde{A}$, $\tilde{K}$ and $\tilde{M}$ act in the Hilbert space
 $H=L_2(0,a)\oplus \mathbb{C}$ according to
\begin{gather}
\label{2.20}
\tilde{A}\begin{pmatrix} y \\ c \end{pmatrix}=\begin{pmatrix} - y''+qy \\
y'(a)+\beta y(a) \end{pmatrix},\\
\label{2.21}
D(\tilde L):=D(\tilde{A})=\big\{\begin{pmatrix} y \\ c \end{pmatrix}: 
 y\in W_2^2(0,a), \ y'(0)-\delta y(0)=0, \ c=y(a) \big\},\\
\label{2.22}
\tilde{M}=\begin{pmatrix} I & 0 \\ 0 & 0 \end{pmatrix},
 \quad \tilde K=\begin{pmatrix} 0 & 0 \\ 0 &  \alpha  \end{pmatrix}.
\end{gather}
Since \eqref{2.11} establishes an isomorphism between $L_2(0,l)$ and 
$L_2(0,a)$ and the corresponding Sobolev spaces of order $2$,
it is clear that $\tilde{M}$, $\tilde{K}$ and $\tilde{A}$ have the same 
properties as $M$, $K$, $A$, respectively, stated in Proposition \ref{prop2.1}. 
In particular, $\tilde A$ is self-adjoint and satisfies $\tilde{A}\gg0$.

\section{Robin-Regge problem and its spectrum}
\subsection{Robin-Regge problem}

In Section 2 we have seen that the string problem \eqref{2.4}--\eqref{2.6} 
can be transformed into a problem of the form
\eqref{2.12}--\eqref{2.14}. In this section we will investigate this latter problem:
\begin{gather}
\label{3.1} -y''+q(x)y=\lambda^2y, \quad x\in[0,a]\\
\label{3.2}
 y'(0)-\delta y(0)=0,\\
\label{3.3}
y'(a)+(i\alpha\lambda+\beta) y(a)=0,
\end{gather}
where we admit general real-valued functions $q\in L_2(0,a)$ and real numbers 
$\alpha >0$, $\beta $ and $\delta $.
We call this problem a Robin-Regge boundary-value problem since the boundary
 condition \eqref{3.2}
is also called a boundary condition of Robin type, see \cite[p.{} 70]{GA}, 
 and   the boundary condition  \eqref{3.3} is of (generalized) Regge type, 
see \cite[(2.1.5)]{MPb}.

 This problem has a representation as a quadratic operator pencil $\tilde L$
as in  \eqref{eq:Ltilde} with the operators $\tilde A$, $\tilde M$, $\tilde K$ 
given by \eqref{2.20}--\eqref{2.22}.
The proof of Proposition \ref{prop2.1}, with $\tilde \beta =\beta $ and except 
for the last paragraph concerning positivity of $A$,
also applies to these operators. Hence we have

\begin{proposition} \label{prop3.1}
The operators $\tilde A$, $\tilde K$ and $\tilde M$ are self-adjoint, 
$\tilde M\ge 0$ and $\tilde K\ge 0$ are bounded, $\tilde M+\tilde K\gg 0$, 
and $\tilde A$ is bounded below and has a
compact resolvent.
\end{proposition}

Therefore \cite[Lemmas 1.1.11 and 1.2.1 and Theorem 1.3.3]{MPb} give the following
result.

\begin{proposition}\label{propsp1}
The spectrum of the pencil $\tilde L$ consists of isolated eigenvalues of finite 
algebraic multiplicity  located in the closed upper half-plane and on the 
imaginary axis, at most finitely many eigenvalues lie on the negative imaginary 
axis, and the total algebraic multiplicity of the eigenvalues
on the negative imaginary axis equals the total multiplicity of the negative 
eigenvalues of $\tilde A$.
The eigenvalues, with multiplicity, are symmetric with respect to the imaginary axis.
\end{proposition}

\subsection{Characteristic function}\label{subs:3.2}

We denote by $\mathcal{L}^{a}$ the Paley-Wiener class of
entire functions of exponential type $\leq a$ whose restrictions to 
the real axis belong to $L_2(-\infty,\infty)$.

Every solution of \eqref{3.1} which satisfies the Robin
 condition \eqref{3.2} is a multiple of the solution $c(\lambda ,\delta ,\cdot )$ 
of \eqref{3.1}  which satisfies the initial conditions $c(\lambda ,\delta ,0)=1$, 
$c'(\lambda ,\delta ,0)=\delta $.

 \begin{proposition}\label{propr}
 The function $c$ has the representation
\begin{equation}
 \begin{aligned}
 c(\lambda,\delta,x)
&=\cos\lambda x+\int_0^xB(x,t,\delta)\cos\lambda t\,dt\\
&=\cos \lambda x+B(x,x,\delta)\frac{\sin\lambda x}{\lambda}
 -\int_0^xB_{t}(x,t,\delta)\frac{\sin\lambda t}{\lambda}\,dt,
\end{aligned}\label{3.4}
\end{equation}
where
\begin{gather*}
B(x,t,\delta)=\delta+\tilde{K}(x,t)+\tilde{K}(x,-t)
+\delta\int_t^x[\tilde K(x,\xi)-\tilde K(x,-\xi)]d\xi, \\
B_{t}(x,t,\delta)=\frac{\partial B(x,t,\delta)}{\partial t}\,,
\end{gather*}
and $\tilde{K}(x,t)$ is the unique solution of the integral equation
\[
\tilde{K}(x,t)=\frac{1}{2}\int_0^{\frac{x+t}{2}}q(\zeta)d\zeta
+\int_0^{\frac{x+t}{2}}\int_0^{\frac{x-t}{2}}
q(\zeta+\xi)\tilde{K}(\zeta+\xi,\zeta-\xi)\,d\xi\,d\zeta
\]
in the region $\{(x,t)\in [0,a]\times [0,a]:|t|\le x\}$.
The solution $\tilde{K}(x,t)$ possesses first order partial derivatives
each belonging to $L_2\left(0,a\right)$ as a function
of each of its variables. Moreover,
\[
B(x,x,\delta)=\delta+\frac{1}{2}\int_0^xq(t)\,dt.
\]
\end{proposition}

\begin{proof}
Let $s(\lambda ,\cdot )$ and $c(\lambda,0 ,\cdot )$ be the solutions of 
the initial value problems for
\eqref{3.1} with the initial conditions $s(\lambda ,0)=0$, $s'(\lambda ,0)=1$ and
$c(\lambda,0 ,0)=1$, $c'(\lambda,0 ,0)=0$, respectively. Then we obviously have
\begin{equation}\label{eq:r}
c(\lambda,\delta,x)=c(\lambda,0,x)+\delta s(\lambda,x).
\end{equation}
 From \cite[Theorem 12.2.9]{MPb}  we know the representations
\begin{gather*}
s(\lambda,x)=\frac{\sin\lambda x}{\lambda}+\int_0^xK(x,t)\frac{\sin\lambda t}
{\lambda}\,dt,\\
c(\lambda,0,x)=\cos\lambda x +\int_0^xB(x,t,0)\cos\lambda t\,dt,
\end{gather*}
where $K(x,t)=\tilde{K}(x,t)-\tilde{K}(x,-t)$,
$B(x,t,0)=\tilde{K}(x,t)+\tilde{K}(x,-t)$ and
\[
B(x,x,0)=K(x,x)=\tilde{K}(x,x)=\frac{1}{2}\int_0^xq(t)\,dt.
\]
Substituting these representations into
\eqref{eq:r} and straightforward integration by parts completes the proof.
\end{proof}

It is clear that
\begin{equation}\label{3.5}
c'(\lambda,\delta, x)=-\lambda\sin\lambda x+B(x,x,\delta)\cos \lambda x
+\int_0^xB_x(x,t,\delta)\cos\lambda t \,dt.
\end{equation}
Using \eqref{3.4} and \eqref{3.5} we obtain
\begin{gather}
\label{3.6}
c(\lambda,\delta,a)=\cos\lambda a+B\frac{\sin\lambda a}{\lambda}
 + \frac{\psi(\lambda)}{\lambda},\\
\label{3.7}
 c'(\lambda,\delta,a)=-\lambda\sin\lambda a+B\cos\lambda a
 +\tilde{\psi}(\lambda),
\end{gather}
where $B:=B(a,a, \delta)$,  $\psi\in \mathcal{L}^{a}$,
 $\tilde{\psi}\in \mathcal{L}^{a}$.
 Moreover,
$\psi(0)=0$ since $c(\cdot ,\delta ,a)$ is an entire function.
Here we have used that by Paley-Wiener's theorem
$\mathcal{L}^a$-functions are the Fourier transformations of square summable
functions supported on $[0,a]$. 
It is clear that $\lambda \in \mathbb{C}$ belongs to
the spectrum of \eqref{3.1}--\eqref{3.3} if and only if
 $c(\lambda ,\delta ,\cdot )$
satisfies the boundary condition \eqref{3.3}, that is, if and only if
\begin{equation} \label{3.8}
\phi(\lambda,\delta):=c'(\lambda,\delta,a)+(i\alpha\lambda+\beta)c(\lambda,\delta,a)
\end{equation}
satisfies $\phi (\lambda ,\delta )=0$.
 The function $\phi (\cdot ,\delta )$ is called the characteristic function of 
\eqref{3.1}--\eqref{3.3}.
It is well known, see e.\,g.{} \cite[Section 6.3]{MM}, 
 that the algebraic multiplicity of the eigenvalue
$\lambda $ equals the multiplicity of the zero $\lambda $ of 
$\phi (\cdot ,\delta )$.

\subsection{Spectrum in the lower half-plane}

\begin{theorem}\label{thmgRHB}
The eigenvalues of the Robin-Regge problem \eqref{3.1}--\eqref{3.3} have the  
following properties:
\begin{itemize}
\item[(1)]  Only a finite number of the eigenvalues lie in the closed lower
 half-plane.

\item[(2)] All nonzero eigenvalues in the closed lower half-plane lie on the 
negative imaginary semiaxis and are
simple. If their number $\kappa$ is positive,  they will be uniquely indexed as
$\lambda_{-j}=-i|\lambda_{-j}|$, $j=1,\dots,\kappa$,  satisfying
$|\lambda_{-j}|<|\lambda_{-(j+1)}|$, $j=1,\dots,\kappa-1$.

\item[(3)] If $\kappa>0$, then the numbers
$i|\lambda_{-j}|$, $j=1,\dots,\kappa$, are not eigenvalues.

\item[(4)] If  $\kappa\ge2$, then  in each of the
intervals $(i|\lambda_{-j}|,i|\lambda_{-(j+1)}|)$, $j=1,\dots,\kappa-1$,  
the number of eigenvalues, counted with multiplicity, is odd.

\item[(5)] If $\kappa >0$ and $0$ is not an eigenvalue, then the interval 
$(0,i|\lambda _{-1}|)$
contains an even number of eigenvalues, counted with multiplicity, or does not
contain any eigenvalues.

\item[(6)] If $\kappa >0$ and $0$ is an eigenvalue, then $0$ is a simple
eigenvalue and the interval $(0,i|\lambda _{-1}|)$ contains an odd number  
of eigenvalues, counted with multiplicity.

\item[(7)] If $\alpha \neq  1$, then the Robin-Regge problem has infinitely
many eigenvalues, which lie in a horizontal strip of the
complex plane.
\end{itemize}
\end{theorem}

\begin{proof}
First we are going to show that there are no eigenvalues $\lambda $ with 
corresponding eigenvectors  of the form $(y,0)^{\rm T}\in D(\tilde A)$. 
Indeed, otherwise the definition of
$D(\tilde A)$ in \eqref{3.23} would imply $y(a)=c=0$, and \eqref{3.3} then 
would give $y'(a)=0$.
Hence $y$ would be a solution of the initial value problem \eqref{3.1},
 $y(a)=y'(a)=0$, which has only the  trivial solution, i.\,e., $y=0$. 
 But this is impossible since $(y,0)^{\rm T}$ was supposed to be an eigenvector.

In view of \cite[Theorem 1.5.6, 1, and Remark 1.5.8]{MPb},
there are no nonzero eigenvalues of type I, i.\,e., eigenvalues $\lambda $ 
such that also $-\lambda $ is an eigenvalue. Therefore all nonzero eigenvalues
are of type II according to \cite[Definition 1.5.2]{MPb}. 
By \cite[Remark 1.5.8]{MPb} there are no nonzero real
eigenvalues of type II, and it follows that there are no nonzero real 
eigenvalues at all.
Therefore statements (1) through (5) follow from \cite[Theorem 1.5.7]{MPb}.

To prove statement (6) we have to consider the case that $0$ is an eigenvalue. 
Then the  first paragraph of this proof shows that 
$N(\tilde A)\not\subset N(\tilde K)$, and $0$ is not an eigenvalue of type II
 by \cite[Remark 1.5.3, 2]{MPb}.   Furthermore,
$N(\tilde M)\cap N(\tilde A)\subset N(M)\cap D(\tilde A)=\{0\}$, and
the geometric multiplicity of the eigenvalue is $1$ by the first paragraph of 
this proof. Next we
are going to show that the eigenvalue is simple, that is, there are no 
associated vectors. By proof of contradiction, assume there are 
$Y_0=(y_0,c_0)^{\rm T}\neq 0$ and
$Y_1=(y_1,c_1)^{\rm T}$ in $D(A)$ such that
\[\tilde A Y_0=0\quad\text{and}\quad -i\tilde K Y_0-\tilde AY_1=0.\]
We recall that $c_0\not=0$. It follows that
\[ -y_j''+qy_j=0,\ y_j'(0)-\delta y_j(0)=0,\quad j=0,1,\]
so that $y_1$ is a multiple of $y_0$, say $y_1=dy_0$. But then also $Y_1=dY_0$, and
\[
(0,c_0)^{\rm T}=\frac 1\alpha \tilde K Y_0=i\frac 1\alpha \tilde AY_1=
 i\frac d\alpha \tilde AY_0=0,
\]
 which is impossible. Hence statement (6) also follows from \cite[Theorem 1.5.7]{MPb}.

We have established in the previous subsection that the spectrum coincides 
 with the zeros of  $\phi (\cdot ,\delta )$, and an application of 
\cite[Lemma 7.1.3]{MPb} shows that statement (7) holds.
\end{proof}

\subsection{Spectral asymptotics}

 For $\alpha \ne1$ the Robin-Regge problem has infinitely many eigenvalues,
  and their the asymptotic representation is given in the following theorem.

\begin{theorem}\label{thm3.3}
\textup{(1)} If $\alpha\in (0,1)$ then the eigenvalues
$(\xi_k)_{k=-\infty, \ k\not=0}^{\infty}$ of problem  \eqref{3.1}--\eqref{3.3} 
can  be indexed such that $\xi_{-k}=-\overline{\xi_k}$ for all not pure imaginary $\xi_k$ and behave
asymptotically as follows:
\begin{equation}
\label{3.11}
\xi_k=\frac{\pi (k-1)}{a}+\frac{i}{2a}\ln\Big(\frac{\alpha+1}{1-\alpha}\Big)
+\frac{P}{k}+\frac{\gamma_k}{k},   \quad k\in\mathbb{N},
\end{equation}
where
\begin{equation}\label{3.11a}
P=\frac{1}{\pi}\Big(B-\frac{\beta}{\alpha^2-1}\Big),
\end{equation}
and where $(\gamma_k)_{k=1}^{\infty}\in l_2$.\\
\textup{(2)} If $\alpha\in (1,\infty)$ then the eigenvalues
$(\xi_k)^{\infty}_{k=-\infty}$ of problem  \eqref{3.1}--\eqref{3.3} can 
 be indexed such that $\xi_{-k}=-\overline{\xi_k}$ for all not pure imaginary $\xi_k$ and behave
asymptotically as follows:
\begin{equation}
\label{3.12}
\xi_k=\frac{\pi \left(k-\frac{1}{2}\right)}{a}+\frac{i}{2a}
\ln\left(\frac{\alpha+1}{\alpha-1}\right)
+\frac{P}{k}+\frac{\tilde{\gamma}_k}{k},  \quad  k\in\mathbb{N},
\end{equation}
where $P$ is given by \eqref{3.11a} and where 
$(\tilde{\gamma}_k)^{\infty}_{k=1}\in l_2$.
\end{theorem}

\begin{proof}
Using \eqref{3.6}, \eqref{3.7} and \eqref{3.8} we obtain that the characteristic 
function $\phi $ of \eqref{3.1}--\eqref{3.3}
given in \eqref{3.8} satisfies
\begin{equation}
\label{3.9}
-i\phi(\lambda,\delta)=\lambda[\alpha \cos\lambda a+i\sin\lambda a]
-i(B+\beta)\cos\lambda a+\alpha B\sin\lambda a+\hat{\psi}(\lambda),
\end{equation}
where $\hat{\phi}\in \mathcal{L}^a$. The right hand side of \eqref{3.9} is 
of the form \cite[(7.1.4)]{MPb} with $M=\alpha B$ and $N=B+\beta $. 
Then the number $P$ defined in \cite[Lemma 7.1.3]{MPb} becomes
\begin{equation}\label{eq:P}
P=\frac{N-\alpha M}{\pi(1-\alpha ^2)}= \frac{B+\beta -\alpha^2 B}{\pi(1-\alpha ^2)}
  =\frac B\pi + \frac{\beta}{\pi(1-\alpha ^2)} \,,
\end{equation}
and \eqref{3.11} and \eqref{3.12} follow from \cite[Lemma 7.1.3]{MPb}.
\end{proof}

\section{Inverse problem}\label{s:inv}

\begin{definition}\label{defnmore} \rm
(1) The function $\theta $ is said to be a Nevanlinna function if: 
\begin{itemize}
\item[(i)] $\theta $ is analytic in the half-planes $\operatorname{Im}\lambda>0$ and
$\operatorname{Im}\lambda<0$;
\item[(ii)] $\theta (\overline{\lambda})=\overline{\theta (\lambda)}$
if $\operatorname{Im}\lambda\ne0$;  
\item[(iii)] $\operatorname{Im}\lambda  \operatorname{Im} \theta (\lambda)\ge 0$ for
$\operatorname{Im}\lambda\ne0$.
\end{itemize}

\noindent(2) The class ${\mathcal N}^{\text{\rm ep}}$ of essentially positive 
Nevanlinna functions is the set of all
Nevanlinna functions  which are analytic in $\mathbb{C}\setminus  [0,\infty)$  
with the possible exception of finitely many poles. 

\noindent(3) The class ${\mathcal N}_+^{\text{\rm ep}}$ is  the set of all 
functions $\theta \in{\mathcal N}^{\text{\rm ep}}$ such that for
some $\gamma\in \mathbb{R}$ we have $\theta (\lambda )>0$ for all 
$\lambda \in(-\infty,\gamma)$.
\end{definition}

\begin{definition}\label{def2.1} \rm
Let $P$ and $Q$ be entire functions with no common zeros such that $P$ and 
$Q$ are real on the real axis and such that
$\frac{Q}{P}$ belongs to ${\mathcal N}^{\text{\rm ep}}_{+}$. Then the
function $\omega $ defined by $\omega (\lambda )=P(\lambda^2)+i\lambda Q(\lambda^2)$ 
is said to belong to the class of symmetric shifted Hermite-Biehler functions.
The class of all symmetric shifted Hermite-Biehler functions is denoted by SSHB.
If the number of negative zeros of $P$
is $\kappa$, then we say that $\omega $
belongs to SSHB$_{\kappa}$.
\end{definition}

\begin{definition}\label{def4.1} \rm
  An entire function $\omega$ of positive
exponential type is said to be a sine type function if
\begin{itemize}
\item[(i)] there is $h>0$ such that all  zeros of $\omega $ lie in the strip
 $\{\lambda \in \mathbb{C}:\left|\operatorname{Im} \lambda \right| <h\}$,

\item[(ii)] there are $h_{1}\in \mathbb{R}$ and positive numbers $m<M$ such that
$m\leq | \omega ( \lambda ) | \leq M$  holds for $\lambda \in \mathbb{C}$ 
with $\operatorname{Im}\lambda =h_{1}$,

\item[(iii)] the type of $\omega$ in the lower half-plane
coincides with that in the upper half-plane.
\end{itemize}
\end{definition}

\begin{lemma}\label{lem4.4} 
Let $(\xi_k)_{k=-\infty}^{\infty}$ be a sequence of complex numbers  
satisfying the following conditions:
\begin{itemize}
\item[(i)]  $\xi_{-k}=-\overline{\xi_k}$ for all not pure imaginary $\xi_k$,
\item[(ii)] the sequence satisfies the asymptotic condition
 \begin{equation} \label{3.10}
\xi_k=\frac{\pi( k-1/2)}{a}+ig +\frac{h}{k}+
\frac{\gamma_k}{k},\quad k\in \mathbb{N},
\end{equation}
where $a,g>0$, $h\in\mathbb{R}$, $(\gamma_k)_{k=1}^{\infty}\in l_2$.
\end{itemize}
Then there exist   $M_1\in \mathbb{R}$, $M_2\in \mathbb{R}$   and 
$\psi\in \mathcal{L}^a$ such that the sequence of  zeros, counted 
with multiplicity, of the entire function
\begin{equation}\label{eq:phi1}
\phi(\lambda)=-\lambda \sin\lambda a+i \lambda \coth(ga)\cos\lambda a
+M_1\cos\lambda a +iM_2\sin\lambda a +\psi(\lambda)
\end{equation}
coincides with the given sequence $(\xi_k)_{k=-\infty}^{\infty}$,
where
\begin{equation}\label{eq:M1M2a}
\frac{1}{\pi (\coth^2(ga)-1)}(\coth(ga) M_2-M_1)=h.
\end{equation}
Furthermore, there exists a number $C>0$ such that
\begin{equation}
\label{4.1a}
\phi(\lambda)=C\prod_{k=-m}^{m}(i(\lambda -\xi _k ))
\prod_{k=m+1}^{\infty}\Big(1-\frac{\lambda}{\xi_k}\Big) 
\Big(1-\frac{\lambda}{\xi_{-k}}\Big) ,
\end{equation}
where we have assumed that the pure imaginary terms (including $0$)
 of the sequence carry the indices $-m,\dots,m$.
\end{lemma}

\begin{proof}
This lemma easily follows from the first part of the proof of 
\cite[Theorem 8.1.4]{MPb}.
Defining
\[
\chi(\lambda)=\tilde C\prod_{\substack{k=-m\\k\not=0}}^{m}(i(\lambda -\xi _k ))
\prod_{k=m+1}^{\infty}\Big(1-\frac{\lambda}{\xi_k}\Big) 
\Big(1-\frac{\lambda}{\xi_{-k}}\Big) 
\]
with $\tilde C\in \mathbb{C}\setminus \{0\}$, this function with $m=1$ 
is of the form \cite[(8.1.20)]{MPb}.
It is easy to see that, for  a suitably chosen $\tilde C$, also in case $m>1$ 
the representation
\begin{equation}\label{eq:chi1a}
\chi(\lambda)=\cos\lambda a+i\alpha \sin\lambda a
+ \frac{\pi h-T\alpha }{\lambda}\sin\lambda a
+\frac{i(T-\pi h\alpha )}{\lambda} \cos\lambda a
+ \frac{\tilde{\psi}(\lambda)}{\lambda}\,,
\end{equation}
in  \cite[(8.1.28)]{MPb} holds, where
 $\alpha = \tanh(ga)$, $T\in \mathbb{R}$ and $\tilde{\psi}\in \mathcal{L}^a$.
Setting
\begin{equation}\label{eq:phi1a}
 \phi (\lambda )=i\coth(ga)(\lambda -\xi _0)\chi (\lambda )
\end{equation}
it is now clear that $\phi $ is of the form \eqref{4.1a} with 
$C\in \mathbb{C}\setminus \{0\}$, and a
straightforward expansion of $\phi$ proves \eqref{eq:phi1}. 
To verify that $C>0$ we observe that the leading term
$i \lambda \coth(ga)\cos\lambda a$ of the representation \eqref{eq:phi1} 
is negative on the positive imaginary axis.
\end{proof}

\begin{lemma}\label{lem4.4a} 
Let $(\xi_k)_{k=-\infty,k\not=0}^{\infty}$ be  a sequence
of complex numbers  satisfying the following conditions:
\begin{itemize}
\item[(i)]  $\xi_{-k}=-\overline{\xi_k}$ for all not pure imaginary $\xi_k$,
\item[(ii)]  the sequence satisfies the asymptotics
 \begin{equation}
\label{3.10a}
\xi_k=\frac{\pi( k-1)}{a}+ig
+\frac{h}{k}+
\frac{\gamma_k}{k},\quad k\in \mathbb{N},
\end{equation}
where $a,g>0$, $h\in\mathbb{R}$, $\{\gamma_k\}_{-\infty}^{\infty}\in l_2$.
\end{itemize}
Then there exist   $M_1\in \mathbb{R}$, $M_2\in \mathbb{R}$   and 
$\psi\in \mathcal{L}^a$ such that the sequence of  zeros, counted with 
multiplicity, of the entire function
\begin{equation}\label{eq:phi2}
\phi(\lambda)=-\lambda\sin\lambda a+i\lambda \tanh(ga)\cos\lambda a
+M_1\cos\lambda a+iM_2\sin\lambda a+\psi(\lambda)
\end{equation}
coincides with the given sequence  $(\xi_k)_{k=-\infty,k\not=0}^{\infty}$,
where
\begin{equation}\label{eq:M1M2b}
\frac{1}{\pi (\tanh^2(ga)-1)}(\tanh(ga) M_2-M_1)=h.
\end{equation}
Furthermore, there exists  $C>0$ such that
\begin{equation}
\label{eq:4.1a}
\phi(\lambda)=C\prod_{k=1}^{m}(i(\lambda -\xi _k))(i(\lambda -\xi _{-k}))
\prod_{k=m+1}^{\infty}\Big(1-\frac{\lambda}{\xi_k}\Big) 
\Big(1-\frac{\lambda}{\xi_{-k}}\Big) ,
\end{equation}
where we have assumed that the pure imaginary terms (including $0$), if any,
 of the sequence carry the indices $\pm1,\dots,\pm m$.
\end{lemma}

\begin{proof}
Define the sequence $(\zeta _k)_{k=-\infty }^\infty $ by 
$\zeta _k=\xi _k+\frac\pi{2a}$ for $k>0$ and
$\zeta _k=\xi _{k-1}+\frac\pi{2a}$ for $k\le0$. For $k\in \mathbb{N}$ we 
therefore have
 \[
\zeta_k=\frac{\pi( k-1/2)}{a}+ig
+\frac{h}{k}+
\frac{\gamma_k}{k},\quad k\not=0,
\]
and, for $k\ge m+1$,
\begin{align*}
\overline{\zeta _{-k}} +\zeta _k&=\overline{\xi _{-k-1}} +\xi _k+\frac\pi a
=-\xi _{k+1}+\xi _k+\frac\pi a\\
&= \frac hk-\frac h{k+1}+\frac {\gamma _k}k-\frac{\gamma _{k+1}}{k+1}
=\frac{\tilde \gamma _k}k,
\end{align*}
where $(\tilde \gamma _k)_{k=m+1}^\infty \in l_2$.
From the proof of \cite[Theorem 8.1.4]{MPb} we see that this condition 
is sufficient to arrive at a slightly  weakened form of the statement of Lemma
\ref{lem4.4}.
 Indeed, the relation \cite[(8.1.23)]{MPb} between the numbers $\zeta _k$, 
$k\not=0$, and the eigenvalues
 of the auxiliary problem \cite[(8.1.21)--(8.1.23)]{MPb}, 
is still valid, and \eqref{eq:chi1a} becomes
 \begin{equation}\label{eq:chi2}
\tilde\chi(\lambda)=\cos\lambda a+i\alpha \sin\lambda a
+ \frac{\pi h-T\alpha }{\lambda}\sin\lambda a
+\frac{i(T-\pi h\alpha )}{\lambda} \cos\lambda a
+ \frac{\tilde{\psi}(\lambda)}{\lambda}\,,
\end{equation}
where $T$ is a complex number which is not necessarily real. 
Then it follows, see \eqref{eq:phi1a}, that the
sequence $(\zeta _k)_{k=-\infty }^\infty$ is the sequence of the zeros 
of the function $\tilde \phi$ defined by
\begin{equation}\label{eq:phi2a}
\tilde \phi (\lambda )=(\lambda -\zeta _0)\tilde\chi (\lambda ).
\end{equation}
  A straightforward calculation shows that
\begin{equation*}
\tilde\phi(\lambda)=\lambda \cos\lambda a+i \lambda \tanh(ga)\sin\lambda a
+\tilde M_1\sin\lambda a -i\tilde M_2\cos\lambda a +\hat\psi(\lambda),
\end{equation*}
where $\tilde M_1$ and $\tilde M_2$ are complex numbers satisfying
\[
\frac{1}{\pi (\tanh^2(ga)-1)}(\tanh(ga) \tilde M_2-\tilde M_1)=h.
\]
By definition of the sequence $(\zeta _k)_{k=-\infty }^\infty $ it follows 
that the sequence
$(\xi_k)_{k=-\infty,k\not=0}^{\infty}$ is the sequence of the zeros of 
the function $\phi $ defined by
$\phi (\lambda )=\tilde \phi (\lambda +\frac\pi{2 a})$.
Therefore the representations \eqref{eq:phi2} and \eqref{eq:M1M2b} 
hold with complex numbers
$M_1$ and $M_2$.   Since $\tilde \chi $ is a sine type function,  \eqref{eq:4.1a} 
follows from
\cite[Lemma  12.2.29]{MPb} with some nonzero complex number $C$. 
The infinite product in \eqref{eq:4.1a} is positive for sufficiently large 
$\lambda $ on the positive imaginary axis, whereas
 in \eqref{eq:phi2} the leading term on the positive imaginary axis,
$-\lambda \sin \lambda a$, is positive there. It follows that $C>0$, and 
then \eqref{eq:4.1a} shows that
$\phi(-\overline{\lambda })=\overline{\phi (\lambda )} $ for all 
$\lambda \in \mathbb{C}$. In particular, $\phi $ is real
on the imaginary axis, which finally gives that $M_1$ and $M_2$ are real.
\end{proof}

\begin{theorem}\label{thm4.5}
Consider a sequence of complex numbers which is indexed as
\begin{itemize}
\item[(i)] $(\xi_k)_{k=-\infty}^{\infty}$, or

\item[(ii)] $(\xi_k)_{-\infty,k\neq 0}^{\infty}$
\end{itemize}
and which satisfies the following conditions:
\begin{itemize}
\item[(1)]  $\xi_{-k}=-\overline{\xi_k}$ for all not pure imaginary $\xi_k$. 
\item[(2)] The numbers $\xi _k$ which lie in the closed lower half-plane and 
on the positive imaginary semiaxis satisfy the properties 
(1)--(6) of Theorem \ref{thmgRHB}, and if $\kappa >0$, then the interval
$(i|\lambda _\kappa |,i\infty )$ contains at least one of the numbers $\xi _k$.
 
\item[(3)] In case (i), the sequence satisfies the asymptotic condition
\begin{equation}
\label{3.17i}
\xi_k=\frac{\pi (k-1/2)}{a}+ig
+\frac{h}{k}+\frac{\gamma_k}{k},\quad k\in \mathbb{N},
\end{equation}
and in case (ii), the sequence satisfies the asymptotic condition
\begin{equation}
\label{3.17ii}
\xi_k=\frac{\pi (k-1)}{a}+ig
+\frac{h}{k}+\frac{\gamma_k}{k},\quad k\in \mathbb{N},
\end{equation}
where $a,g>0$, $h\in\mathbb{R}$ and $(\gamma_k)_{k=1}^{\infty}\in l_2$.
\end{itemize}
Then there exists a unique collection of a real-valued function 
$q\in L_2(0,a)$, $\alpha\in (0,1)\cup(1,\infty)$, $\beta\in \mathbb{R}$
and $\delta\in \mathbb{R}$
such that the spectrum of
problem \eqref{3.1}--\eqref{3.3} coincides with $(\xi_k)$. Here 
$\alpha \in(1,\infty )$ in case (i)
and $\alpha \in (0,1)$ in case (ii).
\end{theorem}


\begin{proof}
We know from Lemmas \ref{lem4.4} and \ref{lem4.4a} that there is an entire 
function $\phi $ with the representation
\begin{equation}\label{eq:phi0}
\phi(\lambda)=-\lambda \sin\lambda a+i \lambda \alpha \cos\lambda a+M_1\cos\lambda a +iM_2\sin\lambda a +\psi(\lambda)
\end{equation}
 whose zeros coincide with the given sequence  $(\xi_k)$, where
\begin{equation}\label{3.18}
\alpha =\coth (ga) \text{ in case (i) and }\alpha =\tanh (ga) \text{ in case (ii),}
\end{equation}
where $M_1$, $M_2$ are real numbers with
\[
\frac{1}{\pi (\alpha ^2-1)}(\alpha M_2-M_1)=h,
\]
and where $\psi\in \mathcal{L}^a$. In view of the assumptions made in part 2, 
we conclude from \cite[Theorem 5.2.16]{MPb} that $\phi $ is of SSHB class. 
Hence there are entire functions $\tilde P$ and $\tilde Q$
which are real on the real axis and which do not have common zeros so that
\[
\phi (\lambda )=\tilde P(\lambda ^2)+i\lambda \tilde Q(\lambda ^2)
\]
and $\frac{\tilde Q}{\tilde P}$ is of class  ${\mathcal N}_+^{\rm ep}$.
We define
\begin{gather}
\label{3.19}
\tilde{\beta}=\frac{\alpha M_1-M_2}{\alpha ^2},\\
\label{3.21}
Q(\lambda)=\frac{1}{\alpha}\tilde Q(\lambda ),\\
\label{3.20}
P(\lambda)=\tilde P(\lambda )-\tilde{\beta} \tilde Q(\lambda).
\end{gather}
From \cite[Lemma 5.1.22]{MPb}  we conclude that
\begin{align*}
\frac QP=\frac 1{\alpha }\,\frac{\tilde Q}{\tilde P-\tilde \beta \tilde Q}
=\frac 1{ \alpha }\,\frac1{\frac {\tilde P}{\tilde Q}-\tilde\beta }
\end{align*}
is a meromorphic Nevanlinna function. Since $\tilde P$ and $\tilde Q$ do not 
have common zeros, also $P$ and $Q$ do not have common zeros, and since zeros 
and poles of meromorphic Nevanlinna functions
are real, simple and interlace, see \cite[Lemma 11.1.3]{MPb}, it follows
 that the zeros of $P$ and $Q$ are real, simple and  interlace.
From \eqref{eq:phi0}, \eqref{3.21} and \eqref{3.20} we obtain
\begin{gather} \label{3.22}
Q(\lambda^2)=\cos\lambda a +\frac{M_2}\alpha \frac{\sin\lambda a}{\lambda}
+\frac{\psi_1(\lambda)}{\lambda}\,,\\
\label{3.23}
P(\lambda^2)=-\lambda\sin\lambda a +\frac{M_2}\alpha \cos\lambda a+\psi_2(\lambda)\,,
\end{gather}
where $\psi_1$ and $\psi _2$ belong to the Paley-Wiener class $\mathcal {L}^a$. 
In \eqref{3.23} we have used that
\[
M_1-\tilde \beta \alpha =\frac{M_2}\alpha 
\]
in view of \eqref{3.19}.

The functions $P$ and $Q$ satisfy the assumptions of \cite[Theorem 2.6]{MP3}.
 Hence there are entire functions $S_0$ and $S_1$ which are real on the real axis,
 whose zeros are real, simple
and interlace, and which are  of the form
\begin{gather}\label{eq:s0a}
S_0(\lambda^2 )= \frac{\sin\lambda a}{\lambda}-A\frac{\cos\lambda a}{\lambda ^2}
+\frac{\hat\psi_3(\lambda)}{\lambda^2}\,,\\
S_1(\lambda^2 )= \cos\lambda a +A \frac{\sin\lambda a}{\lambda}
+\frac{\hat\psi_4(\lambda)}{\lambda }\,,\label{eq:s1a}
\end{gather}
where $A\in \mathbb{R}$ and $\hat\psi_3,\hat\psi_4\in \mathcal L^a$, such that
\begin{equation}\label{eq:PQs}
 QS_1-PS_0=1.
 \end{equation}
In view of \cite[Lemma 12.3.3 and Theorem 12.6.2]{MPb}, there is a unique 
function $q\in L_2(0,a)$
such that $S_0(\lambda ^2)=s(\lambda ,a)$ and $S_1(\lambda ^2)=s'(\lambda ,a)$, 
where $s(\lambda ,\cdot )$ is the solution of the initial value problem \eqref{3.1},
 $s(\lambda ,0)=0$, $s'(\lambda ,0)=1$.

Let $c(\lambda,0 ,\cdot )$ be the solution of the initial value problem \eqref{3.1},
 $c(\lambda,0 ,0)=1$, $c'(\lambda,0 ,0)=0$. By the Lagrange identity, 
the Wronskian
 \begin{equation*}
c(\lambda,0 ,x)s'(\lambda,x )-c'(\lambda,0,x)s(\lambda ,x), \quad 
\lambda \in \mathbb{C},
\end{equation*}
is independent of $x$, and its value at $x=0$ is $1$.
 Together with \eqref{eq:PQs} we conclude that
  \begin{equation}\label{eq:wrdiff}
(Q(\lambda ^2)-c(\lambda,0 ,a))s'(\lambda,a )-(P(\lambda ^2)
-c'(\lambda,0,a))s(\lambda,a )=0, \quad \lambda \in \mathbb{C}.
\end{equation}
Since the zeros of $s(\cdot ,a)$ and $s'(\cdot ,a)$ are distinct, we have 
$Q(\lambda ^2)-c(\lambda,0 ,a)=0$ whenever $s(\lambda ,a)=0$. 
Hence there is an entire function $g$ such that
\[
Q(\lambda ^2)-c(\lambda,0 ,a)=g(\lambda )s(\lambda  ,a).
\]
Because $\lambda \mapsto Q(\lambda ^2)-c(\lambda ,0,a)$ belongs to the
 Payley Wiener class $\mathcal {L}^a$ in view of Proposition \ref{propr}
 and because
$\lambda \mapsto \lambda s(\lambda,a )$ is a sine type function of 
exponential type $a$, it follows that
$\lambda ^{-1}g(\lambda )\to 0$ as $|\lambda |\to \infty $ outside a countable 
set of disjoint discs,
see \cite[Remark 11.2.21 and Lemma 12.2.4]{MPb}, and therefore
$\lambda ^{-1}g(\lambda )=o(\lambda )$ as $|\lambda|\to \infty $.
 It follows by Liouville's theorem that $g$ is constant, say 
$g(\lambda )=\delta $ for some $\delta\in \mathbb{C} $
 and all $\lambda \in \mathbb{C}$. Substitution into \eqref{eq:wrdiff} gives
\begin{gather*}
Q(\lambda ^2)=c(\lambda,0 ,a)+\delta  s(\lambda ,a) 
=c(\lambda,\delta  ,a), \quad \lambda \in \mathbb{C},\\
 P(\lambda ^2)=c'(\lambda ,0,a)+\delta   s'(\lambda ,a)
=c'(\lambda ,\delta ,a), \quad \lambda \in \mathbb{C}.
\end{gather*}
Therefore,
\begin{align*}
\phi (\lambda )
&=\tilde P(\lambda ^2)+i\lambda \tilde Q(\lambda ^2)\\
   &=P(\lambda ^2)+\tilde \beta \alpha Q(\lambda ^2)+i\lambda \alpha Q(\lambda ^2)\\
   &=c'(\lambda,\delta  ,a)+(i\lambda \alpha +\beta )c(\lambda,\delta  ,a),
\end{align*}
where $\beta =\tilde \beta \alpha $.

Thus we have shown that $\phi $ is the characteristic function of a problem of 
the form \eqref{3.1}--\eqref{3.3}. Next
we are going to prove the uniqueness of this problem.

From the spectral asymptotics \eqref{3.11}, \eqref{3.12} it follows that 
$\alpha $ is uniquely determined by \eqref{3.17i}
and \eqref{3.17ii}. In view of \eqref{3.9}, the numbers
\begin{equation}\label{eq:delta}
B=\delta+\frac 12\int_0^aq(x)\,dx
\end{equation}
and $B+\beta $, and therefore $\beta $, are uniquely determined.

Finally, let $q_1\in L_2(0,a)$ and $\delta _1\in \mathbb{R}$ be arbitrary such 
that the eigenvalue problem
\eqref{3.1}--\eqref{3.3} with $q_1$ and $\delta _1$ has the given sequence 
$(\xi _k)$ as eigenvalues. Then its characteristic function is
\begin{equation*}
\phi_1(\lambda,\delta_1):=c_1'(\lambda,\delta_1,a)+(i\alpha\lambda
+\beta)c_1(\lambda,\delta_1,a),
\end{equation*}
see \eqref{3.8}, where $c_1(\lambda ,\delta _1,\cdot )$ is the solution of 
\eqref{3.1} with potential $q_1$ and with the initial conditions 
$c_1(\lambda ,\delta _1,0)=1$, $c_1'(\lambda ,\delta _1,0)=\delta _1$.
The functions $\lambda \mapsto (\lambda -\xi _1)^{-1}\phi _1(\lambda ,\delta _1)$ 
and $\lambda \mapsto (\lambda -\xi _1)^{-1}\phi (\lambda )$  are sine type 
functions whose zeros coincide and
which have the same leading terms. Therefore they are equal, see 
\cite[Lemma 11.2.29]{MPb}, and hence $\phi _1(\cdot ,\delta _1)=\phi $.
Since $c$, $c_1$, $c'$, $c_1'$ are even functions in the variable $\lambda $, 
it follows that
\[ 
c_1(\lambda,\delta_1,a)=c(\lambda,\delta,a),\quad
 c_1'(\lambda,\delta_1,a)=c'(\lambda,\delta,a).
\]
Because the Wronskian of \eqref{3.1} is constant, \eqref{eq:r} shows that
\begin{gather*}
c_1(\lambda,\delta_1 ,a)s_1'(\lambda,a )-c_1'(\lambda,\delta_1 ,a)s(\lambda ,a)=1,\\
c(\lambda,\delta ,a)s'(\lambda,a )-c'(\lambda,\delta ,a)s(\lambda ,a)=1,
\end{gather*}
 and therefore
\[
c(\lambda,\delta ,a)(s_1'(\lambda,a )-s'(\lambda,a ))
-c'(\lambda,\delta ,a)(s_1(\lambda,a )-s(\lambda ,a))=0.
\]
Since $c(\cdot ,\delta ,a)$ and $c'(\cdot ,\delta ,a)$ do not have common zeros,
\[
\chi(\lambda ):=\frac{s_1(\lambda,a )-s(\lambda ,a)}{c(\lambda,\delta ,a)}
=\frac{s_1'(\lambda,a )-s'(\lambda,a )}{c'(\lambda,\delta ,a)}
\]
defines an entire function $\chi$. But $c(\cdot ,\delta ,a)$ is a sine 
type function of type $a$, whereas
$s_1(\cdot ,a)-s(\cdot ,a)\in \mathcal{L}^a$. 
Therefore $\chi=0$, see \cite[Remark 2.3]{MP3}, i.\,e., 
$s_1(\lambda ,a)=s(\lambda ,a)$ and
 $s_1'(\lambda ,a)=s'(\lambda ,a)$. The uniqueness statement in 
\cite[Theorem 12.6.2]{MPb} now shows that $q_1=q$, and finally
 $\delta _1=\delta $ by \eqref{eq:delta}.
\end{proof}

\begin{corollary}\label{coro4.7} 
Consider a sequence of complex numbers which is indexed as
\begin{itemize}
\item[(i)] $(\xi_k)_{k=-\infty}^{\infty}$, or
\item[(ii)] $(\xi_k)_{-\infty,k\neq 0}^{\infty}$
\end{itemize}
and which satisfies the following conditions:
\begin{itemize}
\item[(1)]  $\xi_{-k}=-\overline{\xi_k}$ for all not pure imaginary $\xi_k$. 
\item[(2)] All numbers numbers $\xi _k$ lie in the open upper half-plane.
\item[(3)] In case (i), the sequence satisfies the asymptotic condition
\begin{equation}
\label{3.17ib}
\xi_k=\frac{\pi (k-1/2)}{a}+ig
+\frac{h}{k}+\frac{\gamma_k}{k},\quad k\in \mathbb{N},
\end{equation}
and in case (ii), the sequence satisfies the asymptotics
\begin{equation}
\label{3.17iib}
\xi_k=\frac{\pi (k-1)}{a}+ig
+\frac{h}{k}+\frac{\gamma_k}{k},\quad k\in \mathbb{N},
\end{equation}
where $a,g>0$, $h\in\mathbb{R}$ and $(\gamma_k)_{k=1}^{\infty}\in l_2$.
\end{itemize}
Then for each $l>0$ there exists a unique collection of a real-valued function 
$\rho\gg0$ in $W_2^2(0,a)$ and positive real numbers $p$ and $\nu $
such that the spectrum of
problem \eqref{2.4}--\eqref{2.6} coincides with $(\xi_k)$.
\end{corollary}

\begin{proof}
The assumptions of Theorem \ref{thm4.5} are satisfied, and hence there are
unique real numbers $a>0$,
$\alpha \in (0,1)\cup (1,\infty )$, $\delta $ and $\beta $ and a unique real 
function $q\in L_2(0,a)$ such that $(\xi _k)$ is the spectrum of the 
corresponding eigenvalue problem \eqref{3.1}--\eqref{3.3}. 
As in the proof of \cite[Therorem 8.3.4]{MPb} we can construct a unique function 
$\rho \in W_2^2(0,l)$ which satisfies \eqref{2.15}--\eqref{2.17}
and which is positive on $(0,l]$. The only necessary modification in the 
reasoning is that the Robin boundary condition does not give $\rho (0)\neq 0$ 
directly. Rather,
the reference function $\psi$ considered in the proof of \cite[Therorem 8.3.4]{MPb} 
has to satisfy the Robin condition $\psi'(0)-\delta \psi(0)=0$ here, and this 
gives $\psi(0)\neq 0$ because $\psi\not=0$. Then the Sturm comparison theorem,
see \cite[Theorem 13.1]{Wei} and the continuity of the Pr\"ufer angle show that 
also $\rho (0)\not=0$.

The numbers $\nu >0$ and $p\in \mathbb{R}$ are uniquely determined by \eqref{2.18} 
and \eqref{2.19}. To show that $p>0$ assume that $p\le0$. 
Since $(\lambda _k)$ is the spectrum of the operator pencil $L$ defined 
in \eqref{eq:2.7a}, it follows from \cite[Theorem 1.3]{MPb} that $A\ge 0$ and 
therefore $A\gg0$ since the invertibility of $L(0)$ implies that
$A$ is invertible. Now choose a real function $v\in W_2^2(0,l)$ and 
$\eta \in (0,\frac l2)$ such that $v(0)=2$, $v'(0)=pv(0)$,
$|v'(x)| \le 2|p|$ for $x\in [0,\eta  ]$, $v'(x)=0$ for $x\in(\eta ,l)$ and 
$v(l)\ge 1$.
For example, for $x\in [0,\eta ]$ we can take $v(x)=2+2\frac p\theta \sin (\theta x)$
 with $\theta =\frac{3\pi}{2\eta}$.
Then $Y=(v,v(l))^{\rm T}$ satisfies
$\|Y\|^2\ge1$  and, see \eqref{eq:Apos},
\[
(AY,Y)=\int_0^a|v'(x)|^2\,dx+p|v(0)|^2\le 4\eta |p|^2+4p\le 4\eta |p|^2.
\]
Since $\eta \in (0,\frac l2)$ is arbitrary, this contradicts $A\gg 0$. 
Hence $p>0$ follows.
\end{proof}

\section{Inverse problems by parts of spectra}

Consider the two problems
\begin{gather}
\label{4.1} -y''+q(x)y=\lambda^2y, \quad x\in[0,a],\\
\label{4.2}
 y'(0)-\delta_j y(0)=0,  \quad j=1,2,\\
\label{4.3}
y'(a)+(i\alpha\lambda+\beta) y(a)=0,
\end{gather}
with $\delta_j\in \mathbb{R}$. The corresponding characteristic functions are
\begin{equation}\label{4.4}
\phi(\lambda,\delta_j)=c'(\lambda,\delta_j,a)+(i\alpha\lambda+\beta)
c(\lambda,\delta_j,a),
\end{equation}
see \eqref{3.8}.


\begin{proposition}\label{prop5.1} 
For all $\lambda\in\mathbb{C}$ and all $\delta _1,\delta _2\in\mathbb{R}$
 the following identity holds
\begin{equation}\label{4.5}
\phi(\lambda,\delta_1)\phi(-\lambda,\delta_2)
-\phi(-\lambda,\delta_1)\phi(\lambda,\delta_2)
=2i\alpha\lambda(\delta_2-\delta_1).
\end{equation}
\end{proposition}

\begin{proof}
 Using \eqref{4.4} and taking into account that $c(\lambda,\delta_j,a)$ and 
$c'(\lambda,\delta_j,a)$ are even functions of $\lambda$ we obtain
\begin{equation} \label{4.6}
\begin{aligned}
&\phi(\lambda,\delta_1)\phi(-\lambda,\delta_2)-\phi(-\lambda,\delta_1)
\phi(\lambda,\delta_2)\\
&= 2i\alpha\lambda [c(\lambda, \delta_1,a)c'(\lambda,\delta_2,a)
 -c'(\lambda,\delta_1,a)c(\lambda,\delta_2,a)].
\end{aligned}
\end{equation}
Now using the  Lagrange identity and the initial conditions 
$c(\lambda ,\delta _j,0)=1$, and $c'(\lambda ,\delta _j,0)=1$,
 see Section \ref{s:inv}, we obtain
\begin{align*}
&c(\lambda, \delta_1,a)c'(\lambda,\delta_2,a)-c'(\lambda,\delta_1,a)
 c(\lambda,\delta_2,a) \\
&=c(\lambda, \delta_1,0)c'(\lambda,\delta_2,0)-c'(\lambda,\delta_1,0)
c(\lambda,\delta_2,a)\\&=\delta_2-\delta_1,
\end{align*}
 which proves \eqref{4.5}.
 \end{proof}

\begin{definition}\label{defPhia} \rm
Let  $\alpha >0$.

\noindent(1) For $M_1,M_2\in \mathbb{R}$ let $\Phi _{\alpha,M_1,M_2}$ 
be  the class of entire functions $\phi $  which satisfy the following conditions:
\begin{itemize}
\item[(i)] if $\lambda \in \mathbb{C}\setminus \{0\}$ and $\phi (\lambda ) =0$, 
then $\phi (-\lambda )\not=0$, 
\item[(ii)] $0$ is at most a simple zero of  $\phi$,
\item[(iii)] $\phi(-\overline{\lambda})=\overline{\phi(\lambda)}$ for all
 $\lambda \in \mathbb{C}$,
\item[(iv)] there is $\tau\in \mathcal{L}^a$ such that
\begin{equation}\label{4.7}
\phi(\lambda)=-\lambda\sin\lambda a+i\lambda \alpha\cos\lambda a
+M_1\cos\lambda a+iM_2\sin\lambda a+\tau(\lambda).
\end{equation}
\end{itemize}

\noindent(2) Define
\[
\Phi _\alpha =\cup_{M_1,M_2\in \mathbb{R}}\Phi _{\alpha ,M_1,M_2}.
\]
\end{definition}

\begin{theorem}\label{thm5.3}
Let $\kappa \in \mathbb{N}$, $N_2=\{\pm 1, \pm 2, \dots,\pm \kappa\}$,
$N_1=\{0,\pm (\kappa +1),\pm (\kappa +2),\dots\}$.
 Let   $(\lambda_k)_{k=-\infty}^{\infty}$  be a sequence of complex numbers such that
\begin{itemize}
\item[(1)] $\operatorname{Im} \lambda_k>0$ for all $k\in \mathbb{Z}$,

\item[(2)] $\lambda_{-k}=-\overline{\lambda_k}$ for all not pure imaginary 
$\lambda_{-k}$,

\item[(3)] $\lambda_k\not=\lambda_j$  for all $k\in N_1$ and $j\in N_2$.
\end{itemize}
Let the sequence  $(\lambda_k)$  satisfy \eqref{3.10}
where  $g$ and $h$ are real constants, $(\beta_k)\in l_2$.

 Then for any given $\Delta<0$ there exist a real potential $q\in L_2(0,a)$, 
real numbers $\alpha>1$, $\beta$  and $\delta_1$  such that  
$(\lambda_{k})_{ k\in N_1}$  are eigenvalues  of
problem \eqref{4.1}--\eqref{4.3} with $j=1$, and  $(\lambda_{k})_{ k\in N_2}$  
are eigenvalues  of
problem \eqref{4.1}--\eqref{4.3} with  $j=2$ where $\delta _2=\delta _1+\Delta$.
\end{theorem}

\begin{proof}
The sequence $(\lambda _k)_{k=-\infty }^\infty $ satisfies the assumptions of 
Lemma \ref{lem4.4}, and therefore the zeros of the entire function $\phi $ 
defined in \eqref{eq:phi1} coincide with the given sequence $(\lambda _k)$. 
With $\alpha =\coth(ga)>1$ it follows that $\phi \in \Phi _\alpha $. 
In view of  \cite[Theorems 3.5 and 3.8]{MPb} there are
functions $\tilde \phi $ and $\tilde X\in \Phi _\alpha $ such that
\begin{equation}\label{eq:phiXtilde}
\tilde \phi (\lambda )\tilde X(-\lambda )-\tilde \phi (-\lambda )\tilde X(\lambda )
=2i\alpha \lambda \Delta , \quad \lambda \in \mathbb{C},
\end{equation}
such that the $\lambda _k$ with $k\in N_1$ are zeros of $\tilde \phi $, 
the $\lambda _k$ with $k\in N_2$ are zeros of $\tilde X$,
and  such that all zeros of $\tilde \phi $ lie in the closed upper half-plane.  
From the proof of Theorem \ref{thm3.3} we know that the sequence $(\xi _k) $ of the
zeros of $\tilde \phi $ can be indexed in such a way that they satisfy asymptotic
condition of the form \eqref{3.17i}, and
from the definition of $\Phi _\alpha $ it follows that $\xi _{-k}=-\overline{\xi _k}$ 
for not pure imaginary $\xi _k$ with a suitable indexing. Hence $(\xi _k)$ 
satisfies conditions (3) and (1) of Theorem  \ref{thm4.5}. From the definition of
$\Phi _\alpha $ we know that $0$ can be at most a simple zero of $\tilde \phi $ 
and that $\tilde \phi $ has no nonzero real zeros. Therefore
$(\xi _k)$ also satisfies condition (2) of Theorem \ref{thm4.5}, and hence
Theorem \ref{thm4.5} shows that there are a real function $q\in L_2(0,a)$, 
$\beta \in \mathbb{R}$ and $\delta _1\in \mathbb{R}$ such that the spectrum of
\eqref{4.1}--\eqref{4.3} for $j=1$ coincides with $(\xi _k)$.
 With $\delta _2=\delta _1+\Delta $, Proposition
\ref{prop5.1} shows that the characteristic function $\phi (\cdot ,\delta _2)$ 
of \eqref{4.1}--\eqref{4.3} with $j=2$ satisfies
\[
\tilde \phi (\lambda )\phi (-\lambda ,\delta _2)
-\tilde \phi (-\lambda )\phi (\lambda ,\delta _2)=2ia\lambda (\delta _2-\delta _1)
=2i\alpha \lambda \Delta .
\]
In view of \cite[Theorem 3.2]{MPb}, for any given $\tilde \phi \in \Phi _\alpha $, 
this equation has a unique solution in $\Phi _\alpha $, and therefore 
$\tilde X=\phi (\cdot ,\delta _2)$.
\end{proof}

\begin{theorem}\label{thm5.4}
Let $\kappa \in \mathbb{N}$, $N_2=\pm 1, \pm 2, \dots,\pm \kappa$,
$N_1=\{\pm (\kappa +1),\pm (\kappa +2),\dots\}$.
 Let   $(\lambda_k)_{k-\infty,k\not=0}^{\infty}$  be a sequence of complex
 numbers such that
\begin{itemize}
\item[(1)] $\operatorname{Im} \lambda_k>0$ for all $k\in \mathbb{Z}\setminus \{0\}$,

\item[(2)] $\lambda_{-k}=-\overline{\lambda_k}$ for all not pure imaginary 
$\lambda_{-k}$,

\item[(3)] $\lambda_k\not=\lambda_j$  for all $k\in N_1$ and $j\in N_2$.
\end{itemize}
Let the sequence  $(\lambda_k)$  satisfy \eqref{3.10a}
where  $g$ and $h$ are real constants, $(\beta_k)\in l_2$.

 Then for any given $\Delta\not=0$ there exist a real potential $q\in L_2(0,a)$, 
real numbers $\alpha\in(0,1)$, $\beta$  and $\delta_1$  such that  
$(\lambda_{k})_{ k\in N_1}$  are eigenvalues  of
problem \eqref{4.1}--\eqref{4.3} with $j=1$, and  $(\lambda_{k})_{ k\in N_2}$  
are eigenvalues  of
problem \eqref{4.1}--\eqref{4.3} with $j=2$ where $\delta _2=\delta _1+\Delta$.
\end{theorem}

\begin{proof}
The sequence $(\lambda _k)_{k=-\infty ,k\not=0}^\infty $ satisfies the assumptions
 of Lemma \ref{lem4.4a}, and therefore the zeros of the entire function 
$\phi $ defined in \eqref{eq:phi2} coincide with the given sequence 
$(\lambda _k)$. With $\alpha =\tanh(ga)\in(0,1)$ it follows that 
$\phi \in \Phi _\alpha $.
 If $\Delta <0$, we can proceed as in the proof of Theorem \ref{thm5.3}.

  If $\Delta >0$, it follows from
   \cite[Theorems 3.5 and 3.8]{MPb} that there are
functions $\tilde \phi $ and $\tilde X\in \Phi _\alpha $ such that 
\eqref{eq:phiXtilde} holds, such that the $\lambda _k$ with $k\in N_1$ 
are zeros of $\tilde \phi $, the $\lambda _k$ with $k\in N_2$ are zeros of
$\tilde X$, and  such that all zeros of $\tilde X$ lie in the closed upper 
half-plane. Proceeding as in the proof of Theorem \ref{thm5.3} with
$\tilde \phi $ and $\tilde X$ interchanged completes this proof.
\end{proof}

\begin{remark} \rm
(1) For $\rho \in W_2^2(0,l)$ it was shown in Section 2 that the string 
problem \eqref{2.4}--\eqref{2.6} can be transformed into a Robin-Regge 
eigenvalue problem \eqref{3.1}--\eqref{3.3}. Hence the eigenvalues of  
\eqref{2.4}--\eqref{2.6} lie in the open upper half-plane by
 Proposition \ref{prop2.2} 
and satisfy the asymptotics \eqref{3.11} or \eqref{3.12} if
$\rho ^{-\frac12}(s(a))\nu $ is less than $1$ or larger than $1$, respectively.

 (2) If in the proofs of Theorems \ref{thm5.3} and \ref{thm5.4} both
 $\tilde \phi $ and $\tilde X$ have no zeros in the open lower half-plane, 
then the potentials $q$ and the numbers $\alpha $, $\beta $, 
$\delta _1$ and $\delta _2$ in Theorems
\ref{thm5.3} and \ref{thm5.4} are unique, see \cite[Corollary 3.7]{MPb}.
\end{remark}

\subsection*{Acknowledgments}
M.  M\"oller was supported by a grant from the National Research
Foundation of South Africa (Grant No. 80956).



\begin{thebibliography}{99}

\bibitem{B}  Borg, G.;
 Eine Umkehrung der Sturm-Liouvilleschen
Eigenwertaufgabe.  \emph{Acta Math.},  \textbf{78} (1946), 1--96/

\bibitem{CH} Courant, R.;  Hilbert, D.;
\emph{Methods of Mathematical Physics  I}, Interscience Publishers,  New York, 1953

\bibitem{GRS} del Rio, R.; Gesztesy, F.; Simon,  B.;
Inverse spectral analysis with partial information on the
potential. III. Updating  boundary conditions. 
\textit{Internat. Math. Res. Notices}, \textbf{1997}(15) (1997), 751--758.

\bibitem{GS}  Gesztesy, F.; Simon, B.;
 Inverse spectral analysis with partial information on the
potential. II. The case of discrete  spectrum, 
\emph{Trans. Amer. Math. Soc.}, \textbf{352}  (2000), 2765--2789.

\bibitem{GA} Gustafson, K.;  Abe, T.;
 The third boundary condition--was it Robin's? 
\emph{Math. Intelligencer}, \textbf{20} (1998), 63-71.

\bibitem{Ha}  Hald, O. H.;
 Inverse eigenvalue problem for the mantle, \emph{Geophys. J. R. Astr. Soc.},
\textbf {62} (1980),  41--48.

\bibitem{HL}  Hochstadt, H.;  Lieberman, B.;
An inverse Sturm-Liouville problem with mixed given data,
 \emph{SIAM J. Appl.Math.}, \textbf {34} (1978),  676--680.

\bibitem{H1} Horvath, M.;
 On the inverse spectral theory of Schr\"odinger and Dirac operators.
\emph{Trans. Amer. Math. Soc.}, {\bf 353}(10) (2001),  4155--4171.

\bibitem{H2} Horvath, M.;
 Inverse spectral problems and closed exponential systems.
\emph{Ann. of Math. (2)},  {\bf 162} (2005), 885--918 

\bibitem{KN1} Krein, M. G.; Nudelman. A.  A.;
Direct and inverse problems for frequencies of boundary dissipation 
of a nonuniform string (Russian),
\emph{Dokl. Akad. Nauk SSSR}, \textbf{247}(5) (1979), 1046--1049.

\bibitem{KN2} Krein, M. G., Nudelman, A.  A.;
 Some spectral properties of a nonhomogeneous string with a dissipative
 boundary condition (Russian),
\emph{J. Operator Theory} \textbf{22} (1989), 369--395.

\bibitem{L} Levinson, N.;
 The inverse Sturm-Liouville problem, \emph{Mat. Tidsskr B}, \textbf{1949}
 (1949), 25--30.

\bibitem{LG} Levitan, B. M.;  Gasymov, M. G.;
 Determination of a differential equation by two  spectra (Russian),
\emph{Uspehi Mat. Nauk},
\textbf{19}(2)(116) (1964), 3--63.

\bibitem{M}  Marchenko, V. A.;
\emph{Sturm-Liouville Operators and Applications},  OT \textbf {22},
Birkh\"auser, Basel,  1986.

\bibitem{MP}  Martinyuk, O.;  Pivovarchik, V.;
 On the Hochstadt-Lieberman theorem,  \emph{Inverse problems}, \textbf{26}(3)
(2010),  035011: 6 pp. 

\bibitem{MM} Mennicken, R.; M\"oller, M.;
Non-self-adjoint Boundary Eigenvalue Problems,
North-Holland Mathematics Studies \textbf{192},   Elsevier,  Amsterdam, 2003

\bibitem{MPb}  M\" oller, M.; Pivovarchik, V.;
\emph{Spectral Theory of Operator Pencils, Hermite-Biehler Functions, and 
Their Applications}, OT \textbf {246}, Birkh\"auser, Cham,  2015

\bibitem{MP3} M\"oller, M.; Pivovarchik, V.;
 Functional equations related  to sine type functions, 
\emph{Complex Anal. Oper. Theory},  (2016) DOI 10.1007/s11785-016-0563-2.

\bibitem{Pi4} Pivovarchik, V.;
  On the Hald-Gesztesy-Simon theorem, \emph{Integral Equations Operator Theory}, 
{\bf  73}(3) (2012),  383--393.

\bibitem{PvM}  Pivovarchik, V.;  van der Mee, C.;
  The inverse generalized Regge problem,
\emph{Inverse Problems}, \textbf{17}(6) (2001), 1831--1845.

 \bibitem{Reg} Regge, T.;
Construction of potentials from resonance parameters,
\emph{Nuovo Cimento, X. Ser.}, \textbf{9} (1958), 491--503. 

\bibitem{Su}  Suzuki, T.;
Inverse problems for heat equations on compact intervals and on circles.
 I, \emph{ J. Math. Soc. Japan}, \textbf{38}(1) (1986),  39--65.

\bibitem{WX}  Wei, G.; Xu, H.-K.;
On the missing eigenvalue problem for an inverse Sturm-Liouville problem, 
\emph{J.  Math. Pures. Appl.},  \textbf{91} (2009),  468--475.

\bibitem{Wei} Weidmann, J.;
\emph{Spectral Theory of Ordinary Differential Operators},
 Lecture Notes in Mathematics \textbf{1258},  Springer-Verlag, Berlin, 1987

\end{thebibliography}

\end{document}

