Gaussian bounds for degenerate parabolic equations