need to write in python spyder 

2. Quantum harmonic oscillator. In classical mechanics, a particle has a definite position x. As you may have heard, this is not true in quantum mechanics, which is a more fundamental theory. For a harmonic oscillator (a mass on a spring) the energy is only allowed to have certain discrete values En where n - 0,1,2... and when the energy is E-En the probability of finding the mass at position x is given b p(x)(H( The first factor e-a is a gaussian (bell-shaped curve), while the sec- ond factor is the square of the "Hermite polynomial" Hn(x), one of many special families of functions discovered by mathematicians in the nineteenth century. For this problem, plot p(x) for the four n values 0, 1, 2, and 50. Please make the four plots subplots on one figure, with a title like n- .. on each subplot. Hints: . The function Hn(x) is available as scyp. special. eval-hermite(n,x) The plot of p(x) for each n should cover the range -max sxs I suggest you make your program work for one n at a time before trying to combine the four n's in one figure