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