ar = c(1.6,-0.64)
ma = c(-1.5,0.56)

y=arima.sim(list(ar=ar,ma=ma),1000)
acf(y,lag.max=10)
points(0:10,ARMAacf(ar=ar,ma=ma,lag.max=10))

pacf(y,lag.max=10)
points(1:10,ARMAacf(ar=ar,ma=ma,lag.max=10,pacf=TRUE))