DandelineSpheresEllipseCylinder::usage="Spheres tangent to a plane and a cylinder, their points of contact with the plane determine the foci of the ellipse in the space formed by the intersection of the plane and the cylinder. References: Wolfram MathWorld: Dandeline Spheres; Personal ";
In[]:=
Manipulatej[x_,y_]:=,-,t;[t_]:=,,t;ζ[t_]:={Cos[t],Sin[t],2};[t_]:={Cos[t],Sin[t],-2};d11[t_]:=NNorm-,0,+γ[t];d22[t_]:=NNorm--,0,-+γ[t];d33[t_]:=NNorm-,0,+[t];d44[t_]:=NNorm--,0,-+[t];d3[t_]:=d11[t]+d22[t];d4[t_]:=d33[t]+d44[t];w=Grid[{{,,2a},{d11[t1],d22[t1],d3[t1]}},FrameAll,ItemSize->15,ItemStyle->15];o=Grid[{{,,2a},{d33[t1],d44[t1],d4[t1]}},FrameAll,ItemSize->15,ItemStyle->15];Showsurf1=ParametricPlot3D[{Cos[u],Sin[u],v},{u,0,2π},{v,-2,2},PlotStyleDirective[Green,Specularity[White,20],Opacity[0.1]],MeshFalse,PerformanceGoal"Quality"],surf2=ParametricPlot3D[{Cos[u]Sin[v],Sin[u]Sin[v],Cos[v]+2},{u,0,2π},{v,0,π},ColorFunction"LightTemperatureMap",PlotStyle->Directive[Opacity[0.7]],MeshFalse,PerformanceGoal"Quality"],surf3=ParametricPlot3D[{Cos[u]Sin[v],Sin[u]Sin[v],Cos[v]-2},{u,0,2π},{v,0,π},ColorFunction"LightTemperatureMap",PlotStyle->Directive[Opacity[0.7]],MeshFalse,PerformanceGoal"Quality"],surf4=Plot3D[j[x,y],{x,-5,5},{y,-5,5},PlotStyleDirective[Yellow,Opacity[0.1]],MeshFalse,PerformanceGoal"Quality"],curv1=ParametricPlot3D[ζ[θ],{θ,0,2π},PerformanceGoal"Quality"],curv2=ParametricPlot3D[[θ],{θ,0,2π},PerformanceGoal"Quality"],curv3=ParametricPlot3D[γ[θ],θ,-[θ],θ,-,0,,0.04,Graphics3DBlack,Ball-,0,-,0.04,Graphics3D[{Red,Ball[N[γ[t1]],0.04]}],Graphics3D[{Black,Ball[{0,0,0},0.04]}],Graphics3D[{Red,Ball[N[[t1]],0.04]}],Graphics3DBlack,Dashed,Thick,LineN[γ[t1]],,0,,Graphics3DBlack,Dashed,Thick,LineN[[t1]],,0,,Graphics3DBlack,Dashed,Thick,LineN[γ[t1]],-,0,-,Graphics3DBlack,Dashed,Thick,LineN[[t1]],-,0,-,AxesLabel->(Style[#,15,Blue]&/@{"X","Y","Z"}),AxesOrigin{0,0,0},MeshFalse,AxesTrue,BoxedFalse,PlotRange{{-2,2},{-2,2},{-3,3}},ViewPoint{1.3,-2.4,4.},ImageSize900,Style["Dandeline Spheres Cylinder Ellipse case:",Bold,Medium],t1,-
3
x;γ[t_]:=t
3
3-
2
t
3
γ
1
t
3
3-
2
t
3
ζ
1
3
2
3
2
3
2
3
2
3
2
3
2
γ
1
3
2
3
2
γ
1
PF
1
PF
2
PF
1
PF
2
ζ
1
3
,t1,ColorFunction->"DarkTerrain",PerformanceGoal"Quality"],curv4=ParametricPlot3D[γ
1
3
,t1,ColorFunction->"DarkTerrain",PerformanceGoal"Quality"],Graphics3DBlack,Ball3
2
3
2
3
2
3
2
γ
1
3
2
3
2
γ
1
3
2
3
2
3
2
3
2
γ
1
3
2
3
2
3
+0.00002,"Value (t)",-3
+0.00001,3
,0.00001,Delimiter,{{w,2,"Vectors radii of Ellipse"}},Delimiter,{{o,2,"Vectors radii of Ellipse"}},ControlPlacementTopOut[]=
