model
				{
					x ~ dunif(-1, 1)
					y ~ dunif(-1, 1)
					O <- 0
					O ~ dbern(constraint)
					constraint <- step(x * x + y * y - 1)
				}