Added P[ln==k|n]

Added function to calculate the probability that the length of the
longest run in k bernoulli trials is equal to k.