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

#The data
list(n=10,x=c(1,2,3,4,5,6,7,8,9,10),y=c(6,7,8,4,9,11,12,14,15,19))

#The initial values
list(sigma=1,alpha=0,beta=1)