I am working on Gibbs sampler for my z which is Normal distribution and my code is as follows:
z
x <- matrix(1:3000,nrow=300,ncol=10) y<-ma