code: def binomial(n,k):
      
      if n==k: return 1
      assert n>k, "Attempting to call binomial(%d,%d)" % (n,k)
      return factorial(n)//(factorial(k)*factorial(n-k))
docs: Binomial coefficient
      >>> binomial(5,2)
      10
      >>> binomial(10,5)
      252