R.DNBINOM(x,n,psuc,give_log)
x: observation
n: the number of trials
psuc: the probability of success in each trial
give_log: if true, log of the result will be returned instead
This function returns the probability density function of the negative binomial distribution.