The solution to this equation with an impulse initial condition is Exp[-x 2 /t] , and with a block from -a to a it is (Erf[(a - x)/ √ t ] + Erf[(a + x)/ √ t ])/a .