n = sp.symbols('n')
F = sp.factorial
f = (F(n) + F(n + 1)) / F(n - 2)
print(f)
f = f.simplify()
print(f)
print(f.subs(n, 6))
f = sp.lambdify(n, f, 'numpy')
a = np.arange(2, 10)
print(f(a))
Эффективное вычисление факториалов с помощью Python и Numpy