ggomnbd function

Gamma-Gompertz/NBD model