Abstract: In heating control method of a heat furnace, three non-linear models, model to calculate furnace temperature based on the fuel flow rate by unsteady heat balance system, model to estimate furnace wall temperature distribution based on the furnace temperature and model to estimate material temperature based on the furnace temperature, are used, linearization is performed by perturbation simulation to vary the flow rate stepwise and the optimization to minimize the fuel flow rate is performed and the optimum furnace temperature per each material is determined, and the mixed combustion ratio of plural fuels obtained by dividing total calorific value of each fuel by total calorific value of all fuels and the furnace temperature setting value are calculated using the optimum furnace temperature per material and then set.