Computes the cdf at j using a recursive (iterative) algorithm using logarithms
the probability
the value to be evaluated
number of successes
the probability, must be in range (0,1)