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.