Re: ContourPlot3D Problem
- To: mathgroup at smc.vnet.net
- Subject: [mg75230] Re: ContourPlot3D Problem
- From: dennis <dwangsness at earthlink.net>
- Date: Sat, 21 Apr 2007 23:09:14 -0400 (EDT)
- References: <f079vc$4a0$1@smc.vnet.net>
Sorry for the previous post being so difficult to read.
I am trying to plot a cone with a spherical cap. The code I am using
is shown below.
Needs["Graphics`ContourPlot3D`"]
halfAng = 10.*Degree;
sMax = 10.;
axisVec = {1., 1., 1}/Sqrt[3.];
{x0, y0, z0} = sMax*(1 + Tan[halfAng]^2)*axisVec;
rCirc = sMax*(Tan[halfAng]/Cos[halfAng]);
g1 =
ContourPlot3D[
If[{x, y, z} . axisVec <= sMax,
{x, y, z} . axisVec -
Sqrt[x^2 + y^2 + z^2]*Cos[halfAng],
(x - x0)^2 + (y - y0)^2 + (z - z0)^2 - rCirc^2],
{x, 0, 10}, {y, 0, 10}, {z, 0, 10}, PlotPoints -> {8, 6},
ViewPoint -> {2.651, -2.103, -0.026}, PlotRange ->
{{0, 12}, {0, 12}, {0, 10}}];
This produces the cone with a spherical cap, but it also produces a
plane perpendicular to the cone axis at the junction of the cone and
sphere. How can I get the cone-sphere without the plane?
It also produces a bunch of CompiledFunction errors which I don't
understand.
Thanks,
Dennis