An improved implementation of the discrete ordinates method for computing neutral particle transport in ducts is presented. The considered one-dimensional model makes use of two basis functions to represent the transverse and azimuthal dependencies of the particle angular flux in the duct. It is shown that if the problem is decomposed into uncollided and collided problems prior to using the discrete ordinates approximation, the number of ordinates necessary to achieve a desired degree of accuracy in the solution can be greatly reduced, especially for long ducts with significant wall absorption. Further savings in computer time can be attained by employing a composite quadrature based on a (nonstandard) half-range quadrature that can be generated in an effective and efficient way with one of the classical methods in the constructive theory of orthogonal polynomials.