How to solve it using Python code

3. (4 points) Write a function called derivatives which takes input parameters f, a, and h (with default value h-0.001) and returns approximations of the derivatives f'(a), f"(a), (a) (as a NumPy array) using the formula 2nhn Use either scipy.misc.factorial or scipy.misc.comb to compute n choose k: C) = k!(n-k