The most probable orientational distribution functions of rod-like polar molecules contained in a noncentrosymmetric uniaxial system are established using the first-rank and third-rank Legendre polynomials, (P1(cos theta)) and (P3(cos theta)) order parameters, and the maximum entropy method. Emphasis is put on the different domains of existence in the ((P1), (P3)) plane for the various shapes of the distributions: it is thus shown that, for any positive (P1(cos theta)) value and for decreasing (P3(cos theta)) values, the distribution function may exhibit either a distorted oblate form with an intense maximum at 0 degrees, or a three-leaved rose curve with maxima at 60 degrees, 180 degrees, and 300 degrees, and finally another markedly oblate shape with a strong maximum at 180 degrees. As an illustrative example, we have considered the azobenzene molecular orientations in an electrically poled p(DR1M) homopolymer thin film after a thermal process and several relaxation periods. We have made use not only of the (P1) and (P3) parameters determined from polarized second-harmonic generation (SHG) measurements, but also of the (P2) values extracted from UV-visible spectra and of the (P4) values adjusted according to the information entropy theory. In such a thin film with very large nonlinear properties (d33 coefficients were varying from 437.0 to 117.0 pm/V at 1064 nm) it is evidenced that a strong polar order is maintained even after a long relaxation period of 42 days. So, the distribution functions demonstrate that the poling treatment was quite efficient and they emphasize the importance in the determination of both couples of odd and even order parameters in such uniaxially oriented optical elements.