#The model:
model{
   for(i in 1:n){
       y[i]~dnorm(mu[i],tau) 
       mu[i] <- alpha + beta*x[i]
   }
   tau<-1/sigma*sigma
   sigma~dunif(0,20)             
   alpha~dnorm(0,0.01)
   beta~dnorm(0,0.01)
}