# Load required library
library(forecast)

# Sales data
sales_data <- data.frame(
  Month = seq(as.Date("2023-01-01"), as.Date("2023-06-01"), by = "months"),
  Sales = c(12000, 15000, 18000, 16000, 20000, 22000)
)

# Convert to time series (monthly data)
sales_ts <- ts(sales_data$Sales, frequency = 12, start = c(2023, 1))

# Fit ARIMA model
arima_model <- auto.arima(sales_ts)

# Forecast next 3 months
forecast_result <- forecast(arima_model, h = 3)

# Display forecast and plot
print(forecast_result)
plot(forecast_result)