In the alternating-direction implicit method (ADIM) and the generalized Douglas(GD)scheme, the time-domain beam propagation equation is discretized with a truncation error two orders lower than that of conventional finite-difference scheme.