In this paper the method of density matrix is used to describe two coupled modes. In the case of inverse coupling the density matrix equation of motion is modified by power coupling coefficients caused by random imperfections. The exact analytical solution is obtained.