Arima function

Fit ARIMA model to univariate time series