param <- c(2.5, 0.5, 5)
maxDens <- dgig(gigMode(param = param), param = param)
gigRange <- gigCalcRange(param = param, tol = 10^(-3) * maxDens)
gigRange
curve(dgig(x, param = param), gigRange[1], gigRange[2])
if (FALSE) gigCalcRange(param = param, tol = 10^(-3), density = FALSE)
Run the code above in your browser using DataLab