The spherical harmonics approximation of lth order, applied to the transport equation in slab geometry, ' leads to 21 + 1 couplod first order differential equations. These may be transformed into (21 + 1)/2 second order differential equations similar, in form, to the few-group diffusion equations, and amenable to solution by wellknown iterative techniques. The double-P/sub 1/ equations of Yvon may be manipulated and solved in the same manner. This article describes an IBM 704 code which makes use of such a method. Some of the results obtained with the code are discussod, and machine times for typical problems are compared with times required to solve the same problem by the discrete ordinate methods. (auth)