<t>Hello,Here is my answer to my problem. Change: intg2 <- integrate(nigpdf, iblower, ibupper)$value .... intg2 <- integrate(nigpdf, iblower, ibupper)$valueTo: intg2 <- integrate(nigpdf, iblower, ibupper, alpha=alpha, beta=beta, delta=delta, mu=0)$value .... intg2 <- integrate(nigpdf, iblower, ibupp...