Solve the One Dimensional Eigenvalue Problem: $u_{tt} - u_{xx} = m^2 u$ 
Solve the One Dimensional Eigenvalue Problem: 
  $$u_{tt} - u_{xx} = m^2 u$$

I have learned how to solve PDEs of the form $u_{tt} - u_{xx} = f(x,t)$ using d'Alembert's solution. However, I have not yet studied PDEs of this form. How does one solve this PDE? 
 A: Suppose we try separation of variables, setting
$u(x, t) = f(t)g(x); \tag 1$
then
$u_{xx} = f(t) g''(x), \tag 2$
$u_{tt} = \ddot f(t) g(x), \tag 3$
whence the equation
$u_{tt} - u_{xx} = m^2 u \tag 4$
becomes
$\ddot f g - f g'' = m^2 fg; \tag 5$
if we assume for the moment that $f(t) \ne 0 \ne g(x)$, then we may divide (5) through by $fg$, yielding
$\dfrac{\ddot f}{f} - \dfrac{g''}{g} = m^2, \tag 6$
or
$\dfrac{\ddot f}{f} - m^2 = \dfrac{g''}{g}. \tag 7$
We can now argue in the typical manner that, since the left-hand side depends only on $t$ but the right-hand side only on $x$, they must each equal the same constant value which, for reasons which will hopefully soon become clear, we denote by $-\lambda^2$:
$\dfrac{\ddot f}{f} - m^2 = -\lambda^2 = \dfrac{g''}{g}, \tag 8$
which leads to
$\ddot f + (\lambda^2 - m^2) f = 0, \tag 9$
and
$g'' + \lambda^2 g = 0.\tag{10}$
(10) has solutions of the general form
$g_\lambda(x) = a_+(\lambda)e^{i\lambda x} + a_-(\lambda) e^{-i\lambda x}, \tag{11}$
where $\lambda$ and the $a_\pm(\lambda) \in \Bbb C$ are determined by boundary/initial conditions placed upon $g(x)$; we can't get much more specific at this point without such information.  For each admissible $\lambda$, we then have a solution for $f(t)$:
$f_\lambda(t)= b_+(\lambda) e^{i\sqrt{\lambda^2 - m^2}t} + b_-(\lambda)e^{-i\sqrt{\lambda^2 - m^2}t}; \tag{12}$
therefore,
$u_\lambda(x, t) = f_\lambda(t)g_\lambda(x)$
$= ( b_+(\lambda) e^{i\sqrt{\lambda^2 - m^2}t} + b_-(\lambda) e^{-i\sqrt{\lambda^2 - m^2}t})(a_+(\lambda)e^{i\lambda x} + a_-(\lambda) e^{-i\lambda x}); \tag{13}$
it is easy to check that
$u_{\lambda tt} - u_{\lambda xx} = m^2 u_\lambda. \tag{14}$
Solutions of the form (13) for different $\lambda$ may be composed by summation or integration over $\lambda$ to yield more complicated solutions. 
Finally, it is of course worth observing that the temporal frequencies $\sqrt{\lambda^2 - m^2}$ may become imaginary if $\lambda^2 < m^2$; the we obtain temporally exponentially growing/decreasing solutions $u(x, t)$.
In the above, I have decided to associate the $m^2$ terms with $f(t)$, but we could just as well have taken
$g'' + (\lambda^2 + m^2)g = 0 \tag{15}$
and obtained a different, though similar, set of solutions to (4).  In fact, if we set
$m^2 = m_t^2 - m_x^2, \tag{16}$
we can write
$\dfrac{\ddot f}{f} - \dfrac{g''}{g} = m_t^2 - m_x^2, \tag{17}$
whence
$\dfrac{\ddot f}{f} - m_t^2 = -\lambda^2 = \dfrac{g''}{g} - m_x^2, \tag{18}$
which leads to a more symmetrical analysis and set of solutions.  
This could go on forever, but I can't, so I'll leave it at that.  I would like to say in closing, however, that all this stuff smells like the Special Theory of Relativity and Particle Physics to my somewhat knowlewdge-able knose.
A: Let $p=mx$ ,
Then $u_{tt}=m^2u_{pp}+m^2u$
Let $q=mt$ ,
Then $u_{qq}=u_{pp}+u$
Similar to Solve initial value problem for $u_{tt} - u_{xx} - u = 0$ using characteristics,
Consider $u(p,a)=f(p)$ and $u_q(p,a)=g(p)$ ,
Let $u(p,q)=\sum\limits_{n=0}^\infty\dfrac{(q-a)^n}{n!}\dfrac{\partial^nu(p,a)}{\partial q^n}$ ,
Then $u(p,q)=\sum\limits_{n=0}^\infty\dfrac{(q-a)^{2n}}{(2n)!}\dfrac{\partial^{2n}u(p,a)}{\partial q^{2n}}+\sum\limits_{n=0}^\infty\dfrac{(q-a)^{2n+1}}{(2n+1)!}\dfrac{\partial^{2n+1}u(p,a)}{\partial q^{2n+1}}$
$u_{qqqq}=u_{ppqq}+u_{qq}=u_{pppp}+u_{pp}+u_{pp}+u=u_{pppp}+2u_{pp}+u$
$u_{qqqqqq}=u_{ppppqq}+2u_{ppqq}+u_{qq}=u_{pppppp}+u_{pppp}+2u_{pppp}+2u_{pp}+u_{pp}+u=u_{pppppp}+3u_{pppp}+3u_{pp}+u$
Similarly, $\dfrac{\partial^{2n}u}{\partial q^{2n}}=\sum\limits_{k=0}^nC_k^n\dfrac{\partial^{2k}u}{\partial p^{2k}}$
$u_{qqq}=u_{ppq}+u_q$
$u_{qqqqq}=u_{ppqqq}+u_{qqq}=u_{ppppq}+u_{ppq}+u_{ppq}+u_q=u_{ppppq}+2u_{ppq}+u_q$
$u_{qqqqqqq}=u_{ppppqqq}+2u_{ppqqq}+u_{qqq}=u_{ppppppq}+u_{ppppq}+2u_{ppppq}+2u_{ppq}+u_{ppq}+u_q=u_{ppppppq}+3u_{ppppq}+3u_{ppq}+u_q$
Similarly, $\dfrac{\partial^{2n+1}u}{\partial q^{2n+1}}=\sum\limits_{k=0}^nC_k^n\dfrac{\partial^{2k+1}u}{\partial p^{2k}\partial q}$
$\therefore u(p,q)=\sum\limits_{n=0}^\infty\sum\limits_{k=0}^n\dfrac{C_k^nf^{(2k)}(p)(q-a)^{2n}}{(2n)!}+\sum\limits_{n=0}^\infty\sum\limits_{k=0}^n\dfrac{C_k^ng^{(2k)}(p)(q-a)^{2n+1}}{(2n+1)!}$
