Numerical methods are devised and tested for solution of the neutron transport equation in one-dimensional plane and two-dimensional R-Z geometry. The numerical methods are devised in a manner such that Richardson extrapolation can be applied to significantly increase the accuracy of the angular fluxes. An analysis of the method is presented for slab geometry transport.