Can You write this in python Please

(a) The Planck's theory of thermal radiation tells us that the amount of thermal energy radiated electromagnetically by a black body per unit area is W=42c23kB4T40ex1x3dx where kB is Boltzman's constant, is Planck's constant divided by 2,c is the speed of light, and T is temperature. We cannot calculate this integral analytically. So, your task is to write a program to evaluate it numerically using the Simpson's method (by implementing your own Simpson's function). Do not worry about the temperature or the other constants for this part, focus on the integral only. Hints: See lecture slide \#156 for calculating integrals with infinite limits. First define a function that takes z and returns x. You may find that z=0 or 1 will give you a hard time when it comes to your Simpson's function. Exclude them by shifting your interval by a very small number (say, 1e12) first to the left and then to the right. One of these options should get you a pretty accurate answer, even with 100 slices (I find 6.49394). The other one should get you a wildly inaccurate number. Explain why. (b) Check your results against the output of scipy's simps function. See slide \#155. (c) Even before Planck came up with this formula, it was known that total energy radiated by a black body per unit area per second has the form W=T4, where is the Stefan-Boltzman constant. Use your value for the integral above to compute a value for the Stefan-Boltzman constant in SI units to three significant figures and compare it to the known value (you can find this online)