The eigenvalue equations of multiple waveguides with step-index profile are derived by using the WKB theory. Phase changes unique to step-index discontinuity areintroduced when applying the WKB connection formula to turning points. The transfer matrix method is employed for the analysis of multiple structure and the derived eigenvalue equation are represented in the recursive form. The results by the WKB are compared with those by the FEM for a three-waveguide coupler.