model
{
   for( i in 1 : 30 ) {
      x[i] ~ dnorm(mu1,sig1)
   }
   for( i in 1 : 30 ) {
      y[i] ~ dnorm(mu2,sig2)
   }
   mu1 ~ dnorm(50,1.0E-6)
   sig1 ~ dgamma(0.001,0.001)
   mu2 ~ dnorm(50,1.0E-6)
   sig2 ~ dgamma(0.001,0.001)
   lbda <- mu1 - mu2
}
