d orbitals

We start by defining the  d orbitals

dxz

r3dxz[theta_ , phi_ ] :=
Abs[ Sin[ theta ] Cos[ theta ] Cos[ phi ] ]

ParametricPlot3D[ {r3dxz[ theta, phi ] Sin[ theta ] Cos[ phi ],
r3dxz[ theta, phi ] Sin[ theta ] Sin[ phi ],
r3dxz[ theta, phi ] Cos[ theta ]}, { theta, 0, Pi },
{ phi, - Pi, Pi }, ViewPoint->{4, 4, 1.530} ,
Boxed -> False, AspectRatio -> Automatic];

dx2y2

r3dx2y2[theta_ , phi_ ] :=
Abs[ Sin[ theta ] Sin[ theta ]  Cos[ 2 phi ] ]

ParametricPlot3D[ {r3dx2y2[ theta , phi ] Sin[ theta ] Cos[ phi ],
r3dx2y2[  theta , phi  ] Sin[ theta ] Sin[ phi ],
r3dx2y2[  theta , phi  ] Cos[ theta ]},
{ theta, 0, Pi }, { phi, - Pi, Pi },
ViewPoint->{4.000, 4.000, 1.530},
Boxed -> False, AspectRatio -> Automatic,
PlotRange -> { {-1,1},{-1,1}, {-.5,.5}}  ];

another view

Show[ %,ViewPoint->{4.000, 4.000, 5.530}];

Created by Mathematica  (September 7, 2006)