The solution to the neutron Boltzmann equation is separated into a straightahead component dominating at high energies and an isotropic component dominating at low energies. The high-energy solution is calculated using HZETRN-05, and the low-energy isotropic component is modeled by two non-coupled integro-differential equations describing both forward and backward neutron propagation. Three different solution methods are then used to solve the equations. The collocation method employs linear I3-splines to transform each equation into a system of ODES; the resulting system is then solved exactly and evaluated using numerical integration techniques. Wilson\u27s method uses a perturbational approach in which a fundamental solution is obtained by solving a simple ODE, a new source term is generated by the fundamental solution, and the collocation method is then used to solve the remaining equation. The fixed-point series method extends Wilson\u27s method by continuing the perturbational procedure until desired convergence criteria are met. In all three cases, the total neutron flux is found by adding the forward and backward components. Comparisons are made between the three methods in one, two and three layer configurations in various space environments and compared to Monte Carlo data where available