We derive a general solution of the heat equation through two modied separation methods.

The obtained solution is expressed as linearly combined kernel solutions in terms of Hermite polynomials, which appears to provide an explanation of non-Gaussian behavior observed in various cases. We also consider a typical boundary condition and construct corresponding solutions. It is revealed the boundary-value problem consisting of the heat-transfer partial differential equation and the boundary condition carries infinitely many solutions.