  Field theoretic renormalization group is applied to the Kraichnan model of a
passive scalar advected by the Gaussian velocity field with the covariance
$<{\bf v}(t,{\bf x}){\bf v}(t',{\bf x})> - <{\bf v}(t,{\bf x}){\bf v}(t',{\bf
x'})> \propto\delta(t-t')|{\bf x}-{\bf x'} |^{\epsilon}$. Inertial-range
anomalous exponents, related to the scaling dimensions of tensor composite
operators built of the scalar gradients, are calculated to the order
$\epsilon^{3}$ of the $\epsilon$ expansion. The nature and the convergence of
the $\epsilon$ expansion in the models of turbulence is are briefly discussed.
