  We present a fast method for numerically solving the inhomogeneous Helmholtz
equation. Our iterative method is based on the Born series, which we modified
to achieve convergence for scattering media of arbitrary size and scattering
strength. Compared to pseudospectral time-domain simulations, our modified Born
approach is two orders of magnitude faster and nine orders of magnitude more
accurate in benchmark tests in 1-dimensional and 2-dimensional systems.
