This paper introduces a numerical method for integration of the linear and nonlinear differential dynamic equation of motion. The variation of acceleration in two time steps is approximated as a combination of both trigonometric cosine and hyperbolic cosine functions with weighted coefficient. From which all necessary formulae are elaborated for the direct integration of the governing equation. A number of linear and nonlinear dynamic problems with various degrees of freedom are analysed using both the suggested method and Newmark method for the comparison. The numerical results show high advantages and effectiveness of the new method.