High-accuracy finite-difference beam-propagation method for cylindrical geometry