  This paper describes the design and work of a set of computer routines
capable for numerical computation of generalized parton distributions (GPDs)
evolution at the leading order. The main intention of this work is to present a
fast-working computer code making possible fitting of GPDs parameters to the
data on hard electron-nucleon scattering.
