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