Abstract
A numerical scheme based on the analytic solution of the discrete-ordinate transport equations in each mesh cell is developed. The driving source for the transport equation is represented by a two-moment preserving, strictly positive, exponential distribution obtained using information theory. This numerical scheme is shown to be accurate for both large and small mesh intervals and is strictly positive for all mesh sizes. The scheme has been coded into the existing ONELD(tm) code and tested. We present numerical results to demonstrate the accuracy and positivity of this scheme