param <- c(0, 1, 5, 3, 1)
maxDens <- dghyp(ghypMode(param = param), param = param)
ghypRange <- ghypCalcRange(param = param, tol = 10^(-3) * maxDens)
ghypRange
curve(dghyp(x, param = param), ghypRange[1], ghypRange[2])
if (FALSE) ghypCalcRange(param = param, tol = 10^(-3), density = FALSE)
Run the code above in your browser using DataLab