model;
{
   pi ~ dbeta(0.01,0.01)
   rho ~ dbeta(0.01,0.01)
   diff <- pi - rho
   p <- step(diff)
   x ~ dbin(pi,m)
   y ~ dbin(rho,n)
}

data;
list(m=20,x=8,n=110,y=48)