  We investigate the evolution of globular clusters using N-body calculations
and anisotropic Fokker-Planck (FP) calculations. The models include a mass
spectrum, mass loss due to stellar evolution, and the tidal field of the parent
galaxy. Recent N-body calculations have revealed a serious discrepancy between
the results of N-body calculations and isotropic FP calculations. The main
reason for the discrepancy is an oversimplified treatment of the tidal field
employed in the isotropic FP models. In this paper we perform a series of
calculations with anisotropic FP models with a better treatment of the tidal
boundary and compare these with N-body calculations. The new tidal boundary
condition in our FP model includes one free parameter. We find that a single
value of this parameter gives satisfactory agreement between the N-body and FP
models over a wide range of initial conditions.
  Using the improved FP model, we carry out an extensive survey of the
evolution of globular clusters over a wide range of initial conditions varying
the slope of the mass function, the central concentration, and the relaxation
time. The evolution of clusters is followed up to the moment of core collapse
or the disruption of the clusters in the tidal field of the parent galaxy. In
general, our model clusters, calculated with the anisotropic FP model with the
improved treatment for the tidal boundary, live longer than isotropic models.
The difference in the lifetime between the isotropic and anisotropic models is
particularly large when the effect of mass loss via stellar evolution is rather
significant. On the other hand the difference is small for relaxation-
dominated clusters which initially have steep mass functions and high central
concentrations.
