IAM_VnkNegBin function

Compute the value V(n,k), needed to caclulate the eppf of a Finite Dirichlet process when the prior on the component-weigts of the mixture is a Dirichlet with parameter gamma (i.e. when unnormailized weights are distributed as Gamma(γ\gamma,1) ) when the prior on the number of componet is Negative Binomial with parameter r and pwith mean is mu =1+ r*p/(1-p) TODO: CHECK THIS FORMULA!!!. See Section 9.1.1 of the Paper Argiento de Iorio 2019 for more details