Abstract: The difficulties inherent in the conventional numerical implementation of the discrete ordinate method (following Chandrasekhar's prescription) for solving the radiative transfer equation are discussed. A matrix formulation ...