If a gas flow rate and a pressure in the furnace are set with reference to data of an interstitial oxygen concentration in a crystal pulled in a previous time, a portion in which a change in oxygen concentration is large can be stabilized to a desired value for certain and further, operating conditions can also be efficiently determined when various kinds of operating conditions are selected in or