Monte Carlo method in nonlinear statistical optics