An engine test method that causes a computer to execute a process comprising: acquiring, a first test pattern in which each of first operation variables that are used for an engine test is changed in time series; inputting, based on the first test pattern, a first operation variable among the first operation variables to a mathematical model that represents a time series response of an engine obtained by inputting the first test pattern as a simulation of the engine test; monitoring (S102, S103), as a first monitoring parameter of engine abnormality, at least one of an air excess ratio, pressure and temperature of an intake manifold, pressure and temperature of an exhaust manifold, and a maximum cylinder pressure rise rate that are obtained by inputting the first operation variable to the mathematical model; holding, when the first monitoring parameter exceeds a first threshold value that is previously set for the first monitoring parameter obtained from the mathematical model, the first operation variable to a last time value that is input to the mathematical model until the first monitoring parameter is less than the first threshold value; repeating the acquiring, the inputting, the monitoring and the holding during a test time; creating, after the test time ends, a history of the first operation variable that is input to the mathematical model in the simulation as a second test pattern; monitoring (S202, S203), as a second monitoring parameter, at least one of the air excess ratio, the pressure and the temperature of the intake manifold, the pressure and the temperature of the exhaust manifold, and the maximum cylinder pressure rise rate that are obtained by inputting a second operation variable to a real engine based on the second test pattern; holding, when the second monitoring parameter exceeds a second threshold value that is previously set for the second monitoring parameter obtained from the real engine, the second operation variable to a last time value that is input to the real engine until the second monitoring parameter is less than the second threshold value; and acquiring, time series data of the second operation variable and, from the real engine a controlled variable that is obtained by inputting the time series data of the second operation variable to the real engine; and repeating the monitoring, the holding and the acquiring during a test time, wherein the repeating includes storing, during the test time, the input second operation variable and the acquired controlled variable. acquiring, a first test pattern in which each of first operation variables that are used for an engine test is changed in time series; inputting, based on the first test pattern, a first operation variable among the first operation variables to a mathematical model that represents a time series response of an engine obtained by inputting the first test pattern as a simulation of the engine test; monitoring (S102, S103), as a first monitoring parameter of engine abnormality, at least one of an air excess ratio, pressure and temperature of an intake manifold, pressure and temperature of an exhaust manifold, and a maximum cylinder pressure rise rate that are obtained by inputting the first operation variable to the mathematical model; holding, when the first monitoring parameter exceeds a first threshold value that is previously set for the first monitoring parameter obtained from the mathematical model, the first operation variable to a last time value that is input to the mathematical model until the first monitoring parameter is less than the first threshold value; repeating the acquiring, the inputting, the monitoring and the holding during a test time; creating, after the test time ends, a history of the first operation variable that is input to the mathematical model in the simulation as a second test pattern; monitoring (S202, S203), as a second monitoring parameter, at least one of the air excess ratio, the pressure and the temperature of the intake manifold, the pressure and the temperature of the exhaust manifold, and the maximum cylinder pressure rise rate that are obtained by inputting a second operation variable to a real engine based on the second test pattern; holding, when the second monitoring parameter exceeds a second threshold value that is previously set for the second monitoring parameter obtained from the real engine, the second operation variable to a last time value that is input to the real engine until the second monitoring parameter is less than the second threshold value; and acquiring, time series data of the second operation variable and, from the real engine a controlled variable that is obtained by inputting the time series data of the second operation variable to the real engine; and repeating the monitoring, the holding and the acquiring during a test time, wherein the repeating includes storing, during the test time, the input second operation variable and the acquired controlled variable.