For the Gaussian ansatz, partial differentiation gives
Cancel the positive exponential and compare the constant and quadratic powers of in the diffusion equation. The two coupled ordinary differential equations are
The initial profile requires and . The second equation is a separable differential equation:
Substitute this into the first equation and integrate its logarithmic derivative:
Thus
and the solution is
This is Gaussian diffusion of a smooth initial profile: its width grows while its peak decreases. The Gaussian integral gives , independent of time, which checks conservation of the total amount under the diffusion equation.