Another thing you could have done, to avoid so many calculations was to note that solves the heat equation on and noting that where the means simply that the convolution is with respect to the first variable. Now properties of this operation give you the result based on that for . The fact that it satisfies the initial (or boundary, depending on how you view it) condition on is more delicate since you can't just plug that value.