  In this work we construct an approximate time evolution operator for a system
composed by two coupled Jaynes-Cummings Hamiltonians. We express the full time
evolution operator as a product of exponentials and we analyze the validity of
our approximations contrasting our analytical results with those obtained by
purely numerical methods.
