  A Lagrangian method for the numerical simulation of the Kraichnan passive
scalar model is introduced. The method is based on Monte--Carlo simulations of
tracer trajectories, supplemented by a point-splitting procedure for coinciding
points. Clean scaling behavior for scalar structure functions is observed. The
scheme is exploited to investigate the dependence of scalar anomalies on the
scaling exponent $\xi$ of the advecting velocity field. The three-dimensional
fourth-order structure function is specifically considered.
