Monte Carlo simulations have been performed for a discotic liquid crystal composed of Gay-Berne particles. On the basis of these simulations for the nematic phase, a subset of the spherical harmonic expansion coefficients of the direct pair correlation function (DPCF) were determined from the pair distribution function (PDF) by solving the Ornstein-Zernike (OZ) equation, This was achieved by generalizing the Wiener-Hopf factorization scheme for the numerical solution of the OZ equation. Only the expansion coefficients g(l1,l2,l)(r) (l(alpha)less than or equal to 4) of the PDF in the laboratory frame were used when solving the OZ equation; this means that the DPCF so obtained is equivalent to that for a nematic in which the director is randomly distributed. From the DPCF, the scaled Oseen-Zocher-Frank, elastic constants K-11*, K-22*, and K-33*, as well as the surface constant K-13*, have been calculated from the subset of expansion coefficients. Generally, we find that K-33*