In this paper, a novel numerical approach to calculate properties of the integrated optical waveguides is introduced. Based on the calculated propagation constants and optical power transmissions along taper and Y-branch optical waveguides, a briefly review of the approach is given. We compared the novel approach with widely used FD-BPM quantitatively, from which advantages and limitations of the approach are derived.