SBC.scala

// From Distribution.scala
def gamma(z: Real): Real = {
  val w = (z + 1) + (Real.one / ((12 * (z + 1)) - (Real.one / (10 * (z + 1)))))
  (Real(Math.PI * 2).log / 2) - ((z + 1).log / 2) + ((z + 1) * (w.log - 1)) - z.log
} 

def main(args: Array[String]): Unit = {
    val _ =
      SBC(Uniform(0.5, 1)) { x =>
        Beta(x, x)
      }.prepare(HMC(5), 1000, 1000)
        .animate(3)
 }