Two functions for Generalized SARIMA time series simulation