The numerical solution of the Helmholtz equation with nonlocal radiation boundary conditions is considered. We model sound wave propagation in a multilayered piecewise homogeneous medium. A fourth-order accurate collocated, boundary element method is used, where the solution inside the domain is computed through a representation integral. The method is shown to have the expected fourth-order co...