For any function you write, it is a good idea to test it carefully to ensure it works correctly. The basis of one possible test for your Planck function is the well-known Stefan-Boltzmann radiation law which states that

₀ B(T) d = T⁴/

where is the Stefan-Boltzmann constant above. Use the numpy.trapz function to verify that your Planck function reproduces this law. Make a table or plot in Python to demonstrate this