  In this paper, we consider the regularity of solutions to elliptic equation
with Grushin's operator. By using the Feynman-Kac formula, we first get the
expression of heat kernel, and then by using the properties of heat kernel, the
optimal regularity of solutions will be obtained. The novelty of this paper is
that the Grushin's operator is a degenerate operator.
