A quantum statistical procedure for the study of interacting radiation fields is developed and applied to a model laser operating as a linear amplifier. The method utilizes an expansion of the density matrix in terms of a complete set of orthogonal o...