model;
{
   for( i in 1 : n ) {
      Dec[i] ~ dnorm(mu[i],psi)
   }
   for( i in 1 : n ) {
      mu[i] <- alpha + beta * Nov[i]
   }
   alpha ~ dnorm( 0.0,1.0E-6)
   beta ~ dnorm( 0.0,1.0E-6)
   psi ~ dgamma(0.001,0.001)
   phi <- 1/psi
}


data;
list(n=10,Nov=c(23.9,43.3,36.3,40.6,57.0,52.5,46.1,142.0,112.6,23.7),
          Dec=c(41.0,52.0,18.7,55.0,40.0,29.2,51.0,17.6,46.6,57.0))