// prototype for random Poisson function
int ranpois (double,int&);