In order to achieve high short-term prediction accuracy of ionospheric TEC, first, we transform a seasonal time series for ionospheric Total Electron Content (TEC) into a stationary time series by seasonal differences and regular differences with a full consideration of the Multiplicative Seasonal model. Next, we use the Autoregressive Integrated Moving Average (ARIMA) model taken from time series analysis theory for modeling the stationary TEC values to predict the TEC series. Using TEC data from 2008 to 2012 provided by the Center for Orbit Determination in Europe (CODE) as sample data, we analyzed the precision of this method for prediction of ionospheric TEC values which vary from high to low latitudes during both quiet and active ionospheric periods. The effect of the TEC sample’s length on the predicted accuracy is analyzed, too. Results from numerical experiments show that during the ionospheric quiet period the average relative prediction accuracy for a six day time span reaches up to 83.3% with average prediction residual errors of about 0.18±1.9TECu. During ionospheric active periods it changes to 86.6% with an average prediction residual error of about 0.69±2.6TECu. For the quiet periods, above 90% of predicted residual is less than ±3TECu while during active periods, it is only about 81%. The two periods show that that the higher the latitude, the higher the absolute precision, and the lower the predicted relative accuracy. In addition, the results show that prediction accuracy will improve with an increase of the TEC sample sequences length, but it will gradually reduce if the length exceeds the optimal length, about 30 days. On the other hand, with the same TEC sample, as the predicted days increase, the predictive accuracy decreases. Athough the predictive accuracy is not apparent at the beginning, it will be significantly reduced after 30 days.