def f(x):
    return x ** 3 + 2 * x + 5

def f_prim(a_):
    h = 0.001
    return (f(a_ + h) - f(a_)) / h

a = 4
print(f"f'({a}) = {round(f_prim(a), 1)}")