
\documentclass[twoside]{article}
\usepackage{amsfonts, amsmath} % used for R in Real numbers
\pagestyle{myheadings}

\markboth{\hfil A selfadjoint hyperbolic boundary-value problem \hfil EJDE/Conf/10}
{EJDE/Conf/10 \hfil Nezam Iraniparast \hfil}

\begin{document}
\setcounter{page}{153}
\title{\vspace{-1in}\parbox{\linewidth}{\footnotesize\noindent
Fifth Mississippi State Conference on Differential Equations and
Computational Simulations, \newline
Electronic Journal of Differential Equations,
Conference 10, 2003, pp 153--161. \newline
http://ejde.math.swt.edu or http://ejde.math.unt.edu
\newline ftp  ejde.math.swt.edu (login: ftp)}
 \vspace{\bigskipamount} \\
%
  A selfadjoint hyperbolic boundary-value problem
%
\thanks{ {\em Mathematics Subject Classifications:} 35L05, 35L20, 35P99.
\hfil\break\indent
{\em Key words:} Characteristics, eigenvalues,
 eigenfunctions, Green's function, \hfil\break\indent
 Fredholm alternative.
\hfil\break\indent
\copyright 2003 Southwest Texas State University. \hfil\break\indent
Published February 28, 2003.} }

\date{}
\author{Nezam Iraniparast}
\maketitle

\begin{abstract}
 We consider the eigenvalue wave equation
 $$u_{tt} - u_{ss} = \lambda pu,$$
 subject to $ u(s,0) = 0$,
 where $u\in\mathbb{R}$, is a function of $(s, t) \in \mathbb{R}^2$,
 with $t\ge 0$.
 In the characteristic triangle $T =\{(s,t):0\leq t\leq 1, t\leq s\leq 2-t\}$
 we impose a boundary condition along characteristics so that
 $$
 \alpha u(t,t)-\beta \frac{\partial u}{\partial n_1}(t,t) = \alpha u(1+t,1-t)
 +\beta\frac{\partial u}{\partial n_2}(1+t,1-t),\quad 0\leq t\leq1.
 $$
 The parameters $\alpha$ and $\beta$ are arbitrary except for the condition that
 they are not both zero. The two vectors $n_1$ and $n_2$ are the exterior unit normals
 to the characteristic boundaries and $\frac{\partial u}{\partial n_1}$,
 $\frac{\partial u}{\partial n_2}$ are the normal derivatives in those directions. When
 $p\equiv 1$ we will show that the above characteristic boundary value  problem
 has real, discrete eigenvalues and corresponding eigenfunctions that are complete and
 orthogonal in $L_2(T)$. We will also investigate the case where $p\geq 0$ is an arbitrary
 continuous function in $T$.
\end{abstract}

\numberwithin{equation}{section}
\newtheorem{theorem}{Theorem}[section]

\section{Introduction}

 Consider the wave equation
\begin{equation}
 u_{tt}-u_{ss}=\lambda pu,\quad (s,t) \in
T:=\{(s,t):0 \le t\le 1, t\le s\le 2-t\}, \label{e1}
\end{equation}
where $\lambda$ is a parameter, and $p\geq 0$ is a continuous function in $T$.
We impose the boundary conditions,
\begin{gather}
u(s,0)=0,\quad 0\le s\le 2,  \label{e2}\\
\alpha u(t,t)-\beta \frac{\partial u}{\partial n_1}(t,t) = \alpha
u(1+t,1-t)
+ \beta\frac{\partial u}{\partial n_2}(1+t,1-t),\;\; 0\leq t\leq1, \label{e3}
\end{gather}
where the parameters $\alpha$ and $\beta$ are arbitrary with
$\alpha^2 +\beta^2 \neq 0$. The two vectors $n_1$ and $n_2$ are the exterior unit normals
to the characteristic boundaries and $\frac{\partial u}{\partial n_1}$,
$\frac{\partial u}{\partial n_2}$ are the normal derivatives in those directions. Our
problem is a generalization of the problems studied by Kalmenov [6],
and Kreith \cite{kr:symmetric}, where they consider the boundary conditions,
\begin{gather*}
u(s,0)=0=u(1,1),\quad 0\le s \le 2\\
u(t,t)=u(1+t,1-t), \quad 0\le t \le 1.
\end{gather*}
We note two facts. First, if we only prescribe the
values of $u$ along the characteristics, say, $u(t,t)=f_1(t)$ and $u(1+t,1-t)=f_2(t)$ then
we have a classical characteristic initial value problem,(see, e.g., Garabedian [1]) and
equation \eqref{e1} will have a solution for all values of $\lambda$. However, the conditions \eqref{e2},
\eqref{e3} provide a boundary value problem with spectral properties. Second, if we set $\beta =
0$ in condition \eqref{e3} we have the case of Kreith \cite{kr:symmetric} . If in addition, we set
$p\equiv 1$ we have the case of Kalmenov [6]. As a start we will let
$p\equiv 1$ and use the technique of Kalmenov [6] to show that the equation,
$u_{tt}-u_{ss}=\lambda u,\quad (s,t)
\in T$ subject to conditions \eqref{e2}, \eqref{e3}, is selfadjoint. Then we will study equation \eqref{e1}
subject to conditions \eqref{e2}, \eqref{e3} by converting the problem into a nonhomogeneous eigenvalue
integral equation and using the method of the Fredholm alternative. In this case we will
assume that u along the characteristics is given. This will not weaken the problem because
we still require that condition \eqref{e2}  be satisfied. As a result we will see that not all
values of $\lambda$ will produce a solution.

We would like to add here that the investigation into the spectral properties of the
characteristic initial value problems for the wave equation has been conducted in
different directions by several authors. In this context, beside the above mentioned
references, the works of Haws [2], Kreith
\cite{kr:Mixed, kr:Establishing}
and the author \cite{Ir:A Method, Ir:A Boundary, Ir:A CIV} are the ones most closely
related to the present work.

\section{The selfadjoint problem}

Consider the equation
\begin{equation} \label{e4}
u_{tt}-u_{ss}=\lambda u,\quad (s,t) \in T,
\end{equation}
subject to the conditions \eqref{e2}, \eqref{e3}. Extend $u$ as an odd function of $t$ and
define,
\[
\tilde{u}(s,t)=\begin{cases}
 u(s,t)  &\mbox{if } t>0,\quad(s,t)\in T\\
 -u(s,-t) &\mbox{if } t<0,\quad(s,-t)\in T
\end{cases}
\]
and reflect $T$ in $t=0$ axis to obtain the rectangle
\begin{equation}
R = \{(s, t): |t| \leq 1,\; |t|\leq s \leq 2-|t|\}. \label{e5}
\end{equation}
Now $\tilde{u}$ must satisfy the following conditions:
\begin{gather}
\tilde{u}_{tt}-\tilde{u}_{ss}=\lambda \tilde{u},\quad (s,t) \in R, \label{e6}\\
\alpha \tilde{u}(t,t)-\beta \frac{\partial \tilde{u}}{\partial n_1}(t,t) = \alpha
\tilde{u}(1+t,1-t)
+ \beta\frac{\partial \tilde{u}}{\partial n_2}(1+t,1-t),\quad |t| \leq1, \label{e7}\\
\tilde{u}(s,t) = -\tilde{u}(s,-t), \quad (s, t)\in R. \label{e8}
\end{gather}
Using the change of variables, $ x = s - t$, $y = s + t$ and denoting $\tilde{u}(s,t) =
\tilde{u}(\frac{x+y}{2}, \frac{y-x}{2})$ by $\tilde{U}(x,y)$, the rectangle R maps into $S =
\{(x,y): 0\leq x\leq 2, 0\leq y\leq 2\}$, and we have,
\begin{gather}
-4\tilde{U}_{x,y}=\lambda\tilde{U},\quad (x,y)\in S,\label{e9}\\
\alpha \tilde{U}(0,y) + \beta \tilde{U}_x(0,y)=\alpha \tilde{U}(y,2) + \beta
\tilde{U}_y(y,2),\label{e10}\\
\tilde{U}(x, y) = - \tilde{U}(y,x).\label{e11}
\end{gather}
 From this equation, we have $\tilde{U}_x(x,y)=-\tilde{U}_y(y,x)$; therefore,
\begin{gather}
\tilde{U}(0,y) = -\tilde{U}(y, 0)\label{e12}\\
\tilde{U}_x(0,y) = -\tilde{U}_y(y,0). \label{e13}
\end{gather}
These two identities change the boundary condition \eqref{e10} to,
\begin{equation} \label{e14}
-\alpha \tilde{U}(y,0) - \beta \tilde{U}_y(y,0)=\alpha \tilde{U}(y,2) + \beta
\tilde{U}_y(y,2).
\end{equation}
Now let $\tilde{U}(x,y)=\phi(x)\psi(y)$, and plug it into
equation \eqref{e9}, we have,
\begin{equation} \label{e15}
[2i\phi'(x)][2i\psi'(y)]=\lambda\phi(x)\psi(y)\end{equation}
which leads to
$$\frac{2i\phi'(x)}{\phi(x)}=\lambda \frac{\psi(y)}{2i\psi'(y)}=\mu,$$
which in turn leads to two ODE equations:
\begin{gather}
\psi'(y) = -\frac{i\lambda}{2\mu}\psi(y) =-\frac{i\delta}{2}\psi(y),\quad
\delta=\frac{\lambda}{\mu}, \quad \mu\neq 0,\label{e16}\\
\phi'(x) = -\frac{i\mu}{2}\phi(x). \label{e17}
\end{gather}
The effect of separation of variables on the boundary condition \eqref{e14}
will be
\begin{equation}
-\alpha \phi(y)\psi(0) - \beta \phi(y)\psi'(0)=\alpha \phi(y)\psi(2) + \beta
\phi(y)\psi'(2).
\end{equation}
Upon cancelling $\phi(y)$ here, we have the boundary-value problem
\begin{gather}
\psi'(y) = -\frac{i\delta}{2}\psi(y), \quad 0\leq y\leq 2, \label{e19}\\
-\alpha \psi(0) - \beta \psi'(0)=\alpha \psi(2) + \beta \psi'(2). \label{e20}
\end{gather}
We note that from \eqref{e19} we have
\begin{gather}
\psi'(0)=-\frac{i\delta}{2}\psi(0), \label{e21}\\
\psi'(2)=- \frac{i\delta}{2}\psi(2). \label{e22}
\end{gather}
These relations simplify the condition \eqref{e20} to
\begin{equation} \label{e23}
\psi'(0)=- \psi(2),\end{equation}
The selfadjoint boundary-value problem \eqref{e19}, \eqref{e23} is the same as the one in
\cite{Ka:spectrum}. The eigenvalues and eigenfunctions of problem \eqref{e19}, \eqref{e23}
are
\begin{equation} \label{e24}
\delta = (-2n - 1)\pi, \quad \psi_n(y)=\exp(\frac{(2n+1)\pi iy}{2}),\quad
n=0,\pm1,\pm2,\dots
\end{equation}
Using \eqref{e11} differently to write $\tilde{U}_y(x, y) = -\tilde{U}_x(y,x)$,
\begin{gather}
\tilde{U}(y,2) = -\tilde{U}(2, y) \label{e25}\\
\tilde{U}_y(y,2) = -\tilde{U}_x(2,y).\label{e26}
\end{gather}
The substitution of \eqref{e25}, \eqref{e26} into the boundary condition \eqref{e10} yields
\begin{equation} \label{e27}
\alpha \tilde{U}(0,y) + \beta \tilde{U}_x(0,y)=-\alpha \tilde{U}(2,y) - \beta
\tilde{U}_x(2,y).
\end{equation}
If we again use the separation of variables
$\tilde{U}(x,y)=\phi(x)\psi(y)$ in the equation \eqref{e27} and proceed as for the case of the
function $\psi(y)$, we will have
\begin{gather}
\phi'(x) = -\frac{i\mu}{2}\phi(x), \quad 0\leq x\leq 2, \label{e28}\\
-\alpha \phi(0) - \beta \phi'(0)=\alpha \phi(2) + \beta \phi'(2).\label{e29}
\end{gather}
The condition \eqref{e29} again can be replaced with
\begin{equation} \label{e30}
\phi'(0)=- \phi(2),\end{equation}
The eigenvalues and eigenfunctions of problem \eqref{e28}, \eqref{e30} are
\begin{equation} \label{31}
\mu = (-2m - 1)\pi, \quad \phi_m(y)=\exp(\frac{(2m+1)\pi ix}{2}),\quad
m=0,\pm1,\pm2,\dots
\end{equation}
Both sets of eigenfunctions $\{\psi_n(y)\}$ and $\{\phi_m(x)\}$ are complete and
orthogonal in $L_2(0,2)$. The eigenfunctions and the eigenvalues of \eqref{e9}--\eqref{e11}
will be
\begin{equation}
\begin{gathered}
\lambda_{mn} = (2m + 1)(2n+2)\pi^2,\\
\tilde{U}_{mn}(x,y)=\exp(\frac{(2m+1)x}{2}+\frac{(2n+1)y}{2})\pi i,
\quad m,n=0,\pm1,\pm2,\dots
\end{gathered} \label{e32}
\end{equation}
To find the eigenfunctions of the problem \eqref{e1}--\eqref{e3}, we
introduce the functions
\begin{equation} \label{e33}
\begin{aligned}
 u_{mn}(x,y) &= \tilde{U}_{mn}(x,y)-\tilde{U}_{mn}(y,x)\\
 &=\exp(\frac{(2m+1)x}{2}+\frac{(2n+1)y}{2})\pi i\\
&\quad -\exp(\frac{(2m+1)y}{2}+\frac{(2n+1)x}{2})\pi i,
\end{aligned}
\end{equation}
for $m, n=0,\pm1,\pm2,\dots$.
The set $\{u_{mn}(x,y)\}$ is complete and orthogonal in $L_2(T')$ where $T'=\{(x,y):
0\leq y\leq 2, 0\leq x\leq y\}$.

\begin{theorem} \label{thm1}
Problem \eqref{e1}--\eqref{e3} when $p\equiv 1$ is selfadjoint. It has eigenvalues \eqref{e32} and
eigenfunctions \eqref{e33}. The eigenfunctions are complete and orthogonal in
$L_2(T)$.
\end{theorem}

\section{The non-selfadjoint problem}

Now consider the problem \eqref{e1}--\eqref{e3} and  make the change of variables
$x = s - t$ and $y = s + t$, to obtain
\begin{gather}
U_{xy}=\gamma P(x, y) U(x, y),\quad (x,y)\in T', \label{e34}\\
T' = \{(x, y): 0\leq y\leq 2, 0\leq x \leq y\} \label{e35}\\
\alpha U(0,y) + \beta U_x(0,y)=\alpha U(y,2) + \beta U_y(y,2),\quad 0\leq
y\leq2 \label{e36} \\
U(x, x) = 0,\quad 0\leq x\leq 2, \label{e37}
\end{gather}
where $U(x, y) = u(\frac{x+y}{2}, \frac{y-x}{2})$,
$P(x, y) = p(\frac{x+y}{2},\frac{y-x}{2})$, and $\gamma = -\lambda/4$.
Integrating \eqref{e34} in $T'$ from $0$ to $\xi$, we have
\begin{equation} \label{e38}
U_{y}(\xi, y) - U_{y}(0, y) = \gamma \int_{0}^{\xi}PU\,dx.
\end{equation}
The above equation when $\xi=y$ and $y=2$ will be
\begin{equation} \label{e39}
U_{y}(y, 2) - U_{y}(0, 2) = \gamma \int_{0}^{y}P(x,2)U(x,2)\,dx.
\end{equation}
Integrating equation \eqref{e34} in $T'$ from $y$ to $\eta$, we have
\begin{equation*} %\label{e40}
U_{x}(x, \eta) - U_{x}(x, y) = \gamma \int_{y}^{\eta}PU\,dy.
\end{equation*}
When $x=0$ and $\eta = 2$ this equation becomes
\begin{equation} \label{e41}
U_{x}(0, 2) - U_{x}(0, y) = \gamma \int_{y}^{2}P(0, y)U(0, y)\,dy.
\end{equation}
In \eqref{e38} let $\xi$ lie on the line $y=x$ and integrate the
equation from $\xi$ to $\eta$ in $T'$,
\begin{equation} \label{e42}
U(\xi, \eta) - U(\xi, \xi) -U(0,
\eta)+U(0,\xi)=\gamma\int_{\xi}^{\eta}\int_{0}^{\xi}PUdxdy.
\end{equation}
Since $ U(\xi,\xi) = 0$ by the boundary condition \eqref{e37}, we have
\begin{equation} \label{e43}
U(\xi, \eta) = U(0,\eta) - U(0,\xi)
+ \gamma\int_{\xi}^{\eta}\int_{0}^{\xi}PU\,dx\,dy\,.
\end{equation}
 From \eqref{e39} when $y$ is replaced with $\eta$ we have
\begin{equation} \label{e44}
U_{y}(\eta, 2) = U_{y}(0, 2) + \gamma \int_{0}^{\eta}P(x,2)U(x,2)\,dx.
\end{equation}
 From \eqref{e41}, when $y$ is replaced with $\eta$ we have
\begin{equation} \label{e45}
U_{x}(0, \eta) =U_{x}(0, 2) - \gamma \int_{\eta}^{2}P(0, y)U(0, y)\,dy.
\end{equation}
Finally from equation \eqref{e43} we have,
\begin{equation} \label{e46}
U(\eta, 2) = U(0,2) - U(0,\eta) + \gamma\int_{\eta}^{2}\int_{0}^{\eta}
PU\,dx\,dy\,.
\end{equation}
Now, we substitute the right hand side of the equations \eqref{e44}, \eqref{e45}, \eqref{e46}
into the boundary condition \eqref{e36} with $y$ replaced with $\eta$,
\begin{equation*} %\label{e47}
\begin{aligned}
&\alpha U(0,\eta) + \beta(U_{x}(0, 2) - \gamma
\int_{\eta}^{2}P(0,y)U(0,y)dy)\\
&= \alpha(U(0,2) - U(0,\eta) + \gamma\int_{\eta}^{2}\int_{0}^{\eta}PU\,dx\,dy)
 + \beta(U_{y}(0, 2) \\
 &\quad + \gamma \int_{0}^{\eta}P(x,2)U(x,2)dx).
\end{aligned}
\end{equation*}
Placing $\alpha U(0, \eta)$ in right-hand side and combining,
\begin{equation} \label{e48}
\begin{aligned}
&2\alpha U(0,\eta) + \beta(U_{x}(0, 2) + \gamma
\int_{\eta}^{2}P(0,y)U(0,y)dy)\\
&= \alpha(U(0,2)+ \gamma\int_{\eta}^{2}\int_{0}^{\eta}PU\,dx\,dy) +
\beta(U_{y}(0, 2) \\
&\quad + \gamma \int_{0}^{\eta}P(x,2)U(x,2)dx).
\end{aligned}
\end{equation}
Rewrite this  equation with $\eta$ replaced with $\xi$,
\begin{equation} \label{e49}
\begin{aligned}
&2\alpha U(0,\xi) + \beta(U_{x}(0, 2) + \gamma
\int_{\xi}^{2}P(0,y)U(0,y)dy)\\
&=\alpha(U(0,2) + \gamma\int_{\xi}^{2}\int_{0}^{\xi}PU\,dx\,dy) +
\beta(U_{y}(0, 2) \\
&\quad + \gamma \int_{0}^{\xi}P(x,2)U(x,2)dx).
\end{aligned}
\end{equation}
Subtract equation \eqref{e49} from \eqref{e48},
\begin{align}
&2\alpha (U(0,\eta)-U(0,\xi)) + \beta \gamma (
\int_{\eta}^{2}P(0,y)U(0,y)dy - \int_{\xi}^{2}P(0,y)U(0,y)dy) \nonumber\\
&= \alpha \gamma (\int_{\eta}^{2}\int_{0}^{\eta}PUdxdy -
\int_{\xi}^{2}\int_{0}^{\xi}PU\,dx\,dy)  \label{e50} \\
&\quad + \beta \gamma(\int_{0}^{\eta}P(x,2)U(x,2)dx
-\int_{0}^{\xi}P(x,2)U(x,2)dx). \nonumber
\end{align}
Solve for $U(0,\eta)-U(0,\xi)$ in \eqref{e50}, assuming $\alpha\neq 0$, and
substitute in \eqref{e43},
\begin{equation} \label{e51}
\begin{aligned}
 U(\xi, \eta)
 =& -\frac{\beta}{2\alpha}\gamma (
\int_{\xi}^{2}P(0,y)U(0,y)dy - \int_{\eta}^{2}P(0,y)U(0,y)dy)\\
&+\frac{\gamma}{2} (\int_{\eta}^{2}\int_{0}^{\eta}PU\,dx\,dy -
\int_{\xi}^{2}\int_{0}^{\xi}PU\,dx\,dy) \\
&+ \frac{\beta}{2\alpha}\gamma(\int_{0}^{\eta}P(x,2)U(x,2)dx
-\int_{0}^{\xi}P(x,2)U(x,2)dx)\\
&+ \gamma \int_{\xi}^{\eta}\int_{0}^{\xi} PU\,dx\,dy\,.
\end{aligned}
\end{equation}
Rewrite \eqref{e51} in a compact form using the Green's function,
$G(\xi,\eta;x,y)$ described in Figure 1,
\begin{equation} \label{e52}
\begin{aligned}
U(\xi,\eta) =& \gamma (\int\int_{T'}G(\xi,\eta; x,y)PU(x,y)\,dx\,dy \\
&-\int_{0}^{2}g(\xi,\eta,x) (P(0,x)U(0,x)-P(x,2)U(x,2))dx),
\end{aligned}
\end{equation}
where
\[
g(\xi,\eta,x)=\begin{cases}
 0  &\mbox{if } 0\leq x\leq \xi\\
\beta/(2\alpha) &\mbox{if } \xi\leq x\leq \eta\\
0 &\mbox{if } \eta \leq x \leq 2
\end{cases}
\]
Now, let $G$ be the operator,
\begin{equation} \label{e53}
G[U]=\int\int_{T'}G(\xi,\eta; x,y)PU(x,y)\,dx\,dy
\end{equation}
defined on the Hilbert space of weighted square integrable functions
$H = L_2^P(T')$. Assume the function $U$ along the characteristics
$x = 0$ and $y = 0$ in $T'$ is given, and denote,
\begin{equation}
f(\xi, \eta) = \int_{0}^{2}g(\xi,\eta,x) (P(0,x)U(0,x)-P(x,2)U(x,2))dx.
\end{equation}
Also, let $\Gamma=1/\gamma$ when $\gamma\neq 0$. Then, \eqref{e48} can be
written as
\begin{equation} \label{e55}
G[U]=\Gamma U + f.
\end{equation}
The operator $G$ in \eqref{e53}, is the same as the one in
\cite{kr:symmetric}, where it is shown to be selfadjoint in $H$. Denote the normalized
eigenfunctions of $G$ by $E_k$, and the inner product in $H$ by
$\langle .,.\rangle$. Using the standard
Fredholm alternative \cite{st:Green} we have the following theorem.

\begin{figure}[ht]
\begin{center}
\setlength{\unitlength}{8mm}
\begin{picture}(9,7)(1,1)
\put(2,1){\line(7,0){7}}
\put(2,1){\line(0,7){7}}
\put(1.8,1){\line(6,6){5.4}}
\put (9,.6){x}
\put (7,.6){2}
\put (4,.6){w}
\put (6,.6){z}
\put (1.7,8){y}
\put (1.7,6){2}
\put (1.6,3){w}
\put (1.7,5){z}
\put(2,6){\line(5,0){5}}
\put(2,5){\line(4,0){4}}
\put(2,3){\line(2,0){2}}
\put(6,5){\line(0,1){1}}
\put(4,3){\line(0,3){3}}
\put (3,2.5){0}
\put (5,4.5){0}
\put (6.5,5.6){0}
\put (3,5.6){0}
\put (2.8,4.6){1/2}
\put (4.8,5.6){1/2}
\end{picture}
\end{center}
\caption{The Green's function $G(w,z;x,y)$}
\end{figure}



\begin{theorem} \label{thm2}
Let $\Gamma_k$ be the eigenvalues of $GU=\Gamma U$, where by the selfadjointness of $G$,
are real and satisfy $|\Gamma_k| \rightarrow 0$, as $k\rightarrow \infty$ then,
we have\begin{enumerate}

\item $\Gamma \neq \Gamma_k$ for any integer k. A unique solution of \eqref{e55} exists
and is given in the form $ U = -\frac{f}{\Gamma} +
\sum_k\frac{\Gamma_k}{\Gamma(\Gamma_k-\Gamma)}\langle f, E_k\rangle E_k$.

\item $\Gamma = \Gamma_m$, one of the eigenvalues of $G$, and $\Gamma_m$ is
not degenerate. If $\langle f, E_m\rangle \neq 0$ the equation \eqref{e55} has no
solution. If $\langle f, E_m\rangle = 0$ then, \eqref{e55} has
infinitely many solutions
$$ U = -\frac{f}{\Gamma_m} + \sum_{k\neq m}\frac{\Gamma_k}{\Gamma_m
(\Gamma_k-\Gamma_m)}\langle f, E_k\rangle E_k + cE_m,
$$
with $c$ an arbitrary constant.

\item  $\Gamma_m$ is degenerate, $\Gamma_{m_1}=\Gamma_{m_2}=\dots=\Gamma_{m_j}$,
 for successive indices $m_1$, $m_2$,\dots $m_j$, with some $m_i=m$,
$i=1,2,\dots j$, where $j$ is the multiplicity of $\Gamma_m$.
Then, unless $\langle f, E_i\rangle =0$, $i=1, 2,\dots,j$, the equation \eqref{e55}
has no solution. If however, these $j$ solvability conditions are satisfied,
the solution can be represented by,
$$ U = -\frac{f}{\Gamma_m} + \sum_{k\neq m_i}\frac{\Gamma_k}{\Gamma_m
(\Gamma_k-\Gamma_m)}\langle f, E_k\rangle E_k +\sum_{i=1}^{j}c_iE_{m_i},
$$
where $c_i$'s are arbitrary constants.
\end{enumerate}
\end{theorem}

\paragraph{Remark}
When $p\equiv 1$ in problem \eqref{e1}--\eqref{e3}, the problem is selfadjoint if
$\alpha^2+\beta^2 \neq 0$. When $p$ is not necessarily 1, the case
$\alpha\neq 0, \beta=0$, reduces the problem to the selfadjoint problem
of \cite{kr:symmetric}. When $\alpha \neq 0,
\beta\neq 0$, Theorem 2 holds. When $\alpha = 0$, no conclusion can be drawn.

\begin{thebibliography}{00} \frenchspacing
\bibitem{Ga:Partial}
P. Garabedian. {\em Partial differential equations. Second edition}, Chelsea
Publishing Co., New York, 1986. vxii+672 pp. ISBN: 0-8284-00325

\bibitem{Ha:symmetric} L. Haws. {\em Symmetric Green's functions for certain hyperbolic
problems}, Comput. Math. Appl. 21 (1991), no. 5, 65--78.

\bibitem{Ir:A Method} N. Iraniparast. {\em A method of solving a class of CIV boundary
value problems}, Canad. Math. Bull. 35 (1992), no. 3, 371--375.

\bibitem{Ir:A Boundary} N. Iraniparast. {\em A boundary value problem for the wave
equation}. Int. J. Math. Math. Sci. 22 (1999), no. 4, 835--845.

\bibitem{Ir:A CIV} N. Iraniparast. {\em A CIV boundary value problem for the wave
equation}, Appl. Anal. 76 (2000), no. 3-4, 261--271.

\bibitem{Ka:spectrum} T. Kalmenov. {\em On the Spectrum of a Selfadjoint Problem for the
Wave Equation}, Akad. Nauk. Kazakh. SSR, Vestnik 1 (1983), 63-66.

\bibitem{kr:symmetric} K. Kreith. {\em Symmetric Green Functions for a Class of CIV
Boundary Value Problems}, Canadian Math. Bull. 31 (1988), 272-279.

\bibitem{kr:Mixed} K. Kreith. {\em Mixed selfadjoint boundary conditions for the wave
equation}, Differential equations and its applications (Budapest, 1991), 219--226, Colloq.
Math. Soc. Janos Bolyai, 62, North-Holland, Amsterdam, 1991.

\bibitem{kr:Establishing} K. Kreith. {\em Establishing hyperbolic Green's functions via
Leibniz's rule}, SIAM Rev. 33 (1991), no. 1, 101--105.

\bibitem{st:Green} I. Stakgold. {\em Green's Functions and Boundary Value Problems},
John Wiley \& Sons, New York, 1979.
\end{thebibliography}

\noindent\textsc{Nezam Iraniparast}\\
Department of Mathematics\\
Western Kentucky University\\
E-mail: nezam.iraniparast@wku.edu

\end{document}
