integral floor(x)((x! - x Γ(1 + x) log(x) polygamma(0, 1 + x))/(x (x!)^2)) from 0 to 20