 Solution

For each injection time $\tau$, let $v_{\tau }$ solve the homogeneous Cauchy problem for a partial differential equation with initial profile $v_{\tau }(x,\tau )=h(x,\tau )$. Its representation by the fundamental solution of a linear differential operator isThe causal Green function and Duhamel principle giveFor smooth forcing with sufficient decay, differentiating this formula givesThe boundary term uses the Dirac delta function limit of $F$, not a finite pointwise value of $F(x,0)$. The same formula extends to weaker forcing in an appropriate integrable or distribution sense. Physically, each short time interval adds $h(\cdot ,\tau )d\tau$ to the temperature profile; that addition subsequently undergoes homogeneous diffusion equation. The total solution is the superposition of all such diffusing injections.