用R软件做蒙特卡洛模拟

谁知道怎么样用R软件做下面的蒙特卡洛模拟
模拟1000次，一个AR(1)过程用x[t]=u+r(x[t-1]-u)+sqrt(o.5*(1-r*r))*w公式，其中X0=u，u=1.0
分别对于n=20,50,100,150，的样本空间，以及不同的相关系数r的值，r=0（0.1）0.9



各位大神，跪求啊。