set thin=50 by default as auto-correlation is always present
changed demo to parameters which allow them to be easier identified
Pareto/CNBD now also runs on Unix/Mac - replaced the numerical integration in C++ code, which caused R sessions to terminate sometimes, with simply integration by Simpson 3/8 rule