Beam Propagation Method based on the Iterated Crank-Nicolson Scheme for Solving Large-scale Wave Propagation Problems