gradient for ols model with AR1 errors
[grad]=ar1_grad(param,y,x)
* param = parameter vector (k x 1)
* y = dependent variable vector (n x 1)
* x = explanatory variables matrix (n x m)
* grad = the gradient
g=ar1_grad(parm,grocer_y,grocer_x) // taken from olsar1.