changed geometric distribution geo_p calculation. It now returns Nan if sample size is 0, but it gives correct results otherwise