The paper proposes an algorithm of the numerical integration with the modal analysis for computing inelastic seismic responses, and furthermore, the accuracy of the numerical integration with the Newmark’s =1/4 method that is most popular in the earthquake engineering is discussed by comparing with the response computed by the proposed method