Abstract
We have solved a master equation to obtain the diagonal elements of the density matrix for laser light. The master equation represents a generalization of a master equation derived earlier for a single radiation mode interacting with N stationary two-level atoms. The generalization takes into account the pumping scheme which characterizes a three-level laser. For positive temperatures and a lossless cavity, the equilibrium solution of the master equation gives the correct statistical-mechanical description of the atoms and radiation. For negative temperatures, the distribution function for the number of photons in the mode is either exponential or peaked, depending on whether the laser is operating below or above threshold. The calculated intensity fluctuations are in good agreement with semiclassical results for lasers operating slightly above threshold.