model initial2;
{
tau~dgamma(0.001,0.001)
tau2<-tau*0.0001
mu~dnorm(0.0,tau2)
for(i in 1:n) {
y[i] ~dnorm(mu,tau)

}

sigma <-1.0/sqrt(tau)

}

#initial values

list(mu=404.59, tau=0.04)

#data

list(y= c(375, 392, 393, 397, 398, 398, 399, 399, 399, 399, 399, 399, 399, 400, 400, 400, 400, 401,
          401, 401, 401, 401, 401, 401, 401, 401, 401, 401, 401, 402, 402, 402, 402, 402, 402, 402,
          402, 403, 403, 403, 403, 403, 403, 404, 404, 404, 404, 404, 404, 404, 404, 404, 405, 405,
          405, 405, 405, 406, 406, 406, 406, 406, 406, 406, 406, 406, 406, 406, 406, 407, 407, 407,
          407, 407, 407, 407, 407, 408, 408, 408, 408, 408, 409, 409, 409, 409, 409, 410, 410, 410,
          410, 411, 412, 412, 412, 413, 415, 418, 423, 437), n=100)