Gold price fluctuations have a significant impact because gold is a haven asset. When financial markets are volatile, investors tend to turn to safer instruments such as gold, so gold price forecasting becomes important in economic uncertainty. The novelty of this research is the comparative analysis of time series forecasting models using ARIMA and the NNAR methods to predict gold price movements specifically applied to gold price data with non-stationary and non-linear characteristics. The aim is to identify the strengths and limitations of ARIMA and NNAR on such data. ARIMA can only be applied to time series data that are already stationary or have been converted to stationary form through differentiation. However, ARIMA may struggle to capture complex non-linear patterns in non-stationary data. Instead, NNAR can handle non-stationary data more effectively by modeling the complex non-linear relationships between input and output variables. In the NNAR model, the lag values of the time series are used as input variables for the neural network. The dataset used is the closing price of gold with 1449 periods from January 2, 2018, to October 5, 2023. The augmented Dickey-Fuller test dataset obtained a p-value = 0.6746, meaning the data is not stationary. The ARIMA(1, 1, 1) model was selected as the gold price forecasting model and outperformed other candidate ARIMA models based on parameter identification and model diagnosis tests. Model performance is evaluated based on the RMSE and MAE values. In this study, the ARIMA(1, 1, 1) model obtained RMSE = 16.20431 and MAE = 11.13958. The NNAR(1, 10) model produces RMSE = 16.10002 and MAE = 11.09360. Based on the RMSE and MAE values, the NNAR(1, 10) model produces better accuracy than the ARIMA(1, 1, 1) model.