To forecast the ARIMA model, we want a model with all coefficients being known. After parameter estimation, the fitted model is reported as estModel. So the program will work if you replace "model" by "estModel":
[Y,YMSE,V] = forecast(estModel, 48, 'X0', X, 'Y0', Y(100:length(Y)), 'XF', XF);
Best Answer