Finite difference methods for the analysis of integrated optical waveguides

A finite difference method based on a Crank–Nicolson discretisation procedure has been implemented and applied to the analysis of propagation of optical beams in dielectric waveguides. It is based on a new formula which takes into account the correct boundary conditions. Comparisons with the well known BPM technique show a definite improvement in computer time saving, accuracy and ease of use.