Optical bistability due to a two-photon absorption resonance with fluctuations

Using the intensity-dependent complex dielectric function for a two-photon absorption resonance we derive the Langevin equation for the fluctuating light-field in the non-linear resonator. The corresponding Fokker-Planck equation is solved by expanding the distribution function in terms of products of trigonometric functions and generalized Laguerre polynomials. The expansion coefficients are calculated using the method of matrix continued fractions. Numerical results for the stationary case are given.