bgnbd_pmf function

BG/NBD: Probability Mass Function (PMF)