明治大学 先端数理科学研究科 現象数理学専攻 2019年8月実施 微积分
Author
思齐塾 , 祭音Myyura
Description
(1) f ( x , y ) = x y ( x 2 + y 2 − 1 ) f(x,y) = xy(x^2 + y^2 - 1) f ( x , y ) = x y ( x 2 + y 2 − 1 ) とする。
(a) f f f のヤコビ行列とヘッセ行列を求めよ。
(b) f f f の極値をすべて求めよ。
(c) f f f のグラフ z = f ( x , y ) z = f(x,y) z = f ( x , y ) の ( 1 , 1 , 1 ) (1,1,1) ( 1 , 1 , 1 ) における接平面の方程式を求めよ。
(2) B B B を単位球 x 2 + y 2 + z 2 ≤ 1 x^2 + y^2 + z^2 \leq 1 x 2 + y 2 + z 2 ≤ 1 とし, g ( x , y , z ) = x 2 y 2 z 2 g(x,y,z) = x^2y^2z^2 g ( x , y , z ) = x 2 y 2 z 2 とおく。
(a) g g g の B B B における最大値を求めよ。
(b) ∭ B g ( x , y , z ) d x d y d z \iiint_B g(x,y,z) dx dy dz ∭ B g ( x , y , z ) d x d y d z の値を求めよ。
题目描述
(1) 设
f ( x , y ) = x y ( x 2 + y 2 − 1 ) . f(x,y)=xy(x^2+y^2-1). f ( x , y ) = x y ( x 2 + y 2 − 1 ) .
(a) 求 f f f 的 Jacobian 矩阵与 Hessian 矩阵。
(b) 求 f f f 的全部极值。
(c) 求曲面 z = f ( x , y ) z=f(x,y) z = f ( x , y ) 在点 ( 1 , 1 , 1 ) (1,1,1) ( 1 , 1 , 1 ) 处的切平面方程。
(2) 设 B B B 为单位球体
x 2 + y 2 + z 2 ≤ 1 , x^2+y^2+z^2\leq 1, x 2 + y 2 + z 2 ≤ 1 ,
并令
g ( x , y , z ) = x 2 y 2 z 2 . g(x,y,z)=x^2y^2z^2. g ( x , y , z ) = x 2 y 2 z 2 .
(a) 求 g g g 在 B B B 上的最大值。
(b) 计算
∭ B g ( x , y , z ) d x d y d z . \iiint_B g(x,y,z)\,dx\,dy\,dz. ∭ B g ( x , y , z ) d x d y d z .
Kai
(1)
(a)
f ( x , y ) = x y ( x 2 + y 2 − 1 ) = x 3 y + x y 3 − x y f(x,y) = xy(x^2+y^2-1) = x^3y + xy^3 - xy f ( x , y ) = x y ( x 2 + y 2 − 1 ) = x 3 y + x y 3 − x y
∂ f ∂ x = 3 x 2 y + y 3 − y \frac{\partial f}{\partial x} = 3x^2y + y^3 - y ∂ x ∂ f = 3 x 2 y + y 3 − y
∂ f ∂ y = x 3 + 3 x y 2 − x \frac{\partial f}{\partial y} = x^3 + 3xy^2 - x ∂ y ∂ f = x 3 + 3 x y 2 − x
ヤコビ行列は
J f = ( ∂ f ∂ x ∂ f ∂ y ) = ( 3 x 2 y + y 3 − y x 3 + 3 x y 2 − x ) J_f = \begin{pmatrix} \frac{\partial f}{\partial x} & \frac{\partial f}{\partial y} \end{pmatrix} = \begin{pmatrix} 3x^2y + y^3 - y & x^3 + 3xy^2 - x \end{pmatrix} J f = ( ∂ x ∂ f ∂ y ∂ f ) = ( 3 x 2 y + y 3 − y x 3 + 3 x y 2 − x )
ヘッセ行列は
∂ 2 f ∂ x 2 = 6 x y \frac{\partial^2 f}{\partial x^2} = 6xy ∂ x 2 ∂ 2 f = 6 x y
∂ 2 f ∂ x ∂ y = 3 x 2 + 3 y 2 − 1 \frac{\partial^2 f}{\partial x \partial y} = 3x^2 + 3y^2 - 1 ∂ x ∂ y ∂ 2 f = 3 x 2 + 3 y 2 − 1
∂ 2 f ∂ y 2 = 6 x y \frac{\partial^2 f}{\partial y^2} = 6xy ∂ y 2 ∂ 2 f = 6 x y
H f = ( ∂ 2 f ∂ x 2 ∂ 2 f ∂ x ∂ y ∂ 2 f ∂ y ∂ x ∂ 2 f ∂ y 2 ) = ( 6 x y 3 x 2 + 3 y 2 − 1 3 x 2 + 3 y 2 − 1 6 x y ) H_f = \begin{pmatrix} \frac{\partial^2 f}{\partial x^2} & \frac{\partial^2 f}{\partial x \partial y} \\ \frac{\partial^2 f}{\partial y \partial x} & \frac{\partial^2 f}{\partial y^2} \end{pmatrix} = \begin{pmatrix} 6xy & 3x^2 + 3y^2 - 1 \\ 3x^2 + 3y^2 - 1 & 6xy \end{pmatrix} H f = ( ∂ x 2 ∂ 2 f ∂ y ∂ x ∂ 2 f ∂ x ∂ y ∂ 2 f ∂ y 2 ∂ 2 f ) = ( 6 x y 3 x 2 + 3 y 2 − 1 3 x 2 + 3 y 2 − 1 6 x y )
(b)
∂ f ∂ x = 3 x 2 y + y 3 − y = y ( 3 x 2 + y 2 − 1 ) = 0 \frac{\partial f}{\partial x} = 3x^2y + y^3 - y = y(3x^2 + y^2 - 1) = 0 ∂ x ∂ f = 3 x 2 y + y 3 − y = y ( 3 x 2 + y 2 − 1 ) = 0
∂ f ∂ y = x 3 + 3 x y 2 − x = x ( x 2 + 3 y 2 − 1 ) = 0 \frac{\partial f}{\partial y} = x^3 + 3xy^2 - x = x(x^2 + 3y^2 - 1) = 0 ∂ y ∂ f = x 3 + 3 x y 2 − x = x ( x 2 + 3 y 2 − 1 ) = 0
Case 1: x = 0 x=0 x = 0 . Then y ( y 2 − 1 ) = 0 y(y^2 - 1) = 0 y ( y 2 − 1 ) = 0 , so y = 0 , 1 , − 1 y = 0, 1, -1 y = 0 , 1 , − 1 . Critical points: ( 0 , 0 ) , ( 0 , 1 ) , ( 0 , − 1 ) (0,0), (0,1), (0,-1) ( 0 , 0 ) , ( 0 , 1 ) , ( 0 , − 1 ) .
Case 2: y = 0 y=0 y = 0 . Then x ( x 2 − 1 ) = 0 x(x^2 - 1) = 0 x ( x 2 − 1 ) = 0 , so x = 0 , 1 , − 1 x = 0, 1, -1 x = 0 , 1 , − 1 . Critical points: ( 0 , 0 ) , ( 1 , 0 ) , ( − 1 , 0 ) (0,0), (1,0), (-1,0) ( 0 , 0 ) , ( 1 , 0 ) , ( − 1 , 0 ) .
Case 3: x ≠ 0 x \neq 0 x = 0 and y ≠ 0 y \neq 0 y = 0 . Then 3 x 2 + y 2 = 1 3x^2 + y^2 = 1 3 x 2 + y 2 = 1 and x 2 + 3 y 2 = 1 x^2 + 3y^2 = 1 x 2 + 3 y 2 = 1 . Subtracting gives 2 x 2 − 2 y 2 = 0 2x^2 - 2y^2 = 0 2 x 2 − 2 y 2 = 0 , so x 2 = y 2 x^2 = y^2 x 2 = y 2 , and x = ± y x = \pm y x = ± y . Substituting x = y x = y x = y into 3 x 2 + y 2 = 1 3x^2 + y^2 = 1 3 x 2 + y 2 = 1 gives 4 x 2 = 1 4x^2 = 1 4 x 2 = 1 , so x = ± 1 2 x = \pm \frac{1}{2} x = ± 2 1 , and y = ± 1 2 y = \pm \frac{1}{2} y = ± 2 1 . Critical points: ( 1 2 , 1 2 ) , ( 1 2 , − 1 2 ) , ( − 1 2 , 1 2 ) , ( − 1 2 , − 1 2 ) (\frac{1}{2}, \frac{1}{2}), (\frac{1}{2}, -\frac{1}{2}), (-\frac{1}{2}, \frac{1}{2}), (-\frac{1}{2}, -\frac{1}{2}) ( 2 1 , 2 1 ) , ( 2 1 , − 2 1 ) , ( − 2 1 , 2 1 ) , ( − 2 1 , − 2 1 ) .
Critical points: ( 0 , 0 ) , ( 0 , 1 ) , ( 0 , − 1 ) , ( 1 , 0 ) , ( − 1 , 0 ) , ( 1 2 , 1 2 ) , ( 1 2 , − 1 2 ) , ( − 1 2 , 1 2 ) , ( − 1 2 , − 1 2 ) (0,0), (0,1), (0,-1), (1,0), (-1,0), (\frac{1}{2}, \frac{1}{2}), (\frac{1}{2}, -\frac{1}{2}), (-\frac{1}{2}, \frac{1}{2}), (-\frac{1}{2}, -\frac{1}{2}) ( 0 , 0 ) , ( 0 , 1 ) , ( 0 , − 1 ) , ( 1 , 0 ) , ( − 1 , 0 ) , ( 2 1 , 2 1 ) , ( 2 1 , − 2 1 ) , ( − 2 1 , 2 1 ) , ( − 2 1 , − 2 1 ) .
Evaluate f f f at these points:
f ( 0 , 0 ) = 0 f(0,0) = 0 f ( 0 , 0 ) = 0 , f ( 0 , 1 ) = 0 f(0,1) = 0 f ( 0 , 1 ) = 0 , f ( 0 , − 1 ) = 0 f(0,-1) = 0 f ( 0 , − 1 ) = 0 , f ( 1 , 0 ) = 0 f(1,0) = 0 f ( 1 , 0 ) = 0 , f ( − 1 , 0 ) = 0 f(-1,0) = 0 f ( − 1 , 0 ) = 0 .
f ( 1 2 , 1 2 ) = 1 4 ( 1 4 + 1 4 − 1 ) = 1 4 ( − 1 2 ) = − 1 8 f(\frac{1}{2}, \frac{1}{2}) = \frac{1}{4}(\frac{1}{4} + \frac{1}{4} - 1) = \frac{1}{4}(-\frac{1}{2}) = -\frac{1}{8} f ( 2 1 , 2 1 ) = 4 1 ( 4 1 + 4 1 − 1 ) = 4 1 ( − 2 1 ) = − 8 1 .
f ( 1 2 , − 1 2 ) = − 1 4 ( 1 4 + 1 4 − 1 ) = 1 8 f(\frac{1}{2}, -\frac{1}{2}) = -\frac{1}{4}(\frac{1}{4} + \frac{1}{4} - 1) = \frac{1}{8} f ( 2 1 , − 2 1 ) = − 4 1 ( 4 1 + 4 1 − 1 ) = 8 1 .
f ( − 1 2 , 1 2 ) = − 1 4 ( 1 4 + 1 4 − 1 ) = 1 8 f(-\frac{1}{2}, \frac{1}{2}) = -\frac{1}{4}(\frac{1}{4} + \frac{1}{4} - 1) = \frac{1}{8} f ( − 2 1 , 2 1 ) = − 4 1 ( 4 1 + 4 1 − 1 ) = 8 1 .
f ( − 1 2 , − 1 2 ) = 1 4 ( 1 4 + 1 4 − 1 ) = − 1 8 f(-\frac{1}{2}, -\frac{1}{2}) = \frac{1}{4}(\frac{1}{4} + \frac{1}{4} - 1) = -\frac{1}{8} f ( − 2 1 , − 2 1 ) = 4 1 ( 4 1 + 4 1 − 1 ) = − 8 1 .
The Hessian classification is as follows:
points f det H f classification ( 0 , 0 ) , ( 0 , ± 1 ) , ( ± 1 , 0 ) 0 − 1 or − 4 saddle ( 1 2 , 1 2 ) , ( − 1 2 , − 1 2 ) − 1 8 2 local minimum ( 1 2 , − 1 2 ) , ( − 1 2 , 1 2 ) 1 8 2 local maximum \begin{array}{c|c|c|c}
\text{points} & f & \det H_f & \text{classification}\\ \hline
(0,0),\ (0,\pm1),\ (\pm1,0) & 0 & -1\ \text{or}\ -4 & \text{saddle}\\
(\frac12,\frac12),\ (-\frac12,-\frac12) & -\frac18 & 2 & \text{local minimum}\\
(\frac12,-\frac12),\ (-\frac12,\frac12) & \frac18 & 2 & \text{local maximum}
\end{array} points ( 0 , 0 ) , ( 0 , ± 1 ) , ( ± 1 , 0 ) ( 2 1 , 2 1 ) , ( − 2 1 , − 2 1 ) ( 2 1 , − 2 1 ) , ( − 2 1 , 2 1 ) f 0 − 8 1 8 1 det H f − 1 or − 4 2 2 classification saddle local minimum local maximum
Indeed, at the two minimum points the Hessian eigenvalues are 1 , 2 1,2 1 , 2 , while at the two maximum points they are − 1 , − 2 -1,-2 − 1 , − 2 . Thus these four points give all the extrema.
(c)
z = f ( x , y ) = x y ( x 2 + y 2 − 1 ) z = f(x,y) = xy(x^2+y^2-1) z = f ( x , y ) = x y ( x 2 + y 2 − 1 )
接平面の方程式は
z − f ( 1 , 1 ) = f x ( 1 , 1 ) ( x − 1 ) + f y ( 1 , 1 ) ( y − 1 ) z - f(1,1) = f_x(1,1)(x-1) + f_y(1,1)(y-1) z − f ( 1 , 1 ) = f x ( 1 , 1 ) ( x − 1 ) + f y ( 1 , 1 ) ( y − 1 )
f ( 1 , 1 ) = 1 ( 1 + 1 − 1 ) = 1 f(1,1) = 1(1+1-1) = 1 f ( 1 , 1 ) = 1 ( 1 + 1 − 1 ) = 1
f x ( 1 , 1 ) = 3 ( 1 ) 2 ( 1 ) + ( 1 ) 3 − 1 = 3 + 1 − 1 = 3 f_x(1,1) = 3(1)^2(1) + (1)^3 - 1 = 3+1-1 = 3 f x ( 1 , 1 ) = 3 ( 1 ) 2 ( 1 ) + ( 1 ) 3 − 1 = 3 + 1 − 1 = 3
f y ( 1 , 1 ) = ( 1 ) 3 + 3 ( 1 ) ( 1 ) 2 − 1 = 1 + 3 − 1 = 3 f_y(1,1) = (1)^3 + 3(1)(1)^2 - 1 = 1+3-1 = 3 f y ( 1 , 1 ) = ( 1 ) 3 + 3 ( 1 ) ( 1 ) 2 − 1 = 1 + 3 − 1 = 3
Therefore, z − 1 = 3 ( x − 1 ) + 3 ( y − 1 ) z - 1 = 3(x-1) + 3(y-1) z − 1 = 3 ( x − 1 ) + 3 ( y − 1 )
z − 1 = 3 x − 3 + 3 y − 3 z - 1 = 3x - 3 + 3y - 3 z − 1 = 3 x − 3 + 3 y − 3
z = 3 x + 3 y − 5 z = 3x + 3y - 5 z = 3 x + 3 y − 5
(2)
(a) Using spherical coordinates x = r sin θ cos ϕ x = r \sin \theta \cos \phi x = r sin θ cos ϕ , y = r sin θ sin ϕ y = r \sin \theta \sin \phi y = r sin θ sin ϕ , z = r cos θ z = r \cos \theta z = r cos θ , where 0 ≤ r ≤ 1 0 \le r \le 1 0 ≤ r ≤ 1 , 0 ≤ θ ≤ π 0 \le \theta \le \pi 0 ≤ θ ≤ π , 0 ≤ ϕ ≤ 2 π 0 \le \phi \le 2\pi 0 ≤ ϕ ≤ 2 π .
g ( x , y , z ) = x 2 y 2 z 2 = ( r sin θ cos ϕ ) 2 ( r sin θ sin ϕ ) 2 ( r cos θ ) 2 = r 6 sin 4 θ cos 2 θ cos 2 ϕ sin 2 ϕ = r 6 sin 4 θ cos 2 θ ( 1 4 sin 2 2 ϕ ) g(x,y,z) = x^2y^2z^2 = (r \sin \theta \cos \phi)^2(r \sin \theta \sin \phi)^2(r \cos \theta)^2 = r^6 \sin^4 \theta \cos^2 \theta \cos^2 \phi \sin^2 \phi = r^6 \sin^4 \theta \cos^2 \theta (\frac{1}{4} \sin^2 2\phi) g ( x , y , z ) = x 2 y 2 z 2 = ( r sin θ cos ϕ ) 2 ( r sin θ sin ϕ ) 2 ( r cos θ ) 2 = r 6 sin 4 θ cos 2 θ cos 2 ϕ sin 2 ϕ = r 6 sin 4 θ cos 2 θ ( 4 1 sin 2 2 ϕ )
We want to maximize this expression subject to 0 ≤ r ≤ 1 0 \le r \le 1 0 ≤ r ≤ 1 , 0 ≤ θ ≤ π 0 \le \theta \le \pi 0 ≤ θ ≤ π , 0 ≤ ϕ ≤ 2 π 0 \le \phi \le 2\pi 0 ≤ ϕ ≤ 2 π .
Since r ≤ 1 r \le 1 r ≤ 1 , r 6 r^6 r 6 is maximized when r = 1 r = 1 r = 1 .
sin 4 θ cos 2 θ \sin^4 \theta \cos^2 \theta sin 4 θ cos 2 θ is maximized when θ \theta θ is between 0 and π \pi π .
d d θ ( sin 4 θ cos 2 θ ) = 4 sin 3 θ cos 3 θ − 2 sin 5 θ cos θ = 2 sin 3 θ cos θ ( 2 cos 2 θ − sin 2 θ ) = 0 \frac{d}{d\theta} (\sin^4 \theta \cos^2 \theta) = 4\sin^3 \theta \cos^3 \theta - 2\sin^5 \theta \cos \theta = 2 \sin^3 \theta \cos \theta (2 \cos^2 \theta - \sin^2 \theta) = 0 d θ d ( sin 4 θ cos 2 θ ) = 4 sin 3 θ cos 3 θ − 2 sin 5 θ cos θ = 2 sin 3 θ cos θ ( 2 cos 2 θ − sin 2 θ ) = 0 .
So either sin θ = 0 \sin \theta = 0 sin θ = 0 , cos θ = 0 \cos \theta = 0 cos θ = 0 , or 2 cos 2 θ = sin 2 θ 2 \cos^2 \theta = \sin^2 \theta 2 cos 2 θ = sin 2 θ . Since 2 cos 2 θ = sin 2 θ 2 \cos^2 \theta = \sin^2 \theta 2 cos 2 θ = sin 2 θ , 2 cos 2 θ = 1 − cos 2 θ 2 \cos^2 \theta = 1 - \cos^2 \theta 2 cos 2 θ = 1 − cos 2 θ , so 3 cos 2 θ = 1 3 \cos^2 \theta = 1 3 cos 2 θ = 1 , and cos θ = ± 1 3 \cos \theta = \pm \frac{1}{\sqrt{3}} cos θ = ± 3 1 . Then sin θ = ± 2 3 \sin \theta = \pm \sqrt{\frac{2}{3}} sin θ = ± 3 2 . Thus, the maximum occurs when θ = arccos ( 1 3 ) \theta = \arccos (\frac{1}{\sqrt{3}}) θ = arccos ( 3 1 ) and sin θ = 2 3 \sin \theta = \sqrt{\frac{2}{3}} sin θ = 3 2 .
Thus, the maximum value of sin 4 θ cos 2 θ \sin^4 \theta \cos^2 \theta sin 4 θ cos 2 θ is ( 2 3 ) 2 1 3 = 4 27 (\frac{2}{3})^2 \frac{1}{3} = \frac{4}{27} ( 3 2 ) 2 3 1 = 27 4 .
1 4 sin 2 2 ϕ \frac{1}{4} \sin^2 2\phi 4 1 sin 2 2 ϕ is maximized when sin 2 ϕ = ± 1 \sin 2\phi = \pm 1 sin 2 ϕ = ± 1 , which means 2 ϕ = π 2 , 3 π 2 2\phi = \frac{\pi}{2}, \frac{3\pi}{2} 2 ϕ = 2 π , 2 3 π . So ϕ = π 4 , 3 π 4 \phi = \frac{\pi}{4}, \frac{3\pi}{4} ϕ = 4 π , 4 3 π .
The maximum value of 1 4 sin 2 2 ϕ \frac{1}{4} \sin^2 2\phi 4 1 sin 2 2 ϕ is 1 4 \frac{1}{4} 4 1 .
Therefore, the maximum value of g ( x , y , z ) g(x,y,z) g ( x , y , z ) is 1 6 ⋅ 4 27 ⋅ 1 4 = 1 27 1^6 \cdot \frac{4}{27} \cdot \frac{1}{4} = \frac{1}{27} 1 6 ⋅ 27 4 ⋅ 4 1 = 27 1 .
(b)
∭ B g ( x , y , z ) d x d y d z = ∭ B x 2 y 2 z 2 d x d y d z = ∫ 0 1 ∫ 0 π ∫ 0 2 π ( r sin θ cos ϕ ) 2 ( r sin θ sin ϕ ) 2 ( r cos θ ) 2 r 2 sin θ d ϕ d θ d r = ∫ 0 1 r 8 d r ∫ 0 π sin 5 θ cos 2 θ d θ ∫ 0 2 π cos 2 ϕ sin 2 ϕ d ϕ = 1 9 ∫ 0 π sin 5 θ cos 2 θ d θ ∫ 0 2 π 1 4 sin 2 2 ϕ d ϕ = 1 9 ∫ 0 π sin 5 θ cos 2 θ d θ 1 4 ∫ 0 2 π 1 − cos 4 ϕ 2 d ϕ = 1 9 ∫ 0 π sin 5 θ cos 2 θ d θ 1 4 1 2 2 π = π 36 ∫ 0 π sin 5 θ cos 2 θ d θ \iiint_B g(x,y,z) dx dy dz = \iiint_B x^2y^2z^2 dx dy dz = \int_0^1 \int_0^{\pi} \int_0^{2\pi} (r \sin \theta \cos \phi)^2(r \sin \theta \sin \phi)^2(r \cos \theta)^2 r^2 \sin \theta d\phi d\theta dr = \int_0^1 r^8 dr \int_0^{\pi} \sin^5 \theta \cos^2 \theta d\theta \int_0^{2\pi} \cos^2 \phi \sin^2 \phi d\phi = \frac{1}{9} \int_0^{\pi} \sin^5 \theta \cos^2 \theta d\theta \int_0^{2\pi} \frac{1}{4} \sin^2 2\phi d\phi = \frac{1}{9} \int_0^{\pi} \sin^5 \theta \cos^2 \theta d\theta \frac{1}{4} \int_0^{2\pi} \frac{1 - \cos 4\phi}{2} d\phi = \frac{1}{9} \int_0^{\pi} \sin^5 \theta \cos^2 \theta d\theta \frac{1}{4} \frac{1}{2} 2\pi = \frac{\pi}{36} \int_0^{\pi} \sin^5 \theta \cos^2 \theta d\theta ∭ B g ( x , y , z ) d x d y d z = ∭ B x 2 y 2 z 2 d x d y d z = ∫ 0 1 ∫ 0 π ∫ 0 2 π ( r sin θ cos ϕ ) 2 ( r sin θ sin ϕ ) 2 ( r cos θ ) 2 r 2 sin θ d ϕ d θ d r = ∫ 0 1 r 8 d r ∫ 0 π sin 5 θ cos 2 θ d θ ∫ 0 2 π cos 2 ϕ sin 2 ϕ d ϕ = 9 1 ∫ 0 π sin 5 θ cos 2 θ d θ ∫ 0 2 π 4 1 sin 2 2 ϕ d ϕ = 9 1 ∫ 0 π sin 5 θ cos 2 θ d θ 4 1 ∫ 0 2 π 2 1 − c o s 4 ϕ d ϕ = 9 1 ∫ 0 π sin 5 θ cos 2 θ d θ 4 1 2 1 2 π = 36 π ∫ 0 π sin 5 θ cos 2 θ d θ
Using the beta integral,
∫ 0 π sin 5 θ cos 2 θ d θ = 2 ∫ 0 π / 2 sin 5 θ cos 2 θ d θ = B ( 3 , 3 2 ) = Γ ( 3 ) Γ ( 3 / 2 ) Γ ( 9 / 2 ) = 16 105 . \begin{aligned}
\int_0^{\pi}\sin^5\theta\cos^2\theta\,d\theta
&=2\int_0^{\pi/2}\sin^5\theta\cos^2\theta\,d\theta\\
&=B\left(3,\frac32\right)
=\frac{\Gamma(3)\Gamma(3/2)}{\Gamma(9/2)}
=\frac{16}{105}.
\end{aligned} ∫ 0 π sin 5 θ cos 2 θ d θ = 2 ∫ 0 π /2 sin 5 θ cos 2 θ d θ = B ( 3 , 2 3 ) = Γ ( 9/2 ) Γ ( 3 ) Γ ( 3/2 ) = 105 16 .
Therefore,
∭ B x 2 y 2 z 2 d x d y d z = π 36 16 105 = 4 π 945 . \boxed{
\iiint_Bx^2y^2z^2\,dx\,dy\,dz
=\frac{\pi}{36}\frac{16}{105}
=\frac{4\pi}{945}
}. ∭ B x 2 y 2 z 2 d x d y d z = 36 π 105 16 = 945 4 π .