“…The resulting operators not only give evanescent waves a desired damping, but also allow more accurate approximations to the Helmholtz equation than real Padé(m,n) approximant operators. In this work, rather than using real Padé approximant operators as the existing approaches the beam propagation the equation for the analysis of non-paraxial propagation in nonlinear waveguides is based on the modified Padé(1,1) approximant operator as follows: [4] Furthermore, by adding an extra calculation step for dealing with nonlinear effects as well-done for a different method in [5] and using PMLs to absorb incident radiation without any additional parasitic reflections, regardless of wavelength, incident angle or polarization as boundary conditions, the CJI method is thus a promising technique to solve the nonlinear beam propagation equation quickly and effectively.…”