Computational method for intrinsic optical bistability: Steady‐state solutions and their stability

An efficient numerical algorithm is described for solving problems with coupled nonlinear differential equations with boundary conditions applied at opposite ends of the medium. It is applied to a problem of optical bistability with coupled Maxwell equations and heat‐flow equations, and associated boundary conditions. It is shown that the numerical procedure is well suited to perform a stability analysis of the solution.