To provide a technique of obtaining a simulation model with high accuracy.

A first coefficient matrix indicating a relationship between the quantity of state and a load when an analysis object is vibrated in a specific axial direction is derived from the relationship between the quantity of state observed when the analysis object is vibrated in a specific axial direction and a given load. The quantity of state of the analysis object calculated based on the load when the analysis object is vibrated in a composite direction and the first coefficient matrix is removed from the quantity of state of the analysis object observed when the analysis object is vibrated in a composite direction. A second coefficient matrix indicating the relationship between the quantity of state of the analysis object when the analysis object is vibrated in a composite direction and a load is calculated from the relationship between the quantity of state of the removed analysis object and the load given to the analysis object when the analysis object is vibrated in a composite direction. A combination of multiplication of a matrix indicating the quantity of state when the analysis object is vibrated in a specific axial direction by the first coefficient matrix and multiplication of a matrix indicating the quantity of state of the analysis object when the analysis object is vibrated in a composite direction by a second coefficient matrix is derived as a simulation model of the analysis object.