So, I’m trying to compute the mean of a distribution, it goes smoothly in python and the results match the article I’m reading, but it doesn’t work in R, probably because of the negative inputs for the gamma function since b-i is going to be negative most of the time. Is there any package or modification I could do to make it works? The resources that I have for this distribution are mostly in R, so I don’t really need it in python, just saw it was working there.