Create a function that calculates the normalized radial solution. This function should take 3 float arguments and return the value of the normalized radial solution for the specific n, l, and r. The return value should be normalized to a^(-3/2), where a is the Bohr's radius, and rounded to 5 decimal places. My code so far in Python 3.6 ONLY:

import scipy.constants as c import numpy as np

def radial_wave_func(n,l,r): a=c.physical_constants['Bohr radius'][0]