model{
 for (i in 1:n1)
  {
     data1[i] ~ dnorm(mu1,tau)
  }
 for (i in 1:n2)
  {
     data2[i] ~ dnorm(mu2,tau)
  }

 mu1 ~ dnorm(0,0.001)
 mu2 ~ dnorm(0,0.001)
 diff <- mu1 - mu2
 tau ~ dgamma(0.01,0.01)
sigma <- 1/tau
}

list(mu1=0,mu2=0,tau=1.0)


list(n1=20,n2=20)

data1[ ]  data2[ ]
4.81	4.17
4.17	3.05
4.41	5.18
3.59	4.01
5.87	6.11
3.83	4.1
6.03	5.17
4.98	3.57
4.9	5.33
5.75	5.59
5.36	4.66
3.48	5.58
4.69	3.66
4.44	4.5
4.89	3.9
4.71	4.61
5.48	5.62
4.32	4.53
5.15	6.05
6.34	5.14
END