The output of a nuclear reactor is increased by a predetermined magnitude, and the neutron beam is measured as time-series data. The temperature of the moderator in the reactor is acquired as time-series data. Time-series data on the reactivity is acquired from the time-series data on the neutron beam by the reverse dynamic characteristic method with respect to a one-point reactor kinetics equation. Time-series data on the fuel temperature of a predetermined average acquired by using the time-series data on the reactor output and a predetermined dynamic characteristic model is acquired. The reactivity feedback contribution component is determined by using the time-series data on the reactivity and the applied reactivity. The Doppler reactivity coefficient is determined by using the time-series data on the average temperature of the moderator in the reactor, the time-series data on the fuel temperature of the predetermined average, the isothermal temperature reactivity coefficient, and the reactivity feedback contribution component.