  Dynamical interactions that occur between objects in dense stellar systems
are particularly important for the question of formation of X-ray binaries. We
present results of numerical simulations of 70 globular clusters with different
dynamical properties and a total stellar mass of 2*10^7 Msun. We find that in
order to retain enough neutron stars to match observations we must assume that
NSs can be formed via electron-capture supernovae. Our simulations explain the
observed dependence of the number of LMXBs on ``collision number'' as well as
the large scatter observed between different globular clusters. For millisecond
pulsars, we obtain good agreement between our models and the numbers and
characteristics of observed pulsars in the clusters Terzan 5 and 47 Tuc
