expr = (factorial(365, evaluate=False) / factorial(365 - x)) / 365**x
f = lambdify(x, expr)
f(5)
# Expected:
## nan