I tried to simulate 200 random normal numbers for a given normal distribution:
How can i create 200 numbers for x1?
I tried 
covMat=matrix(c(1,0,0,0,0,1,rho.xz,rho.xeps,0,rho.xz, 1, rho.zeps,0, rho.xeps, rho.zeps, 1), nrow=4, ncol=4)
rmvnorm(200,sigma=covMat,mean=c(0,0,0,0))
but that doesn't really work for me.
Any ideas?

 
    