Re: 3D plotting problem

*To*: mathgroup at smc.vnet.net*Subject*: [mg51755] Re: 3D plotting problem*From*: "Steve Luttrell" <steve_usenet at _removemefirst_luttrell.org.uk>*Date*: Mon, 1 Nov 2004 02:52:31 -0500 (EST)*References*: <cm20mr$g60$1@smc.vnet.net>*Sender*: owner-wri-mathgroup at wolfram.com

This is the sort of thing you want to do: <<Graphics`ContourPlot3D` fz[x_,y_,z_,a_]:=z*Exp[-a*Sqrt[x^2+y^2+z^2]] ContourPlot3D[Abs[fz[x,y,z,1]]-0.3,{x,-1.5,1.5},{y,-1.5,1.5},{z,-3,3}, PlotPoints->{7,7}]; Comments: 1. Your notation for the 2p orbital function (the left hand side of your definition) could not possibly have worked in Mathematica so I have rewritten it. 2. The surface in ContourPlot3D is defined as Abs[fz[x,y,z,1]]-0.3 == 0. I have used Abs so you get both of the lobes of the 2p orbital in the same plot. 3. PlotPoints->{7,7}] gives you sufficient resolution for a good quality plot. If the resolution is set too low then the plot might not appear at all! 4. Play around with the various parameters to get whatever plot you actually want. Steve Luttrell "Henning Heiberg-Andersen" <Henning.Heiberg-Andersen at fi.uib.no> wrote in message news:cm20mr$g60$1 at smc.vnet.net... > > Hi, > > I find myself unable to make a 3D plot of the atomic > 2p orbital > > f_z(2p,x,y,z;a)=z*Exp[-a*Sqrt[x^2+y^2+z^2]]. > > As far as I can understand, 'ContourPlot3D' is the > only option, but it gives either an empty plot box or > a flat sheet. > > Can anyone suggest a way out? > > Best regards, > Henning Heiberg-Andersen > > > > > >