source("hdr.txt")
priormean <- 0.2
priorworth <- 9
a0 <- priormean*priorworth
a0 <- round(a0)
b0 <- priorworth - a0
n <- 12
x <- 2
a <- a0 + x
b <- b0 + n - x
cat(hbeta(0.9,a,b),"\n")