  In this paper we have solved the nonlinear Gribov-Levin-Ryskin-Mueller-Qiu
(GLR-MQ) evolution equation for gluon distribution function G(x,Q^2) and
studied the effects of the nonlinear GLR-MQ corrections to the Leading Order
(LO) Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equations.
Here we incorporate a Regge like behaviour of gluon distribution function to
obtain the solution of GLR-MQ evolution equation. We have also investigated the
Q^2-dependence of gluon distribution function from the solution of GLR-MQ
evolution equation. Moreover it is interesting to observe from our results that
nonlinearities increase with decreasing correlation radius (R) between two
interacting gluons. Results also confirm that the steep behavior of gluon
distribution function is observed at R=5 GeV^{-1}, whereas it is lowered at R=2
GeV^{-1} with decreasing x as Q^2 increases. In this work we have also checked
the sensitivity of \lambda_G in our calculations. Our computed results are
compared with those obtained by the global DGLAP fits to the parton
distribution functions viz. GRV, MRST, MSTW and with the EHKQS model.
