model {
  for(i in 1:N) {
    mu[i] <- beta * x[i]
    y[i] ~ dnorm(mu[i],psi)
  } 
  beta ~ dnorm(0.5, 5)
  psi ~ dgamma(3, 1)
} 
