A modular spherical harmonics approach to the neutron transport equation