Two methods for determination of molecular orientation distribution from angular dependence of EPR spectra are compared. One of these methods is based on direct expansion of the orientation distribution function in a series of generalized spherical harmonics, in another method the mean field potential is used. The limitations of these approaches as well as their capabilities and advantages are considered