The alternate-direction implicit finite difference beam propagation method (FD-BPM) is used to analyze the two-dimensional (2D) symmetrical multimode interference (MMI) couplers. The positions of the images at the output plane and the length of multimode waveguide are accurately determined numerically. In order to reduce calculation time, the parallel processing of the arithmetic is implemented by the message passing interface and the simulation is accomplished by eight personal computers.