  Considering well localized atom, the mean field theory (MFT) was applied to
replace the operators by equivalent expectation values. The Rabi model was
reduced to a fourth orders NDE describing atoms position. Solution by the
harmonic balance method (HBM) showed good accuracy and consistency to the
numerical results, which introduces it as a useful tool in the quantum dynamics
studies.
