model{
prec.ybar <- n/sigma2
prec.theta <- 1/tau2
ybar ~ dnorm(theta,prec.ybar)
theta ~ dnorm(mu,prec.theta)
}