ar1_ma0_sim function

Slightly faster arima.sim implementation which assumes AR(1) and ma=0.