  The main features of the electron structure function method for calculations
of the higher order QED radiative effects to polarized deep-inelastic
ep-scattering are presented. A new FORTRAN code ESFRAD based on this method was
developed. A detailed quantitative comparison between the results of ESFRAD and
other methods implemented in the codes POLRAD and RADGEN for calculation of the
higher order radiative corrections is performed.
