M87, a cD galaxy in the Virgo cluster, has 3-10 times larger enclosed mass than the Milky Way at any given galactocentric radius. Thus the globular cluster (GC) system in M87 is expected to have undergone a more significant dynamical evolution than that of the Milky Way if it had started from the same initial mass function (MF) and radial distribution (RD) as the Milky Way. The evolution of MF and RD of the M87 GC system has been calculated using an advanced, realistic Fokker-Planck (FP) model that considers dynamical friction, disk/bulge shocks, and eccentric cluster orbits. We perform hundreds of FP calculations with different initial cluster conditions, and then search a wide parameter space for the best-fit initial GC MF and RD that evolves into the observed present-day GC MF and RD. We also find best-fit initial MFs and RDs for blue and red GC groups, separately.