ÔñéðëÜ ïëïêëçñþìáôá

NB Init. ÌÇÍ ÁÃÃÉÆÅÔÅ !!!

Mathematica 2. Initializations
Book 3

Needs["Graphics`Colors`"]
Off[General::spell];
Off[General::spell1];
Off[ParametricPlot3D::ppcom];
Off[ParametricPlot::ppcom];
Off[Plot::plnr];
SetOptions[Limit,Analytic->True];
Off[Clear::ssym];

CMView = {2.7, 1.6, 1.2};

ThreeAxes[u,v]

ThreeAxes[u_,v_] := 
      Graphics3D[{
            {Blue,Line[{{-u,0,0},{u,0,0}}]},
                    Text["x",{u + v,0,0}],
            {Blue,Line[{{0,-u,0},{0,u,0}}]},
                    Text["y",{0, u + v,0}],
            {Blue,Line[{{0,0,0},{0,0,u}}]},
                    Text["z",{0,0,u + v}]}]; 
ThreeAxes[u_] := ThreeAxes[u,u/8];
ThreeAxes::usage = "ThreeAxes[a,b] makes a standard cartesian axis graphics object
	with x, y, and z running from -a to a, and with axis
	labels b units beyond the tips of the axes.  ThreeAxes[a]
	is ThreeAxes[a,a/8].";

Vector Drawers

BeginPackage["Vectors3D`","Graphics`Colors`"];

(* Last repair, February, 1993,  Bill Davis   *)

Perp::usage = "Perp[expr] returns a single 3-d vector which is 
		perpendicular to the vector, 'expr'.";
		
Cross::usage = "Cross[a,b] returns the cross product of a and b.";

Norm::usage = "Norm[expr] returns the Euclidean length of the
		vector, 'expr'.";
		
Turn::usage = "Turn[expr] returns an orthogonal 3 by 3 matrix 
		which takes {0,0,1} to 'expr'.";
		
Vector::notnum = "The vectors `1` and `2` aren't numeric.";
		
Arrow::usage = "Arrow[expr] produces a 2 or 3 dimensional 
	vector from the origin to expr. Arrow[expr,Tail->tail] gives a 
	vector from tail to expr.  RGBColor and Thickness may be added as 
	options";
	
Vector::usage = "Vector[a,b] gives a vector from a to b. 
			RGBColor and Thickness may be added as options";

Vector2D::usage = "Vector2D[a,b] produces a 2-d graphics object, a 
  		vector from a to b. Vector[a,b,options] produces the 
  		graphics object with the styles such as RGBColor and 
  		Thickness included.";
		
Vector3D::usage = "Vector3D[a,b] produces the 3-d graphics object,
		a vector from a to b.  Vector3D[a,b,options] produces the 
		graphics object with the options such as RGBColor and 
		Thickness included.";
		
Begin["`Private`"];

SetOptions[ParametricPlot,AspectRatio->Automatic];
SetOptions[Plot,AspectRatio->Automatic];
SetOptions[Graphics,AspectRatio->Automatic];

colorQ[a_]:=(Head[a]===RGBColor);
grayQ[a_]:=(Head[a]===GrayLevel);
thickQ[a_]:=(Head[a]===Thickness);
notruleQ[u_] := !(Head[u]===Rule);
ruleQ[u_] := (Head[u]===Rule);

Perp[a_?VectorQ] := If[a[[1]]==0,Return[{1,0,0}],
			 If[a[[2]]==0,Return[{0,1,0}],
			  Return[{-a[[2]],a[[1]],0}]]];
			
				
Norm[a_?VectorQ]:=Sqrt[a.a];

Turn[a_?VectorQ] := Module[{A = a, B},
		B = Perp[A]/Norm[Perp[A]];
		A = a/Norm[a];
		Return[Transpose[{B,A~Cross~B,A}]]
					];
					
corn = Table[N[1/4 {1.3Cos[2 k Pi/7],1.3Sin[2 k Pi/7],-4}],
	{k,0,9}]; 

corners[a_] := N[Table[(a.corn[[j]]),{j,1,10}]]; 

head[A_,a_,b_,hr___] := 
	Module[{HR},
	HR=ToExpression["HeadRatio"]/.Select[hr,ruleQ]/.
		(ToExpression["HeadRatio"]->.17);
	Return[Flatten[
		Table[{b,b+HR Norm[b-a] corners[A][[j-1]],
			b + HR Norm[b-a] corners[A][[j]]},{j,2,9}],1]]];
				

Vector3D[a_,b_,opts___] := 
	Module[{A,B,AA,BB,m,aa,bb,bob,boob,booob},
		aa=N[a];bb=N[b];
		If[!NumberQ[aa.aa+bb.bb],
			Message[Vector::notnum,a,b];Return[Null],					
			If[Chop[aa-bb]=={0.,0.,0.},
				B={bb,bb},
				A = Turn[bb-aa];
				B = head[A,aa,bb,{opts}]];
		bob=Join[Select[{opts},colorQ],Select[{opts},thickQ],
			Select[{opts},grayQ]];
		boob=Join[Select[{opts},colorQ],Select[{opts},grayQ]];
		booob=If[Length[boob]==0,RGBColor[0,0,0],
				boob[[1]]];			
		Return[Graphics3D[{Join[bob,{Line[{aa,bb}]}],
			{SurfaceColor[booob],Polygon[B]}}]]
					]];
									
									

Vector2D[from_,to_,opts___] := Module[{A,B,AA,BB,c,d,a},
			AA=N[to];
			BB=N[from];
			
			If[!NumberQ[AA.AA+BB.BB],
				Message[Vector::notnum,from,to];
				Return[Null],
			A=AA-BB;
			
			a = Join[Select[{opts},colorQ],Select[{opts},
				thickQ],Select[{opts},grayQ]];
			
			If[Chop[AA-BB]=={0.,0.},
				Return[Graphics[{Join[a,{Line[{BB,AA}]}]}]],
			
			B = RotateLeft[A] {-1,1};
			e = ToExpression["HeadRatio"]/.
					Select[{opts},ruleQ]/.
				{ToExpression["HeadRatio"]->.21};
			c = AA - e A + e/3 B;
			d = AA - e A - e/3 B;
						
			Return[
			Graphics[{Join[a,{Line[{BB,AA}]}],
					Join[a,{Polygon[{AA,c,d}]}]}]]]]];
						
					
Vector[from_,to_,opts___]:=
	If[Length[N[from]]==2,Return[Vector2D[from,to,opts]],
		Return[Vector3D[from,to,s]]];
		
Arrow[vec_,opts___]:= 
	Module[{a,b},
				If[Length[N[vec]]==2,
				a=ToExpression["Tail"]/.Select[{opts},ruleQ]/.
					{ToExpression["Tail"]->{0,0}};
				Return[Vector2D[a,a+vec,opts]],
				a=ToExpression["Tail"]/.Select[{opts},ruleQ]/.
					{ToExpression["Tail"]->{0,0,0}};
				Return[Vector3D[a,a+vec,opts]]]
				];
	
	
(*Protect[Perp,Cross,Turn,Vector3D,Vector2D,Vector,Arrow];*)

End[];


EndPackage[];
Null

Spherical Setup

Clear[polar]
polar[{x_,y_,z_}]:=
{Sqrt[x^2+y^2+z^2],ArcTan[y/x],ArcCos[z/Sqrt[x^2+y^2+z^2]]};
Clear[X]
length[X_]:=Sqrt[X.X];
rr=3;
R=5;
scale=0.8;
red=RGBColor[1,0,0];
blue=RGBColor[0,0,1];
color1=RGBColor[0.8,0.2,0.8];
color2=RGBColor[0.2,0.8,0.2];
black=RGBColor[0,0,0];
Clear[axes]
axes:=Graphics3D[{
Line[{X,{1.5rr,0,0}}],
Line[{X,{0,1.5rr,0}}],
Line[{X,{0,0,2rr}}]}
];
Clear[xyplane]
xyplane:=Graphics3D[Polygon[
{
{-2rr,-2rr,0},
{-2rr, 2rr,0},
{ 2rr, 2rr,0},
{ 2rr,-2rr,0}
}]];
X={0,0,0};
origin=Graphics3D[{black,Point[X]}];
Z={2,3,5};
Clear[point]
point:=Graphics3D[{red,PointSize[0.03],Point[scale Z]}];
Z0={0,0,length[0.8 Z]};
Clear[zpoint]
zpoint=Graphics3D[{black,Point[Z0]}];
XY=scale Z {1,1,0};
Clear[xypoint]
xypoint:=Graphics3D[{black,Point[XY]}];
X0={length[scale Z {1,1,0}],0,0};
xpoint=Graphics3D[{black,Point[X0]}];
Y0={0,length[scale Z {1,1,0}],0};
ypoint=Graphics3D[{black,Point[Y0]}];
Clear[ray]
ray:=Graphics3D[{blue,Line[{X,scale Z}]}];
Clear[ray1]
ray1:=Graphics3D[{color1,Line[{scale Z,Z}]}];
Clear[xyray]
xyray:=Graphics3D[{red,Line[{X, Z {1,1,0}}]}];
rays:={ray,xyray};
xyaux=
length[0.8 Z](0.8 Z {1,1,0})/length[0.8 Z {1,1,0}];
Clear[zarc]
zarc=Graphics3D[{blue,Line[
Table[Cos[t]{0,0,length[0.8 Z]} + Sin[t]xyaux,
{t,0,polar[0.8 Z][[3]],.1}]
]}];
Clear[zarc1]
zarc1=Graphics3D[{color1,Line[
Table[Cos[t]{0,0,length[0.8 Z]} + Sin[t]xyaux,
{t,polar[0.8 Z][[3]],1,.1}]
]}];
Clear[zarc2]
zarc2=Graphics3D[{color1,Line[
Table[Cos[t]{0,0,length[0.8 Z]} + Sin[t]xyaux,
{t,-.8,0,.1}]
]}];
Clear[xyarc]
xyarc=Graphics3D[{blue,Line[
Table[Cos[t] X0 + Sin[t] Y0,
{t,0,polar[XY][[2]],.1}]
]}];
Clear[xyarc1]
xyarc1=Graphics3D[{color1,Line[
Table[Cos[t] X0 + Sin[t] Y0,
{t,-.5,0,.1}]
]}];
Clear[xyarc2]
xyarc2=Graphics3D[{color1,Line[
Table[Cos[t] X0 + Sin[t] Y0,
{t,polar[XY][[2]],polar[XY][[2]]+0.3,.1}]
]}];
Clear[stick]
stick=Graphics3D[{red,Line[{scale Z,XY}]}];

xaxislabel=Text["x",{0.363,0.26}];
yaxislabel = Text["y",{0.73, 0.30}];
zaxislabel = Text["z",{0.50, 0.67}];

stickplot =
Show[xyplane,axes,ray,point,Boxed->False,
ViewPoint->{2.747, 1.594, 1.168},
Epilog->{xaxislabel,yaxislabel,zaxislabel},
DisplayFunction->Identity];

rho=Text["r",{0.53, 0.45}];
rlabel = Show[xyplane,axes,ray,point,
Boxed->False,ViewPoint->{2.747, 1.594, 1.168},
Epilog->{rho,xaxislabel,yaxislabel,zaxislabel},
DisplayFunction->Identity];

phi=Text["s",{0.55, 0.57}];
slabel =
Show[xyplane,zpoint,axes,ray,point,zarc,zarc1,zarc2,Boxed->False,ViewPoint->{2.747, 1.594, 1.168},
Epilog->{phi,xaxislabel,yaxislabel,zaxislabel},DisplayFunction->Identity];

theta=Text["t",{0.5, 0.3}];

tlabel=Show[xyplane,axes,xypoint,xpoint,ray,stick,point,xyarc,xyarc1,xyarc2,Boxed->False,ViewPoint->{2.747, 1.594, 1.168},
Epilog->{theta,xaxislabel,yaxislabel,zaxislabel},
DisplayFunction->Identity];

all =
Show[xyplane,axes,ray,stick,point,zpoint,xpoint,xypoint,xyarc,xyarc1,xyarc2,zarc,zarc1,zarc2,Boxed->False,ViewPoint->{2.747, 1.594, 1.168},
Epilog->{theta,phi,rho,xaxislabel,yaxislabel,zaxislabel},DisplayFunction->Identity];

Crisscross

Clear[x]
crisscross = Plot[{Cos[x],Sin[x]},{x,Pi/8,3 Pi/8},
    PlotStyle->{{RGBColor[1,0,0],Thickness[0.02]}},
    Axes->None,AspectRatio->1/3,
    DisplayFunction->Identity];

Local Intializations

Clear[grad, expr, vars, point];
\(\(grad[expr_, \ vars_List, 
          point_List]\  := \ \(D\ [
expr, \ #] &\)\  /@ \ vars\ \  /. \ 
          Thread[vars\  -> point];\)

\)
grad[expr_, vars_] := grad[expr, vars, vars]
grad[expr_] := grad[expr, {x, y, z}, {x, y, z}]

\(
Clear[Curl]\)
\(Curl[expr_List, 
        vars_List] := \ 
{D[expr[\([3]\)], vars[\([2]\)]]\  - \ 
            D[expr[\([2]\)], vars[\([3]\)]], 
\ \(-\ 
              D[expr[\([3]\)], vars[\([1]\)]]\)\  + \ 
            D[expr[\([1]\)], vars[\([3]\)]], 
\ 
          D[expr[\([2]\)], vars[\([1]\)]]\  - \ 
            D[expr[\([1]\)], vars[\([2]\)]]} /; 
\ \(\(Length[expr] == 
          Length[vars]\)\(
\)\(
\)\)\)
Curl[expr_List] := Curl[expr, {x, y, z}];

\(
Clear[Div, expr, vars]\)
\(Div[expr_List, vars_List] := 
      Inner[D, expr, vars]\  /; \ 
\(\(Length[expr] == 
          Length[vars]\)\(
\)\(
\)\)\)
Div[expr_List] := Div[expr, {x, y, z}]

Clear[jacobian, jacobiandet]
\(jacobian[expr_List, vars_List] := 
      Outer[D, \ \ expr\ , 
          vars] /; 
\(\(Length[expr]\)\(==\)\(Length[
            vars]\)\(
\)\)\)
\(jacobiandet[expr_List, vars_List] := 

      Det[jacobian[expr, vars]] /; 
\(\(Length[expr] == 
          Length[vars]\)\(
\)\(\ \)\)\)

Clear[hessian, funs, vars]
\(\(hessian[funs_, \ vars_List, point_List]\  := \ 
        Outer[
D[funs, \ #1, \ #2] &, \ vars, \ vars] /. \ 

          Thread[vars -> point]\ ;\)\)
hessian[funs_,  vars_List] :=  
        hessian[funs, vars, vars];
hessian[funs_] := hessian[funs, {x, y}, {x, y}];
\(hessiandet[funs_, \ vars_List, point_List] := 
      Det[hessian[
funs, vars, point]]\)
\(hessiandet[funs_, \ vars_List] := 
      Det[hessian[
funs, vars, vars]]\)
hessiandet[funs_] := hessiandet[funs, {x, y}, {x, y}]

Clear[box]
\(box[
{x1_, x2_}, {y1_, y2_}, {z1_, z2_}, 

        color_:  RGBColor[1, 0, 0]] := 
      Module\ [
{side1, side2, side3, side4}, 


        side1 = Graphics3D[{color, Thickness[ .02], 

              Line[{{\ x1, y1, z1}, {x1, y2, z1}, 
\ \ \ \ \ \ {x1, y2, 
                    z2}, {x1, y1, z2}, {\ x1, y1, z1}}]}]; 


        side2 = Graphics3D[{color, Thickness[ .02], 

              Line[{{\ x1, y1, z1}, {x2, y1, z1}, 
\ \ \ \ \ \ {x2, y1, 
                    z2}, {x1, y1, z2}, {\ x1, y1, z1}}]}]; 


        side3 = Graphics3D[{color, Thickness[ .02], 

              Line[{{\ x2, y1, z1}, {x2, y2, z1}, 
\ \ \ \ \ \ {x2, y2, 
                    z2}, {x2, y1, z2}, {\ x2, y1, z1}}]}]; 
\ \ \ \ \ \ 

        side4 = Graphics3D[{color, Thickness[ .02], 

              Line[{{\ x1, y2, z1}, {x2, y2, z1}, 
\ \ \ \ \ \ {x2, y2, 
                    z2}, {x1, y2, z2}, {\ x1, y2, z1}}]}]; 
\ \ \ \ \ \ 
\ 
        Return[{side1, side2, side3, side4}]\ \ \ \ \ 

]


\)
Unprotect[SA];
Clear[SA, expr, vars];
\(\(SA[expr_List, vars_List] := 

        Norm[\ Cross[\ 
              D[expr, vars[\([1]\)]\ ]\ \ , 
\ \ \ \ \ \ \ \ \ \ \ \ \ D[\ 
                expr, vars[\([2]\)]]\ ]] /; 
\((Length[expr] == 3 && 
              Length[vars] == 2)\);\)
\)
\(\(Protect[SA];\)

\)
Unprotect[NV];
Clear[NV, expr, vars];
\(\(NV[expr_List, vars_List] := 

        Cross[\ D[expr, vars[\([1]\)]\ ]\ \ , 
\ \ \ \ \ \ \ \ \ \ \ \ \ D[\ 
              expr, vars[\([2]\)]]\ ] /; 
\((Length[expr] == 3 && 
              Length[vars] == 2)\);\)
\)
Protect[NV];

\(\(\( (*\(:\)\(Version : \ Mathematica\ 2.0\)\ *) \)\(
\)\)\)

◆Ðïõ ÷ñçóéìïðïéïýíôáé ôá ôñéðëÜ
      ïëïêëçñþìáôá ;

Áêñéâþò üðùò ôá äéðëÜ ïëïêëçñþìáôá ÷ñçóéìïðïéÞèçêáí ãéá ôïí õðïëïãéóìü åìâáäþí, ôá ôñéðëÜ ïëïêëçñþìáôá ìðïñïýí íá ÷ñçóéìïðïéçèïýí ãéá ôïí õðïëïãéóìü ôïõ 
üãêïõ åíüò 3-äéÜóôáôïõ ÷ùñßïõ
×ñçóéìåýïõí åðßóçò ãéá ôçí åýñåóç ôçò ìÝóçò ôéìÞò ìéáò óõíÜñôçóçò ïñéóìÝíçò óå Ýíá 3-äéÜóôáôï ÷ùñßï êáèþò åðßóçò ôçò ìÜæáò, ôùí óõíôåôáãìÝíùí ôïõ êÝíôñïõ ìÜæáò êáé ôçò ñïðÞò áäñáíåßáò 3-äéÜóôáôùí óùìÜôùí.

◆Ïñéóìüò ôñéðëïý ïëïêëçñþìáôïò

¼ðùò óôá äéðëÜ ïëïêëçñþìáôá, îåêéíÞóáìå ìå ïëïêëÞñùóç ïñèïãùíßùí ðáñáëëçëïãñÜììùí, óôá ôñéðëÜ ïëïêëçñþìáôá ôï áðëïýóôåñï ÷ùñßï åðß ôïõ ïðïßïõ áðïäåéêíýåôáé ç ýðáñîç ïëïêëçñþìáôïò ãéá ìéá óõíå÷Þ óõíÜñôçóç f åßíáé ôï ïñèïãþíéï ðáñáëëçëåðßðåäï, äçëáäÞ, Ýíá óôåñåü ìå  ðëåõñÝò  ðáñÜëëçëåò óôïõò Üîïíåò x, y êáé z.  ÓõãêåêñéììÝíá, 
                     R = [(Subscript[a, 1]), (Subscript[b, 1])]×[(Subscript[a, 2]), 
            (Subscript[b, 2])]×[(Subscript[a, 3]), (Subscript[b, 3])] 
êáé ï üãêïò ôïõ R  åßíáé 
                    Vol(
          R) = ((Subscript[b, 1]) - (Subscript[a, 1])) ((Subscript[b, 2]) - (Subscript[a, 2])) ((Subscript[b, 3]) - (Subscript[a, 3]))
                    
Ìßá äéáìÝñéóç P ôïõ R áðáñôßæåôáé áðü äéáìåñßóåéò (Subscript[P, 1]),(Subscript[P, 2]),(Subscript[P, 3]) ôùí äéáóôçìÜôùí [(Subscript[a, 1]), (Subscript[b, 1])], [(Subscript[a, 2]), (Subscript[b, 2])], [(Subscript[a, 3]), 
        (Subscript[b, 3])] áíôßóôïé÷á êáé ÷ùñßæåé ôï ïñèïãþíéï R óå ìéêñÜ õðï-ïñèïãþíéá ôá ïðïßá ôá óõìâïëßæïõìå ìå ôï ãñÜììá S. ¼ðùò áêñéâþò êÜíáìå óôá äéðëÜ ïëïêëçñþìáôá, èåùñïýìå Üíù êáé êÜôù áèñïßóìáôá ùò åîÞò:
                L(P, f) = (∑_S)(inf f) Vol(S)  
êáé
               U(P, f) = (∑_S)(sup f) Vol(S).

Ìßá óõíÜñôçóç f èá ëÝãåôáé ïëïêëçñþóéìç åðß ôïõ ïñèïãùíßïõ  R  åÜí õðÜñ÷åé ìïíáäéêüò ðñáãìáôéêüò áñéèìüò  Á  ôÝôïéïò þóôå ãéá êÜèå äýï äéáìåñßóåéò  P, Q ôïõ  R  íá éó÷ýåé 
                  L(P, f) ≤ A ≤ U(Q, f)
                  
Èåþñçìá: ¸óôù R ïñèïãþíéï êáé  f : R->R öñáãìÝíç óõíÜñôçóç ç ïðïßá åßíáé óõíå÷Þò óôï R  åêôüò ßóùò áðü óçìåßá ôá ïðïßá áíÞêïõí óå ðåðåñáóìÝíåò ôï ðëÞèïò äéáöïñßóéìåò åðéöÜíåéåò óôï R. Ôüôå ç f åßíáé ïëïêëçñþóéìç.
Óõìâïëéóìüò: ∫∫(∫_R) 
              f =\(TraditionalForm\`∫\(∫\(∫\_R\ \(f(x, y, 
                z)\) dx dy dz\)\)\).

Åðåêôåßíïõìå ôçí ïêëïêëÞñùóç óå öñáãìÝíá ÷ùñßá  Ó  ãåíéêüôåñá ôïõ ïñèïãùíßïõ ðáñáëëçëåðéðÝäïõ  R  ùò åîÞò: èåùñïýìå ðáñáëëçëåðßðåäï  R ðïõ íá ðåñéÝ÷åé ôï Ó êáé èÝôïõìå
             f(x) = 0 ãéá   x ϵ R \ Ó êáé \(TraditionalForm\`∫\(∫\(\(∫\_Ó\)\(\ \)\(f\)\(\ \)\)\)\)=\(TraditionalForm\`∫\(∫\(\(∫\_R\)\(\ \)\(f\)\(\ \)\)\)\)

◆ÏëïêëÞñùóç åðß ïñèïãùíßïõ 
      ðáñáëëçëåðéðÝäïõ (õðïëïãéóìüò)

Ãéá ôïí õðïëïãéóìü ôùí ôñéðëþí ïëïêëçñùìÜôùí ÷ñçóéìïðïéïýìå ôï Èåþñçìá Fubini äéáôõðùìÝíï óå 3 äéáóôÜóåéò. Ìå áõôü ôïí ôñüðï ï õðïëïãéóìüò åíüò ôñéðëïý ïëïêëçñþìáôïò áíÜãåôáé óå äéáäï÷éêÝò ïëïêëçñþóåéò Riemann ðñáãìáôéêþí óõíáñôÞóåùí.

ÐáñÜäåéãìá

Õðïëïãéóìüò ôïõ ôñéðëïý ïëïêëçñþìáôïò
 \(TraditionalForm\`∫\(∫∫\_\(\(\ \)\(B\)\)\)\) (x^2) (y^3) z dx dy dz   
 üðïõ B ðáñáëëçëåðßðåäï ìå  0 ≤ x ≤ 2; -1 ≤ y ≤ 3 ; 4 ≤ z ≤ 7.

×ñçóéìïðïéþíôáò ôï Mathematica Ý÷ïõìå ôçí ðáñáêÜôù ôñéóäéÜóôáôç ãñáöéêÞ ðáñÜóôáóç ôïõ ðáñáëëçëåðéðÝäïõ.

Clear[setup1]
setup1 = 
        Show[box[{0, 2}, {-1, 3}, {4, 7}], Axes -> True, 
          AxesLabel -> {"x", "y", "z"}];

×ñçóéìïðïéþíôáò êáé ðÜëé ôï Mathematica  ìðïñïýìå íá õðïëïãßóïõìå ôï ïëïêëÞñùìá áìÝóùò. Ôï ìüíï ðïõ ÷ñåéÜæåôáé íá êÜíïõìå åßíáé íá êáèïñßóïõìå ôá üñéá ôçò ïëïêëÞñùóçò (ôá ïðïßá ìáò äßíïíôáé).

Clear[x, y, z]
(∫_4^7)(∫_(-1)^3)(∫_0^2) (x^2) (y^3) z \[DifferentialD]x \[DifferentialD]y \[DifferentialD]z
Integrate[
      (x^2) (y^3) z,  {z,  4,  7},  {y,  -1,  3},  {x,  0,  
        2}]

Ï õðïëïãéóìüò ãßíåôáé óå ôñßá âÞìáôá ÷ñçóéìïðïéþíôáò ôï Èåþñçìá Fubini. ÊÜèå âÞìá åßíáé Ýíá áðëü ïëïêëÞñùìá. 
Ðñþôï âÞìá:

Clear[first, second, third, x, y, z, f]
f[x_, y_, z_] = (x^2) (y^3) z;
first[y_, 
        z_] = (∫_0^2) f[x, y, z] \[DifferentialD]x
\(		Integrate[f[x, \ y, \ z], \ {x, \ 0, \ 2}]\)

Äåýôåñï âÞìá,

second[
        z_] = (∫_(-1)^3) first[y, z] \[DifferentialD]y
\(		Integrate[first[y, \ z], \ {y, \ \(-1\), \ 3}]\)

Êáé ôï ôñßôï âÞìá:

third = (∫_4^7) second[z] \[DifferentialD]z
\(		Integrate[second[z], \ {z, \ 4, \ 7}]\)

Èá åîåôÜóïõìå êÜèå Ýíá áðü ôá ôñßá âÞìáôá ÷ùñéóôÜ.
Ôé âñÞêáìå üôáí õðïëïãßóáìå ôï ðñþôï ïëïêëÞñùìá, äçëáäÞ, üôáí õðïëïãßóáìå ôçí óõíÜñôçóç first[y,z] ?

Clear[setup2]
setup2[y0_, z0_] := 
        Graphics3D[{Blue, Thickness[ .02], 
            Line[{{0, y0, z0}, {2, y0, z0}}]}];
firstintegral = Show[setup1, setup2[1, 5]];

AðÜíôçóç: Õðïëïãßóáìå ôï ïëïêëÞñùìá ôçò óõíÜñôçóçò 
                   (x^2) (y^3) z 
åðß  åõèõãñÜììïõ ôìÞìáôïò ðïõ õðïäåéêíýåôáé óôçí ðáñáðÜíù ãñáöéêÞ ðáñÜóôáóç.
Åäþ ôï y êáé ôï z åßíáé óôáèåñÜ, äçëáäÞ èåùñïýìå ôçí
 (x^2) (y^3) z  ùò óõíÜñôçóç ìüíï ôïõ x.

Ôé êÜíáìå üôáí õðïëïãßóáìå ôï äåýôåñï ïëïêëÞñùìá, äçëáäÞ, üôáí õðïëïãßóáìå ôçí óõíÜñôçóç second[z] ?

Clear[setup2, strip]
step =  .4;
setup2[y0_, z0_] := 
        Graphics3D[{Blue, Thickness[ .02], 
            Line[{{0, y0, z0}, {2, y0, z0}}]}];
strip[z0_] := 
        Table[setup2[y0, z0], {y0, -1, 3, step}];
Show[setup1, strip[5]];

AðÜíôçóç: Ïëïêëçñþóáìå åðß ìßáò ëùñßäáò ðïõ õðïäåéêíýåôáé óôçí ðáñáðÜíù ãñáöéêÞ ðáñÜóôáóç. 
Åäþ êñáôÞóáìå ôï z óôáèåñü êáé ç ìáôáâëçôÞ y êéíåßôï áðü -1 Ýùò 3.

Óôçí ðáñáðÜíù ãñáöéêÞ ðáñÜóôáóç ÷ñçóéìïðïéÞóáìå äéÜöïñåò äéáêñéôÝò ôéìÝò ãéÜ  ôï y. ÅÜí ðÜñïõìå ìéêñüôåñï "âÞìá" ãéá ôéò ôéìÝò ôïõ y, ôüôå ç ãñáöéêÞ ðáñÜóôáóç èá öáßíåôáé ùò åðéöÜíåéá.

Clear[setup2, strip]
step =  .1;
setup2[y0_, z0_] := 
        Graphics3D[{Blue, Thickness[ .02], 
            Line[{{0, y0, z0}, {2, y0, z0}}]}];
strip[z0_] := 
        Table[setup2[y0, z0], {y0, -1, 3, step}];
secondintegral = Show[setup1, strip[5]];

Ôé êÜíáìå üôáí õðïëïãßóáìå ôï ôñßôï ïëïêëÞñùìá third ?

Clear[setup3, strips, stepz]
stepz =  .5;
setup3[z0_] := 
        Graphics3D[
          Polygon[{{0, -1, z0}, {2, -1, z0}, {2, 3, z0}, {0, 3, 
                z0}, {0, -1, z0}}]];
strips = Table[setup3[z0], {z0, 4, 7, stepz}];
Show[setup1, strips];

ÁðÜíôçóç: ÁöÞíïíôáò ôçí ìåôáâëçôÞ z íá êéíåßôáé áðü 4 óôï 7, ç ëùñßäá ðïõ åßäáìå ðáñáðÜíù óÜñùóå ïëüêëçñï ôï ðáñáëëçëåëßðåäï.

Èåùñþíôáò ìéêñü "âÞìá" ãéá ôéò ôéìÝò ôïõ z (i.e. stepz=.1) âëÝðïõìå ôé ãßíåôáé ìå ôçí ïëïêëÞñùóç.

Clear[setup3, strips, stepz]
stepz =  .1;
setup3[z0_] := 
        Graphics3D[
          Polygon[{{0, -1, z0}, {2, -1, z0}, {2, 3, z0}, {0, 3, 
                z0}, {0, -1, z0}}]];
strips = Table[setup3[z0], {z0, 4, 7, stepz}];
thirdintegral = Show[setup1, strips];

Ôá ðáñáðÜíù âÞìáôá óå ôáéíßåò ìéêñïý ìÞêïõò.
ÂÞìá 1: y, z óôáèåñÜ

Clear[line, r, step]
step =  .25;
line[x_] := 
      Graphics3D[{Blue, Thickness[ .02], 
          Line[{{0, 1, 5}, {x, 1, 5}}]}]
Table[Show[setup1, line[x]];, {x, 0, 2, step}];

ÂÞìá 2: z óôáèåñü

Clear[setup2, strip, polygon]
step =  .25;
polygon[r_, z0_] := 
        Graphics3D[
          Polygon[{{0, -1, z0}, {2, -1, z0}, {2, -1 + r, 
                z0}, {0, -1 + r, z0}, {0, -1, z0}}]];
Table[Show[setup1, polygon[r, 5]];, {r, 0, 4, 
          step}];

ÂÞìá 3: ôßðïôå óôáèåñü (x, y êáé z åßíáé ìåôáâëçôÝò).

Clear[stepz]
stepz =  .5;
Table[Show[setup1, setup3[z0]];, {z0, 4, 7, 
          stepz}];

◆ÔñéðëÜ ïëïêëçñþìáôá óôåñåþí: 
      êáèïñéóìüò ïñßùí ïëïêëÞñùóçò  
      êáé õðïëïãéóìïß

ÏëïêëÞñùóç åðß óôåñåïý ðïõ öñÜóóåôáé áðü åðéöÜíåéåò 
z= (Subscript[g, 1])[x,y] êáé z = (Subscript[g, 2])[x,y],  
üðïõ  x êáé y áíÞêïõí óå êÜðïéï ÷ùñßï ôïõ xy-åðéðÝäïõ.

ÐáñÜäåéãìá 
Õðïëïãßæïõìå ôï ôñéðëü ïëïêëÞñùìá 
                  \(TraditionalForm\`∫\(∫∫\_\(\(\ \)\(R\)\)\)\)  (x^2) (y^2) z dz dx dy, 
 üðïõ R åßíáé ôï ÷ùñßï ðïõ öñÜóóåôáé áðü ôéò åðéöÜíåéåò
 (Subscript[g, 1])[x,y]= 6 -((x^2)+ (y^2)) êáé (Subscript[g, 2])[x,y] = ((x^2)+ (y^2))-2.

Åßíáé ðÜíôá ÷ñÞóéìï íá Ý÷ïõìå óôï ìõáëü ìáò ìßá åéêüíá ôïõ óôåñåïý åðß ôïõ ïðïßïõ ïëïêëçñþíïõìå. 
Ôï Mathematica åßíáé ÷ñÞóéìï åñãáëåßï ãéá ãñáöéêÝò ðáñáóôÜóåéò.

Óôçí ðáñáêÜôù ãñáöéêÞ ðáñÜóôáóç ðÞñáìå Ýíáí êýêëï áêôßíáò 1,5 (ç áêôßíá äéáëÝ÷èçêå ôõ÷áßá, áëëÜ ç åðéëïãÞ ó÷Þìáôïò êýêëïõ äåí åßíáé ôõ÷áßá) êáé æçôÜìå áðü ôï  Mathematica íá ó÷åäéÜóåé ôá ôìÞìáôá ôùí åðéöáíåéþí  (Subscript[g, 1])[x,y] êáé (Subscript[g, 2])[x,y] ðïõ âñßóêïíôáé ðÜíù áðü áõôü ôï êýêëï.

Clear[x, y, z, g1, g2, r, t]
g1[x_, y_] = 6 - ((x^2) + (y^2));
g2[x_, y_] = ((x^2) + (y^2)) - 2;
x[r_, t_] = r Cos[t];
y[r_, t_] = r Sin[t];
ParametricPlot3D[
        Evaluate[{{x[r, t], y[r, t], g1[x[r, t], y[r, t]]}, {x[r, t], 
              y[r, t], g2[x[r, t], y[r, t]]}}], {t, 0, 2 π}, {r, 0, 
          1.5}];

ÊÜíïõìå ôçí ßäéá ãñáöéêÞ ðáñÜóôáóç åðéëÝãïíôáò ôçí áêôßíá  0≤r≤2 (áõôÞ ç åðéëïãÞ áêôßíáò äåí åßíáé ôõ÷áßá äéüôé 
6 -((x^2)+ (y^2)) = ((x^2)+ (y^2))-2 ⟹ 8 = 2 ((x^2)+ (y^2))  ⟹  4 = (x^2)+ (y^2) 
ðïõ åßíáé êýêëïò áêôßíáò 2).

Clear[x, y, z, g1, g2, r, t]
g1[x_, y_] = 6 - ((x^2) + (y^2));
g2[x_, y_] = ((x^2) + (y^2)) - 2;
x[r_, t_] = r Cos[t];
y[r_, t_] = r Sin[t];
ParametricPlot3D[{ {x[r, t], y[r, t], 
            g1[x[r, t], y[r, t]]}, {x[r, t], y[r, t], 
            g2[x[r, t], y[r, t]]}}, {t, 0, 2 π}, {r, 0, 
          2}];

Ïé äýï åðéöÜíåéò ôÝìíïíôáé üôáí r = 2.  Áõôü óçìáßíåé üôé ïé óõíôåôáãìÝíåò x êáé y (ôùí óçìåßùí ôïõ óôåñåïý) ðñÝðåé íá âñßóêïíôáé ìÝóá óôïí êýêëï êÝíôñïõ 0 êáé áêôßíáò 2 óôï xy-åðßðåäï.

Èåùñïýìå ôá x êáé y óôáèåñÜ.  Ôüôå ôï z äéáôñÝ÷åé ôï äéÜóôçìá áðü ôï (Subscript[g, 1])[x,y] óôï (Subscript[g, 2])[x,y].  Ïëïêëçñþíïíôáò ùò ðñïò z, ôï ôñéðëü ïëïêëÞñùìá áíÜãåôáé óå Ýíá äéðëü ïëïêëÞñùìá åðß ôïõ êýêëïõ áêôßíáò 2 óôï xy-åðßðåäï.

Ï õðïëïãéóìüò ôïõ ïëïêëçñþìáôïò (ìå ÷ñÞóç ôïõ
Mathematica ãéá ôïõò õðïëïãéóìïýò). 

Ðñþôá ïëïêëçñþíïõìå ùò ðñïò z ìå üñéá ïëïêëÞñùóçò  (Subscript[g, 1]) [x, y] and (Subscript[g, 2]) [x, y].

Clear[first, second, third, x, y, z, f]
f[x_, y_, z_] = (x^2) (y^2) z;
first[y_, x_] = 
      Expand[(∫_(g2[x, y])^(g1[x, y]))f[x, y, 
              z] \[DifferentialD]z]

Óôçí óõíÝ÷åéá ïëïêëçñþíïõìå óôïí êýêëï óôï xy-åðßðåäï.  ÊñáôÜìå ôï x óôáèåñü êáé ïëïêëçñþíïõìå ùò ðñïò y ìå üñéá ïëïêëÞñùóçò -(Sqrt[4 - (x^2)]) êáé (Sqrt[4 - (x^2)]).

second[
        x_] = (∫_(-(Sqrt[4 - (x^2)]))^((Sqrt[4 - (x^2)])))first[y, 
            x] \[DifferentialD]y

ÔÝëïò ïëïêëçñþíïõìå ùò ðñïò x.

third = (∫_(-2)^2) second[x] \[DifferentialD]x

ÅëÝã÷ïõìå åáí ç åöáñìïãÞ ìáò äïõëåýåé óùóôÜ:

(∫_(-2)^2)((∫_(-(Sqrt[4 - (x^2)]))^((Sqrt[4 - (x^2)])))((∫_(g2[x, y])^(g1[x, y]))f[x, y, 
                      z] \[DifferentialD]z) \[DifferentialD]y) \[DifferentialD]x

ÏëïêëÞñùóç åðß óôåñåïý ðïõ öñÜóóåôáé áðü åðéöÜíåéåò 
y = (Subscript[g, 1])[x,z] êáé y = (Subscript[g, 2])[x,z],  
üðïõ  x êáé z áíÞêïõí óå êÜðïéï ÷ùñßï ôïõ xz-åðéðÝäïõ.

Õðïëïãßæïõìå ôï ôñéðëü ïëïêëÞñùìá 
                   ∫∫(∫_R) (x^4) y (z^8) dz dx dy
 üðïõ R åßíáé ôï ÷ùñßï ðïõ öñÜóóåôáé áðü ôéò åðéöÜíåéåò
y= ((x^2)+(z^2)) êáé y = 18 - ((x^2)+(z^2)).

Âñßóêïõìå ôçí ôïìÞ ôùí äýï åðéöáíåéþí üðùò ðñïçãïõìÝíùò, äçëáäÞ,  ëýíïíôáò ôçí åîßóùóç
          ((x^2)+(z^2)) = 18 - ((x^2)+(z^2)) 
 ç ïðïßá ìáò äßíåé (x^2)+(z^2) = 9.  
Áõôü óçìáßíåé üôé ïé óõíôåôáãìÝíåò x êáé z (ôùí óçìåßùí ôïõ óôåñåïý) ðñÝðåé íá âñßóêïíôáé ìÝóá óôïí êýêëï êÝíôñïõ 0 êáé áêôßíáò 3 óôï xz-åðßðåäï.

ÊÜíïõìå ôçí ãñáöéêÞ ðáñÜóôáóç ìå åðéëïãÞ áêôßíáò (óôéò ðáñáìåôñéêÝò åîéóþóåéò) r=2.

Clear[x, y, z, g1, g2, r, t]
g1[x_, z_] = (x^2) + (z^2);
g2[x_, z_] = 18 - ((x^2) + (z^2));
x[r_, t_] = r Cos[t];
z[r_, t_] = r Sin[t];
ParametricPlot3D[
        Evaluate[{{x[r, t], g1[x[r, t], z[r, t]], z[r, t]}, {x[r, t], 
              g2[x[r, t], z[r, t]], z[r, t]}}], {t, 0, 2 π}, {r, 0, 2}, 
        AxesLabel -> {"x", "y", "z"}];

ÊÜíïõìå ôþñá ôçí ãñáöéêÞ ðáñÜóôáóç ìå áêôßíá óôéò ðáñáìåôñéêÝò åîéóþóåéò r=3 ðïõ âñÞêáìå üôé õëïðïéåß ôçí ôïìÞ ôùí äýï åðéöáíåéþí.

Clear[x, y, z, g1, g2, r, t]
g1[x_, z_] = (x^2) + (z^2);
g2[x_, z_] = 18 - ((x^2) + (z^2));
x[r_, t_] = r Cos[t];
z[r_, t_] = r Sin[t];
ParametricPlot3D[
        Evaluate[{{x[r, t], g1[x[r, t], z[r, t]], z[r, t]}, {x[r, t], 
              g2[x[r, t], z[r, t]], z[r, t]}}], {t, 0, 2 π}, {r, 0, 3}, 
        AxesLabel -> {"x", "y", "z"}];

Ï õðïëïãéóìüò ôïõ ïëïêëçñþìáôïò (ìå ÷ñÞóç ôïõ
Mathematica). 

Ðñþôá ïëïêëçñþíïõìå ùò ðñïò y ìå üñéá ïëïêëÞñùóçò  (Subscript[g, 1]) [x, z] êáé (Subscript[g, 2]) [x, z].

Clear[first, second, third, x, y, z, f]
f[x_, y_, z_] = (x^4) y (z^8);
first[x_, z_] = 
      Expand[(∫_(g1[x, z])^(g2[x, z]))f[x, y, 
              z] \[DifferentialD]y]

Óôçí óõíÝ÷åéá äïõëåýïõìå óôï xz-åðßðåäï, êñáôþíôáò ôï x óôáèåñü êáé ïëïêëçñþíïíôáò ùò ðñïò z.

second[
        x_] = (∫_(-(Sqrt[9 - (x^2)]))^((Sqrt[9 - (x^2)])))first[x, 
            z] \[DifferentialD]z

ÔÝëïò ïëïêëçñþíïõìå ùò ðñïò x.

third = (∫_(-3)^3) second[x] \[DifferentialD]x

Ôá âÞìáôá óôçí ðáñáðÜíù ïëïêëÞñùóç Ý÷ïõí êáèáñÜ åðåîçãçìáôéêü ñüëï. ÅÜí èÝëïõìå, ïëüêëçñïò ï ðáñáðÜíù õðïëïãéóìüò ìðïñåß íá ãßíåé ìå ÷ñÞóç ôïõ  Mathematica  áð' åõèåßáò.

(∫_(-3)^3)((∫_(-(Sqrt[9 - (x^2)]))^((Sqrt[9 - (x^2)])))((∫_(g1[x, z])^(g2[x, z]))f[x, y, 
                      z] \[DifferentialD]y) \[DifferentialD]z) \[DifferentialD]x

ÏëïêëÞñùóç åðß óôåñåïý ðïõ ðåñéãñÜöåôáé áðü áíéóüôçôåò 
a ≤ x ≤ b, (Subscript[g, 1])[x] ≤ y ≤ (Subscript[g, 2])[x] êáé
(Subscript[h, 1])[x,y] ≤ z ≤ (Subscript[h, 2])[x,y]

Ç ïëïêëÞñùóç åðß åíüò óôåñåïý Á ôï ïðïßï áðïôåëåßôáé áðü üëá ôá (x,y,z)  ãéá ôá ïðïßá éó÷ýïõí ïé ðáñáðÜíù áíéóüôçôåò ãéá êÜðïéåò äéáöïñßóéìåò óõíáñôÞóåéò 
 
 (Subscript[g, 1]), (Subscript[g, 2]): [a,b]->R ìå (Subscript[g, 1])(x) ≤ (Subscript[g, 2])(x) êáé 
 (Subscript[h, 1]), (Subscript[h, 2]): {(x,y) : a ≤ x ≤ b êáé (Subscript[g, 1])(x) ≤ y ≤ (Subscript[g, 2])(x)}->R ìå 
                     (Subscript[h, 1])(x,y) ≤ (Subscript[h, 2])(x,y)
 
 áíÜãåôáé óôïí õðïëïãéóìü ôïõ ðáñáêÜôù ïëïêëçñþìáôïò
 
        \(TraditionalForm\`∫\(∫∫\_\(\(\ \)\(Á\)\)\)\) f = (∫_a^b)((∫_((Subscript[g, 1])[x])^((Subscript[g, 2])[x]))((∫_((Subscript[h, 1])[x, y])^((Subscript[h, 2])[x, y]))f(
                          x, y, z) dz)dy)dx
 
 ÁõôÞ ç ðåñßðôùóç åìðåñéÝ÷åôáé óôéò ðñïçãïýìåíåò

ÏëïêëÞñùóç åðß óôåñåïý ðïõ ðåñéãñÜöåôáé áðü äýï áíéóüôçôåò: åýñåóç üãêïõ ôïõ óôåñåïý ðïõ êáèïñßæåôáé áðü ôçí ôïìÞ ôùí óöáéñþí (x^2) + (y^2) + (z^2) = (4^2) êáé (((x - 3))^2) + (y^2) + (z^2) = (2^2).

Áíáæçôïýìå ôïí üãêï ôïõ óôåñåïý Á ôï ïðïßï áðïôåëåßôáé áðü üëá ôá (x,y,z)  ãéá ôá ïðïßá éó÷ýïõí ïé áíéóüôçôåò 
   (x^2) + (y^2) + (z^2) ≤ (4^2) êáé (((x - 3))^2) + (y^2) + (z^2) ≤ (2^2)
Õðïëïãßæïõìå ôï ÷ùñßï R ôïõ xy-åðéðÝäïõ ðïõ ðåñéÝ÷åôáé óôåñåü Á:
{({{(x^2) + (y^2) = (4^2)}, {(((x - 3))^2) + (y^2) = (2^2)}})}⟹ ⋯ ⟹ {({{x = 21/6}, {y = ±(Sqrt[(4^2) - (((21/6))^2)])}})}
Óõíåðþò ôï ÷ùñßï R áðïôåëåßôáé áðü ôá ÷ùñßá  R = (Subscript[R, 1]) ⋃ (Subscript[R, 2]) üðïõ 
(Subscript[R, 1]) =  {(x, 
                y) │  1 ≤ x ≤ 
            21/6 ,  -(Sqrt[(2^2) - (((x - 3))^2)]) ≤ 
            y ≤ (Sqrt[(2^2) - (((x - 3))^2)])}
  êáé 
(Subscript[R, 2]) =  {(x, y) │ 21/
                  6 ≤ x ≤ 
              4 ,  -(Sqrt[(4^2) - (x^2)]) ≤ 
              y ≤ (Sqrt[(4^2) - (x^2)])}.

Ç ãñáöéêÞ ðáñÜóôáóç ôçò ïñéæüíôéáò ôïìÞò óôï ýøïò z = 0 (äçë. ôïõ xy-åðéðÝäïõ)

Clear[v, w, s1, s2,  c1,  c2,  t1,  t2,  t3]
s1 = 
        ParametricPlot[{4  Cos[w] ,  4  Sin[w]}, {w, -Pi/2, Pi/2},  
          DisplayFunction  ->  Identity];
s2 = 
        ParametricPlot[{3 + 2  Cos[v] ,  2  Sin[v]}, {v, Pi/3, 7  Pi/4},  
          DisplayFunction  ->  Identity];
Show[s1, s2,  
      DisplayFunction  ->  $DisplayFunction]

Ôï óôåñåü  Á ìðïñïýìå íá ôï "äïýìå" íá ó÷çìáôßæåôáé áöïý ðñïçãïõìÝíùò  ðåñéïñéóôïýìå óå êáôÜëëçëá ôìÞìáôá ôùí åìðëåêïìÝíùí óöáéñþí:

Clear[t, u, sph1, sph2]
sph1 = 
        ParametricPlot3D[{4  Cos[t] Cos[u], 4  Sin[t] Cos[u],  
            4  Sin[u]}, {t, 0, 5  Pi/3}, {u,  -Pi/3,  Pi/3},  
          DisplayFunction  ->  Identity];
sph2 = 
        ParametricPlot3D[{3 + 2  Cos[t] Cos[u], 2  Sin[t] Cos[u],  
            2  Sin[u]}, {t, Pi/6, 6  Pi/4}, {u,  -Pi/2,  Pi/2}, 
          DisplayFunction  ->  Identity];
Show[sph1,  sph2, Axes -> True, 
        AxesLabel -> {"x", "y", "z"}, Boxed -> False,  
        DisplayFunction  -> $DisplayFunction ];

ÔÝëïò ï üãêïò ôïõ óôåñåïý Á äßíåôáé áðü ôï æçôïýìåíï ïëïêëÞñùìá

(∫_1^((21/6)))((∫_(-(Sqrt[(2^2) - (((x - 3))^2)]))^((Sqrt[(2^2) - (((x - 3))^2)])))((∫_(-(Sqrt[(2^2) - (y^2) - (((x - 3))^2)]))^((Sqrt[(2^2) - (y^2) - (((x - 3))^2)]))) 
                    1 \[DifferentialD]z) \[DifferentialD]y) \[DifferentialD]x+(∫_((21/6))^4)((∫_(-(Sqrt[(4^2) - (x^2)]))^((Sqrt[(4^2) - (x^2)])))((∫_(-(Sqrt[(4^2) - (y^2) - (x^2)]))^((Sqrt[(4^2) - (y^2) - (x^2)]))) 1 \[DifferentialD]z) \[DifferentialD]y) \[DifferentialD]x

◆ÁëëáãÞ ôçò óåéñÜò ïëïêëçñþóåùò

Äßíåôáé ôï ôñéðëü ïëïêëÞñùìá  \(∫\_0\%1\(∫\_0\%z\(∫\_0\%y\ \ xy\^2\ z\^3\)\)\) dx dy dz 
êáé Ýóôù üôé åðéèõìïýìå íá áëëÜîïõìå ôçí óåéñÜ ïëïêëÞñùóçò þóôå íá åßíáé dx dz dy.

Õðïëïãßæïõìå ôï ïëïêëÞñùìá (êÜôé ôÝôïéï, åí ãÝíåé, äåí ÷ñåéÜæåôáé)

Clear[x, y, z]
(∫_0^1)(∫_0^z)(∫_0^y) 
              x (y^2) (z^3) \[DifferentialD]x \[DifferentialD]y \[DifferentialD]z

ÊÜíïõìå ôçí ãñáöéêÞ ðáñÜóôáóç ôïõ óôåñåïý åðß ôïõ ïðïßïõ ïëïêëçñþíïõìå (áõôü åßíáé ðÜíôá ÷ñÞóéìï ü÷é üìùò êáé áðáñáßôçôï).

Clear[setup1, setupyz]
setupyz = 
        Graphics3D[{Red, Thickness[ .02], 
            Line[{{0, 0, 0}, {0, 0, 1}, {0, 1, 1}, {0, 0, 0}}]}];
setup1 = 
        Graphics3D[{Blue, Thickness[ .02], 
            Line[{{0, 0, 0}, {1, 1, 1}, {0, 0, 1}}]}];
setup2 = 
        Graphics3D[{Blue, Thickness[ .02], 
            Line[{{1, 1, 1}, {0, 1, 1}}]}];
Show[setupyz, setup1, setup2, Axes -> True, 
        AxesLabel -> {"x", "y", "z"}, 
        Boxed -> False];

Ôï óôåñåü åðß ôïõ ïðïßïõ ïëïêëçñþíïõìå åßíáé ôï ðáñáðÜíù ôåôñÜåäñï ðïõ Ý÷åé âÜóç óôï yz-åðßðåäï (êüêêéíï ôñßãùíï). Ìðïñïýìå, åíáëáêôéêÜ, íá èåùñÞóïõìå üôé ç âÜóç ôïõ âñßóêåôáé óôï åðßðåäï y=x, Þ óôï åðßðåäï y=z.

Tï äïèÝí ïëïêëÞñùìá (∫_0^1)(∫_0^z)(∫_0^y) 
              (y^2) (z^3)dx dy dz õðïëïãßæåôáé ìå ôá åîÞò âÞìáôá (ùò ðñïò ôçí óåéñÜ ïëïêëÞñùóçò dx dy dz).

Clear[strip, setup3, step]
step =  .6;
strip[t_] := 
        Graphics3D[
          Polygon[{{0, 0, t}, {0, t, t}, {t, t, t}, {0, 0, t}}]];
setup3 = strip[step];
line = 
        Graphics3D[{Goldenrod, Thickness[ .018], 
            Line[{{0, (step/2), step}, {(step/2), (step/2), 
                  step}}]}];
point1 = 
        Graphics3D[{Red, 
            Text["{0,y,z}", {0, (step/2), step +  .12}]}];
point2 = 
        Graphics3D[{Red, 
            Text["{y,y,z}", {step -  .12, (step/2) -  .10, 
                step +  .2}]}];
point3 = 
        Graphics3D[{Red, 
            Text["{0,z,z}", { .2, step, step +  .10}]}];
point4 = 
        Graphics3D[{Red, 
            Text["{0,0,z}", { .16, 0, step +  .14}]}];
Show[setupyz, setup1, setup2, setup3, line, point1, point2, point3, 
        point4, Axes -> True, Boxed -> False, 
        AxesLabel -> {"x", "y", "z"}, 
        Boxed -> False];

Ç åãêÜñóéá óôï Üîïíá 0z ôïìÞ (ôï ìðëÝ ôñßãùíï) "óáñþíåé" ôï óôåñåü ìáò ùò åîÞò:

Clear[strip, setup3, step]
step =  .1
strip[t_] := 
        Graphics3D[
          Polygon[{{0, 0, t}, {0, t, t}, {t, t, t}, {0, 0, t}}]];
setup3 = Table[strip[t], {t, step, 1, step}];
Show[setupyz, setup1, setup2, setup3, Axes -> True, 
        AxesLabel -> {"x", "y", "z"}, 
        Boxed -> False];

Clear[strip, setup3, step]
step =  .05
strip[t_] := 
        Graphics3D[
          Polygon[{{0, 0, t}, {0, t, t}, {t, t, t}, {0, 0, t}}]];
setup3 = Table[strip[t], {t, step, 1, step}];
\(\(method1 = 
	
        Show[setupyz, setup1, setup2, setup3, Axes -> True, 
          AxesLabel -> {"x", "y", "z"}, 
          Boxed -> False];\)\)

Ôï ßäéï "óÜñùìá" ôïõ ôåôñáÝäñïõ ìðïñïýìå íá ôï åðéôý÷ïõìå üðùò öáßíåôáé óôï ðáñáêÜôù ó÷Þìá

Clear[strip2, setup4, strip2, step]
step =  .5;
strip2[t_] := 
        Graphics3D[
          Polygon[{{0, t, t}, {0, t, 1}, {t, t, 1}, {t, t, t}, {0, t, 
                t}}]];
setup4 = Table[strip2[t], {t, step, 1, step}];
line = 
        Graphics3D[{Goldenrod, Thickness[ .01], 
            Line[{{0, step, step +  .25}, {step, step, 
                  step +  .25}}]}];
point1 = 
        Graphics3D[{Red, 
            Text["{0,y,z}", {0, step -  .15, step +  .30}]}];
point2 = 
        Graphics3D[{Red, 
            Text["{y,y,z}", {step +  .20, step -  .15, 
                step +  .35}]}];
point3 = 
        Graphics3D[{Red, 
            Text["{0,y,y}", {0, step -  .15, step +  .10}]}];
point4 = 
        Graphics3D[{Red, 
            Text["{0,y,1}", {0, step -  .15, 1.10}]}];
Show[setupyz, setup1, setup2, setup4, line, point1, point2, point3, 
        point4, Axes -> True, 
        AxesLabel -> {"x", "y", "z"}, 
        Boxed -> False];

Ôþñá, ç åãêÜñóéá óôïí Üîïíá 0y ôïìÞ åßíáé Ýíá ïñèïãþíéï ðáñáëëçëüãñáììï.

Clear[strip2, setup4, strip2, step]
step =  .1;
strip2[t_] := 
        Graphics3D[
          Polygon[{{0, t, t}, {0, t, 1}, {t, t, 1}, {t, t, t}, {0, t, 
                t}}]];
setup4 = Table[strip2[t], {t, step, 1, step}];
Show[setupyz, setup1, setup2, setup4, Axes -> True, 
        AxesLabel -> {"x", "y", "z"}, 
        Boxed -> False];

Clear[strip2, setup4, strip2, step]
step =  .05;
strip2[t_] := 
        Graphics3D[
          Polygon[{{0, t, t}, {0, t, 1}, {t, t, 1}, {t, t, t}, {0, t, 
                t}}]];
setup4 = Table[strip2[t], {t, step, 1, step}];
method2 = 
        Show[setupyz, setup1, setup2, setup4, Axes -> True, 
          AxesLabel -> {"x", "y", "z"}, 
          Boxed -> False];

Óôï äåýôåñï "óÜñùìá" ïëïêëçñþóáìå ðñþôá ùò ðñïò x (êßôñéíç ãñáììÞ). Óôï äåýôåñï âÞìá êñáôÜìå ôï y óôáèåñü êáé ïëïêëçñþíïõìå ùò ðñïò z  ìå üñéá ôï y êáé ôï 1. ÔÝëïò ïëïêëçñþíïõìå ùò ðñïò y áðü ôï 0 Ýùò ôï 1. 
ÔåëéêÜ, åÜí áëëÜîïõìå ôçí óåéñÜ ïëïêëÞñùóçò áðü 
åßíáé dx dy dz óå dx dz dy ôï ïëïêëÞñùìá 

\(∫\_0\%1\(∫\_0\%z\(∫\_0\%y\ \ xy\^2\ z\^3\)\)\) dx dy dz  ãßíåôáé    (∫_0^1)(∫_y^1)(∫_0^y)x (y^2) (z^3)dx dz dy. 
ÃéÜ ôïõ ëüãïõ ôï áëçèÝò :

Clear[x, y, z]
(∫_0^1)(∫_y^1)(∫_0^y) 
                x (y^2) (z^3) \[DifferentialD]x \[DifferentialD]z \[DifferentialD]y  ==  (∫_0^1)(∫_0^z)(∫_0^y) 
                x (y^2) (z^3) \[DifferentialD]x \[DifferentialD]y \[DifferentialD]z

◆ÊÝíôñï ìÜæáò êáé ñïðÞ áäñáíåßáò

ÅÜí D åßíáé Ýíá 3-äéÜóôáôï óôåñåü êáé  δ(x,y,z) ìßá óõíå÷Þò óõíÜñôçóç ðïõ äßíåé ôçí ðõêíüôçôá êáé åßíáé ïñéóìÝíç ãéÜ êÜèå óçìåßï (x,y,z) óôï D, ôüôå Ý÷ïõìå ôïõò ðáñáêÜôù ôýðïõò:

ÌÜæá:		M = \(TraditionalForm\`∫\(∫∫\_\(\(\ \)\(D\)\)\)\)  δ(x,y,z) dV
   
   
ÊÝíôñï ìÜæáò:
	
          x*=((∫∫(∫_( D))x  δ(x, y, z)dV)/(\(\(M\)\(\ \)\(
\)\)))
          y*=((∫∫(∫_( D))y δ(x, y, z)dV    )/(M   )) 
      
          z* = ((∫∫(∫_( D)) z  δ(x, y, z)dV    )/(M   ))
	
	
ÑïðÞ áäñáíåßáò ùò ðñïò ôïõò Üîïíåò :

(Subscript[I, x])=∫∫(∫_( D))  ((y^2) + (z^2)) δ(x,y,z) dV

(Subscript[I, y])=∫∫(∫_( D))  ((x^2) + (z^2)) δ(x,y,z) dV

(Subscript[I, z])=∫∫(∫_( D))  ((x^2) + (y^2)) δ(x,y,z) dV	

	
ÑïðÞ áäñáíåßáò ùò ðñïò ôá åðßðåäá :

(Subscript[I, xy])=∫∫(∫_( D))  (z^2) δ(x,y,z) dV


(Subscript[I, yz])=∫∫(∫_( D))  (x^2) δ(x,y,z) dV


(Subscript[I, xz])=∫∫(∫_( D))  (y^2) δ(x,y,z) dV

ÐáñÜäåéãìá 1

Õðïëïãßóôå ôéò óõíôåôáãìÝíåò ôïõ êÝíôñïõ ìÜæáò ôïõ óôåñåïý R ðïõ öñÜóóåôáé áðü ôéò åðéöÜíåéåò 
          y= ((x^2)+(z^2)) êáé y = 18 - ((x^2)+(z^2)), 
åÜí ç ðõêíüôçôá óå êÜèå óçìåßï äßíåôáé áðü ôçí óõíÜñôçóç  δ(x,y,z) = (x^2) + (y^2)+ (z^2)

Tï óôåñåü R åßíáé ôï ßäéï óôåñåü ìå áõôü ðïõ ÷ñçóéìïðïéÞóáìå óôá ðáñáðÜíù ðáñáäåßãìáôá ïëïêëÞñùóçò.

Clear[δ, x, y, z, g1, 
      g2]
g1[x_, z_] = (x^2) + (z^2);
g2[x_, z_] = 
        18 - ((x^2) + (z^2));
δ[x_, y_, z_] = 
        (x^2) + (y^2) + (z^2);
ÌÜæá=(∫_(-3)^3)(∫_(-(Sqrt[9 - (x^2)]))^((Sqrt[9 - (x^2)])))(∫_(g1[x, z])^(g2[x, z]))δ[x, y, 
                z] \[DifferentialD]y \[DifferentialD]z \[DifferentialD]x

xÊÝíôñïõÌÜæáò=(((∫_(-3)^3)(∫_(-(Sqrt[9 - (x^2)]))^((Sqrt[9 - (x^2)])))(∫_(g1[x, z])^(g2[x, z]))x δ[
                  x, y, z] \[DifferentialD]y \[DifferentialD]z \[DifferentialD]x)/(ÌÜæá))

yÊÝíôñïõÌÜæáò=(((∫_(-3)^3)(∫_(-(Sqrt[9 - (x^2)]))^((Sqrt[9 - (x^2)])))(∫_(g1[x, z])^(g2[x, z])) y δ[x, y, 
                  z] \[DifferentialD]y \[DifferentialD]z \[DifferentialD]x)/(ÌÜæá))

zÊÝíôñïõÌÜæáò=(((∫_(-3)^3)(∫_(-(Sqrt[9 - (x^2)]))^((Sqrt[9 - (x^2)])))(∫_(g1[x, z])^(g2[x, z])) z δ[x, y, 
                  z] \[DifferentialD]y \[DifferentialD]z \[DifferentialD]x)/(ÌÜæá))

ÐáñÜäåéãìá 2

Õðïëïãßóôå ôçí ñïðÞ áäñáíåßáò ôïõ óôåñåïý R ãýñù áðü ôïí Üîïíá 0y ðïõ öñÜóóåôáé áðü ôéò åðéöÜíåéåò 
          y= ((x^2)+(z^2)) êáé y = 18 - ((x^2)+(z^2)), 
åÜí ç ðõêíüôçôá óå êÜèå óçìåßï äßíåôáé áðü ôçí óõíÜñôçóç  δ(x,y,z) = (x^2) + (y^2)+ (z^2)

Clear[δ, x, y, z, g1, 
      g2]
g1[x_, z_] = (x^2) + (z^2);
g2[x_, z_] = 
        18 - ((x^2) + (z^2));
δ[x_, y_, z_] = 
        (x^2) + (y^2) + (z^2);
yÑïðÞÁäñáíåßáò=(∫_(-3)^3)(∫_(-(Sqrt[9 - (x^2)]))^((Sqrt[9 - (x^2)])))(∫_(g1[x, z])^(g2[x, z]))((x^2) + (z^2)) δ[x, y, 
                z] \[DifferentialD]y \[DifferentialD]z \[DifferentialD]x

ÐáñÜäåéãìá 3

Õðïëïãßóôå ôçí ñïðÞ áäñáíåßáò ôïõ óôåñåïý R, ãýñù áðü ôï åðßðåäï xy, ðïõ öñÜóóåôáé áðü ôéò åðéöÜíåéåò 
          y= ((x^2)+(z^2)) êáé y = 18 - ((x^2)+(z^2)), 
åÜí ç ðõêíüôçôá óå êÜèå óçìåßï äßíåôáé áðü ôçí óõíÜñôçóç  δ(x,y,z) = (x^2) + (y^2)+ (z^2)

Clear[δ, x, y, z, g1, g2]
g1[x_, z_] = (x^2) + (z^2);
g2[x_, z_] = 18 - ((x^2) + (z^2));
δ[x_, y_, z_] = (x^2) + (y^2) + (z^2);
xyÑïðÞÁäñáíåßáò = (∫_(-3)^3)(∫_(-(Sqrt[9 - (x^2)]))^((Sqrt[9 - (x^2)])))(∫_(g1[x, z])^(g2[x, z]))(z^2) δ[x, y, 
                z] \[DifferentialD]y \[DifferentialD]z \[DifferentialD]x