A numerical study of the small-x behavior of deep inelastic structure functions in QCD

A numerical study of the Gribov-Levin-Ryskin equation for deep inelastic structure functions in the small-x regions is performed. In the first part we study a simplified version which illustrates general features of the solution and also allows to investigate certain analytic properties. In the second part we formulate a more realistic version of the GLR equation and compute its solution, for different input distributions and for different values of the strength of the nonlinearity in the GLR equation. One of the conclusions is that, provided the input distribution is not too flat and the nonlinear coupling is not too small, deviations from the standard evolution may be visible at HERA