  We establish two-sided Gaussian bounds for fundamental solutions of general
non-divergence form parabolic operators with H\"older continuous coefficients.
The result we obtain is essentially based on parametrix method.
