Solution (source code)

= Solution

Let there be $J$ interior points. The <Dirichlet discrete Laplacian> has positive values
$$
\eta_j=4\sin^2\left(\frac{j\pi}{2(J+1)}\right)
$$
for the eigenvalues of its negative, and the amplification eigenvalues are $(1+\mu\eta_j)^{-1}$. Therefore the <Backward Euler diffusion stability on a finite Dirichlet interval> is
$$
\boxed{\mu\geq0
\quad\text{or}\quad
\mu\leq-\frac1{2\sin^2(\pi/(2(J+1)))}}.
$$
If, as usual, a Courant number is restricted to nonnegative values, this again reduces to all $\mu\geq0$. The additional negative branch is a finite-grid artefact and disappears to $-\infty$ as the mesh is refined.