MathGroup Archive 2007

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: ContourPlot3D Problem

  • To: mathgroup at smc.vnet.net
  • Subject: [mg75241] Re: ContourPlot3D Problem
  • From: Roger Bagula <rlbagula at sbcglobal.net>
  • Date: Sat, 21 Apr 2007 23:14:53 -0400 (EDT)
  • References: <f079vc$4a0$1@smc.vnet.net>

dennis wrote:

>I am trying to plot a cone with a spherical cap.  The code I am using
>is shown below.
>
>=CE=B8 = 10.*Degree;
>sMax = 10;
>a = {1., 1., 1}/Sqrt[3.];
>{x0, y0, z0} = sMax*(1 + Tan[=CE=B8]^2)*a;
>rCirc = sMax*(Tan[=CE=B8]/Cos[=CE=B8]);
>g1 = ContourPlot3D[
>If[{x, y, z} . a <= sMax, {x, y, z} . a -
>      Sqrt[x^2 + y^2 + z^2]*Cos[=CE=B8], (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?
>
>Thanks,
>Dennis
>
>
>  
>
Your code doesn't work here either.
Try plotting it as g1 and g2:
g1-> cone
{x, y, z} . a -
      Sqrt[x^2 + y^2 + z^2]*Cos[ang]
g2->sphere
(x - x0)^2 + (y - y0)^2 + (z - z0)^2 - rCirc^2
Show[{g1,g2}]
Doing this isolates the plane to the cut off cone section.
You have to specify that:
 Sqrt[x^2 + y^2 + z^2]*Cos[ang]<>10
My code works worse than yours but elementates the plane 
((0.5773502691896258` x + 0.5773502691896258` y + 0.5773502691896258` 
z-10==0)

\!\(\*
  StyleBox[\(g1\  = \
    ContourPlot3D[If[0.5773502691896258`\ x + 0.5773502691896258`\
          y + 0.5773502691896258`\ z > 10, \n103.10912041257637`\
\[InvisibleSpace] - 11.906015685221298`\ x +
              x\^2 - 11.906015685221298`\ y + y\^2 - 
11.906015685221298`\ z +
           z\^2, 0], \ {x, \ 0, \ 10}, \ {y, \ 0, \ 10}, \ {z, \ 0, \ 
10}, \ \
PlotPoints\  -> \ {8, \n6}, \n\ \ \ \ ViewPoint\  -> \ {
            2.651, \ \(-2.103\), \ \(-0.026\)}, \ PlotRange\  -> \n\ \ \ 
\ \ \
{{0, \ 12}, \ {0, \ 12}, \ {0, \ 10}}]\),
    FontFamily->"Lucida Grande",
    FontSize->13]\)

\!\(\*
  StyleBox[\(g2\  = \
    ContourPlot3D[If[\((
      0.5773502691896258`\ x + 0.5773502691896258`\
              y + 0.5773502691896258`\ z > 10)\) && 
\((0.984807753012208`\ \@\
\(\(x\^2\)\(+\)\(y\^2\)\(+\)\(z\^2\)\(\ \ \ \ \ \ \ \ \)\))\) <> 10, \n0, \
0.5773502691896258`\ x +
              0.5773502691896258`\ y + 0.5773502691896258`\ z - \
0.984807753012208`\ \@\(x\^2 + y\^2 + z\^2\)], \ {x, \ 0, \ 10}, \ {y, \ 
0, \ \
10}, \ {z, \ 0, \ 10}, \ PlotPoints\  -> \ {8, \n6}, \n\ \ \ \ 
ViewPoint\  -> \
\ {2.651, \ \(-2.103\), \ \(-0.026\)}, \ PlotRange\  -> \n\ \ \ \ \ {{0, \ \
12}, \ {0, \ 12}, \ {0, \ 10}}]\),
    FontFamily->"Lucida Grande",
    FontSize->13]\)

Show[{g1, g2}]

The plane shows up when:
0.984807753012208*Sqrt[x^2+y^2+z^2]==0


  • Prev by Date: Re: set versus set-delay
  • Next by Date: Re: set versus set-delay
  • Previous by thread: Re: ContourPlot3D Problem
  • Next by thread: Complex bessel function