model{
for(i in 1:5){
sample[i,1:5] ~ dmnorm(mean0[],tau[,]) }
for(a in 1:5){
sigma[a,a]<-1}
sigma[1,2]~dnorm(0,1)I(-1,1)
sigma[2,1]<-sigma[1,2]
sigma[1,3]~dnorm(0,1)I(-1,1)
sigma[3,1]<-sigma[1,3]
sigma[1,4]~dnorm(0,1)I(-1,1)
sigma[4,1]<-sigma[1,4]
sigma[1,5]~dnorm(0,1)I(-1,1)
sigma[5,1]<-sigma[1,5]
sigma[2,3]~dnorm(0,1)I(-1,1)
sigma[3,2]<-sigma[2,3]
sigma[2,4]~dnorm(0,1)I(-1,1)
sigma[4,2]<-sigma[2,4]
sigma[2,5]~dnorm(0,1)I(-1,1)
sigma[5,2]<-sigma[2,5]
sigma[3,4]~dnorm(0,1)I(-1,1)
sigma[4,3]<-sigma[3,4]
sigma[3,5]~dnorm(0,1)I(-1,1)
sigma[5,3]<-sigma[3,5]
sigma[4,5]~dnorm(0,1)I(-1,1)
sigma[5,4]<-sigma[4,5]
tau[1:5,1:5]<-inverse(sigma[,])
}