model {
Z ~ dnorm(0, 1)
Y <- pow(2*Z + 1, 3)
P10 <- step(Y - 10)
}

node mean sd MC error 2.5% median 97.5% start sample
P10 0.2797 0.4488 0.0014 0.0 0.0 1.0 1 100000
Y 12.83 39.78 0.1329 -25.59 0.9584 118.5 1 100000
Z -0.005727 0.9987 0.003153 -1.973 -0.007038 1.956 1 100000