[en] A method is described in which the reflector regions in neutron diffusion codes may be approximately represented by the use of a matrix boundary conditions relating currents and fluxes at the core/reflector interface. The method is incorporated into the EROS computer code on the IBM 370 computer. (author)