model {
  h~dbin(pi, n)  #data
  pi~dbeta(a, b) #prior
  a <- 3.26 
  b <- 7.19
  pi2 <- pow(pi,2)
  }