Computes the probability mass function at j using a recursive (iterative) algorithm using logarithms
cumulative probability of j
the value to evaluate
number of trials
success probability