A high precise integration scheme for analysis of dynamical response of the structures with elastic or rigid modes is developed in the paper. The scheme can be used in the excitations with linear as well as polynomial forms. Moreover, if the nonlinear excitation can be presented in the polynomial series, the dynamic response of structure is also obtained by the proposed method, and the computed precision can be adjusted by changing the order of the polynomial.