#The model:
model{
   for(i in 1:n){
 y[i]~dnorm(mu[i],taue)
      mu[i] <- alpha + beta*x[i]

   }
   
     alpha~dnorm(alphamean,taualpha)
   beta~dnorm(betamean,taubeta)
   
alphamean ~ dnorm(0,0.0001)
betamean ~ dnorm(0,0.0001)

taualpha ~ dgamma(1,0.001)
taubeta ~ dgamma(1,0.001)
taue ~ dgamma(1,0.001)

 
}