  We numerically study the effects of high gluon density at small x on the
evolution of gluon distribution function in both hadrons and nuclei. Using a
newly derived, Wilson renormalization group-based evolution equation which
includes n to 1 gluon ladder fusion, we find significant reduction in nuclear
gluon distribution function for large nuclei at zero impact parameter at
energies relevant for RHIC and LHC experiments.
