[en] A numerical method has been developed for obtaining of the Legendre polynomials expansion coefficients of the anisotropy density function of the elastic neutron scattering. The method uses natural cubic splines as a basis for representing of the anisotropy density function and the Legendre polynomials. The procedure is programmed and incorporated into the SUPERTOG package for preparation of multigroup neutron data