CylinderCircumscribedSphere::usage="The axis of a cylinder is the line:ℒ: {h, k, l} + t {A, B, C1}. The family of spheres that has its center on that axis is:(x - (h + t A) + {y - (k + t B) + {z - (l + t C1) = , with A, B ≠ 0 and r > 0 , because the radius of each element of the family is r.For t ϵ .The cylinder is circumscribed to the family of spheres. That is, the cylinder is the envelope of the family of spheres.Reference: Personal. ";
2
}
2
}
2
}
2
r
Example 1:
In[]:=
Showh=0;k=0;l=0;A=7;B=-1;C1=5;t1=1;r=3;Solve[++==,(++-)==0,{x,y}]x105t-7z--15t+z-7105t-7z+-15t+z+7105t-7z--15t+z-7+14xy+74-70xz+10yz+500;{Graphics3D[{Opacity[0.8],Red,Specularity[White,40],Sphere[{h+t1A,k+t1B,l+t1C1},r]}],ParametricPlot3D[{h,k,l}+t2{A,B,C1},{t2,-8,8},PlotStyle{DarkGreen,Thick}],ContourPlot3D[-675+26+14xy+74-70xz+10yz+500,{x,-15,15},{y,-15,15},{z,-15,15},PerformanceGoal->"Quality",PlotPoints->60,MeshFalse,ContourStyle{Opacity[0.2],Green,Specularity[White,40]}]},AxesLabel->(Style[#,15,Blue]&/@{"X","Y","Z"}),AxesOrigin{0,0,0},AxesTrue,BoxedFalse,BoxRatiosAutomatic,PlotRange{{-15,15},{-15,15},{-15,15}},ImageSize900,ViewPoint{1.3,-2.4,2.}
2
(x-(h+tA))
2
(y-(k+tB))
2
(z-(l+tC1))
2
r
∂
t
2
(x-(h+tA))
2
(y-(k+tB))
2
(z-(l+tC1))
2
r
1
10
3
6-25+10tz-
,y2
t
2
z
1
10
3
6-25+10tz-
,x2
t
2
z
1
10
3
6-25+10tz-
,y2
t
2
z
1
10
3
6-25+10tz-
;Eliminatex==2
t
2
z
1
10
3
6-25+10tz-
,y==2
t
2
z
1
10
3
6-25+10tz-
,z==z,t-675+262
t
2
z
2
x
2
y
2
z
2
x
2
y
2
z
Out[]=
Example 2:
In[]:=
Showh=2;k=-5;l=8;A=1;B=1;C1=1;r=4;t1=9;Solve[++==,(++-)==0,{x,y}]x12+3t-z--2+3t-z+12+3t-z+-2+3t-z-12+3t-z--2+3t-z++20y-xy+-19z-xz-yz+0;Graphics3D[{Opacity[0.8],Red,Specularity[White,40],Sphere[{h+t1A,k+t1B,l+t1C1},r]}],Graphics3DOpacity[0.8],Red,Specularity[White,40],Sphereh+A,k+B,l+C1,r,ParametricPlot3D[{h,k,l}+t2{A,B,C1},{t2,-8,20},PlotStyle{DarkGreen,Thick}],ContourPlot3D[103-x++20y-xy+-19z-xz-yz+0,{x,-30,30},{y,-30,20},{z,0,25},PerformanceGoal->"Quality",PlotPoints->60,MeshFalse,ContourStyle{Opacity[0.2],Green,Specularity[White,40]}],AxesLabel->(Style[#,15,Blue]&/@{"X","Y","Z"}),AxesOrigin{0,0,0},AxesTrue,BoxedFalse,BoxRatiosAutomatic,PlotRangeAll,ImageSize900,ViewPoint{4,-2.4,4.}
2
(x-(h+tA))
2
(y-(k+tB))
2
(z-(l+tC1))
2
r
∂
t
2
(x-(h+tA))
2
(y-(k+tB))
2
(z-(l+tC1))
2
r
1
2
-160-48t-3+48z+6tz-3
,y2
t
2
z
1
2
-160-48t-3+48z+6tz-3
,x2
t
2
z
1
2
-160-48t-3+48z+6tz-3
,y2
t
2
z
1
2
-160-48t-3+48z+6tz-3
;Eliminatex==2
t
2
z
1
2
-160-48t-3+48z+6tz-3
,y==2
t
2
z
1
2
-160-48t-3+48z+6tz-3
,z==z,t103-x+2
t
2
z
2
x
2
y
2
z
t1
9
t1
9
t1
9
2
x
2
y
2
z
Out[]=