A known analytical physical quantity of an analysis subject and an element level consistent mass matrix of each element are obtained by acquiring units. A bubble function is integrated for each element, and the element level diagonal mass matrix of each element is computed by substituting the value that is integrated for the element level consistent mass matrix of each element. The diagonal mass matrix for the entire analysis region is computed by the summation (superposition) of the element level diagonal mass matrices of each element and inverse matrix thereof is computed. A motion of the analysis subject is analyzed based on the known analytical physical quantity, the diagonal mass matrix for the entire analysis region, and the inverse matrix thereof.