This paper presents a precise and stable time integration method for dynamic analysis of vibration or multibody systems. A total system is divided into several subsystems and their responses are calculated separately, while the coupling effect is treated equivalently as constant force during time steps. By using iterative procedure to improve equivalent coupling forces, a precise and stable solution is obtained. Some examples such as a seismic response and multibody analyses were carried out to demonstrate its usefulness.