Applying the Milburn equation to describe intrinsic decoherence, we study the interaction of threecoupled quantum harmonic oscillators or quantized fields. We give an explicit solution for the complete equation, i.e., beyond the usual second order approximation used to arrive to the Lindblad form. Then we calculate the expectation value of the number operator of each oscillator or mode for one of the modes given in an initial coherent state.