import sympy as sym

x = sym.symbols('x')
f = x**3 - 2*x - 1

dfdx = sym.diff(f, x, 1)

print("f'(x) =", dfdx)