  We present a simple analytical method to solve master equations for finite
temperatures and any initial conditions, which consists in the expansion of the
density operator into normal modes. These modes and the expansion coefficients
are obtained algebraically by using ladder superoperators. This algebraic
technique is successful in cases in which the Liouville superoperator is
quadratic in the creation and annihilation operators.
