A solution of the neutron transport equation using spherical harmonics

A solution of the neutron transport equation is obtained by expanding the flux as a series. Preliminary investigations in one dimension indicated that the first-order differential equations resulting for the unknown coefficients or moments could be solved by eliminating terms with odd L (L = order of Legendre polynomial) to give a second-order system. FORTRAN subroutines have been written to calculate the necessary coefficients and specify the relevant differentials. A finite-difference or finite-element approximation can then be used. (U.K.)