A demand forecasting procedure based on ARIMA (p, d, q) model is developed. This procedure consists of 5 steps. They are (1) determination of (p, d, q) (2) estimation of initial parameter values (3) parameter estimation (4) test of the estimated model, and (5) adaptive forecasting. An algorithm based on the modified Gauss-Newton method is proposed, as a parameter estimation method.