n <- 100
x <- rgamma(n,3,1)
target <- sum(log(x))/n

alpha <- seq(0.1,10,0.1)

plot(alpha,digamma(alpha))
abline(a=target,b=0)