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.