l=0 H-Atom Wave Functions Set Up of Contour Plots:
(Most output statements have been deleted from this presentation.)

(* H-Atom Wave-Functions: CONTOUR PLOTS *)

(* First, set up all wave-functions through n=3 *)

a0=0.5292

R=Sqrt[x^2+y^2+z^2]/a0

C1=Sqrt[1/(Pi*(a0^3))]

Y1s=C1*Exp[-R]

y=0
z=0
Plot[Evaluate[Y1s^2],{x,0,10}]

-Graphics-

Clear[y]
Clear[z]

Y2s=C1*(Sqrt[2]/8)*(2-R)*Exp[-R/2]

y=0
z=0
Plot[Evaluate[Y2s^2],{x,0,10}]

-Graphics-

Clear[y]
Clear[z]

Y2pz=C1*(Sqrt[2]/8)*R*Exp[-R/2]*(z/Sqrt[x^2+y^2+z^2])

x=0
y=0
Plot[Evaluate[Y2pz^2],{z,0,10}]

-Graphics-

Clear[x]
Clear[y]

Y2px=(C1/8)*R*Exp[-R/2]*(Sqrt[x^2+y^2]/Sqrt[x^2+y^2+z^2])*
(x/Sqrt[x^2+y^2])

General::spell1:
Possible spelling error: new symbol name "Y2px"
is similar to existing symbol "Y2pz".

y=0
z=0
Plot[Evaluate[Y2px^2],{x,0,10}]

-Graphics-

Clear[y]
Clear[z]

Y2py=(C1/8)*R*Exp[-R/2]*(Sqrt[x^2+y^2]/Sqrt[x^2+y^2+z^2])*
(y/Sqrt[x^2+y^2])

General::spell:
Possible spelling error: new symbol name "Y2py"
is similar to existing symbols {Y2px, Y2pz}.

Y3s=(C1*Sqrt[3]/243)*(27-18R+2R^2)*Exp[-R/3]

y=0
z=0
Plot[Evaluate[Y3s^2],{x,0,10}]

-Graphics-

Clear[y]
Clear[z]

Y3pz=(C1*Sqrt[2]/81)*R*(6-R)*Exp[-R/3]*(z/Sqrt[x^2+y^2+z^2])

x=0
y=0
Plot[Evaluate[Y3pz^2],{z,0,10}]

-Graphics-

Clear[x]
Clear[y]

Y3px=(C1/81)*R*(6-R)*Exp[-R/3]*
(Sqrt[x^2+y^2]/Sqrt[x^2+y^2+z^2])*(x/Sqrt[x^2+y^2])

General::spell1:
Possible spelling error: new symbol name "Y3px"
is similar to existing symbol "Y3pz".

Y3py=(C1/81)*R*(6-R)*Exp[-R/3]*
(Sqrt[x^2+y^2]/Sqrt[x^2+y^2+z^2])*(y/Sqrt[x^2+y^2])

Y3dz2=(C1*Sqrt[6]/486)*(R^2)*Exp[-R/3]*
((3*(z^2)/(x^2+y^2+z^2))-1)

(* Mathematica can't deal with 1/R^2 in expression
above, therefore rewriting Y3dz2 *)

Y3dz2=(C1*Sqrt[6]/486)*Exp[-R/3]*
(2*z^2-x^2-y^2)

x=0
y=0
Plot[Evaluate[Y3dz2^2],{z,0,10}]

-Graphics-

Clear[x]
Clear[y]

Y3dxz=(C1/81)*(R^2)*Exp[-R/3]*
(Sqrt[x^2+y^2]/Sqrt[x^2+y^2+z^2])*
(z/Sqrt[x^2+y^2+z^2])*(x/Sqrt[x^2+y^2])

y=0
Plot3D[Evaluate[Y3dxz^2],{x,-10,10},{z,-10,10}]

-SurfaceGraphics-

Clear[y]

Y3dyz=(C1/81)*(R^2)*Exp[-R/3]*
(Sqrt[x^2+y^2]/Sqrt[x^2+y^2+z^2])*
(z/Sqrt[x^2+y^2+z^2])*(y/Sqrt[x^2+y^2])

General::spell1:
Possible spelling error: new symbol name "Y3dyz"
is similar to existing symbol "Y3dxz".

Y3dx2y2=(C1/162)*(R^2)*Exp[-R/3]*
((x^2+y^2)/(x^2+y^2+z^2))*
((x^2-y^2)/(x^2+y^2))

z=0
Plot3D[Evaluate[Y3dx2y2^2],{x,-10,10},{y,-10,10}]

-SurfaceGraphics-

Clear[z]

Y3dxy=(C1/162)*(R^2)*Exp[-R/3]*
((x^2+y^2)/(x^2+y^2+z^2))*
(2x*y/(x^2+y^2))

General::spell1:
Possible spelling error: new symbol name "Y3dxy"
is similar to existing symbol "Y3dxz".

H-atmWF's;SetUp.contourPlot3D.JKB

---

This ends part 2 of the L=0 H-Atom Wave Functions Mathematica files. You can now go on to the 3rd and final file, l=0 H-Atom Wave Functions Contour Plots in 3-D, or return to the notes.

---