An algorithm for solving nonlinear Heisenberg operator equations by extending the solution into series is presented. it is used for solving quantum nonlinear equations. The method allows fast analysis of radiation statistics if the Hamiltonian of the process is known. The method can be used for solving a set of Heisenberg equations in quantum optics and other fields of physics.