Solve this problem by using python. I need code.

Write a function that can return each of the lint three spherical Bessel functions j_n(x): j_0 (x) = sin x/x j_1 (x) = sin x/x^2 - cos x/x j_1 (x) = (3/z^2 - 1) sin x/x - 3 cos x/x^2 Your function should take as argument a NumPy array x and the order n, and should return an array of the designated order n spherical Bessel function. Take care to make sure that your functions behave properly at x = 0. Demonstrate the use of your function by writing a Python routine that plots the three Bend functions for 0 lessthanorequalto x lessthanorequalto 20. Your plot should look like the one below. Something to think about: You might note that f_1(x) can be written in terms of j_0(x), can be written in terms of j_1(x) and j_0(x). Can you take advantage of this to write a more efficient function for the calculations of j_1(x) and j_2(x)