A finite difference scheme based on the polynomial interpolation is constructed to solve the quasi-vector equations for optical waveguides with step-index profiles. The discontinuities of the normal components of the electric field across abrupt diel