We have obtained a set of coupled differential equations from the continuous limit of the transfer matrix method. Decoupling such a set of equations yields an extension to the Wentzel-Kramers-Brillouin (WKB) approximation for the Schrödinger equation with a position-dependent effective mass (PDEM). In the classically allowed region, the decoupling is to ignore the reflection resulting from the variations of both the potential and effective mass. By considering an infinite-well example with a PDEM, it is shown that the extended WKB approximation can provide not only an estimation of the eigenenergies, but also an analytic form for the approximate wavefunctions.