Uloha je spocitat integral 1 na mnozine M = {(x,y,z) | (x^2/4 + y^2/9 + z^2/16)^3 <= z*x^2 }
Nejprve se pouziji kulove souradnice(nebo tak nejak se jmenuji :) )
x = 2r sina sinb
y = 3r sina cosb
z = 4r cosa
Koeficienty 2,3,4 jsou tam proto, aby se zkratily ty zlomky...
0 < r < oo
0 < a < PI
0 < b < 2PI
Jakobian = -24r^2 sinb
Podle vety o substituci:
Int f na mnozine A = Int (f o g)*|Jakobian g| na mnozine g^-1(A)
poznamky:
a) g musi byt regularni a prosta.
b) g^-1 tim minim inverzni funkci
V nasem pripade je g rovna funkci (r,b,a) -> (r sina sinb, r sina cosb, r cosa), je regularni a prosta.
(x^2/4 + y^2/9 + z^2/16)^3 <= z*x^2
Dosadim substituci za x,y,z ->
r^6 <= 16r^3 sin^2 (a) sin^2(b) cos(a)
r^3 <= 16 sin^2 (a) sin^2(b) cos(a)
Vzhledem k tomu ze r musi byt > 0, tak i cos(a) musi byt >0 => 0 < a < pi/2
Takze nase nova uloha zni:
Int (|-24*r^2 sina|) na mnozine {(r,b,a) | r^3 <= 16 sin^2 (a) sin^2(b) cos(a) }
0<r<= treti_odmocnina(16 sin^2 (a) sin^2(b) cos(a))
0<b<2PI
0<a<PI/2
Ted se veme Fubiniho veta:
PI/2 2PI treti_odmocnina(16 sin^2 (a) sin^2(b) cos(a))
Int Int Int (24r^2 sina) dr db da
0 0 0
Mno to uz preci umime spocitat :) Mathematica mi to spocetla jako 32PI.
Myslim ze az na pripadne preklepy by to melo byt sprave.