Prove the Lyapunov stability in $L^2([0, \ell])$ of the diffusion problem governed by $ \frac{\partial u}{\partial t} − a^2\frac{\partial^2 u}{\partial x^2}= f(x, t), 0 \leq x \leq \ell, t \geq 0 $ subject to Neumann boundary conditions $u_x(0, t) = 0, u_x(\ell, t) = 0$, with respect to variations in the initial condition, $u(x, 0) = \varphi(x)$.
I'm unsure about how to start this problem. I know i need to find a $C>0$ such that it satisfies $$\int_0^\ell [u_1(x, t) − u_2(x, t)]^2{\rm d}x ≤ C\int_0^\ell[\varphi_1(x) − \varphi_2(x)]^2{\rm d}x$$
but I am unsure how to do this.