We develop a full-vectorial higher-order finitedifference formulation to model optical waveguides. Our proposed scheme yields improved convergence and field calculation results. Applications to beam propagation analysis of a photonic crystal coupler is demonstrated.