Milne’s half-space problem: A numerical solution of the related integral equation

The one speed, isotropic scattering, planar Milne’s problem in neutron transport has been a test bed for many methods. We further explore here a numerical solution of the related integral equation. We find that with the use of singularity subtraction, high order quadratures, some interval transformations, and very high precision in all computations, it is possible to obtain results of benchmark accuracy in a straightforward manner.