  The saturation model with DGLAP evolution is shown to give good description
of the production of the charm and beauty quarks in deep inelastic scattering.
The modifications of saturation properties caused by the presence of heavy
quarks are also discussed.
