Please only use python to finish this problem,Its better add comments for me to learn,thank you very much.

reference [1] Real spherical harmonics, in Table of spherical harmonics - Wikipedia. https://en.wikipedia.org/wiki/Table_of_spherical_harmonics [2] The Radial Wavefunction solutions. https://quantummechanics.ucsd.edu/ph130a/130_notesode233.html

The eigenfunction of an electron in a hydrogen atom is nlm(r,,)=Rnl(r)Ylm(,) which is complex. In many applications, real wavefunctions are used instead, they can be formed by taking linear combinations of the eigenfunctions. 2. (a) Following ref. and other references on atomic orbitals, write down px as linear combinations of spherical harmonics. (For example, dx,y=i21(Y22+Y22).) What is the quantum number l for px and dxy orbitals, respectively? 2. (b) After finding out s,p,d,f correspond to which values of l, and how to form real wavefunctions with different m. Plot the atomic orbitals 2px and 3dx,y. The numbers 2 and 3 are the principal quantum number n. See ref. 2. Plot the intersection of the atomic wavefunction with the xy plane (z=0). Use pseudocolor plot with colormap "bwr". Choose your own range and figure size. Use appropriate units. The constant terms can be ignored and the magnitude of the wavefunction is not important here. Your plots should look like fig. 2 below. FIG. 2 Pseudocolorplots of the atomic orbitals 2px and 3dxy