Semianalytical calculations of fusion neutron spectra in case of the presence of substantial suprathermal tails in fuel nuclei velocity distributions have been implemented using two different approaches. Numerical techniques are described. Results for a variety of isotropic and anisotropic distributions of suprathermal fuel nuclei are presented. Verification of calculated fusion product spectra has been performed using several methods.