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