model
	{
		for( i in 1 : N ) {
			Y[i] ~ dnorm(mu, prec)
		}
		mu ~ dnorm (0, 0.0001)
		prec ~ dgamma (0.001, 0.001)
		s2 <- 1/prec
	}