pnbd_expectation function

Pareto/NBD: Unconditional Expectation