Matlab implementation of high accuracy Fast-Fourier-transform-based simulations of atmospheric phase screens