For a prototype Hamiltonian describing a driven, dissipative bosonic mode in a Kerr-like medium, the exact matrix elements of the reduced density matrix are obtained from a generating function in terms of the normal characteristic functions. The approach is based on the Heisenberg equations and operator calculus. The special and limiting cases are discussed.