Code for prompt numerical computation of the leading order GPD evolution

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.