AreaCylinderSpherePlaneXY::usage="The surface area between the cylinder (x-5
2
)
+
2
y
= 25 and the sphere
2
x
+
2
y
+
2
z
= 100,​ above the XY plane, it will be obtained. Reference: Personal.";
In[]:=
Show[{​​ParametricPlot3D[{10Cos[u]Sin[v],10Sin[u]Sin[v],10Cos[v]},{u,0,2π},{v,0,π},PlotStyle->Directive[Opacity[0.7],Specularity[White,30]],MeshFalse,ColorFunction->"BlueGreenYellow",PerformanceGoal->"Quality",PlotPoints->40],ParametricPlot3D[{5Cos[t]+5,5Sin[t],v},{t,0,2π},{v,-10,10},PlotStyle->Directive[Opacity[0.7],Specularity[White,30]],MeshFalse,ColorFunction->"BlueGreenYellow",PerformanceGoal->"Quality",PlotPoints->40],ContourPlot3D[z==0,{x,-10,10},{y,-10,10},{z,-5,5},MeshNone,ContourStyle->{Green,Opacity[0.8]}]},​​AxesLabel->(Style[#,15,Blue]&/@{"X","Y","Z"}),AxesOrigin{0,0,0},AxesTrue,BoxedFalse,BoxRatiosAutomatic,PlotRangeAll,ImageSize700,ViewPoint{1.3,-2.4,1}]
Out[]=
Development for surface area:
Method 1:
Viviani’s vault:
In[]:=
​​f[x_,y_]:=
100-
2
x
-
2
y
;​​
f
x
[x_,y_]=D[f[x,y],x];​​
f
y
[x_,y_]=D[f[x,y],y];​​
1+
2
f
x
[x,y]
+
2
f
y
[x,y]
//FullSimplify
10
-
1
-100+
2
x
+
2
y
;
In[]:=
x=rCos[θ];y=rSin[θ];​​
2
x
+
2
y
-10x==0;​​
2
r
2
Cos[θ]
+
2
r
2
Sin[θ]
-10rCos[θ]==0;​​r(r-10Cos[θ])==0;​​0<=r<=10Cos[θ];​​-
π
2
<=θ<=
π
2
;​​
A
1
=2(10)
π
2
∫
0
10Cos[θ]
∫
0
-
1
-100+
2
r
rrθ
100(-2+π);
Cylinder base:
A
2
=π
2
(5)
=25π;
Portion of the cylinder:
In[]:=
Solve[{
2
(x-5)
+
2
y
==25,
2
x
+
2
y
+
2
z
==100},{y,z}]
y-
10x-
2
x
,z-
10
10-x
,y-
10x-
2
x
,z
10
10-x
,y
10x-
2
x
,z-
10
10-x
,y
10x-
2
x
,z
10
10-x
;
In[]:=
g[x_,z_]:=
25-
2
(x-5)
;​​
g
x
[x_,z_]=D[g[x,z],x];​​
g
z
[x_,z_]=D[g[x,z],z];​​
1+
2
g
x
[x,z]
+
2
g
z
[x,z]
//FullSimplify​​
5
1
10x-
2
x
;
By symmetry:
In[]:=
A
3
=2(5)
10
∫
0
10
10-x
∫
0
1
10x-
2
x
zx
200;
Total area:
A
T
=100(-2+π)+25π+200;
125π;
Method 2:
In[]:=
ℜ=ImplicitRegion[
2
(x-5)
+
2
y
<=25&&
2
x
+
2
y
+
2
z
<=100&&z>=0,{x,y,z}];​​SurfaceArea[ℜ]
125π;
Graph:
In[]:=
Solve[{
2
(x-5)
+
2
y
==25,
2
x
+
2
y
+
2
z
==100},{x,y}]
x
1
10
(100-
2
z
),y-
1
10
100
2
z
-
4
z
,x
1
10
(100-
2
z
),y
1
10
100
2
z
-
4
z
;
In[]:=
Reduce[100
2
z
-
4
z
>=0,z]
-10≤z≤10;
Only for z ≥ 0:
0≤z≤10;
In[]:=
gr=Showν[t_]:=
1
10
(100-
2
t
),-
1
10
100
2
t
-
4
t
,t;
ν
1
[t_]:=
1
10
(100-
2
t
),
1
10
100
2
t
-
4
t
,t;​​{RegionPlot3D[
2
(x-5)
+
2
y
<=25&&
2
x
+
2
y
+
2
z
<=100&&z>=0,{x,0,10},{y,-5,5},{z,0,10},Mesh->None,PlotStyle->Directive[Green,Opacity[0.3],Specularity[White,40]],PerformanceGoal->"Quality",PlotPoints120,BoundaryStyleDirective[Black,Thick]],ParametricPlot3D[{ν[t],
ν
1
[t]},{t,0,10},PlotStyleBlack]/.Line[pts_,rest___]:>Tube[pts,0.026,rest]},AxesLabel->(Style[#,15,Blue]&/@{"X","Y","Z"}),Axes->True,AxesOrigin{0,0,0},Boxed->False,BoxRatios->Automatic,ViewPoint->{1.3,-2.4,2.},ImageSize700,PlotRange{{0,10},{-5,5},{0,10}}
Out[]=