Missing data in univariate time series are a recurring problem causing bias and leading to inefficient analyses. Most existing statistical methods which address the missingness problem do not consider the characteristics of the time series when imputing the missing values and, most of all, do not allow the imputation in a univariate time series context. Moreover, just a few methods can be applied to all missing data patterns. Finally, no intuitive procedure addressing the missingness obstacle exists in the literature. In this work of investigation, an algorithm having the aim of filling in these gaps is presented; its main purpose is to find a procedure that gives reliable imputations of the missing values, i.e. not far from the true ones. To this aim, the reliability and robustness of the algorithm have been tested through the simulations campaigns approach. Its innovative feature is the combination of the ARMA models, used to impute the missing values through a forecast and a backcast approach, and the Expectation-Maximization algorithm, used to achieve the parameters convergence. This approach was evaluated through the RMSE and the MAPE metrics, which showed that the algorithm can be used in almost every model setting among the tested ones, with a good reliability. However, one of the main limitations of the introduced procedure is that the nonconvergence of the algorithm could bring to biased imputations. The algorithm can be applied step by step by a common analyst, in a more intuitive way than the majority of other existing approaches.