acfARp function

Compute autocovariances of an AR(p) process