model{
   for(i in 1:n){
      y[i]~dnorm(mu[i],phi) 
      mu[i] <- alpha + beta*x[i]
   }
   phi~dgamma(0.001,0.001)
   alpha~dnorm(0,0.001)
   beta~dnorm(0,0.001)
}