Two-sided Gaussian bounds for fundamental solutions of non-divergence form parabolic operators with H\"older continuous coefficients

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.