F binomial_coeff(n, k)
   V result = 1
   L(i) 1..k
      result = result * (n - i + 1) / i
   R result

print(binomial_coeff(5, 3))