A fourth-order finite-difference formula is applied to bidirectional beam-propagation method based on scattering operators. The method is shown to provide significant improvement in accuracy and ...