 Solution

The two homogeneous Dirichlet boundary conditions on $x$ force negative separation constants. The initial profile is an even function, so only the even sine-family eigenfunctions are needed:They obey $X_m(\pm 1)=0$ and ${\int }_{-1}^1{X}_m{X}_{n}dx={\delta }_{mn}$. Odd eigenfunctions have zero projection onto the even initial data. Hence the initial coefficients areThe desired eigenfunction expansion is thereforeAt $y=0$ this is the Fourier series of $1-x^2$ on the interval; its coefficients are absolutely summable. Each term satisfies the partial differential equation and side conditions. For $y>0$ the additional Gaussian decay in $m$ justifies termwise differentiation. Since $k_m^2+1>0$ and $e^{\cos y-1}\le 1$, the whole series decays uniformly as $y\to \infty$.