Univariate Weibull Count Probability with gamma heterogeneity
dWeibullgammaCount_mat( x, shape, r, alpha, time = 1, logFlag = FALSE, jmax = 100L )
r, alpha
double gamma parameters