The behavior of gluon distributions in the small- x region are semi-numerically predicted by using the nonlinear Gribov–Levin–Ryskin and Mueller–Qiu (GLR–MQ) evolution equation with the nonlinear shadowing term incorporated. It is interesting to observe that gluon recombination at the twist-4 level lowers the rapid growth of gluon densities with decreasing x compared to the linear DGLAP evolution. The nonlinear gluon density functions are calculated and compared with different parameterizations as well as models. Results show that our predictions for the gluon density function are comparable with different parameterizations. We make a deliberate attempt to explore the effect of nonlinear or shadowing corrections and we observe that nonlinea...