Calculational procedures have been implemented for solving importance and higher harmonic neutronics problems. Solutions are obtained routinely to support analysis of reactor core performance, treating up to three space coordinates with the multigroup diffusion theory approximation to neutron transport. The techniques used and some of the calculational difficulties are discussed.