On the spectrum of the pencil of high order differential operators with almost periodic coefficients

Abstract

In this paper, the spectrum and the resolvent of the operator $L_{\lambda}$ which is generated by the differential expression $\ell_{\lambda }(y)=y^{(m)}+\sum_{\gamma=1}^{m} ( \sum_{k=0}^{\gamma }\lambda^{k}p_{\gamma k}(x) ) y^{(m-\gamma)}$ has been investigated in the space $L_{2}(\mathbb{R})$. Here the coefficients $p_{\gamma k}(x)=\sum_{n=1}^{\infty }p_{\gamma kn}e^{i\alpha_{n}x}$, $k=0,1,\ldots,\gamma-1$; $p_{\gamma\gamma }(x)=p_{\gamma\gamma}$, $\gamma=1,2,\ldots,m$, are constants, $p_{mm}\neq0$ and $p_{\gamma k}^{(\nu)}(x)$, $\nu=0,1,2,\ldots,m-\gamma$, are Bohr almost-periodic functions whose Fourier series are absolutely convergent. The sequence of Fourier exponents of coefficients (these are positive) has a unique limit point at +∞. It has been shown that if the polynomial $\phi (z)=z^{m}+p_{11}z^{m-1}+p_{22}z^{m-2}+\cdots+p_{m-1,m-1}z+p_{mm}$ has the simple roots $\omega_{1},\omega_{2},\ldots,\omega_{m}$ (or one multiple root $\omega_{0}$), then the spectrum of operator $L_{\lambda}$ is pure continuous and consists of lines $\operatorname{Re}(\lambda\omega _{k})=0$, $k=1,2,\ldots,m$ (or of line $\operatorname{Re}(\lambda\omega_{0})=0$). Moreover, a countable set of spectral singularities on the continuous spectrum can exist which coincides with numbers of the form $\lambda=0$, $\lambda _{sjn}=i\alpha_{n} ( \omega_{j}-\omega_{s} ) ^{-1}$, $n\in \mathbb{N}$, $s,j=1,2,\ldots,m$, $j\neq s$. If $\phi(z)=(z-\omega_{0})^{m}$, then the spectral singularity does not exist. The resolvent $L_{\lambda}^{-1}$ is an integral operator in $L_{2}(\mathbb{R})$ with the kernel of Karleman type for any $\lambda\in\rho (L_{\lambda})$.

1 Introduction

In this study, we investigate the spectrum and the resolvent of the maximal differential operator $L_{\lambda}$ which is generated by the linear differential expression

$$ \ell_{\lambda}(y)=y^{(m)}+\sum_{\gamma=1}^{m}p_{\gamma}(x, \lambda )y^{(m-\gamma)} $$

in the space $L_{2}(\mathbb{R})$, where λ is a complex parameter,

$$ \begin{aligned} &p_{\gamma}(x,\lambda)=\sum_{k=0}^{\gamma} \lambda^{k}p_{\gamma k}(x), \qquad p_{\gamma\gamma}(x)=p_{\gamma\gamma}, \\ &p_{\gamma k}(x)=\sum_{n=1}^{\infty}p_{\gamma kn}e^{i\alpha_{n}x}, \quad \gamma =1,2,\ldots,m,\ k=0,1,\ldots,\gamma-1, \end{aligned} $$

(1)

with $p_{\gamma\gamma},p_{\gamma kn}\in \mathbb{C}$, $p_{mm}\neq0$ and the condition

$$ \sum_{\gamma=1}^{m}\sum _{k=0}^{\gamma-1}\sum_{n=1}^{\infty} \alpha _{n}^{m-\gamma} \vert p_{\gamma kn}\vert < + \infty $$

(2)

is satisfied. Here $(\alpha_{n})_{n=1}^{\infty}$ is an increasing sequence of positive numbers with $\alpha_{n}\rightarrow+\infty$ and the set $\{ \alpha_{n}:n\in \mathbb{N} \}$ is an additive semigroup.

The operator $L_{\lambda}$ is defined in the domain

$$\begin{aligned} D(L_{\lambda}) =&\bigl\{ y(x)|y(x),y^{\prime }(x), \ldots,y^{(m-1)}(x)\in \operatorname{AC}[a,b]\text{ for all }a,b\in \mathbb{R}, \\ & y(x),\ell_{\lambda}(y)\in L_{2}( \mathbb{R})\bigr\} . \end{aligned}$$

If at least one of the functions $p_{\gamma k}(x)$, $\gamma =1,2,\ldots,m$, $k=0,1,\ldots,\gamma-1$, is not zero, then the operator $L_{\lambda}$ is non-self-adjoint for each $\lambda\in\mathbb{C}$.

Let $\mathit{AP}^{+}$ be a class of Bohr almost-periodic functions $q(x)=\sum_{n=1}^{\infty}q_{n}e^{i\alpha_{n}x}$, where $\Vert q\Vert =\sum_{n=1}^{\infty} \vert q_{n}\vert <+\infty$. In the case $\alpha_{n}=n$, $n\in \mathbb{N}$, we denote this class by $Q^{+}$. It is clear that $\mathit{AP}^{+}$ is a normed space and (2) means that $p_{\gamma k}^{(\nu)}(x)\in \mathit{AP}^{+}$ for $\gamma =1,2,\ldots,m$, $k=0,1,\ldots,\gamma-1$, $\nu=0,1,\ldots,m-\gamma$.

Under the assumed conditions, coefficients $p_{\gamma}(x,\lambda)$ can be represented as

$$\begin{aligned} p_{\gamma}(x,\lambda) =&p_{\gamma\gamma}\lambda^{\gamma }+\sum _{k=0}^{\gamma-1}\lambda^{k}\sum _{n=1}^{\infty }p_{\gamma kn}e^{i\alpha_{n}x}=p_{\gamma\gamma} \lambda^{\gamma }+\sum_{n=1}^{\infty} \Biggl( \sum_{k=0}^{\gamma-1}\lambda ^{k}p_{\gamma kn} \Biggr) e^{i\alpha_{n}x} \\ =&p_{\gamma\gamma}\lambda^{\gamma}+\sum_{n=1}^{\infty} \tilde{p}_{\gamma n}(\lambda)e^{i\alpha_{n}x},\quad \gamma=1,2, \ldots,m. \end{aligned}$$

Here, $\tilde{p}_{\gamma n}(\lambda)$ is an algebraic polynomial whose degree does not exceed $\gamma-1$. Moreover, according to (2) the series $\sum_{n=1}^{\infty} \vert \tilde{p}_{\gamma n}{ (\lambda)}\vert$ is majorized in every compact set $S\subseteq \mathbb{C}$, i.e. for $\tilde{p}_{\gamma n}=\sup_{\lambda\in S}\vert \tilde{p}_{\gamma n}{ (\lambda)}\vert$, $\gamma =1,2,\ldots,m$, $n\in \mathbb{N}$, the series $\sum_{n=1}^{\infty}\tilde{p}_{\gamma n}$ converges.

Let $\omega_{k}$, $k=1,2,\ldots,m$, denote the roots of the characteristic polynomial

$$ \phi(z)=z^{m}+p_{11}z^{m-1}+p_{22}z^{m-2}+ \cdots+p_{m-1,m-1}z+p_{mm} $$

corresponding to the linear differential expression $\ell_{\lambda}(y)$, numbered as $0\leq\arg\omega_{1}\leq\arg\omega_{2}\leq\cdots\leq \arg \omega_{m}<2\pi$ and $\lambda_{sjn}=i\alpha_{n} ( \omega _{j}-\omega_{s} ) ^{-1}$ for $n\in \mathbb{N}$, $j,s=1,2,\ldots,m$, $j\neq s$.

Let $\Lambda= \{ \lambda_{sjn}:s,j=1,2,\ldots,m,\ j\neq s, n\in \mathbb{N} \}$ and $\Lambda_{0}=\Lambda\cup \{ 0 \}$. It is obvious that the roots are different from zero according to the condition $p_{mm}\neq0$. Below we shall assume that these roots are different or all coincide and any three of these roots are not on the same line in the complex plane. Under these conditions, for each constant s the numbers $\lambda_{sjn}=i\alpha_{n} ( \omega_{j}-\omega_{s} ) ^{-1}$, $n\in \mathbb{N}$, $j=1,2,\ldots,m$, $j\neq s$, are located on the $m-1$ rays from the origin. Moreover, since $\lambda_{sjn}=-\lambda_{jsn}$, the set Λ is symmetric with respect to the origin.

The lines $l_{k}= \{ \lambda:\lambda\in \mathbb{C}, \operatorname{Re}(\lambda\omega_{k})=0 \}$, $k=1,2,\ldots,m$, divide the complex λ-plane into $2m_{0}$ open sectors $S_{k}$, $k=1,2,\ldots,2m_{0}$. Let us assume that beginning from sector $S_{1}$ whose closure contains positive numbers, these sectors are numbered $S_{1},S_{2},\ldots,S_{2m_{0}}$ ($m_{0}\leq m$) counterclockwise successively. It is clear that if there are different roots $\omega_{k}$, $\omega_{j}$ such that $\omega_{k}/\omega_{j}\in \mathbb{R}$, then the lines $l_{k}$ and $l_{j}$ coincide. Therefore, the number of sectors $S_{k}$ may be less than 2m. In the case $\phi(z)=(z-\omega _{0})^{m}$ the line $l_{0}= \{ \lambda:\lambda\in \mathbb{C}, \operatorname{Re}\lambda\omega_{0}=0 \}$ divides the complex λ-plane into two half planes $S^{+}= \{ \lambda:\lambda\in \mathbb{C}, \operatorname{Re}\lambda\omega_{0}>0 \}$, $S^{-}= \{ \lambda :\lambda\in\mathbb{C}, \operatorname{Re}\lambda\omega_{0}<0 \}$. In the sequel we shall see that the resolvent set of the operator $L_{\lambda }$ consists of the above defined sectors.

The interest of the investigation of the spectral properties of the differential operators with coefficients belonging to class $\mathit{AP}^{+}$ has been increased after the study of [1]. In [1] the differential operator $L(y)=-y^{\prime\prime}+q(x)y$ with periodic potential $q(x)\in Q^{+}$ (case $m=2$) has been investigated in the space $L_{2}( \mathbb{R})$. In this study the spectral data $\{s_{n}\}_{n\in \mathbb{N}}$ has been determined and sufficient conditions have been obtained for solvability of the inverse problem according to the spectral data. Afterward, in [2] were found the necessary and sufficient conditions for a set $\{s_{n}\}_{n\in\mathbb{N}}$ to be the spectral data of the operator $L(y)=-y^{\prime\prime }+q(x)y$ with periodic potential $q(x)\in L_{2} ( 0,2\pi )$. In [3], the results of [2] are generalized for almost-periodic potential $q(x)$ having only positive Fourier exponents. The spectral properties of ordinary differential operators of high order with coefficients from $\mathit{AP}^{+}$ have been investigated in [4, 5]. The spectrum and the resolvent of the bundle of 2n order and second order differential operators with coefficients from $\mathit{AP}^{+}$ and from $Q^{+}$ have been examined in [6–9] and in [10], respectively. In all of these studies the examined operators in the space $L_{2}(\mathbb{R})$ have a pure continuous spectrum which consists of a half-line or a union of lines passing from the origin. Moreover, there may be at most a countable number of spectral singularities on the continuous spectrum of the examined operators.

In the present paper the spectrum and the resolvent of the class of a pencil of m order differential operators, with coefficients from $\mathit{AP}^{+}$, have been investigated under more general conditions. It has been proved that the operator $L_{\lambda}$ has a pure continuous spectrum. If the characteristic polynomial $\phi(z)$ has only simple roots $\omega _{k}$, $k=1,2,\ldots,m$, the continuous spectrum consists of the lines $\operatorname{Re}(\lambda\omega_{k})=0$, $k=1,2,\ldots,m$. Moreover, there may be spectral singularities (in the sense of Naimark [11]) on the continuous spectrum which coincide with numbers of the form $\lambda =0$, $\lambda_{sjn}=i\alpha _{n}( \omega _{j}-\omega _{s}) ^{-1}$, $s,j=1,2,\ldots,m$, $j\neq s$, $n\in \mathbb{N}$. If the characteristic polynomial $\phi (z)$ has one multiple root $\omega _{0}$ then the continuous spectrum consists of the line $\operatorname{Re} (\lambda \omega _{0})=0$ and a spectral singularity does not exist. The resolvent operator $L_{\lambda}^{-1}$ is an integral operator in $L_{2}( \mathbb{R})$ with the kernel of Karleman type for any $\lambda\in\rho (L_{\lambda})$. Under weakened conditions, the obtained results of this paper generalize all results of [5, 10] and some parts of results of [7–9].

2 Floquet solutions of the equation $\ell_{\lambda}(y)=0$

Here, we will show the existence of the Floquet solutions of the equation $\ell_{\lambda}(y)=0$, which plays an important role in the investigation of the spectrum of the operator $L_{\lambda}$. If the characteristic polynomial has more than one multiple root, then there may arise various cases to obtain the fundamental system of solutions. Below, we consider the cases when there exist simple roots or one multiple root.

Case I. The characteristic polynomial $\phi(z)$ has different simple roots $\omega_{1},\omega_{2},\ldots,\omega_{m}$.

Theorem 1

If (1), (2) hold and ω is any root of $\phi(z)=0$ then for each $\lambda\in \mathbb{C}$, $\lambda\neq i\alpha_{n} ( \omega_{j}-\omega ) ^{-1}$, $j=1,2,\ldots,m$, $\omega_{j}\neq\omega$, $n\in \mathbb{N}$, the differential equation

$$ y^{(m)}+\sum_{\gamma=1}^{m}p_{\gamma}(x, \lambda)y^{(m-\gamma)}=0 $$

(3)

has the solution

$$ f(x,\lambda)=e^{\omega\lambda x} \Biggl( 1+\sum_{n=1}^{\infty }U_{n}( \lambda )e^{i\alpha_{n}x} \Biggr) , $$

(4)

where

$$ U_{n}(\lambda)=U_{n}+ \sum_{k=1}^{n} \sum_{\substack{ j=1 \\ \omega_{j}\neq\omega}}^{m}\frac{U_{jkn}}{ [ i\alpha_{k}+(\omega -\omega_{j})\lambda ] },\quad \forall n\in \mathbb{N}, $$

with $U_{n},U_{jkn}\in \mathbb{C}$ and series

$$ \sum_{n=1}^{\infty}\alpha_{n}^{\gamma} \bigl\vert U_{n}(\lambda )\bigr\vert ,\quad \gamma=0,1, \ldots,m, $$

(5)

is majorized in each compact set $S\subseteq \mathbb{C}$ which does not contain the numbers $\lambda=i\alpha _{n} ( \omega_{j}-\omega ) ^{-1}$ for $j=1,2,\ldots,m$, $\omega _{j}\neq\omega$, $n\in \mathbb{N}$.

Proof

Let ω be any root of the characteristic polynomial $\phi(\omega)$. If we assume the existence of the solution of equation (3) represented as (4) with convergent series (5), then we can find the derivatives of $f(x,\lambda)$ with respect to x as

$$ f^{(\gamma)}(x,\lambda)=e^{\omega\lambda x} \Biggl( (\omega\lambda )^{\gamma}+\sum_{n=1}^{\infty}(i \alpha_{n}+\omega\lambda)^{\gamma }U_{n}( \lambda)e^{i\alpha_{n}x} \Biggr) ,\quad \gamma=0,1,\ldots,m. $$

(6)

If we substitute these derivatives in (3) and divide both sides by $e^{\omega \lambda x}$, then we get

$$\begin{aligned}& (\omega\lambda)^{m}+\sum_{n=1}^{\infty} (i\alpha_{n} +\omega\lambda)^{m}{U}_{n} ( \lambda)e^{i\alpha_{n}x}+\sum_{\gamma=1}^{m}(\omega \lambda)^{m-\gamma }p_{\gamma\gamma}\lambda^{\gamma} \\& \quad {}+\sum _{\gamma=1}^{m}p_{\gamma \gamma}\lambda^{\gamma}\sum _{n=1}^{\infty}(i\alpha_{n}+\omega \lambda)^{m-\gamma}U_{n}(\lambda)e^{i\alpha_{n}x} \\& \quad {}+\sum_{\gamma=1}^{m}(\omega \lambda)^{m-\gamma }\sum_{n=1}^{\infty} \tilde{p}_{\gamma n}(\lambda)e^{i\alpha _{n}x} \\& \quad {} +\sum _{\gamma=1}^{m}\sum_{n=1}^{\infty} \tilde{p}_{\gamma n}(\lambda)e^{i\alpha_{n}x}\sum _{n=1}^{\infty}(i\alpha _{n}+\omega \lambda)^{m-\gamma}U_{n}(\lambda)e^{i\alpha_{n}x}=0 \end{aligned}$$

or

$$\begin{aligned}& {\lambda}^{m}\phi(\omega)+\sum_{n=1}^{\infty} \Biggl[ (i\alpha_{n}+\omega\lambda)^{m}+\sum _{\gamma=1}^{m}p_{\gamma \gamma}\lambda^{\gamma}(i \alpha_{n}+\omega\lambda)^{m-\gamma } \Biggr] U_{n}( \lambda)e^{i\alpha_{n}x} \\& \quad {}+\sum_{n=1}^{\infty }\sum _{\gamma=1}^{m}(\omega\lambda)^{m-\gamma} \tilde{p}_{\gamma n}(\lambda)e^{i\alpha_{n}x} \\& \quad {}+\sum_{\gamma=1}^{m}\sum _{n=2}^{\infty} \biggl( \sum _{\alpha_{r}+\alpha_{s}=\alpha_{n}}\tilde{p}_{\gamma s}(\lambda) (i \alpha_{r}+\omega\lambda)^{m-\gamma}U_{r}(\lambda ) \biggr) e^{i\alpha_{n}x}=0 . \end{aligned}$$

Taking into account the uniqueness theorem for almost-periodic functions we have

$$\begin{aligned}& \Biggl[ (i\alpha_{n}+\omega\lambda)^{m}+\sum _{\gamma =1}^{m}p_{\gamma\gamma}\lambda^{\gamma}(i \alpha_{n}+\omega\lambda )^{m-\gamma} \Biggr] U_{n}( \lambda)+\sum_{\gamma=1}^{m}(\omega \lambda)^{m-\gamma}\tilde{p}_{\gamma n}(\lambda) \\& \quad {} +\sum_{\gamma=1}^{m} \sum _{\alpha_{r}+\alpha _{s}=\alpha_{n}}\tilde{p}_{\gamma s}(\lambda) (i\alpha _{r}+\omega \lambda)^{m-\gamma}U_{r}(\lambda)=0 , \quad n\in \mathbb{N} . \end{aligned}$$

Using the expansion

$$\begin{aligned}& (i\alpha_{n}+\omega\lambda)^{m}+\sum _{\gamma=1}^{m}p_{\gamma \gamma}\lambda^{\gamma}(i \alpha_{n}+\omega\lambda)^{m-\gamma} \\& \quad =\lambda^{m} \Biggl[ \biggl(\frac{i\alpha_{n}}{\lambda}+\omega \biggr)^{m}+\sum_{\gamma =1}^{m}p_{\gamma\gamma} \biggl(\frac{i\alpha_{n}}{\lambda}+\omega \biggr)^{m-\gamma} \Biggr]= \lambda^{m}\phi\biggl(\frac{i\alpha_{n}}{\lambda}+\omega\biggr) \\& \quad = \bigl[ i\alpha _{n}+\lambda ( \omega-\omega_{1} ) \bigr] \cdot \bigl[ i\alpha _{n}+\lambda ( \omega- \omega_{2} ) \bigr] \cdot\ldots\cdot \bigl[ i\alpha_{n}+ \lambda ( \omega-\omega_{m} ) \bigr] , \end{aligned}$$

we obtain

$$ U_{n}(\lambda)=- \frac{\sum_{\gamma=1}^{m}(\lambda { \omega})^{m-\gamma}\tilde{p}_{\gamma n}(\lambda )+\sum_{\gamma=1}^{m} \sum_{\alpha_{r}+\alpha _{s}=\alpha_{n}}(i\alpha_{r}+{ \omega}\lambda)^{m-\gamma}\tilde{p}_{\gamma s}(\lambda){U_{r}(\lambda)}}{[i\alpha _{n}+\lambda(\omega-\omega_{1})]\cdot[ i\alpha_{n}+\lambda (\omega-\omega_{2})]\cdot\ldots\cdot[ i\alpha_{n}+\lambda (\omega -\omega_{m})]} $$

(7)

for $\lambda\in \mathbb{C}$, $\lambda\neq i\alpha_{n} ( \omega_{j}-\omega ) ^{-1}$, $j=1,2,\ldots,m$, $\omega_{j}\neq\omega$, $n\in \mathbb{N}$.

On the contrary, if $\{ { U}_{n}{ (\lambda)} \}$ satisfies the system of equations (7) and the series (5) converges, then it can be shown that $f(x,\lambda)$ determined by (4) is a solution of (3). Therefore, the solvability of (7) and the convergence of the series (5) are sufficient to prove the theorem.

From (7), $\{ U_{n}(\lambda) \}$ is determined by the recurrent manner uniquely. It is possible to see that $U_{n}(\lambda)$ is the rational function which can have simple poles $\lambda=i\alpha _{k}(\omega _{j}-\omega)^{-1}$, $j=1,2,\ldots,m$, $\omega_{j}\neq\omega$, $k=1,2,\ldots,n$, and therefore it can be uniquely written as

$$ U_{n}(\lambda)=U_{n}+ \sum_{k=1}^{n} \sum_{\substack{ j=1 \\ \omega_{j}\neq\omega}}^{m}\frac{U_{jkn}}{ [ i\alpha_{k}+(\omega -\omega_{j})\lambda ] },\quad \forall n\in \mathbb{N}, $$

where $U_{n},U_{jkn}\in \mathbb{C}$. Let $S\subseteq \mathbb{C}$ be a compact set which does not contain the points $\lambda=i\alpha _{n}(\omega_{j}-\omega)^{-1}$ for $j=1,2,\ldots,m$, $\omega_{j}\neq \omega$, $n\in\mathbb{N}$. Let us show that the series (5) is majorized in S for $\{ { U}_{n}{ (\lambda)} \}$ which is determined from (7).

It is obvious that there exist $c_{0}>0$, $q>1$ such that

$$ c_{0}\alpha_{n}^{m}\leq\bigl\vert \bigl[ i \alpha_{n}+(\omega-\omega _{1})\lambda\bigr]\cdot\bigl[ i \alpha_{n}+(\omega-\omega_{2})\lambda \bigr]\cdot \ldots \cdot\bigl[ i\alpha_{n}+(\omega-\omega_{m})\lambda\bigr]\bigr\vert $$

and $\vert \omega\lambda \vert \leq q$ for $\forall n\in \mathbb{N}$, $\forall\lambda\in S$. Then from (7), we have

$$\begin{aligned} c_{0}\alpha_{n}^{m}\bigl\vert U_{n}(\lambda)\bigr\vert \leq& \sum_{\gamma=1}^{m} \vert { \omega}\lambda \vert ^{m-\gamma}\bigl\vert \tilde{p}_{\gamma n}( \lambda)\bigr\vert +\sum_{\gamma=1}^{m} \sum _{\alpha_{r}+\alpha _{s}=\alpha_{n}}\bigl\vert U_{r}(\lambda)\bigr\vert \cdot \vert i{{\alpha _{r}+\omega\lambda}} \vert ^{m-\gamma}\bigl\vert \tilde{p}_{\gamma s}(\lambda)\bigr\vert \\ \leq&\sum_{\gamma=1}^{m}q^{m-\gamma} \tilde{p}_{\gamma n}+\sum_{\gamma=1}^{m} \sum _{\alpha_{r}+\alpha _{s}=\alpha_{n}}\bigl\vert U_{r}(\lambda)\bigr\vert \bigl(\vert {{\alpha _{r}}}\vert +q\bigr)^{m-\gamma} \tilde{p}_{\gamma s} \\ \leq& q^{m-1}\sum_{\gamma=1}^{m} \tilde{p}_{\gamma n}+\sum_{\gamma=1}^{m} \sum _{\alpha_{r}+\alpha _{s}=\alpha_{n}}\bigl\vert U_{r}(\lambda)\bigr\vert \alpha _{r}^{m-\gamma }\biggl(1+\frac{q}{{{\alpha}_{1}}} \biggr)^{m-\gamma}\tilde{p}_{\gamma s} \\ \leq& q^{m-1}\sum_{\gamma=1}^{m} \tilde{p}_{\gamma n}+\biggl(1+\frac{q}{{{\alpha}_{1}}}\biggr)^{m-1} \sum _{\alpha_{r}+\alpha_{s}=\alpha _{n}}\bigl\vert U_{r}(\lambda)\bigr\vert \alpha _{r}^{m-1}\cdot\sum _{\gamma=1}^{m}\tilde{p}_{\gamma s} \end{aligned}$$

for $\forall\lambda\in S$, $\forall n\in \mathbb{N}$ and $\tilde{p}_{\gamma n}=\sup_{\lambda\in S}\vert \tilde{p}_{\gamma n}{ (\lambda)}\vert$, $\gamma=1,2,\ldots,m$.

Let

$$\begin{aligned}& A_{n}= \frac{q^{m-1}}{c_{0}}\sum_{\gamma=1}^{m}\tilde{p}_{\gamma n} , \\& B_{n}= \frac{1}{c_{0}} \biggl( 1+\frac{q}{\alpha_{1}} \biggr) ^{m-1}\sum_{\gamma=1}^{m} \tilde{p}_{\gamma n} . \end{aligned}$$

Then from the last inequality we obtain

$$ \alpha_{n}^{m}\bigl\vert U_{n}(\lambda)\bigr\vert \leq A_{n}+\sum_{\alpha_{r}+\alpha_{s}=\alpha_{n}} \alpha _{r}^{m-1}\bigl\vert U_{r}(\lambda)\bigr\vert B_{s} ,\quad n\in \mathbb{N}. $$

If $u_{n}= \sup_{\lambda\in S}\vert U_{n}(\lambda )\vert$, then we have

$$ \alpha_{n}^{m}u_{n}\leq A_{n}+ \sum _{\alpha_{r}+\alpha _{s}=\alpha_{n}}\alpha_{r}^{m-1}u_{r}B_{s} , \quad n\in \mathbb{N}$$

or

$$\begin{aligned} \sum_{n=1}^{t}\alpha_{n}^{m}u_{n} \leq& \sum_{n=1}^{t}A_{n}+\sum _{n=1}^{t}\sum_{\alpha_{r}+\alpha_{s}=\alpha _{n}} \alpha_{r}^{m-1}u_{r}B_{s} \\ \leq& A+\sum_{r=1}^{t-1}\alpha _{r}^{m-1}u_{r}\sum_{s=1}^{t-1}B_{s} \leq A+B\sum_{n=1}^{t-1}\alpha _{n}^{m-1}u_{n},\quad n\in \mathbb{N}. \end{aligned}$$

From (2) it is clear that $A= \sum_{n=1}^{\infty}A_{n}<+\infty$ and $B= \sum_{n=1}^{\infty}B_{n}<+\infty$.

Therefore, for all $t\in \mathbb{N}$, $\sum_{n=1}^{t}\alpha_{n}^{m}u_{n}\leq A+B\sum_{n=1}^{t-1}\alpha_{n}^{m-1}u_{n}$ is satisfied. By using this inequality and $\alpha_{n}\rightarrow+\infty$, we can easily show the convergence of the series $\sum_{n=1}^{+\infty}\alpha _{n}^{m}u_{n}$ according to the lemma in [12] (see [12], pp.21-22). In this case, $\sum_{n=1}^{+\infty}\alpha_{n}^{m}| U_{n}(\lambda)|$ is a majorized series in S. According to the Weierstrass theorem, the series (5) is uniform convergent in S. Since $S\subseteq \mathbb{C}$ is an arbitrarily chosen compact set, the series (5) is convergent for all $\lambda\neq i\alpha_{n}(\omega_{j}-\omega)^{-1}$, $j=1,2,\ldots,m$, $\omega_{j}\neq\omega$, and $n\in \mathbb{N}$. Thus $f(x,\lambda)$ is a solution of equation (3). The theorem is proved. □

It is clear that $\lambda_{jn}=i\alpha_{n}(\omega_{j}-\omega)^{-1}$ may be a singular point of $f(x,\lambda)$ for any $j=1,2,\ldots,m$, $\omega _{j}\neq\omega$, $n\in \mathbb{N}$. Actually, according to Theorem 1, the functional series in the representation

$$ \bigl[ i\alpha_{n}+(\omega-\omega_{j}) \bigr] f(x, \lambda)=e^{\omega \lambda x} \Biggl( 1+\sum_{r=1}^{\infty} \bigl[ i\alpha_{n}+(\omega -\omega _{j})\lambda \bigr] U_{r}(\lambda)e^{i\alpha_{r}x} \Biggr) $$

and the obtained series by m times term by term differentiation are absolutely and uniformly convergent with respect to λ in the closed disk with a small radius centered in point $\lambda_{jn}$. Therefore, the finite limits

$$ \lim_{\lambda\rightarrow\lambda_{jn}} \bigl[ i\alpha _{n}+(\omega- \omega_{j})\lambda \bigr] \frac{\partial ^{s}f(x,\lambda)}{\partial x^{s}}=e^{\omega\lambda_{jn}x}{\sum _{r=n}^{\infty }U_{jnr}(\omega \lambda_{jn}+i\alpha_{r})^{s}}e^{i\alpha_{r}x}, \quad s=0,1,2,\ldots,m, $$

exist. Moreover, the series ${\sum_{r=n}^{\infty}} \vert {U_{jnr}}\vert {\alpha_{r}^{m}}$ is convergent. If this limit is not zero, then it means that the point $\lambda_{jn}$ is a simple pole of the functions $\frac{\partial^{s}f(x,\lambda)}{\partial x^{s}}$, $s=0,1,\ldots,m$.

Corollary 1

For $\forall x\in R$, the functions $\frac {\partial ^{s}f(x,\lambda)}{\partial x^{s}}$, $s=0,1,\ldots,m$, are meromorphic functions with respect to λ and they can have only simple poles ${\lambda_{jn}}=i\alpha_{n}(\omega_{j}-\omega)^{-1}$, $j=1,2,\ldots,m$, $\omega_{j}\neq\omega$, $n\in\mathbb{N}$. Moreover, these functions are also continuous functions of the pair $(x,\lambda)$ for all $(x,\lambda)\in \mathbb{R}\times\mathbb{C}$, $\lambda\neq i\alpha_{n}(\omega_{j}-\omega)^{-1}$, $j=1,2,\ldots,m$, $\omega_{j}\neq\omega$, $n\in \mathbb{N}$.

Corollary 2

For $\forall\lambda\neq\lambda _{sjn}$, $s,j=1,2,\ldots,m$, $j\neq s$, $n\in \mathbb{N}$, equation (3) has the Floquet solutions represented as

$$ f_{s}(x,\lambda)=e^{\omega_{s}\lambda x} \Biggl( {1+\sum _{n=1}^{\infty}U_{n}^{(s)}( \lambda)e^{i\alpha_{n}x}} \Biggr) , $$

where

$$ U_{n}^{(s)}(\lambda)=U_{n}^{(s)}+ \sum _{k=1}^{n}\sum _{\substack{ j=1 \\ j\neq s}}^{m}\frac{U_{jkn}^{(s)}}{i\alpha_{k}+(\omega _{s}-\omega _{j})\lambda}. $$

The Wronskian of the functions $f_{1}(x,\lambda ), f_{2}(x,\lambda),\ldots, f_{m}(x,\lambda)$ for $\lambda\in\mathbb{C}\backslash\Lambda$ is found as

$$ W[f_{1},f_{2},\ldots,f_{m}]=W(x,\lambda)= \lambda^{\frac{m(m-1)}{2}}W_{m}e^{-\lambda p_{11}x-\sum_{n=1}^{\infty}\frac{p_{10n}}{i\alpha _{n}}e^{i\alpha_{n}x}}. $$

(8)

Here $W_{m}$ is the Vandermonde determinant of the numbers $\omega_{1},\omega_{2},\ldots,\omega_{m}$. The functions $f_{1}(x,\lambda), f_{2}(x,\lambda),\ldots, f_{m}(x,\lambda)$ form the fundamental system of solutions of equation (3) in the interval $( -\infty,+\infty )$ for $\forall\lambda\in \mathbb{C}\backslash\Lambda_{0}$.

It is clear that the existence of the solutions $f_{1}(x,\lambda), f_{2}(x,\lambda),\ldots, f_{m}(x,\lambda)$ follows from Theorem 1 for $\omega=\omega_{s}$, $s=1,2,\ldots,m$. Since every $q(x)\in \mathit{AP}^{+}$ can be extended to the upper semi-plane as an analytic function of x and $\lim_{\operatorname{Im}x\rightarrow+\infty}q(x)=0$, by passing to the limit as $\operatorname{Im}x\rightarrow+\infty$ on both sides of the equation

$$ W(x,\lambda)e^{\lambda p_{11}x}=W(0,\lambda )e^{-\int_{0}^{x}p_{10}(t)\, dt} $$

we get equation (8). From (8) it follows that if $\lambda\neq0$, then $W(x,\lambda)\neq0$. Thus the system of functions is independent.

Note that the solutions of type $f_{s}(x,\lambda )$, $s=1,2,\ldots,m$, are obtained in [6–9] under the different conditions and in various forms of the representation.

According to Corollary 1, it is obvious that the function

$$ f_{sjn}(x)= {\lim_{\lambda\rightarrow\lambda _{sjn}}f_{s}(x,\lambda)} \bigl[ {i\alpha_{n}+(\omega_{s}-\omega _{j}) \lambda} \bigr] {=e^{{\omega_{s}}\lambda_{sjn}x}\sum_{k=n}^{\infty }U_{jnk}^{(s)}e^{i\alpha_{k}x}} $$

is a solution of equation (3) for $\lambda=\lambda_{sjn}$, where the series ${\sum_{k=n}^{\infty}\alpha_{k}^{m}}\vert {U_{jnk}^{(s)}} \vert$ is convergent. As in [4], writing the equation which is satisfied by the coefficients $U_{jnk}^{(s)}$ (see [4], p.778), it is seen easily that $U_{jnk}^{(s)}=0$ for every $k\geq n$ if $U_{jnn}^{(s)}=0$. Therefore, $f_{sjn}(x)\equiv0$ if and only if $U_{jnn}^{(s)}=0$. In this case, $f_{s}(x,\lambda)$ is regular at point $\lambda=\lambda_{sjn}$ and $f_{s}(x,\lambda_{sjn})$ is a solution of equation (3). It can be shown that the functions $f_{sjn}(x)$ and $f_{j}(x,\lambda_{sjn})$ are linearly dependent. Moreover, $f_{sjn}(x)=U_{jnn}^{(s)}f_{j}(x,\lambda_{sjn})$ for any $s=1,2,\ldots,m$, $s\neq j=1,2,\ldots,m$, $n\in \mathbb{N}$, is valid which is important for establishing the fundamental system of solutions of equation (3) for $\lambda=\lambda_{sjn}$.

Let s, j, n be fixed and $f_{s}(x,\lambda), f_{s_{1}}(x,\lambda ),f_{s_{2}}(x,\lambda),\ldots,f_{s_{\mu}}(x,\lambda)$ be all functions which have a pole at the point $\lambda=\lambda_{sjn}$. It is only possible when the equality $\lambda_{sjn}=\lambda_{s_{\beta}j_{\beta}{ n}_{\beta }}$, $\beta=1,2,\ldots,\mu$, is valid for some different indices $n_{1},n_{2},\ldots,n_{\mu}\in \mathbb{N}$, and $1\leq j_{1},j_{2},\ldots,j_{\mu}\leq m$. Then all other functions $f_{j}(x,\lambda), f_{j_{1}}(x,\lambda), f_{j_{2}}(x,\lambda ),\ldots,f_{j_{\mu}}(x,\lambda), f_{j_{\mu+1}}(x,\lambda ),\ldots,f_{j_{\nu }}(x,\lambda)$, $\mu+\nu+2=m$, are regular at the point $\lambda =\lambda _{sjn}$.

If we define the functions $f_{kjn}(x)= {\lim_{\lambda \rightarrow\lambda_{sjn}}f_{k}(x,\lambda)} [ {i\alpha _{n}+(\omega _{s}-\omega_{j})\lambda} ]$, $k=s_{1},s_{2},\ldots,s_{\mu}$, then it is obvious that the functions

$$ \begin{aligned} &f_{sjn}(x),\qquad f_{s_{1}jn}(x),\qquad f_{s_{2}jn}(x), \qquad \ldots, \\ &f_{s_{\mu }jn}(x),\qquad f_{k}(x, \lambda_{sjn}),\quad k=j,j_{1},j_{2}, \ldots,j_{\nu},\end{aligned} $$

(9)

are solutions of equation (3) for $\lambda=\lambda_{sjn}$ and the functions of this system are linear dependent in $( -\infty ,+\infty )$, since their Wronskian is equal to zero. Moreover, any three of the numbers $\operatorname{Re}(\lambda _{sjn}\omega _{k})$, $k=1,2,\ldots,m$, can not be equal and there are some equal pairs between them. These equal pairs are $\operatorname{Re}(\lambda _{sjn}\omega _{s})=\operatorname{Re}(\lambda _{sjn}\omega _{j})$, $\operatorname{Re}(\lambda _{sjn}\omega _{s_{\beta }})= \operatorname{Re}(\lambda _{sjn}\omega _{j_{\beta }})$, $\beta =1,2,\ldots,\mu$. Then taking into our account the behaviors as $x\rightarrow \pm \infty$ of the functions belonging to the system (9), as it is shown in [12] (see pp.43-45), we have the existence of some constants $b_{k}$, $k=s$, $s_{1},s_{2},\ldots,s_{\mu }$ such that

$$ f_{sjn}(x)=b_{s}f_{j}(x,\lambda _{sjn}),\qquad f_{s_{\beta }jn}(x)=b_{s_{\beta }}f_{j_{\beta }}(x,\lambda _{sjn}),\quad \beta =1,2,\ldots,\mu . $$

From the equality $f_{sjn}(x)=b_{s}f_{j}(x,\lambda_{sjn})$, according to the uniqueness theorem for almost-periodic functions, it is seen that $b_{s}=U_{jnn}^{(s)}$. Using the equalities

$$ f_{sjn}(x)=b_{s}f_{j}(x,\lambda_{sjn}), \qquad f_{s_{\beta }jn}(x)=b_{s_{\beta}}f_{j_{\beta}}(x, \lambda_{sjn}),\quad \beta =1,2,\ldots,\mu, $$

the system of linearly independent solutions of equation (3) corresponding to $\lambda=\lambda_{sjn}$ can be established.

Since the functions $f_{k}(x,\lambda)$, $k=j,j_{1},j_{2},\ldots,j_{\mu }$, are regular at $\lambda=\lambda_{sjn}$, the functions

$$\begin{aligned}& \tilde{f}_{sjn}(x)= \lim_{\lambda\rightarrow\lambda _{sjn}} \biggl( f_{s}(x,\lambda)- {\frac{b_{s}f_{j}(x,\lambda)}{i\alpha_{n}+(\omega_{s}-\omega_{j})\lambda}} \biggr) , \\& \tilde{f}_{s_{\beta}jn}(x)= \lim_{\lambda\rightarrow \lambda_{sjn}} \biggl( f_{s_{\beta}}(x,\lambda)- {\frac{b_{s_{\beta}}f_{j_{\beta}}(x,\lambda)}{i\alpha_{n}+(\omega _{s}-\omega _{j})\lambda}} \biggr) ,\quad \beta=1,2, \ldots,\mu, \end{aligned}$$

are also solutions of equation (3) corresponding to $\lambda =\lambda_{sjn}$. According to the expressions of the functions $f_{k}(x,\lambda )$, $k=s, s_{1},s_{2},\ldots,s_{\mu},j,j_{1},j_{2},\ldots,j_{\mu}$, we conclude that $\tilde{f}_{kjn}(x)=e^{i\omega_{k}\lambda_{sjn}x}(\psi _{kjn}(x)+x\phi _{kjn}(x))$, where $\psi_{kjn}(x)$ and $\phi_{kjn}(x)$ are Bohr almost-periodic functions for $k=s,s_{1},s_{2},\ldots,s_{\mu}$. From the explicit form of the functions $\tilde{f}_{sjn}(x),\tilde{f}_{s_{1}jn}(x), \tilde{f}_{s_{2}jn}(x),\ldots, \tilde{f}_{s_{\mu}jn}(x), f_{k}(x,\lambda_{sjn})$, $k=j,j_{1},j_{2},\ldots,j_{\nu}$, it is seen that these functions are linearly independent on $(-\infty,+\infty)$. Therefore, these functions form a fundamental system of solutions of equation (3) for $\lambda=\lambda_{sjn}$.

Now let us construct the linearly independent solutions of equation (3) for $\lambda=0$.

Note that, since the Wronskian of the solutions $f_{s}(x,\lambda)$, $s=1,2,\ldots,m$, is equal to zero for $\lambda=0$, they are linearly dependent. Linearly independent solutions of equation (3) corresponding to $\lambda =0$ are established according to Theorem 1. It is clear that solutions of the equation

$$ y^{(m)}+\sum_{\gamma=1}^{m}p_{\gamma0}(x)y^{(m-\gamma)}= \lambda ^{m}y $$

(10)

corresponding to $\lambda=0$ are also solutions of equation (3) for $\lambda=0$. By Theorem 1, equation (10) has the solution

$$ \tilde{f}(x,\lambda)=e^{\lambda x} \Biggl( {1+\sum _{n=1}^{\infty}}\widetilde{{U}} {_{n}( \lambda)e^{i\alpha_{n}x}} \Biggr), $$

which is analytic with respect to λ in some small neighborhood of $\lambda=0$. By putting $\tilde{f}(x,\lambda)$ in (10) and by differentiating equation (10) with respect to λ, it is sure that functions $\tilde{f}_{s}(x)= \frac{\partial^{s}\tilde{f}(x,\lambda)}{\partial\lambda^{s}}|_{\lambda=0}$, $s=0,1,\ldots,m-1$, are also solutions of (10) and (3) corresponding to $\lambda =0$. We can see easily that $\tilde{f}_{0}(x)=\alpha_{00}(x)$, $\tilde {f}_{1}(x)=x\alpha _{11}(x)+\alpha_{10}(x)$, … , $\tilde{f}_{m-1}(x)=x^{m-1}\alpha _{m-1,m-1}(x)+x^{m-2}\alpha_{m-1,m-2}(x)+\cdots+\alpha_{m-1,0}(x)$, where $\alpha_{sj}(x)$, $s=0,1,\ldots,m-1$, $j=0,1,\ldots,s$, are Bohr almost-periodic functions and $\alpha_{ss}(x)$, $s=0,1,\ldots,m-1$, are nonzero. The linear independence of $\tilde{f}_{s}(x)$, $s=0,1,\ldots,m-1$, in $(-\infty ,+\infty)$ is seen from their open form.

Case II. The characteristic polynomial has a unique multiple root $\omega_{0}$, i.e. $\phi(z)=(z-\omega_{0})^{m}$.

In this case, to find the particular solutions of equation (3) we will use the following theorem.

Theorem 2

If the characteristic polynomial has a unique multiple root $\omega_{0}$, then for each function $g(x,\lambda )=e^{\omega_{0}\lambda x}\sum_{n=1}^{\infty}g_{n}(\lambda )e^{i\alpha_{n}x}$ such that $g_{n}(\lambda)$, $n\in \mathbb{N}$, are polynomials whose degree does not exceed $n(m-1)$ and the series $\sum_{n=1}^{\infty} \vert g_{n}(\lambda )\vert$ is majorized in any compact set $S\subseteq \mathbb{C}$, the equation

$$ y^{(m)}+\sum_{\gamma=1}^{m}p_{\gamma}(x, \lambda)y^{(m-\gamma )}=g(x,\lambda) $$

(11)

has a solution

$$ h(x,\lambda)=e^{\omega_{0}\lambda x} \Biggl( {1+\sum _{n=1}^{\infty}h_{n}(\lambda)e^{i\alpha_{n}x}} \Biggr) $$

(12)

in $(-\infty,+\infty)$ for every $\lambda\in \mathbb{C}$. Here the coefficients $h_{n}(\lambda)$, $n\in \mathbb{N}$, are polynomials whose degrees do not exceed $n(m-1)$, and the series $\sum_{n=1}^{\infty} \vert {h_{n}(\lambda)} \vert \alpha _{n}^{m}$ is majorized in each compact set $S\subseteq \mathbb{C}$.

Proof

If we substitute the function (12) in (11), to find the coefficients sequence $\{ {h_{n}(\lambda)} \}$ as in the proof of Theorem 1, we obtain a system of equations,

$$\begin{aligned}& \Biggl[ (i\alpha_{n}+\omega\lambda)^{m}+\sum _{\gamma =1}^{m}p_{\gamma\gamma}\lambda^{\gamma}(i \alpha_{n}+\omega\lambda )^{m-\gamma} \Biggr] h_{n}( \lambda)+\sum_{\gamma=1}^{m}(\omega \lambda)^{m-\gamma}\tilde{p}_{\gamma n}(\lambda) \\ & \quad {}+\sum_{\gamma=1}^{m} \sum _{\alpha_{r}+\alpha _{s}=\alpha_{n}}\tilde{p}_{\gamma s}(\lambda) (i\alpha _{r}+\omega \lambda)^{m-\gamma}h_{r}( \lambda)=g_{n}(\lambda),\quad n\in \mathbb{N}, \lambda\in \mathbb{C} \end{aligned}$$

or

$$\begin{aligned}& ( i\alpha_{n} ) ^{m}h_{n}(\lambda)+\sum _{\gamma =1}^{m}(\omega\lambda)^{m-\gamma} \tilde{p}_{\gamma n}(\lambda) \\ & \quad {}+\sum_{\gamma=1}^{m} \sum _{\alpha_{r}+\alpha _{s}=\alpha_{n}}\tilde{p}_{\gamma s}(\lambda) (i\alpha _{r}+\omega \lambda)^{m-\gamma}h_{r}( \lambda)=g_{n}(\lambda), \quad n\in \mathbb{N}, \lambda\in \mathbb{C}, \end{aligned}$$

or

$$\begin{aligned}& h_{n}(\lambda)=- \frac{\sum_{\gamma=1}^{m}(\lambda { \omega})^{m-\gamma}\tilde{p}_{\gamma n}(\lambda )+\sum_{\gamma=1}^{m} {\sum_{\alpha_{r}+\alpha _{s}=\alpha_{n}}(i\alpha_{r}+{ \omega}\lambda)}^{m-\gamma}\tilde{p}_{\gamma s}(\lambda){h_{r}(\lambda)-}g_{n}(\lambda)}{ ( i\alpha_{n} ) ^{m}}, \\ & \quad n\in \mathbb{N}, \lambda \in \mathbb{C}. \end{aligned}$$

(13)

The coefficients ${h_{n}(\lambda)}$ are found uniquely from equation (13). In fact, the degree of the polynomial $h_{1}(\lambda)=-\frac{\sum_{\gamma =1}^{m}(\lambda{ \omega})^{m-\gamma}\tilde{p}_{\gamma 1}(\lambda)-g_{1}(\lambda)}{ ( i\alpha_{1} ) ^{m}}$ for $n=1$ does not exceed $m-1$. Subsequently, for $n=2$ the degree of the polynomial $h_{2}(\lambda)$ does not exceed $2 ( m-1 )$ and, for each n, $h_{n}(\lambda)$ is found as a polynomial whose degree does not exceed $n ( m-1 )$. If for the obtained coefficients $h_{n}(\lambda )$, $n\in\mathbb{N}$, $h_{n}= \sup_{\lambda\in S}\vert h_{n}(\lambda )\vert$, then convergence of the series $\sum_{n=1}^{+\infty }\alpha_{n}^{m}h_{n}$ and so majorization of the series $\sum_{n=1}^{+\infty}\alpha_{n}^{m}\vert h_{n}(\lambda) \vert$ in the set $S\subseteq \mathbb{C}$ easily can be shown as in the proof of Theorem 1. Therefore, for each $\lambda\in S$ the function $h(x,\lambda)=e^{\omega\lambda x} ( {1+\sum_{n=1}^{\infty}h_{n}(\lambda)e^{i\alpha_{n}x}} )$ is the solution of equation (11) in $(-\infty,+\infty)$. The theorem is proved. □

Corollary 3

If the characteristic polynomial $\phi(z)$ has a unique multiple root $\omega_{0}$, then equation (3) has a solution $\hat{f}(x,\lambda)=e^{\omega\lambda x}q(x,\lambda )$ in $(-\infty,+\infty)$ for every $\lambda\in \mathbb{C}$. Here $q(x,\lambda)=1+\sum_{n=1}^{\infty}q_{n}(\lambda )e^{i\alpha_{n}x}$ is a Bohr almost-periodic function. The $q_{n}(\lambda)$, $n\in \mathbb{N}$, are polynomials whose degree does not exceed $n(m-1)$, the series $\sum_{n=1}^{+\infty}\alpha_{n}^{m}\vert q_{n}(\lambda)\vert$ is majorized in each compact set $S\subseteq\mathbb{C}$. $\hat{f}(x,\lambda),\frac{\partial\hat {f}(x,\lambda)}{\partial x},\ldots,\frac{\partial^{m}\hat{f}(x,\lambda)}{\partial x^{m}}$ are continuous functions in $\mathbb{R}\times \mathbb{C}$ with respect to the ordered pair $(x,\lambda)$ and they are an entire function of λ.

To prove Corollary 3, it is enough to take $g(x,\lambda )=0$ in Theorem 2 and to see m times differentiability term by term of the series in the expression of the obtained solution $\hat{f}(x,\lambda )$ with respect to x. Here the obtained series are uniformly convergent in each bounded set of the ordered pairs $(x,\lambda )$, therefore functions $\hat{f}(x,\lambda ),\frac{\partial \hat{f}(x,\lambda )}{\partial x},\ldots,\frac{\partial ^{m}\hat{f}(x,\lambda )}{\partial x^{m}}$ are continuous functions of the ordered pairs $(x,\lambda )$ and they are entire functions of λ.

Theorem 3

If $\phi(z)=(z-\omega_{0})^{m}$, then equation (3) has Floquet solutions in the interval $(-\infty ,+\infty)$ as

$$\begin{aligned}& \hat{f}_{1}(x,\lambda) = e^{\omega\lambda x}q_{1}(x,\lambda ),\qquad \hat{f}_{2}(x,\lambda)=e^{\omega\lambda x} \bigl[ xq_{1}(x,\lambda )+q_{2}(x,\lambda) \bigr] ,\qquad \ldots, \\& \hat{f}_{m}(x,\lambda) = e^{\omega\lambda x} \biggl[ \frac {x^{m-1}}{(m-1)!}q_{1}(x,\lambda)+\cdots+xq_{m-1}(x, \lambda)+q_{m}(x,\lambda ) \biggr], \end{aligned}$$

where the functions $q_{1}(x,\lambda), q_{2}(x,\lambda ),\ldots,q_{m}(x,\lambda)$ are almost-periodic function:

$$ q_{s}(x,\lambda)=1+ \sum_{n=1}^{\infty}q_{sn}( \lambda )e^{i\alpha_{n}x},\quad s=1,2,\ldots,m $$

for $\forall\lambda\in \mathbb{C}$. Here $q_{sn}(\lambda)$, $s=1,2,\ldots,m$, $n\in\mathbb{N}$, are polynomials whose degrees do not exceed $n(m-1)$ and the series $\sum_{n=1}^{\infty}\alpha_{n}^{m}\vert q_{sn}(\lambda )\vert$ are majorized in each compact set $S\subseteq \mathbb{C}$.

Proof

When $\phi(z)= ( z-\omega_{0} ) ^{m}$, equation (3) has a solution $\hat{f}_{1}(x,\lambda)=e^{\omega \lambda x}q_{1}(x,\lambda)$ according to Corollary 3. In order to obtain other solutions which form a fundamental system of solutions of equation (3) together with $\hat{f}_{1}(x,\lambda)$, let us use the properties of the linear differential operator $L:C^{m}( \mathbb{R} )\rightarrow C( \mathbb{R})$, which is defined as

$$ L(y)=p_{0}(x)y^{(m)}+p_{1}(x)y^{(m-1)}+p_{2}(x)y^{(m-2)}+ \cdots+p_{m-1}(x)y^{ \prime}+p_{m}(x)y, $$

(14)

where $p_{j}(x)\in C( \mathbb{R})$, $j=0,1,\ldots,m$. Let us define the operators $L^{(k)}(y)=\sum_{ \gamma=0}^{m-k}A_{m}^{k}p_{\gamma}(x) y^{(m-\gamma-k)}$, $L^{(k)}:C^{m}(\mathbb{R} )\rightarrow C( \mathbb{R})$, $k=1,2,\ldots,m$. Here, $A_{m}^{k}=m(m-1)\cdot\ldots\cdot(m-k+1)$, $k=1,2,\ldots,m$, and $A_{m}^{0}=1$. For any system of functions $y_{1}(x),y_{2}(x),\ldots,y_{m}(x)\in C^{m}( \mathbb{R})$ it is not difficult to show that the identities

$$\begin{aligned}& L(xy_{1}+y_{2})=L(y_{2})+L^{(1)}(y_{1})+xL(y_{1}), \\& L\biggl(\frac {x^{2}}{2!} y_{1}+xy_{2}+y_{3} \biggr)= L(y_{3})+L^{(1)}(y_{2})+ \frac{1}{2!}L^{(2)}(y_{1})+x \bigl[ L(y_{2})+L^{(1)}(y_{1}) \bigr] + \frac{x^{2}}{2!}L(y_{1}), \\& \ldots, \\& L\biggl(\frac{x^{s-1}}{(s-1)!}y_{1}+\frac{x^{s-2}}{(s-2)!}y_{2}+ \cdots+\frac {x}{1!}y_{s-1}+y_{s}\biggr) \\& \quad =L(y_{s})+L^{(1)}(y_{s-1})+ \frac{1}{2!}L^{(2)}(y_{s-2})+\cdots+\frac{1}{(s-1)!}L^{(s-1)}(y_{1}) \\& \qquad {}+x \biggl[ L(y_{s-1})+\frac{1}{1!}L^{(1)}(y_{s-2})+ \cdots+\frac{1}{(s-2)!}L^{(s-2)}(y_{1}) \biggr] +\cdots \\& \qquad {}+\frac{x^{s-2}}{(s-2)!} \bigl[ L(y_{2})+L^{(1)}(y_{1}) \bigr] +\frac {x^{s-1}}{(s-1)!}L(y_{1}),\quad s=2,3,\ldots,m, \forall x\in \mathbb{R}, \end{aligned}$$

hold. Therefore, when the equations

$$\begin{aligned}& L(y_{1})=0, \\& L(y_{2})+L^{(1)}(y_{1})=0, \\& L(y_{3})+L^{(1)}(y_{2})+\frac{1}{2!}L^{(2)}(y_{1})=0, \\& \ldots, \\& L(y_{m})+L^{(1)}(y_{m-1})+\frac{1}{2!}L^{(2)}(y_{m-2})+ \cdots+\frac {1}{(m-1)!}L^{(m-1)}(y_{1})=0, \end{aligned}$$

(15)

are satisfied, the functions

$$ \tilde{y}_{1}=y_{1},\qquad \tilde{y}_{j}= \frac{x^{j-1}}{(j-1)!} y_{1}+\frac{x^{j-2}}{(j-2)!}y_{2}+\cdots+ \frac{x}{1!}y_{j-1}+y_{j},\quad j=2,3, \ldots,m, $$

(16)

are solutions of the equation $L(y)=0$.

Let us show the existence of functions $y_{s}=e^{\omega\lambda x}q_{s}(x,\lambda)$, $s=1,2,\ldots,m$, satisfying the system of equations (15) for the operators $L=L_{\lambda}$, $L^{(k)}=L_{\lambda}^{(k)}$, $k=1,2,\ldots,m$. If we set in these equations $L=L_{\lambda}$ and $y_{1}=\hat{f}_{1}(x,\lambda)$, the solution $y_{2}=e^{\omega\lambda x}q_{2}(x,\lambda )$, which satisfies the equation

$$ L_{\lambda}(y_{2})+L_{\lambda}^{(1)}(y_{1})=0, $$

exists according to Theorem 2 for $g(x,\lambda)=-L_{\lambda}^{(1)}(y_{1})$. It is not difficult to verify that the conditions of Theorem 2 are satisfied. In the same manner, when the functions $y_{s}=e^{\omega\lambda x}p_{s}(x,\lambda)$, $s=1,2,\ldots,k-1$, were found, the existence of the function $y_{k}=e^{\omega\lambda x}q_{k}(x,\lambda)$ which satisfies the equation

$$ L_{\lambda}(y_{k})=-L_{\lambda}^{(1)}(y_{k-1})- \frac{1}{2!}L_{\lambda }^{(2)}(y_{k-2})-\cdots- \frac{1}{(k-1)!}L_{\lambda}^{(k-1)}(y_{1}),\quad k=2,3, \ldots,m, $$

is obtained according to Theorem 2 by induction for $g(x,\lambda )=-L_{\lambda}^{(1)}(y_{k-1})-\frac{1}{2!}L_{\lambda }^{(2)}(y_{k-2})-\cdots-\frac{1}{(k-1)!}L_{\lambda}^{(k-1)}(y_{1})$. Consequently, according to (15), (16) the functions

$$\begin{aligned}& \hat{f}_{1}(x,\lambda) = e^{\omega\lambda x}q_{1}(x,\lambda ), \qquad \hat{f}_{2}(x,\lambda)=e^{\omega\lambda x} \bigl[ xq_{1}(x,\lambda )+q_{2}(x,\lambda) \bigr] , \\& \ldots, \\& \hat{f}_{m}(x,\lambda) = e^{\omega\lambda x}\biggl[\frac {x^{m-1}}{(m-1)!}q_{1}(x, \lambda)+\cdots+xq_{m-1}(x,\lambda)+q_{m}(x,\lambda)\biggr] \end{aligned}$$

are solutions of equation (3) in $( -\infty,+\infty )$ for $\lambda\in \mathbb{C}$. The theorem is proved. □

Note that the solutions of type $\hat{f}_{s}(x,\lambda )$, $s=1,2,\ldots,m$, are obtained in [7] under the different conditions and in various form of the representation.

Corollary 4

When $\phi(z)= ( z-\omega_{0} ) ^{m}$, for each $\lambda\in \mathbb{C}$, $x\in \mathbb{R}$ Wronskian of functions $\hat{f}_{1}(x,\lambda),\hat {f}_{2}(x,\lambda),\ldots,\hat{f}_{m}(x,\lambda)$ is found as

$$ \widehat{W}(x,\lambda)=e^{m\omega_{0}\lambda x-\sum_{n=1}^{\infty}\frac{p_{10n}}{i\alpha_{n}}e^{i\alpha_{n}x}}\neq0 $$

(17)

and hence for each $\lambda\in \mathbb{C}$, the functions $\hat{f}_{1}(x,\lambda),\hat {f}_{2}(x,\lambda ),\ldots,\hat{f}_{m}(x,\lambda)$ form the fundamental system of solutions of equation (3) in the interval $( -\infty ,+\infty )$.

3 The spectrum and resolvent of the operator $L_{\lambda}$

Here we investigate the structure of the spectrum of the operator $L_{\lambda}$ and the resolvent operator $L_{\lambda}^{-1}$.

Theorem 4

The operator $L_{\lambda}$ does not have eigenvalues, i.e. $\sigma_{p}(L_{\lambda})=\varnothing$.

Proof

Let us show that equation $L_{\lambda}y=0$ has only a trivial solution which belongs to $L_{2}( \mathbb{R})$ for $\forall\lambda\in \mathbb{C}$. In the case of simple roots of characteristic polynomial for $\lambda\neq0$, $\lambda\neq\lambda_{sjn}$, $s,j=1,2,\ldots,m$, $j\neq s$, $n\in \mathbb{N}$, it follows from the properties of the solutions $f_{k}(x,\lambda)$, $k=1,2,\ldots,m$. Really, the solution $y(x,\lambda)=c_{1}f_{1}(x,\lambda )+c_{2}f_{2}(x,\lambda)+\cdots+c_{m}f_{m}(x,\lambda)$ is in $L_{2}(\mathbb{R})$ if and only if $c_{1}=c_{2}=\cdots=c_{m}=0$. If we take linearly independent solutions of (3) according to $\lambda=0$ or $\lambda=\lambda_{sjn}$, then a similar result is also valid. Hence $\sigma_{p}(L_{\lambda })=\varnothing$. The theorem is proved. □

Theorem 5

The residual spectrum of the operator $L_{\lambda}$ is an empty set, i.e. $\sigma_{r}(L_{\lambda })=\varnothing$.

Proof

Since $\sigma _{p}(L_{\lambda })=\varnothing$ and for every $\lambda \in \mathbb{C}$ the operator $L_{\lambda }$ is one to one, $\lambda \in \sigma _{r}(L_{\lambda })$ if and only if the range $R(L_{\lambda })$ is not dense in $L_{2}( \mathbb{R})$. It means the equation $L_{\lambda}^{\ast}(z)=0$ has a nontrivial solution $z(x,\lambda)\in L_{2}( \mathbb{R})$, in other words $\overline{z(x,\lambda)}$ satisfies the conjugate equation

$$ (-1)^{m}z^{(m)}+\sum_{\gamma=1}^{m}(-1)^{(m-\gamma)} \bigl[ p_{\gamma}(x,\lambda)z \bigr] ^{(m-\gamma)}=0. $$

(18)

Since (18) can be written as

$$ z^{(m)}+\sum_{\gamma=1}^{m}p_{\gamma}^{\ast}(x, \lambda )z^{(m-\gamma)}=0 $$

(19)

which is in type of equation (3), equation (18) does not have a nontrivial solution which belongs to $L_{2}( \mathbb{R})$. Therefore $\overline{z(x,\lambda)}\equiv0$ and $\sigma _{p}(L_{\lambda }^{\ast})=\varnothing$, which means $\overline{R(L_{\lambda})}=L_{2}( \mathbb{R})$, or $\sigma_{r}(L_{\lambda})=\varnothing$. The theorem is proved. □

From Theorem 4 and Theorem 5 it follows that $\sigma(L_{\lambda })=\sigma _{c}(L_{\lambda})$ and $L_{\lambda}^{-1}$ is defined in a dense set in $L_{2}(\mathbb{R})$ for each$\lambda\in \mathbb{C}$.

In order to find $L_{\lambda}^{-1}$and the resolvent set $\rho (L_{\lambda })$, let us investigate the existence of the solution $y(x,\lambda)\in L_{2}( \mathbb{R})$ of the equation

$$ y^{(m)}+\sum_{\gamma=1}^{m}p_{\gamma}(x, \lambda)y^{(m-\gamma )}=f(x) $$

(20)

when $f(x)\in L_{2}( \mathbb{R})$.

Let us consider cases I and II separately.

In case I if we apply the Lagrange method by using the Floquet solutions $f_{s}(x,\lambda)$, $s=1,2,\ldots,m$, of equation (3) and linearly independent solutions (see [13], pp.208-210)

$$ z_{s}(x,\lambda)=(-1)^{m+s}\frac{W [ f_{1},f_{2},\ldots,f_{s-1},f_{s+1},\ldots,f_{m} ] (x,\lambda)}{W [ f_{1},f_{2},\ldots,f_{m} ] (x,\lambda)}, \quad s=1,2, \ldots,m, $$

of (18) for $\lambda\in \mathbb{C}\backslash\Lambda_{0}$, then we find the solution of (20) as

$$ y(x,\lambda)= \int_{-\infty}^{+\infty}G(x,t,\lambda)f(t)\, dt, $$

where the expression of $G(x,t,\lambda)$ can be written explicitly via $f_{s}(x,\lambda)$, $z_{s}(x,\lambda)$, $s=1,2,\ldots,m$. Using the properties of the functions $f_{s}(x,\lambda)$ we can show that

$$\begin{aligned}& W [ f_{1},f_{2},\ldots,f_{s-1},f_{s+1}, \ldots,f_{m} ] (x,\lambda) \\& \quad = \left \vert \textstyle\begin{array}{@{}c@{\quad}c@{\quad}c@{\quad}c@{\quad}c@{\quad}c@{}} f_{1}(x,\lambda) & \cdots& f_{s-1}(x,\lambda) & f_{s+1}(x,\lambda) & \cdots& f_{m}(x,\lambda) \\ f_{1}^{\prime}(x,\lambda) & \cdots& f_{s-1}^{\prime}(x,\lambda) & f_{s+1}^{\prime}(x,\lambda) & \cdots& f_{m}^{\prime}(x,\lambda) \\ f_{1}^{\prime\prime}(x,\lambda) & \cdots& f_{s-1}^{\prime\prime }(x,\lambda) & f_{s+1}^{\prime\prime}(x,\lambda) & \cdots& f_{m}^{\prime\prime}(x,\lambda) \\ \cdots& \cdots& \cdots& \cdots& \cdots& \cdots \\ f_{1}^{(m-2)}(x,\lambda) & \cdots& f_{s-1}^{(m-2)}(x,\lambda) & f_{s+1}^{(m-2)}(x,\lambda) & \cdots& f_{m}^{(m-2)}(x,\lambda)\end{array}\displaystyle \right \vert \\& \quad =e^{-(p_{11}+\omega_{s})\lambda x} \Biggl( A^{(s)}(\lambda)+{\sum _{n=1}^{\infty}A_{n}^{(s)}( \lambda)e^{i\alpha_{n}x}} \Biggr), \end{aligned}$$

where $A^{(s)}(\lambda)$, ${A_{n}^{(s)}(\lambda)}$ are complex valued functions of λ for which the series ${\sum_{n=1}^{\infty }}\vert {A_{n}^{(s)}(\lambda)}\vert$ is convergent and

$$\begin{aligned} (-1)^{m+s}A^{(s)}(\lambda) =&(-1)^{m+s}\left \vert \textstyle\begin{array}{@{}c@{\quad}c@{\quad}c@{\quad}c@{\quad}c@{\quad}c@{}} 1 & \cdots& 1 & 1 & \cdots& 1 \\ \omega_{1}\lambda& \cdots& \omega_{s-1}\lambda& \omega _{s+1}\lambda& \cdots& \omega_{m}\lambda \\ (\omega_{1}\lambda)^{2} & \cdots& (\omega_{s-1}\lambda)^{2} & (\omega _{s+1}\lambda)^{2} & \cdots& (\omega_{m}\lambda)^{2} \\ \cdots& \cdots& \cdots& \cdots& \cdots& \cdots \\ (\omega_{1}\lambda)^{m-2} & \cdots& (\omega_{s-1}\lambda)^{m-2} & (\omega_{s+1}\lambda)^{m-2} & \cdots& (\omega_{m}\lambda)^{m-2}\end{array}\displaystyle \right \vert \\ =&(-1)^{m+s}\lambda^{\frac{(m-2)(m-1)}{2}}\left \vert \textstyle\begin{array}{@{}c@{\quad}c@{\quad}c@{\quad}c@{\quad}c@{\quad}c@{}} 1 & \cdots& 1 & 1 & \cdots& 1 \\ \omega_{1} & \cdots& \omega_{s-1} & \omega_{s+1} & \cdots& \omega_{m} \\ \omega_{1}^{2} & \cdots& \omega_{s-1}^{2} & \omega_{s+1}^{2} & \cdots& \omega_{m}^{2} \\ \cdots& \cdots& \cdots& \cdots& \cdots& \cdots \\ \omega_{1}^{m-2} & \cdots& \omega_{s-1}^{m-2} & \omega_{s+1}^{m-2} & \cdots& \omega_{m}^{m-2}\end{array}\displaystyle \right \vert \\ =&\lambda^{\frac{(m-2)(m-1)}{2}}W_{ms}, \end{aligned}$$

where $W_{ms}$, $s=1,2,\ldots,m$, are cofactors of elements of mth row of the Vandermonde determinant $W_{m}$.

Thus we find that

$$\begin{aligned} z_{s}(x,\lambda) =&(-1)^{m+s}\frac{e^{-(p_{11}+\omega_{s})\lambda x} ( A^{(s)}(\lambda)+{\sum_{n=1}^{\infty}A_{n}^{(s)}(\lambda)e^{i\alpha _{n}x}} ) }{\lambda^{\frac{m(m-1)}{2}}W_{m}}e^{p_{11}\lambda x+\sum_{n=1}^{\infty}\frac{p_{10n}}{i\alpha_{n}}e^{i\alpha_{n}x}} \\ =&(-1)^{m+s}\frac{e^{-i\omega_{s}\lambda x} ( A^{(s)}(\lambda)+{\sum_{n=1}^{\infty}A_{n}^{(s)}(\lambda)e^{i\alpha_{n}x}} ) }{\lambda ^{\frac{m(m-1)}{2}}W_{m}}e^{\sum_{n=1}^{\infty}\frac{p_{10n}}{i\alpha_{n}}e^{i\alpha_{n}x}},\quad s=1,2,\ldots,m. \end{aligned}$$

On the other hand, it can be seen that equation (18) which is equivalent to (19) is in the type of equation (3) and, moreover, the characteristic polynomial corresponding to the expression

$$ \ell_{\lambda}^{\ast}(\lambda)=z^{(m)}+\sum _{\gamma=1}^{m}p_{\gamma }^{\ast}(x, \lambda)z^{(m-\gamma)} $$

is in the form of $\phi^{\ast}(\lambda)=(-1)^{m}\phi(-\lambda)$. Then, if $\omega_{1},\omega_{2},\ldots,\omega_{m}$ are roots of $\phi (\lambda)$, then $-\omega_{1},-\omega_{2},\ldots,-\omega_{m}$ will be roots of $\phi ^{\ast}(\lambda)$. According to Theorem 1, for all $\lambda\in \mathbb{C}\backslash\Lambda_{0}$ equation (19) has solutions as

$$ \varphi_{s}(x,\lambda)=e^{-\omega_{s}\lambda x} \Biggl( {1+\sum _{n=1}^{\infty}V_{n}^{(s)}( \lambda)e^{i\alpha_{n}x}} \Biggr) , \quad s=1,2,\ldots,m, $$

where

$$ V_{n}^{(s)}(\lambda)=V_{n}^{(s)}+ \sum _{k=1}^{n}\sum _{\substack{ j=1 \\ j\neq s}}^{m}\frac{V_{jkn}^{(s)}}{i\alpha_{k}-(\omega _{s}-\omega _{j})\lambda},\quad n\in \mathbb{N}. $$

Taking into consideration that for every point $\lambda\neq\lambda _{sjn}$, $s,j=1,2,\ldots,m$, $j\neq s$, each function $z_{s}(x,\lambda)$ is a linear combination of the functions $\varphi_{1}(x,\lambda),\varphi _{2}(x,\lambda),\ldots,\varphi_{m}(x,\lambda)$ and by virtue of the behavior of the function $z_{s}(x,\lambda)$ as $x\rightarrow\pm\infty$, it is possible to show that for some constants $C_{s}(\lambda)$ the equality

$$ z_{s}(x,\lambda)=C_{s}(\lambda)\varphi_{s}(x, \lambda),\quad s=1,2,\ldots,m, $$

is satisfied.

This means that for every $x\in\mathbb{R}$ the equality

$$\begin{aligned}& (-1)^{m+s}\frac{e^{-\omega_{s}\lambda x} ( A^{(s)}(\lambda)+{\sum_{n=1}^{\infty}A_{n}^{(s)}(\lambda)e^{i\alpha_{n}x}} ) }{\lambda ^{\frac{m(m-1)}{2}}W_{m}}e^{\sum_{n=1}^{\infty}\frac{p_{10n}}{i\alpha_{n}}e^{i\alpha_{n}x}} \\& \quad =C_{s}(\lambda) \varphi_{s}(x,\lambda),\quad s=1,2,\ldots,m, \end{aligned}$$

holds and hence dividing this equality by $e^{-\omega_{s}\lambda x}$ we obtain

$$\begin{aligned}& (-1)^{m+s}\frac{ ( A^{(s)}(\lambda)+{\sum_{n=1}^{\infty }A_{n}^{(s)}(\lambda)e^{i\alpha_{n}x}} ) }{\lambda^{\frac {m(m-1)}{2}}W_{m}}e^{\sum_{n=1}^{\infty}\frac{p_{10n}}{i\alpha_{n}}e^{i\alpha _{n}x}} \\& \quad =C_{s}(\lambda) \Biggl( {1+\sum_{n=1}^{\infty }V_{n}^{(s)}( \lambda)e^{i\alpha_{n}x}} \Biggr) , \quad s=1,2,\ldots,m. \end{aligned}$$

According to the uniqueness theorem of analytic functions, this equality is satisfied for the analytic continuations of the functions on the semi-plane $\operatorname{Im} x\geq0$ with respect to x, which are on both sides of the given equality.

Thus, taking the limit as $\operatorname{Im} x\rightarrow+\infty$ we have

$$ C_{s}(\lambda)=(-1)^{m+s}\frac{A^{(s)}(\lambda)}{\lambda^{\frac {m(m-1)}{2}}W_{m}}= \frac{\lambda^{\frac{(m-2)(m-1)}{2}}W_{ms}}{\lambda^{\frac {m(m-1)}{2}}W_{m}}=\frac{W_{ms}}{\lambda^{m-1}W_{m}}, \quad s=1,2,\ldots,m. $$

Therefore, the equalities

$$ z_{s}(x,\lambda)=\frac{W_{ms}}{\lambda^{m-1}W_{m}}\varphi _{s}(x, \lambda),\quad s=1,2,\ldots,m, $$

are valid.

For all $j=1,2,\ldots,m$, $k=1,2,\ldots,2m_{0}$, and $\forall\lambda\in S_{k}$ the value of $\operatorname{Re}(\lambda\omega_{j})$ has a constant sign. Hence, there are $M_{k}\subseteq \{ 1,2,\ldots,2m \}$ and $M_{k}^{\prime }= \{ 1,2,\ldots,2m \} \backslash M_{k}$ such that if $j\in M_{k}$, $\forall\lambda\in S_{k}$ then $\operatorname{Re}(\lambda\omega _{j})<0$, if $j\in M_{k}^{\prime}$ then $\operatorname{Re}(\lambda\omega_{j})>0$. Thus, for any $k=1,2,\ldots,2m_{0}$, $\forall a\in \mathbb{R}$, and $\forall\lambda\in S_{k}\backslash\Lambda$, if $l\in M_{k}$ then $f_{l}(x,\lambda)\in L_{2}(a,+\infty)$, $f_{l}(x,\lambda)\notin L_{2}(-\infty,a)$, and if $l\in M_{k}^{\prime}$ then $f_{l}(x,\lambda )\in L_{2}(-\infty,a)$, $f_{l}(x,\lambda)\notin L_{2}(a,+\infty)$. In the same manner, for $\forall a\in \mathbb{R}$, $\forall\lambda\in S_{k}\backslash\Lambda$ and $l\in M_{k}$ then $\varphi_{l}(x,\lambda)\notin L_{2}(a,+\infty)$, $\varphi _{l}(x,\lambda )\in L_{2}(-\infty,a)$, if $l\in M_{k}^{\prime}$ then $\varphi _{l}(x,\lambda)\notin L_{2}(-\infty,a)$, $\varphi_{l}(x,\lambda)\in L_{2}(a,+\infty)$. Taking into account all these properties, for $k=1,2,\ldots,2m_{0}$, $\forall\lambda\in S_{k}\backslash\Lambda$, the kernel $G(x,t,\lambda)$ can be written as

$$ G(x,t,\lambda)= \frac{1}{\lambda^{m-1}W_{m}}\left \{ \textstyle\begin{array}{l@{\quad}l} \sum_{l\in M_{k}}W_{ml}f_{l}(x,\lambda)\varphi_{l}(t,\lambda), & t\leq x, \\ {-}\sum_{l\in M_{k}^{\prime}}W_{ml}f_{l}(x,\lambda)\varphi _{l}(t,\lambda), & t>x.\end{array}\displaystyle \right . $$

(21)

From the expression of the functions $f_{s}(x,\lambda)$ and $\varphi _{s}(x,\lambda)$ and from (21) it follows that for every $x,t\in R$ and $\lambda\in S_{k}\backslash\Lambda$, $k=1,2,\ldots,2m_{0}$,

$$ \bigl\vert G(x,t,\lambda)\bigr\vert \leq C(\lambda)e^{-\tau \vert x-t\vert }, $$

(22)

where $C(\lambda)>0$, $\tau=\min \{ \vert \operatorname {Re}(\lambda \omega_{s})\vert :s=1,2,\ldots,m \}$. From (22) we have

$$ \int_{-\infty}^{+\infty}\bigl\vert G(x,t,\lambda)\bigr\vert ^{2}\, dt< +\infty\quad \text{and}\quad \int_{-\infty}^{+\infty}\bigl\vert G(x,t,\lambda)\bigr\vert ^{2}\, dx< +\infty. $$

Using (22) it can be proved by the standard method (see [11], pp.302-304) that the operator

$$ L_{\lambda}^{-1}f(x)= \int_{-\infty}^{+\infty}G(x,t,\lambda)f(t)\, dt $$

as $L_{\lambda}^{-1}:L_{2}( \mathbb{R} )\rightarrow L_{2}( \mathbb{R})$ is bounded for $\lambda\in S_{1}\cup S_{2}\cup\cdots\cup S_{2m_{0}}$ (it is $\lambda\in\rho(L_{\lambda})$). In the case $\lambda\in \bigcup_{k=1}^{m}l_{k}$ the operator $L_{\lambda}^{-1}$ is a closed operator defined on a dense proper subset $R(L_{\lambda})$ of $L_{2}( \mathbb{R})$ and so $L_{\lambda}^{-1}$ is an unbounded operator which means $\lambda\in \sigma_{c}(L_{\lambda})$. On the other hand, since the functions $f_{l}(x,\lambda)$ and $\varphi_{l}(t,\lambda)$ do not have the same poles the points $\lambda=\lambda_{sjn}$ may be simple poles of $G(x,t,\lambda)$. If any $\lambda=\lambda_{sjn}$ belong to any set $S_{k}$ then these points can be only eigenvalues of the operator $L_{\lambda}$. Since $L_{\lambda}$ does not have eigenvalue, there is no singularity of operator $L_{\lambda}^{-1}$ at these points. Therefore $\lambda_{sjn}\in\rho (L_{\lambda})$ and $G(x,t,\lambda)$ is regular at these points too. So the lines $\operatorname{Re}(\lambda\omega_{s})=0$, $s=1,2,\ldots,m$, consist of a continuous spectrum of $L_{\lambda}$, i.e. $\sigma _{c}(L_{\lambda })=\bigcup_{k=1}^{m}l_{k}$. Then the resolvent set of operator $L_{\lambda}$ is $\rho(L_{\lambda})=S_{1}\cup S_{2}\cup\cdots\cup S_{2m_{0}}$. An analytic continuation of the kernel $G(x,t,\lambda)$, with respect to λ, may have some poles at the points $\lambda\in\Lambda _{0}$, which belong to $\sigma_{c}(L_{\lambda})$. These poles are called spectral singularities (in the sense of [11], p.306) of the operator $L_{\lambda }$.

Similarly, in case II the kernel $G(x,t,\lambda)$ can be written as

$$ G(x,t,\lambda)=\left \{ \textstyle\begin{array}{l@{\quad}l} \sum_{k=1}^{m}\hat{f}_{k}(x,\lambda)\hat{\varphi}_{k}(t,\lambda), & t< x, \\ 0, & t\geq x \end{array}\displaystyle \right . $$

for $\lambda\in S^{-}$ or as

$$ G(x,t,\lambda)=\left \{ \textstyle\begin{array}{l@{\quad}l} 0, & t< x, \\ {-}\sum_{k=1}^{m}\hat{f}_{k}(x,\lambda)\hat{\varphi}_{k}(t,\lambda), & t\geq x \end{array}\displaystyle \right . $$

for $\lambda\in S^{+}$, where $\hat{f}_{k}(x,\lambda)$, $k=1,2,\ldots,m$, are solutions of equation (3) and $\hat{\varphi }_{k}(x,\lambda)$ are solutions of (11) according to Theorem 4. In order to have $y(x,\lambda )=\int_{-\infty}^{+\infty}G(x,t,\lambda)f(t)\, dt\in L_{2}(\mathbb{R})$, for each $\lambda\in \mathbb{C}$, $\operatorname{Re}(\omega_{0}\lambda)\neq0$ and $f(x)\in L_{2}( \mathbb{R})$ we use

$$ \sum_{k=1}^{m}\hat{f}_{k}(x, \lambda)\hat{\varphi}_{k}(t,\lambda)=\frac{H(x,t,\lambda)}{W [ \hat{f}_{1},\hat {f}_{2},\ldots,\hat{f}_{m} ] (t)}. $$

Here

$$\begin{aligned}& H(x,t,\lambda) \\& \quad =\left \vert \textstyle\begin{array}{@{}c@{\quad}c@{\quad}c@{\quad}c@{}} \hat{f}_{1}(t,\lambda) & \hat{f}_{2}(t,\lambda) & \cdots & \hat{f}_{m}(t,\lambda) \\ \hat{f}_{1}^{\prime}(t,\lambda) & \hat{f}_{2}^{\prime }(x,\lambda) & \cdots & \hat{f}_{m}^{\prime}(t,\lambda) \\ \cdots & \cdots & \cdots & \cdots \\ \hat{f}_{1}^{(m-2)}(t,\lambda) & \hat {f}_{2}^{(m-2)}(x,\lambda) & \cdots & \hat{f}_{m}^{(m-2)}(t,\lambda) \\ \hat{f}_{1}(x,\lambda) & \hat{f}_{2}(x,\lambda) & \cdots & \hat{f}_{m}(x,\lambda) \end{array}\displaystyle \right \vert \\& \quad =\left \vert \textstyle\begin{array}{@{}c@{\quad}c@{\quad}c@{\quad}c@{}} e^{\omega\lambda t}q_{1}(t,\lambda) & e^{\omega\lambda t} [ tq_{1}(t,\lambda)+q_{2}(t,\lambda) ] & \cdots & e^{\omega \lambda t}\sum_{j=1}^{m}\frac{t^{m-j}}{(m-j)!}q_{j}(t,\lambda) \\ ( e^{\omega\lambda t}q_{1}(t,\lambda) ) ^{\prime} & ( e^{\omega\lambda t} [ tq_{1}(t,\lambda)+q_{2}(t,\lambda) ] ) ^{\prime} & \cdots & ( e^{\omega\lambda t}\sum_{j=1}^{m}\frac{t^{m-j}}{(m-j)!}q_{j}(t,\lambda) ) ^{\prime} \\ \cdots & \cdots & \cdots & \cdots \\ ( e^{\omega\lambda t}q_{1}(t,\lambda) ) ^{(m-2)} & ( e^{\omega\lambda t} [ tq_{1}(t,\lambda)+q_{2}(t,\lambda) ] ) ^{(m-2)} & \cdots & ( e^{\omega\lambda t}\sum_{j=1}^{m}\frac{t^{m-j}}{(m-j)!}q_{j}(t,\lambda) ) ^{(m-2)} \\ e^{\omega\lambda x}q_{1}(x,\lambda) & e^{\omega\lambda x} [ xq_{1}(x,\lambda)+q_{2}(x,\lambda) ] & \cdots & e^{\omega \lambda x}\sum_{j=1}^{m}\frac{x^{m-j}}{(m-j)!}q_{j}(x,\lambda) \end{array}\displaystyle \right \vert . \end{aligned}$$

When we make elementary operations on rows and columns of this determinant, we can transform this determinant as follows:

$$\begin{aligned}& H(x,t,\lambda) \\& \quad =e^{\omega\lambda [ (m-1)t+x ] }\left \vert \textstyle\begin{array}{@{}c@{\quad}c@{\quad}c@{\quad}c@{}} q_{1}(t,\lambda) & tq_{1}(t,\lambda)+q_{2}(t,\lambda) & \cdots & \sum_{j=1}^{m}\frac{t^{m-j}}{(m-j)!}q_{j}(t,\lambda) \\ q_{1}^{\prime}(t,\lambda) & [ tq_{1}(t,\lambda )+q_{2}(t,\lambda) ] ^{\prime} & \cdots & ( \sum_{j=1}^{m}\frac{t^{m-j}}{(m-j)!}q_{j}(t,\lambda) ) ^{\prime} \\ \cdots & \cdots & \cdots & \cdots \\ q_{1}^{(m-2)}(t,\lambda) & [ tq_{1}(t,\lambda )+q_{2}(t,\lambda) ] ^{(m-2)} & \cdots & ( \sum_{j=1}^{m}\frac{t^{m-j}}{(m-j)!}q_{j}(t,\lambda) ) ^{(m-2)} \\ q_{1}(x,\lambda) & xq_{1}(x,\lambda)+q_{2}(x,\lambda) & \cdots & \sum_{j=1}^{m}\frac{x^{m-j}}{(m-j)!}q_{j}(x,\lambda) \end{array}\displaystyle \right \vert \\& \quad =e^{\omega\lambda [ (m-1)t+x ] }\left \vert \textstyle\begin{array}{@{}c@{\quad}c@{\quad}c@{\quad}c@{}} q_{1}(t,\lambda) & q_{2}(t,\lambda) & \cdots & q_{m}(t,\lambda) \\ q_{1}^{\prime}(t,\lambda) & q_{2}^{\prime}(t,\lambda )+q_{1}(t,\lambda ) & \cdots & q_{m}^{\prime}(t,\lambda)+q_{m-1}(t,\lambda) \\ \cdots & \cdots & \cdots & \cdots \\ q_{1}^{(m-2)}(t,\lambda) & q_{2}^{(m-2)}(t,\lambda )+(m-2)q_{1}(t,\lambda) & \cdots & \sum_{j=0}^{m-2}C_{m-2}^{j}q_{m-j}^{(m-j-2)}(t,\lambda) \\ q_{1}(x,\lambda) & (x-t)q_{1}(x,\lambda)+q_{2}(x,\lambda) & \cdots & \sum_{j=1}^{m}\frac{(x-t)^{m-j}}{(m-j)!}q_{j}(x,\lambda) \end{array}\displaystyle \right \vert \\& \quad = e^{\omega\lambda [ (m-1)t+x ] } \biggl( c_{1}(x,t,\lambda )+(x-t)c_{2}(x,t, \lambda)+\cdots+\frac {(x-t)^{m-1}}{(m-1)!}c_{m}(x,t,\lambda ) \biggr) . \end{aligned}$$

(23)

Here, $c_{i}(x,t,\lambda)$, $i=1,2,\ldots,m$, are continuous and bounded functions of $(x,t)$ in $\mathbb{R}^{2}$ for each constant λ. According to equations (17), (23), for $\tau=\operatorname{Re}(\omega_{0}\lambda)$ and some constant $C^{\prime}(\lambda)>0$, we obtain the inequality

$$ \bigl\vert G(x,t,\lambda)\bigr\vert \leq C^{\prime}( \lambda)e^{-\tau \vert x-t\vert }\bigl(1+\vert x-t\vert \bigr)^{m-1},\quad \forall x,t\in R, \forall\lambda\in \mathbb{C} . $$

(24)

From (24) we have

$$ \int_{-\infty}^{+\infty}\bigl\vert G(x,t,\lambda)\bigr\vert ^{2}\, dt< +\infty\quad \text{and}\quad \int_{-\infty}^{+\infty}\bigl\vert G(x,t,\lambda)\bigr\vert ^{2}\, dx< +\infty. $$

By considering inequality (24) it can be proved by the standard method (see [11], pp.302-304) that the operator

$$ L_{\lambda}^{-1}f(x)= \int_{-\infty}^{+\infty}G(x,t,\lambda)f(t)\, dt $$

as $L_{\lambda}^{-1}:L_{2}( \mathbb{R} )\rightarrow L_{2}( \mathbb{R})$ is bounded for $\lambda\in \mathbb{C}$, $\operatorname{Re}(\lambda\omega_{0})\neq0$ (it is $\lambda\in\rho(L_{\lambda })$). If $\operatorname{Re}(\lambda\omega_{0})=0$, then, as in case I, the operator $L_{\lambda}^{-1}$ is a closed operator defined on a dense proper subset $R(L_{\lambda})$ of $L_{2}( \mathbb{R})$ and so $L_{\lambda}^{-1}$ is an unbounded operator, which means $\lambda\in\sigma_{c}(L_{\lambda})$. From the expression of $G(x,t,\lambda)$ we see that $G(x,t,\lambda)$ is a holomorphic function in $S^{+}$ and $S^{-}$. Analytic continuation of the function $G(x,t,\lambda)$ with respect to λ out of sectors $S^{+}$ and $S^{-}$ does not have a singularity on the line $\operatorname{Re}(\lambda\omega_{0})=0$. Therefore, the line $\operatorname {Re}(\lambda\omega _{0})=0$ consists of a continuous spectrum of $L_{\lambda}$, i.e. $\sigma _{c}(L_{\lambda})=l_{0}$ and does not have a spectral singularity of this operator. Then the resolvent set of the operator $L_{\lambda}$ is $\rho (L_{\lambda})=S^{+}\cup S^{-}$. Thus, the following theorem is true.

Theorem 6

The operator $L_{\lambda}$ has a pure continuous spectrum $\sigma_{c}(L_{\lambda})$. If the characteristic polynomial $\phi (z)$ has simple roots $\omega_{s}$, $s=1,2,\ldots,m$, then it is made up of lines $\operatorname{Re}(\lambda\omega_{s})=0$, $s=1,2,\ldots,m$. The countable set of simple spectral singularities may exist at the points $\lambda=\lambda _{sjn}\in \sigma_{c}(L_{\lambda})$ and the spectral singularity degree, which does not exceed $m-1$, may exist at the point $\lambda=0$. If the characteristic polynomial $\phi(z)$ has a unique multiple root $\omega_{0}$, then a continuous spectrum consists of the line $\operatorname{Re}(\lambda\omega_{0})=0$ and spectral singularities of $L_{\lambda}$ do not exist. For any $\lambda \in\rho(L_{\lambda})$ the resolvent $L_{\lambda}^{-1}$ is an integral operator with a kernel of Karleman type.