pnbd function

Pareto/NBD models