model {
theta[1] ~ dnorm(0, 1)
theta[2] ~ dnorm(0, 1)
z <- 1
z ~ dbern(constraint)
constraint <- step(theta[1]*theta[1] + theta[2]*theta[2] - 1)
}

[example-9_7_2-doughnut0]