Kaksinkertaista integraalia voidaan käyttää kolmiulotteisen kappaleen tilavuuden määrittämiseen. Kuvan esimerkissä pinnan
f
(
x
,
y
)
=
10
−
x
2
−
y
2
8
{\textstyle f(x,y)=10-{\frac {x^{2}-y^{2}}{8}}}
alle jäävä tilavuus saadaan integroimalla funktio
f
{\textstyle f}
suorakaiteen muotoisen pohjan pinta-alan yli.
Integrointi moniulotteisessa avaruudessa tarkoittaa kahden tai useamman muuttujan reaaliarvoisten funktioiden määrätyn integraalin selvittämistä. Moniulotteisen avaruuden määrättyjä integraaleja kutsutaan moninkertaisiksi integraaleiksi , ja niiden määrittäminen perustuu yhden muuttujan määrättyjen integraalien määrittämiseen toistuvasti.[ 1]
Kaksiulotteisen avaruuden
R
2
{\textstyle \mathbb {R} ^{2}}
osajoukoissa määriteltyjen funktioiden
f
(
x
,
y
)
{\textstyle f(x,y)}
integraaleja kutsutaan kaksinkertaisiksi integraaleiksi [ 1] ja kolmiulotteisen avaruuden
R
3
{\textstyle \mathbb {R} ^{3}}
osajoukoissa määriteltyjen funktioiden
f
(
x
,
y
,
z
)
{\textstyle f(x,y,z)}
integraaleja vastaavasti kolminkertaisiksi integraaleiksi .[ 2] Yksinkertaisin moniulotteisessa avaruudessa integroinnin sovellus on kolmiulotteisen alueen tilavuuden määrittäminen kaksinkertaisen integraalin avulla.[ 1]
Tavanomaista, välillä
[
a
,
b
]
⊂
R
{\textstyle [a,b]\subset \mathbb {R} }
määritellyn, yhden muuttujan funktion määrättyä integraalia merkitään yleisesti:
∫
a
b
f
(
x
)
d
x
{\displaystyle \int _{a}^{b}f(x)\,\mathrm {d} x}
.
Moniulotteisten avaruuksien integraaleissa merkinnät ovat samankaltaisia. Toisaalta eri lähteissä moninkertaisia integraaleja merkitään eri tavoin. Kahden muuttujan funktion
f
(
x
,
y
)
{\textstyle f(x,y)}
integraalia yli tason
D
⊂
R
2
{\textstyle D\subset \mathbb {R} ^{2}}
(merkinnästä
R
2
{\textstyle \mathbb {R} ^{2}}
ks. karteesinen tulo ) merkitään kaksinkertaisella integraalimerkillä:
∬
D
f
(
x
,
y
)
d
A
{\displaystyle \iint _{D}f(x,y)\,\mathrm {d} A}
.[ 1]
Vastaavasti kolmen muuttujan funktion
f
(
x
,
y
,
z
)
{\textstyle f(x,y,z)}
integraalia yli avaruuden
R
⊂
R
3
{\textstyle R\subset \mathbb {R} ^{3}}
merkitään kolminkertaisella integraalimerkillä:
∭
R
f
(
x
,
y
,
z
)
d
V
{\displaystyle \iiint _{R}f(x,y,z)\,\mathrm {d} V}
.[ 1]
Useamman kuin kolmen muuttujan funktion integraaleille toistuvia integraalimerkkejä ei enää tilan säästämisen vuoksi käytetä (vrt. moninkertaisten derivaattojen merkintä ). Sen sijaan avaruudessa
A
⊂
R
n
{\textstyle A\subset \mathbb {R} ^{n}}
, missä
n
∈
N
{\textstyle n\in \mathbb {N} }
,
n
≥
3
{\textstyle n\geq 3}
, määritellyn funktion
f
(
x
1
,
x
2
,
…
,
x
n
)
{\textstyle f(x_{1},x_{2},\dotsc ,x_{n})}
integraalissa käytetään eri lyhennysmerkintöjä:
∫
…
∫
A
f
(
x
1
,
x
2
,
…
,
x
n
)
d
V
{\displaystyle \int \dots \int _{A}f(x_{1},x_{2},\dotsc ,x_{n})\,\mathrm {d} V}
,[ 3]
∫
…
∫
A
f
(
x
1
,
x
2
,
…
,
x
n
)
d
x
1
…
d
x
n
{\displaystyle \int \dots \int _{A}f(x_{1},x_{2},\dotsc ,x_{n})\,\mathrm {d} x_{1}\dotsc \mathrm {d} x_{n}}
.[ 3]
Joskus kaksin-, kolmin- tai useammankertaista integraalia merkitään myös yksinkertaisella integraalimerkillä (jolloin integraalissa esiintyvän avaruuden ulottuvuus on selvitettävä kontekstista ):
∫
A
f
(
x
1
,
…
,
x
n
)
d
x
1
…
d
x
n
{\displaystyle \int _{A}f(x_{1},\dotsc ,x_{n})\,\mathrm {d} x_{1}\dotsc \mathrm {d} x_{n}}
,
∫
A
f
(
x
)
d
V
{\displaystyle \int _{A}f(\mathbb {x} )\,\mathrm {d} V}
tai
∫
A
f
(
x
)
d
x
{\displaystyle \int _{A}f(\mathbb {x} )\,\mathrm {d} \mathbb {x} }
, [ 3] [ 4]
missä kahdessa viimeisessä käytetään merkintää
f
(
x
)
=
f
(
x
1
,
x
2
,
…
,
x
n
)
{\textstyle f(\mathbb {x} )=f(x_{1},x_{2},\dotsc ,x_{n})}
. Kummallakin merkintätavalla on omat hyvät ja huonot puolensa, ja kumpaankin käytetään yleisesti.
Kaksinkertaisen integraalin määritelmä esitetään (yksinkertaisuutensa ja helpommin hahmotettavuutensa takia) siten, kuin se on esitetty lähteessä [ 1] . Useamman muuttujan funktioiden integrointi määritellään yleisemmin, kuten se on tehty lähteessä [ 4] . Nämä kaksi määritelmää eivät ole ristiriidassa keskenään.
Suorakulmion
D
{\textstyle D}
(tummennettu alue) ositus pienemmiksi suorakulmioiksi
R
i
j
{\textstyle R_{ij}}
, missä
1
≤
i
≤
m
{\textstyle 1\leq i\leq m}
ja
1
≤
j
≤
n
{\textstyle 1\leq j\leq n}
. Olkoon
D
⊂
R
2
{\textstyle D\subset \mathbb {R} ^{2}}
suljettu suorakulmio , jonka sivut ovat koordinaattiakselien suuntaiset ja
f
:
D
→
R
{\textstyle f:D\to \mathbb {R} }
rajoitettu kuvaus . Jos
D
{\textstyle D}
:hen kuuluvat kaikki pisteet
(
x
,
y
)
{\textstyle (x,y)}
siten, että
x
∈
[
a
,
b
]
{\textstyle x\in [a,b]}
ja
y
∈
[
c
,
d
]
{\textstyle y\in [c,d]}
, niin
D
{\textstyle D}
:hen voidaan määritellä ositus
P
{\textstyle P}
siten, että
a
=
x
0
<
x
1
<
x
2
<
⋯
<
x
m
−
1
<
x
m
=
b
c
=
y
0
<
y
1
<
y
2
<
⋯
<
y
n
−
1
<
y
n
=
d
.
{\displaystyle {\begin{aligned}a&=x_{0}<x_{1}<x_{2}<\dots <x_{m-1}<x_{m}=b\\c&=y_{0}<y_{1}<y_{2}<\dots <y_{n-1}<y_{n}=d.\end{aligned}}}
Ositus
P
{\textstyle P}
koostuu
m
⋅
n
{\textstyle m\cdot n}
:stä suorakulmiosta
R
i
j
{\textstyle R_{ij}}
, missä
1
≤
i
≤
m
{\textstyle 1\leq i\leq m}
ja
1
≤
j
≤
n
{\textstyle 1\leq j\leq n}
, jotka edelleen koostuvat pisteistä
(
x
,
y
)
{\textstyle (x,y)}
, missä
x
∈
[
x
i
−
1
,
x
i
]
{\textstyle x\in [x_{i-1},x_{i}]}
ja
y
∈
[
y
j
−
1
,
y
j
]
{\textstyle y\in [y_{j-1},y_{j}]}
. Suorakulmion
R
i
j
{\textstyle R_{ij}}
pinta-ala on
Δ
A
i
j
=
Δ
x
i
Δ
y
j
=
(
x
i
−
x
i
−
1
)
(
y
j
−
y
j
−
1
)
{\displaystyle \Delta A_{ij}=\Delta x_{i}\Delta y_{j}=(x_{i}-x_{i-1})(y_{j}-y_{j-1})}
ja sen lävistäjän pituus on Pythagoraan lauseen nojalla
d
i
a
m
(
R
i
j
)
=
(
Δ
x
i
)
2
+
(
Δ
y
j
)
2
=
(
x
i
−
x
i
−
1
)
2
+
(
y
j
−
y
j
−
1
)
2
{\displaystyle \mathrm {diam} (R_{ij})={\sqrt {(\Delta x_{i})^{2}+(\Delta y_{j})^{2}}}={\sqrt {(x_{i}-x_{i-1})^{2}+(y_{j}-y_{j-1})^{2}}}}
.
Määritellään osituksen
P
{\textstyle P}
suurimman suorakulmion lävistäjän pituus osituksen normiksi :
‖
P
‖
=
max
1
≤
i
≤
m
,
1
≤
j
≤
n
d
i
a
m
(
R
i
j
)
{\displaystyle \lVert P\rVert =\max _{1\leq i\leq m,~1\leq j\leq n}\mathrm {diam} (R_{ij})}
.
Valitaan jokaisesta suorakulmiosta
R
i
j
{\textstyle R_{ij}}
mielivaltainen piste
(
x
i
j
∗
,
y
i
j
∗
)
{\textstyle (x_{ij}^{*},y_{ij}^{*})}
ja muodostetaan niiden avulla kaksiulotteinen Riemannin summa :
S
P
(
f
)
=
∑
i
=
1
m
∑
j
=
1
n
f
(
x
i
j
∗
,
y
i
j
∗
)
Δ
A
i
j
{\displaystyle S_{P}(f)=\sum _{i=1}^{m}\sum _{j=1}^{n}f(x_{ij}^{*},y_{ij}^{*})\Delta A_{ij}}
.
Suorakulmio
R
i
j
{\textstyle R_{ij}}
muodostaa suorakulmaisen laatikon pohjan. Laatikon korkeus on
z
=
f
(
x
i
j
∗
,
y
i
j
∗
)
{\textstyle z=f(x_{ij}^{*},y_{ij}^{*})}
. Riemannin summa on kaikkien näiden laatikoiden tilavuuksien summa.
Riemannin summan jokainen termi kuvaa sellaisen suorakulmaisen laatikon tilavuutta, jonka pohja on suorakulmio
R
i
j
{\textstyle R_{ij}}
ja korkeus on funktion
f
{\textstyle f}
arvo pisteessä
(
x
i
j
∗
,
y
i
j
∗
)
{\textstyle (x_{ij}^{*},y_{ij}^{*})}
(siis mikäli
f
(
x
i
j
∗
,
y
i
j
∗
)
≥
0
{\textstyle f(x_{ij}^{*},y_{ij}^{*})\geq 0}
). Näin ollen positiivisille funktioille Riemannin summa approksimoi suorakulmion
D
{\textstyle D}
ja funktion
f
{\textstyle f}
kuvaajan väliin jäävän avaruuden tilavuutta .[ 1] Kun laatikoiden lukumäärää kasvatetaan ja vastaavasti pinta-aloja
Δ
A
i
j
{\textstyle \Delta A_{ij}}
pienennetään (jolloin jaon normi
‖
P
‖
{\textstyle \lVert P\rVert }
pienee), muuttuu approksimaatio yhä tarkemmaksi. Kaksinkertaisen integraalin määritelmä perustuukin siihen, että
‖
P
‖
→
0
{\textstyle \lVert P\rVert \to 0}
.[ 1] Lopullinen määritelmä kuuluu:
Funktio
f
:
D
→
R
{\textstyle f:D\to \mathbb {R} }
on integroituva suorakulmion
D
{\textstyle D}
yli ja sillä on kaksinkertainen integraali
I
=:
∬
D
f
(
x
,
y
)
d
A
{\displaystyle I=:\iint _{D}f(x,y)\,\mathrm {d} A}
,
jos kaikilla
ε
>
0
{\textstyle \varepsilon >0}
on olemassa (
ε
{\textstyle \varepsilon }
:sta riippuva) luku
δ
>
0
{\textstyle \delta >0}
siten, että ehto
|
S
P
(
f
)
−
I
|
<
ε
{\displaystyle \left|S_{P}(f)-I\right|<\varepsilon }
pätee kaikille
D
{\textstyle D}
:n osituksille
P
{\textstyle P}
, joille
‖
P
‖
<
δ
{\textstyle \lVert P\rVert <\delta }
sekä kaikkien osajoukkojen
R
i
j
{\textstyle R_{ij}}
kaikissa pisteissä
(
x
i
j
∗
,
y
i
j
∗
)
{\textstyle (x_{ij}^{*},y_{ij}^{*})}
.[ 1]
Merkintä
d
A
{\textstyle \mathrm {d} A}
tarkoittaa differentiaalista pinta-alaelementtiä, joka on eräs esitystapa pinta-alan
Δ
A
=
Δ
x
Δ
y
{\textstyle \Delta A=\Delta x\Delta y}
raja-arvolle . Pinta-alaelementti voidaan kirjoittaa muodossa
d
A
=
d
x
d
y
=
d
y
d
x
{\displaystyle \mathrm {d} A=\mathrm {d} x\,\mathrm {d} y=\mathrm {d} y\,\mathrm {d} x}
,
mistä on erityistä hyötyä myöhemmin, kun integraaleja lasketaan iteroimalla.[ 1]
Olkoon
n
∈
N
{\textstyle n\in \mathbb {N} }
,
n
≥
3
{\textstyle n\geq 3}
. Joukko
I
=
[
a
1
,
b
1
]
×
[
a
2
,
b
2
]
×
…
×
[
a
n
,
b
n
]
{\textstyle I=[a_{1},b_{1}]\times [a_{2},b_{2}]\times \dotsc \times [a_{n},b_{n}]}
on avaruuden
R
n
{\textstyle \mathbb {R} ^{n}}
kompakti väli [ 4] (nk. hypersuorakulmio ). Olkoon lisäksi
f
:
I
→
R
{\textstyle f:I\to \mathbb {R} }
rajoitettu funktio. Jos
1
≤
i
≤
n
{\textstyle 1\leq i\leq n}
, niin välin
[
a
i
,
b
i
]
{\textstyle [a_{i},b_{i}]}
jako
P
i
{\textstyle P_{i}}
on:
a
i
=
t
0
≤
t
1
≤
…
≤
t
k
i
=
b
i
{\displaystyle a_{i}=t_{0}\leq t_{1}\leq \dotsc \leq t_{k_{i}}=b_{i}}
.
Ts. väli
[
a
i
,
b
i
]
{\textstyle [a_{i},b_{i}]}
voidaan esittää
k
i
{\textstyle k_{i}}
:n eri osavälin (joista osa voi olla myös surkastunut pisteeksi )
[
t
i
−
1
,
t
i
]
{\textstyle [t_{i-1},t_{i}]}
unioinina . Jos jokaiselle
1
≤
i
≤
n
{\textstyle 1\leq i\leq n}
on olemassa välin
[
a
i
,
b
i
]
{\textstyle [a_{i},b_{i}]}
jako
P
i
{\textstyle P_{i}}
, niin joukon
I
{\textstyle I}
jako on:
P
=
P
1
×
P
2
×
…
×
P
n
{\displaystyle P=P_{1}\times P_{2}\times \dotsc \times P_{n}}
.
Jako
P
{\textstyle P}
määrää
κ
=
k
1
⋅
k
2
⋅
…
⋅
k
n
{\textstyle \kappa =k_{1}\cdot k_{2}\cdot \dotsc \cdot k_{n}}
kappaletta kompakteja jakovälejä
I
1
,
I
2
,
…
,
I
κ
{\textstyle I_{1},I_{2},\dots ,I_{\kappa }}
, jotka osittavat välin
I
{\textstyle I}
. Merkitään mielivaltaista joukon
I
{\textstyle I}
pistettä
(
x
1
,
x
2
,
…
,
x
n
)
=
x
{\textstyle (x_{1},x_{2},\dotsc ,x_{n})=\mathbb {x} }
sekä kaikilla
j
∈
{
1
,
2
,
…
,
κ
}
{\textstyle j\in \{1,2,\dotsc ,\kappa \}}
:
g
j
=
inf
x
∈
I
j
f
(
x
)
{\displaystyle g_{j}=\inf _{\mathbb {x} \in I_{j}}f(\mathbb {x} )}
ja
G
j
=
sup
x
∈
I
j
f
(
x
)
{\displaystyle G_{j}=\sup _{\mathbb {x} \in I_{j}}f(\mathbb {x} )}
(ks. infimum ja supremum ). Määritellään näiden avulla jakoon
P
{\textstyle P}
liittyvä alasumma kaavalla
s
P
(
f
)
=
∑
j
=
1
κ
[
g
j
⋅
∏
i
=
1
n
(
b
i
−
a
i
)
]
=
∑
j
=
1
κ
[
g
j
⋅
(
b
1
−
a
1
)
⋅
(
b
2
−
a
2
)
⋅
…
⋅
(
b
n
−
a
n
)
]
{\displaystyle s_{P}(f)=\sum _{j=1}^{\kappa }\left[g_{j}\cdot \prod _{i=1}^{n}(b_{i}-a_{i})\right]=\sum _{j=1}^{\kappa }\left[g_{j}\cdot (b_{1}-a_{1})\cdot (b_{2}-a_{2})\cdot \dotsc \cdot (b_{n}-a_{n})\right]}
ja vastaavasti jakoon
P
{\textstyle P}
liittyvä yläsumma kaavalla
S
P
(
f
)
=
∑
j
=
1
κ
[
G
j
⋅
∏
i
=
1
n
(
b
i
−
a
i
)
]
=
∑
j
=
1
κ
[
G
j
⋅
(
b
1
−
a
1
)
⋅
(
b
2
−
a
2
)
⋅
…
⋅
(
b
n
−
a
n
)
]
{\displaystyle S_{P}(f)=\sum _{j=1}^{\kappa }\left[G_{j}\cdot \prod _{i=1}^{n}(b_{i}-a_{i})\right]=\sum _{j=1}^{\kappa }\left[G_{j}\cdot (b_{1}-a_{1})\cdot (b_{2}-a_{2})\cdot \dotsc \cdot (b_{n}-a_{n})\right]}
.
Joukolle
I
{\textstyle I}
voidaan asettaa useita eri ala- ja yläsummia riippuen siitä, millä tavalla (ja kuinka monilla osaväleillä) jako
P
{\textstyle P}
tehdään. Jos jakopisteitä lisätään, ei yläsumma kasva eikä alasumma pienene.[ 4] Näin ollen, jos
P
{\textstyle P}
on välin
I
{\textstyle I}
jako,
P
′
{\textstyle P'}
on
P
{\textstyle P}
:n alajako (jako, joka sisältää (vähintään) kaikki jaon
P
{\textstyle P}
jakopisteet) ja
P
″
{\textstyle P''}
on kummankin edellisen jaon alajako, niin:
s
P
(
f
)
≤
s
P
″
(
f
)
≤
S
P
″
(
f
)
≤
S
P
′
(
f
)
{\displaystyle s_{P}(f)\leq s_{P''}(f)\leq S_{P''}(f)\leq S_{P'}(f)}
[ 4]
Tästä seuraa edelleen, että alasummien joukon supremum on korkeintaan yläsummien joukon infimum. Näin ollen, jos
P
{\textstyle P}
on joukon
I
{\textstyle I}
mielivaltainen jako, niin aina
sup
P
s
P
(
f
)
≤
inf
P
S
P
(
f
)
{\displaystyle \sup _{P}s_{P}(f)\leq \inf _{P}S_{P}(f)}
.
Tällöin funktio
f
:
I
→
R
{\textstyle f:I\to \mathbb {R} }
on (Riemann- )integroituva (joukossa
I
{\textstyle I}
tai joukon
I
{\textstyle I}
yli), jos
sup
P
s
P
(
f
)
=
inf
P
S
P
(
f
)
{\displaystyle \sup _{P}s_{P}(f)=\inf _{P}S_{P}(f)}
.[ 4]
Tämä määritelmä yhtyy kaksinkertaisen integraalin määritelmään seuraavasti: Olkoon
I
=
[
a
1
,
b
1
]
×
[
a
2
,
b
2
]
×
…
×
[
a
n
,
b
n
]
⊂
R
n
{\textstyle I=[a_{1},b_{1}]\times [a_{2},b_{2}]\times \dotsc \times [a_{n},b_{n}]\subset \mathbb {R} ^{n}}
,
n
≥
2
{\textstyle n\geq 2}
ja
f
:
I
→
R
{\textstyle f:I\to \mathbb {R} }
on rajoitettu funktio. Tällöin funktio
f
{\textstyle f}
on integroituva, jos ja vain, jos kaikilla
ε
>
0
{\textstyle \varepsilon >0}
on olemassa välin
I
{\textstyle I}
jako
P
{\textstyle P}
siten, että
S
P
(
f
)
−
s
P
(
f
)
<
ε
{\displaystyle S_{P}(f)-s_{P}(f)<\varepsilon }
.[ 4]
Integrointialue
A
{\textstyle A}
on suorakulmion
I
{\textstyle I}
osajoukko .
Yleisesti ottaen joukko, jonka yli integrointi suoritetaan, ei välttämättä ole suorakulmio (tai hypersuorakulmio). Tällöin integrointialuetta voidaan yksinkertaistaa laajentamalla sitä siten, että uusi integrointialue on sellainen (hyper-)suorakulmio, joka sisältää alkuperäisen integrointialueen sekä määrittelemällä integroitava funktio paloittain . Oletetaan nyt, että
I
=
[
a
1
,
b
1
]
×
[
a
2
,
b
2
]
×
…
×
[
a
n
,
b
n
]
⊂
R
n
{\textstyle I=[a_{1},b_{1}]\times [a_{2},b_{2}]\times \dotsc \times [a_{n},b_{n}]\subset \mathbb {R} ^{n}}
on hypersuorakulmio. Olkoon lisäksi integrointialue
A
⊂
I
{\textstyle A\subset I}
ja funktio
f
:
A
→
R
{\textstyle f:A\to \mathbb {R} }
rajoitettu . Määritellään funktio
f
^
:
I
→
R
{\textstyle {\hat {f}}:I\to \mathbb {R} }
siten, että:
f
^
(
x
)
=
{
f
(
x
)
,
jos
x
∈
A
0
,
jos
x
∈
I
∖
A
.
{\displaystyle {\hat {f}}(\mathbb {x} )={\begin{cases}f(\mathbb {x} ),&{\text{ jos }}\mathbb {x} \in A\\0,&{\text{ jos }}\mathbb {x} \in I\,\backslash \,A.\end{cases}}}
Jos
f
^
{\textstyle {\hat {f}}}
on integroituva yli
I
{\textstyle I}
:n, niin myös funktio
f
{\textstyle f}
on integroituva yli
A
{\textstyle A}
:n ja:
∫
A
f
(
x
)
d
V
=
∫
I
f
^
(
x
)
d
V
{\displaystyle \int _{A}f(\mathbb {x} )\,\mathrm {d} V=\int _{I}{\hat {f}}(\mathbb {x} )\,\mathrm {d} V}
.[ 1] [ 4]
Integraalin arvo ei riipu joukon
I
{\textstyle I}
valinnasta, kunhan
A
⊂
I
{\textstyle A\subset I}
.[ 4]
Olkoon seuraavassa
n
∈
N
{\textstyle n\in \mathbb {N} }
,
A
⊂
R
n
{\textstyle A\subset \mathbb {R} ^{n}}
,
x
=
(
x
1
,
x
2
,
…
,
x
n
)
{\textstyle \mathbb {x} =(x_{1},x_{2},\dotsc ,x_{n})}
, funktiot
f
{\textstyle f}
ja
g
{\textstyle g}
integroituvia funktioita joukossa
A
{\textstyle A}
sekä
a
∈
R
{\textstyle a\in \mathbb {R} }
ja
b
∈
R
{\textstyle b\in \mathbb {R} }
vakioita . Seuraavat ominaisuudet pätevät moniulotteisille integraaleille:
∫
A
f
(
x
)
d
V
=
∫
i
n
t
(
A
)
f
(
x
)
d
V
{\displaystyle \int _{A}f(\mathbb {x} )\,\mathrm {d} V=\int _{\mathrm {int} (A)}f(\mathbb {x} )\,\mathrm {d} V}
(merkintä
i
n
t
(
A
)
{\textstyle \mathrm {int} (A)}
tarkoittaa joukon
A
{\textstyle A}
sisäpisteiden joukkoa).[ 1]
∫
A
1
d
V
{\displaystyle \int _{A}1\,\mathrm {d} V}
= integrointialueen
A
{\textstyle A}
tilavuus (pinta-ala , jos
n
=
2
{\textstyle n=2}
ja pituus , jos
n
=
1
{\textstyle n=1}
).[ 1]
∫
A
f
(
x
)
d
V
=
0
{\displaystyle \int _{A}f(\mathbb {x} )\,\mathrm {d} V=0}
, jos
i
n
t
(
A
)
=
∅
{\textstyle \mathrm {int} (A)=\varnothing }
(joukolla
A
{\textstyle A}
ei ole sisäpisteitä ).[ 1]
Funktio
a
⋅
f
+
b
⋅
g
{\textstyle a\cdot f+b\cdot g}
on integroituva joukossa
A
{\textstyle A}
ja
∫
A
(
a
f
(
x
)
+
b
g
(
x
)
)
d
V
=
a
∫
A
f
(
x
)
d
V
+
b
∫
A
g
(
x
)
d
V
{\displaystyle \int _{A}\left(af(\mathbb {x} )+bg(\mathbb {x} )\right)\,\mathrm {d} V=a\int _{A}f(\mathbb {x} )\,\mathrm {d} V+b\int _{A}g(\mathbb {x} )\,\mathrm {d} V}
.[ 1] [ 4]
Jos
f
(
x
)
≤
g
(
x
)
{\textstyle f(\mathbb {x} )\leq g(\mathbb {x} )}
kaikilla
x
∈
A
{\textstyle \mathbb {x} \in A}
, niin
∫
A
f
(
x
)
d
V
≤
∫
A
g
(
x
)
d
V
{\displaystyle \int _{A}f(\mathbb {x} )\,\mathrm {d} V\leq \int _{A}g(\mathbb {x} )\,\mathrm {d} V}
.[ 1] [ 4]
Kolmioepäyhtälö :
|
∫
A
f
(
x
)
d
V
|
≤
∫
A
|
f
(
x
)
|
d
V
{\displaystyle \left|\int _{A}f(\mathbb {x} )\,\mathrm {d} V\right|\leq \int _{A}\left|f(\mathbb {x} )\right|\,\mathrm {d} V}
.[ 1] [ 4]
Jos joukot
A
1
,
A
2
,
…
,
A
k
{\textstyle A_{1},A_{2},\dotsc ,A_{k}}
osittavat joukon
A
{\textstyle A}
, niin
∫
A
f
(
x
)
d
V
=
∑
i
=
1
k
∫
A
i
f
(
x
)
d
V
{\displaystyle \int _{A}f(\mathbb {x} )\,\mathrm {d} V=\sum _{i=1}^{k}\int _{A_{i}}f(\mathbb {x} )\,\mathrm {d} V}
.[ 1]
Kaksin- ja kolminkertaisia integraaleja voidaan käyttää erityisesti pinta-alojen ja tilavuuksien määrittämiseen. Tätä ominaisuutta voidaan käyttää myös toisin päin: jos tiedetään integrointialueen pinta-ala tai tilavuus, voidaan integrointia helpottaa tai jopa jättää kokonaan pois. Vastaavasti voidaan käyttää integroitavan funktion symmetriaominaisuuksia, kuten parillisuutta ja parittomuutta .
Olkoon integrointialue
D
⊂
R
2
{\textstyle D\subset \mathbb {R} ^{2}}
suorakulmio
a
≤
x
≤
b
{\textstyle a\leq x\leq b}
,
c
≤
y
≤
d
{\textstyle c\leq y\leq d}
. Tällöin
∬
D
3
d
A
=
3
∬
D
d
A
=
3
⋅
(
D
:n pinta-ala
)
=
3
(
b
−
a
)
(
d
−
c
)
.
{\displaystyle {\begin{aligned}\iint _{D}3\,\mathrm {d} A&=3\iint _{D}\mathrm {d} A=3\cdot (D{\text{:n pinta-ala}})\\&=3(b-a)(d-c).\end{aligned}}}
Määritetään integraalin
∬
x
2
+
y
2
≤
1
(
sin
x
+
y
3
+
4
)
d
A
{\textstyle \iint _{x^{2}+y^{2}\leq 1}(\sin x+y^{3}+4)\,\mathrm {d} A}
arvo. Ensimmäiseksi havaitaan, että integrointialue on
x
y
{\textstyle xy}
-tason origokeskinen yksikkökiekko (symmetrinen origon suhteen). Integraali voidaan jakaa kolmeksi eri integraaliksi:
∬
x
2
+
y
2
≤
1
(
sin
x
+
y
3
+
4
)
d
A
=
∬
x
2
+
y
2
≤
1
sin
x
d
A
+
∬
x
2
+
y
2
≤
1
y
3
d
A
+
∬
x
2
+
y
2
≤
1
4
d
A
=:
I
1
+
I
2
+
I
3
{\displaystyle {\begin{aligned}\iint _{x^{2}+y^{2}\leq 1}(\sin x+y^{3}+4)\,\mathrm {d} A&=\iint _{x^{2}+y^{2}\leq 1}\sin x\,\mathrm {d} A+\iint _{x^{2}+y^{2}\leq 1}y^{3}\,\mathrm {d} A+\iint _{x^{2}+y^{2}\leq 1}4\,\mathrm {d} A\\&=:I_{1}+I_{2}+I_{3}\end{aligned}}}
Koska integrointialue on origon suhteen symmetrinen ja funktio
f
(
x
,
y
)
=
sin
x
{\textstyle f(x,y)=\sin x}
on pariton funktio , niin kuvaajan ja
x
y
{\textstyle xy}
-tason väliin jäävä tilavuus on yhtä suuri alueissa
x
<
0
{\textstyle x<0}
ja
x
>
0
{\textstyle x>0}
. Parittomuudesta johtuen nämä tilavuudet kumoavat toisensa, jolloin
I
1
=
0
{\textstyle I_{1}=0}
. Vastaavasti, funktio
g
(
x
,
y
)
=
y
3
{\textstyle g(x,y)=y^{3}}
on myös pariton, joten
I
2
=
0
{\textstyle I_{2}=0}
. Näin ollen
I
1
+
I
2
+
I
3
=
I
3
=
∬
x
2
+
y
2
≤
1
4
d
A
=
4
⋅
(
ympyrän
x
2
+
y
2
≤
1
pinta-ala
)
=
4
π
.
{\displaystyle {\begin{aligned}I_{1}+I_{2}+I_{3}&=I_{3}=\iint _{x^{2}+y^{2}\leq 1}4\,\mathrm {d} A\\&=4\cdot ({\text{ympyrän }}x^{2}+y^{2}\leq 1{\text{ pinta-ala}})\\&=4\pi .\end{aligned}}}
Alue
D
{\textstyle D}
on säännöllinen
y
{\textstyle y}
-suunnassa, muttei
x
{\textstyle x}
-suunnassa. Alue
E
{\textstyle E}
on säännöllinen
x
{\textstyle x}
-suunnassa, muttei
y
{\textstyle y}
-suunnassa.
Käytetään tässä yksinkertaisuuden vuoksi esimerkkinä integrointia kaksiulotteisessa avaruudessa. Sanotaan, että integrointialue
D
⊂
R
2
{\textstyle D\subset \mathbb {R} ^{2}}
on
säännöllinen
y
{\textstyle y}
-suunnassa , jos sitä rajoittavat suorat
x
=
a
{\textstyle x=a}
ja
x
=
b
{\textstyle x=b}
, sekä mikä tahansa
y
{\textstyle y}
-akselin suuntainen suora leikkaa alueen
D
{\textstyle D}
reunan korkeintaan kahdesti.[ 5]
säännöllinen
x
{\textstyle x}
-suunnassa , jos sitä rajoittavat suorat
y
=
c
{\textstyle y=c}
ja
y
=
d
{\textstyle y=d}
, sekä mikä tahansa
x
{\textstyle x}
-akselin suuntainen suora leikkaa alueen
D
{\textstyle D}
reunan korkeintaan kahdesti.[ 5]
Oletetaan esimerkiksi, että integrointialue
D
{\textstyle D}
on säännöllinen
y
{\textstyle y}
-suunnassa ja sitä rajoittavat suorat
x
=
a
{\textstyle x=a}
ja
x
=
b
{\textstyle x=b}
sekä käyrät
y
=
α
(
x
)
{\textstyle y=\alpha (x)}
ja
y
=
β
(
x
)
{\textstyle y=\beta (x)}
siten, että
α
(
x
)
≤
β
(
x
)
{\textstyle \alpha (x)\leq \beta (x)}
kaikilla
x
∈
[
a
,
b
]
{\textstyle x\in [a,b]}
. Funktion
z
=
f
(
x
,
y
)
{\textstyle z=f(x,y)}
kuvaajan ja
x
y
{\textstyle xy}
-tason väliin jäävä avaruus voidaan ''viipaloida''
y
z
{\textstyle yz}
-tason suuntaisilla tasoilla (joissa siis
x
=
vakio
{\textstyle x={\text{vakio}}}
) pitkin
x
{\textstyle x}
-akselia. Tällaisen tason pinta-ala saadaan yksiulotteisella määrätyllä integraalilla:
Integroitaessa
y
{\textstyle y}
-suunnassa säännöllisen alueen yli funktion
z
=
f
(
x
,
y
)
{\textstyle z=f(x,y)}
kuvaajan ja
x
y
{\textstyle xy}
-tason väliin jäävä alue ''viipaloidaan''
x
{\textstyle x}
-akselia vastaan kohtisuorasti.
A
(
x
)
=
∫
α
(
x
)
β
(
x
)
f
(
x
,
y
)
d
y
{\displaystyle A(x)=\int _{\alpha (x)}^{\beta (x)}f(x,y)\,\mathrm {d} y}
.[ 5]
Tällöin kaksinertainen integraali
∬
D
f
(
x
,
y
)
d
A
{\textstyle \iint _{D}f(x,y)\,\mathrm {d} A}
saadaan summaamalla differentiaalisen paksujen viipaleiden (pinta-ala
A
(
x
)
{\textstyle A(x)}
ja paksuus
d
x
{\textstyle \mathrm {d} x}
) tilavuudet suorien
x
=
a
{\textstyle x=a}
ja
x
=
b
{\textstyle x=b}
välillä:
∬
D
f
(
x
,
y
)
d
A
=
∫
a
b
A
(
x
)
d
x
=
∫
a
b
(
∫
α
(
x
)
β
(
x
)
f
(
x
,
y
)
d
y
)
d
x
{\displaystyle \iint _{D}f(x,y)\,\mathrm {d} A=\int _{a}^{b}A(x)\,\mathrm {d} x=\int _{a}^{b}\left(\int _{\alpha (x)}^{\beta (x)}f(x,y)\,\mathrm {d} y\right)\,\mathrm {d} x}
.[ 5]
Samaa integraalia voidaan merkitä myös:
∫
a
b
∫
α
(
x
)
β
(
x
)
f
(
x
,
y
)
d
x
d
y
{\displaystyle \int _{a}^{b}\int _{\alpha (x)}^{\beta (x)}f(x,y)\,\mathrm {d} x\,\mathrm {d} y}
tai
∫
a
b
d
x
∫
α
(
x
)
β
(
x
)
f
(
x
,
y
)
d
y
{\displaystyle \int _{a}^{b}\mathrm {d} x\int _{\alpha (x)}^{\beta (x)}f(x,y)\,\mathrm {d} y}
.[ 5]
Näin päästään kaksinkertaisten integraalien määrittämiseen iteroimalla:
Jos
D
⊂
R
2
{\displaystyle D\subset \mathbb {R} ^{2}}
on rajoitettu ja säännöllinen
y
{\textstyle y}
-suunnassa (
a
≤
x
≤
b
{\displaystyle a\leq x\leq b}
,
α
(
x
)
≤
y
≤
β
(
x
)
{\displaystyle \alpha (x)\leq y\leq \beta (x)}
) sekä
f
:
D
→
R
{\displaystyle f:D\to \mathbb {R} }
on jatkuva funktio, niin
∬
D
f
(
x
,
y
)
d
A
=
∫
a
b
d
x
∫
α
(
x
)
β
(
x
)
f
(
x
,
y
)
d
y
{\displaystyle \iint _{D}f(x,y)\,\mathrm {d} A=\int _{a}^{b}\mathrm {d} x\int _{\alpha (x)}^{\beta (x)}f(x,y)\,\mathrm {d} y}
.[ 5]
Jos
D
⊂
R
2
{\displaystyle D\subset \mathbb {R} ^{2}}
on rajoitettu ja säännöllinen
x
{\textstyle x}
-suunnassa (
c
≤
y
≤
d
{\displaystyle c\leq y\leq d}
,
γ
(
y
)
≤
x
≤
δ
(
y
)
{\displaystyle \gamma (y)\leq x\leq \delta (y)}
) sekä
f
:
D
→
R
{\displaystyle f:D\to \mathbb {R} }
on jatkuva funktio, niin
∬
D
f
(
x
,
y
)
d
A
=
∫
c
d
d
y
∫
γ
(
y
)
δ
(
y
)
f
(
x
,
y
)
d
x
{\displaystyle \iint _{D}f(x,y)\,\mathrm {d} A=\int _{c}^{d}\mathrm {d} y\int _{\gamma (y)}^{\delta (y)}f(x,y)\,\mathrm {d} x}
.[ 5]
D
⊂
R
2
{\textstyle D\subset \mathbb {R} ^{2}}
on neliö , jonka määrittelevät rajat
0
≤
x
≤
1
{\textstyle 0\leq x\leq 1}
ja
1
≤
y
≤
2
{\textstyle 1\leq y\leq 2}
. Määritetään neliön
D
{\textstyle D}
ja tason
z
=
4
−
x
−
y
{\textstyle z=4-x-y}
väliin jäävän avaruuden tilavuus. Integrointialue
D
{\textstyle D}
on säännöllinen sekä
x
{\textstyle x}
- että
y
{\textstyle y}
-suunnissa, joten iterointi voidaan suorittaa kummassa järjestyksessä tahansa.
Integroidaan ensin
x
{\textstyle x}
- ja sitten
y
{\textstyle y}
-suunnassa. Integroitaessa
x
{\textstyle x}
-suunnassa muuttujaa
y
{\textstyle y}
kohdellaan kuten vakiota . Tason ja neliön
D
{\textstyle D}
rajoittaman kappaleen tilavuus on:
V
=
∬
D
(
4
−
x
−
y
)
d
A
=
∫
y
=
1
2
d
y
∫
x
=
0
1
(
4
−
x
−
y
)
d
x
=
∫
y
=
1
2
d
y
|
x
=
0
1
(
4
x
−
1
2
x
2
−
x
y
)
=
∫
1
2
(
7
2
−
y
)
d
y
=
|
1
2
(
7
2
y
−
1
2
y
2
)
=
2
{\displaystyle {\begin{aligned}V&=\iint _{D}(4-x-y)\,\mathrm {d} A\\&=\int _{y=1}^{2}\mathrm {d} y\int _{x=0}^{1}(4-x-y)\,\mathrm {d} x\\&=\int _{y=1}^{2}\mathrm {d} y\,{\Bigg |}_{x=0}^{1}\left(4x-{\frac {1}{2}}x^{2}-xy\right)\\&=\int _{1}^{2}\left({\frac {7}{2}}-y\right)\,\mathrm {d} y\\&={\Bigg |}_{1}^{2}\left({\frac {7}{2}}y-{\frac {1}{2}}y^{2}\right)=2\end{aligned}}}
Integroidaan ensin
y
{\textstyle y}
- ja sitten
x
{\textstyle x}
-suunnassa. Integroitaessa
y
{\textstyle y}
-suunnassa muuttujaa
x
{\textstyle x}
kohdellaan kuten vakiota. Tason ja neliön
D
{\textstyle D}
rajoittaman kappaleen tilavuus on:
V
=
∬
D
(
4
−
x
−
y
)
d
A
=
∫
x
=
0
1
d
x
∫
y
=
1
2
(
4
−
x
−
y
)
d
y
=
∫
x
=
0
1
d
x
|
y
=
1
2
(
4
y
−
x
y
−
1
2
y
2
)
=
∫
0
1
(
5
2
−
x
)
d
x
=
|
0
1
(
5
2
x
−
1
2
x
2
)
=
2
{\displaystyle {\begin{aligned}V&=\iint _{D}(4-x-y)\,\mathrm {d} A\\&=\int _{x=0}^{1}\mathrm {d} x\int _{y=1}^{2}(4-x-y)\,\mathrm {d} y\\&=\int _{x=0}^{1}\mathrm {d} x\,{\Bigg |}_{y=1}^{2}\left(4y-xy-{\frac {1}{2}}y^{2}\right)\\&=\int _{0}^{1}\left({\frac {5}{2}}-x\right)\,\mathrm {d} x\\&={\Bigg |}_{0}^{1}\left({\frac {5}{2}}x-{\frac {1}{2}}x^{2}\right)=2\end{aligned}}}
Alue
D
{\textstyle D}
on yksinkertainen sekä
x
{\textstyle x}
- että
y
{\textstyle y}
-suunnissa.
Ratkaistaan itegraali
I
=
∬
D
e
y
3
d
A
{\textstyle I=\iint _{D}e^{y^{3}}\,\mathrm {d} A}
, missä aluetta
D
⊂
R
2
{\textstyle D\subset \mathbb {R} ^{2}}
rajaavat suorat
x
=
0
{\textstyle x=0}
ja
y
=
1
{\textstyle y=1}
sekä käyrä
x
=
y
2
{\textstyle x=y^{2}}
(tai
y
=
x
{\textstyle y={\sqrt {x}}}
). Alue
D
{\textstyle D}
on säännöllinen sekä
x
{\textstyle x}
- että
y
{\textstyle y}
-suunnissa. Koska funktiolle
f
:
R
→
R
{\textstyle f:\mathbb {R} \to \mathbb {R} }
,
f
(
y
)
=
e
y
3
{\textstyle f(y)=e^{y^{3}}}
ei voi kirjoittaa antiderivaattaa , täytyy iterointi suorittaa siten, että ensin integroidaan
x
{\textstyle x}
-suunnassa:
I
=
∬
D
e
y
3
d
A
=
∫
y
=
0
1
d
y
∫
x
=
0
y
2
e
y
3
d
x
=
∫
y
=
0
1
d
y
|
x
=
0
y
2
x
e
y
3
=
∫
0
1
y
2
e
y
3
d
y
=
|
0
1
1
3
e
y
3
=
e
−
1
3
{\displaystyle {\begin{aligned}I&=\iint _{D}e^{y^{3}}\,\mathrm {d} A\\&=\int _{y=0}^{1}\mathrm {d} y\int _{x=0}^{y^{2}}e^{y^{3}}\,\mathrm {d} x\\&=\int _{y=0}^{1}\mathrm {d} y\,{\Bigg |}_{x=0}^{y^{2}}xe^{y^{3}}\\&=\int _{0}^{1}y^{2}e^{y^{3}}\,\mathrm {d} y\\&={\Bigg |}_{0}^{1}{\frac {1}{3}}e^{y^{3}}={\frac {e-1}{3}}\end{aligned}}}
Olkoon
A
⊂
R
3
{\textstyle A\subset \mathbb {R} ^{3}}
suorakulmainen laatikko siten, että
0
≤
x
≤
a
{\textstyle 0\leq x\leq a}
,
0
≤
y
≤
b
{\textstyle 0\leq y\leq b}
ja
0
≤
z
≤
c
{\textstyle 0\leq z\leq c}
, missä
a
{\textstyle a}
,
b
{\textstyle b}
ja
c
{\textstyle c}
ovat positiivisia vakioita. Ratkaistaan integraali
I
=
∭
A
(
x
y
2
+
z
3
)
d
V
{\textstyle I=\iiint _{A}(xy^{2}+z^{3})\,\mathrm {d} V}
. Avaruus
A
{\textstyle A}
''viipaloidaan'' nyt (esimerkiksi)
x
y
{\textstyle xy}
-tason suuntaisilla tasoilla, jolloin ensimmäisenä integroidaan muuttujan
z
{\textstyle z}
suhteen. Nämä ''viipaleet'' ovat suorakulmioita, joten integrointi niiden yli voidaan myös suorittaa iteroimalla kummassa järjestyksessä tahansa:
I
=
∭
A
(
x
y
2
+
z
3
)
d
V
=
∫
z
=
0
c
d
z
∫
y
=
0
b
d
y
∫
x
=
0
a
(
x
y
2
+
z
3
)
d
x
=
∫
z
=
0
c
d
z
∫
y
=
0
b
d
y
|
x
=
0
a
(
1
2
x
2
y
2
+
x
z
3
)
=
∫
z
=
0
c
d
z
∫
y
=
0
b
(
1
2
a
2
y
2
+
a
z
3
)
d
y
=
∫
z
=
0
c
d
z
|
y
=
0
b
(
1
6
a
2
y
3
+
a
y
z
3
)
=
∫
0
c
(
1
6
a
2
b
3
+
a
b
z
3
)
d
z
=
|
0
c
(
1
6
a
2
b
3
z
+
1
4
a
b
z
4
)
=
1
6
a
2
b
3
c
+
1
4
a
b
c
4
{\displaystyle {\begin{aligned}I&=\iiint _{A}(xy^{2}+z^{3})\,\mathrm {d} V\\&=\int _{z=0}^{c}\mathrm {d} z\int _{y=0}^{b}\mathrm {d} y\int _{x=0}^{a}(xy^{2}+z^{3})\,\mathrm {d} x\\&=\int _{z=0}^{c}\mathrm {d} z\int _{y=0}^{b}\mathrm {d} y\,{\Bigg |}_{x=0}^{a}\left({\frac {1}{2}}x^{2}y^{2}+xz^{3}\right)\\&=\int _{z=0}^{c}\mathrm {d} z\int _{y=0}^{b}\left({\frac {1}{2}}a^{2}y^{2}+az^{3}\right)\,\mathrm {d} y\\&=\int _{z=0}^{c}\mathrm {d} z\,{\Bigg |}_{y=0}^{b}\left({\frac {1}{6}}a^{2}y^{3}+ayz^{3}\right)\\&=\int _{0}^{c}\left({\frac {1}{6}}a^{2}b^{3}+abz^{3}\right)\,\mathrm {d} z\\&={\Bigg |}_{0}^{c}\left({\frac {1}{6}}a^{2}b^{3}z+{\frac {1}{4}}abz^{4}\right)\\&={\frac {1}{6}}a^{2}b^{3}c+{\frac {1}{4}}abc^{4}\end{aligned}}}
Koordinaatistomuunnoksessa
{
x
=
x
(
u
,
v
)
y
=
y
(
u
,
v
)
{\displaystyle {\begin{cases}x=x(u,v)\\y=y(u,v)\end{cases}}}
suorat
u
=
u
0
{\textstyle u=u_{0}}
ja
v
=
v
0
{\textstyle v=v_{0}}
kuvautuvat käyriksi
{
x
=
x
(
u
0
,
v
)
y
=
y
(
u
0
,
v
)
{\displaystyle {\begin{cases}x=x(u_{0},v)\\y=y(u_{0},v)\end{cases}}}
ja
{
x
=
x
(
u
,
v
0
)
y
=
y
(
u
,
v
0
)
{\displaystyle {\begin{cases}x=x(u,v_{0})\\y=y(u,v_{0})\end{cases}}}
x
y
{\textstyle xy}
-tasossa. Piste
(
u
0
,
v
0
)
{\textstyle (u_{0},v_{0})}
kuvautuu pisteeksi
(
x
0
,
y
0
)
{\textstyle (x_{0},y_{0})}
.
Kuten yksiulotteisessa tapauksessa , myös moniulotteisen integraalin selvittämistä voidaan tuntuvasti helpottaa muuttujanvaihdolla . Käytetään tässä yksinkertaisuuden vuoksi esimerkkinä integrointia kaksiulotteisessa avaruudessa. Joskus integraalin tai integrointialueen kannalta on luontevampaa käyttää karteesisten koordinaattien
x
{\textstyle x}
ja
y
{\textstyle y}
sijaan muita koordinaattijärjestelmiä. Oletetaan, että muuttujat
x
{\textstyle x}
ja
y
{\textstyle y}
voidaan esittää kahden muun muuttujan,
u
{\textstyle u}
ja
v
{\textstyle v}
, funktioina:
{
x
=
x
(
u
,
v
)
y
=
y
(
u
,
v
)
.
{\displaystyle {\begin{cases}x=x(u,v)\\y=y(u,v).\end{cases}}}
Tätä funktioparia sanotaan koordinaatistomuunnokseksi
u
v
{\textstyle uv}
-tason osajoukolta
D
{\textstyle D}
x
y
{\textstyle xy}
-tason osajoukolle
S
{\textstyle S}
. Muunnoksen täytyy olla injektio , jotta integrointi muuttujanvaihdolla onnistuisi.[ 6] Tällöin on olemassa käänteismuunnos
{
u
=
u
(
x
,
y
)
v
=
v
(
x
,
y
)
{\displaystyle {\begin{cases}u=u(x,y)\\v=v(x,y)\end{cases}}}
joukolta
S
{\textstyle S}
joukolle
D
{\textstyle D}
.[ 6] Jos funktioilla
x
(
u
,
v
)
{\textstyle x(u,v)}
ja
y
(
u
,
v
)
{\textstyle y(u,v)}
on olemassa jatkuvat ensimmäisen kertaluvun osittaisderivaatat ja Jacobin determinantti
∂
(
u
,
v
)
∂
(
x
,
y
)
=
|
∂
u
∂
x
∂
u
∂
y
∂
v
∂
x
∂
v
∂
y
|
≠
0
{\displaystyle {\frac {\partial (u,v)}{\partial (x,y)}}={\begin{vmatrix}{\frac {\partial u}{\partial x}}&{\frac {\partial u}{\partial y}}\\{\frac {\partial v}{\partial x}}&{\frac {\partial v}{\partial y}}\end{vmatrix}}\neq 0}
pisteessä
(
u
,
v
)
{\textstyle (u,v)}
, niin myös käänteismuunnos on injektio pisteen
(
u
,
v
)
{\textstyle (u,v)}
ympäristössä . Myös käänteismuunnoksella on olemassa jatkuvat 1. kertaluvun osittaisderivaatat ja nollasta eroava Jacobin determinantti, jolloin
∂
(
u
,
v
)
∂
(
x
,
y
)
=
(
∂
(
x
,
y
)
∂
(
u
,
v
)
)
−
1
{\displaystyle {\frac {\partial (u,v)}{\partial (x,y)}}=\left({\frac {\partial (x,y)}{\partial (u,v)}}\right)^{-1}}
[ 6]
Koordinaatistomuunnoksessa joukolta
S
{\textstyle S}
(
x
y
{\textstyle xy}
-tasolta) joukolle
D
{\textstyle D}
(
u
v
{\textstyle uv}
-tasolle) integroitava funktio
f
:
S
→
R
{\textstyle f:S\to \mathbb {R} }
,
f
(
x
,
y
)
{\textstyle f(x,y)}
muuntuu funktioksi
g
:
D
→
R
{\textstyle g:D\to \mathbb {R} }
,
g
(
u
,
v
)
=
f
(
x
(
u
,
v
)
,
y
(
u
,
v
)
)
{\displaystyle g(u,v)=f\left(x(u,v),y(u,v)\right)}
.
Osoittautuu, että Jacobin determinantin itseisarvo on eri koordinaatistoissa esitettyjen differentiaalisten pinta-alaelementtien suhde :
d
A
=
d
x
d
y
=
|
∂
(
x
,
y
)
∂
(
u
,
v
)
|
d
u
d
v
{\displaystyle \mathrm {d} A=\mathrm {d} x\,\mathrm {d} y=\left|{\frac {\partial (x,y)}{\partial (u,v)}}\right|\,\mathrm {d} u\,\mathrm {d} v}
[ 6]
Tällöin kaksinkertainen integraali voidaan ratkaista muuttujanvaihdolla (eli koordinaatistomuunnoksella):
Olkoon
x
=
x
(
u
,
v
)
{\textstyle x=x(u,v)}
ja
y
=
y
(
u
,
v
)
{\textstyle y=y(u,v)}
injektioita
u
v
{\textstyle uv}
-tason osajoukolta
D
{\textstyle D}
x
y
{\textstyle xy}
-tason osajoukolle
S
{\textstyle S}
. Olkoon funktioilla
x
(
u
,
v
)
{\textstyle x(u,v)}
ja
y
(
u
,
v
)
{\textstyle y(u,v)}
olemassa joukossa
D
{\textstyle D}
jatkuvat 1. kertaluvun osittaisderivaatat muuttujien
u
{\textstyle u}
ja
v
{\textstyle v}
suhteen. Jos funktio
f
(
x
,
y
)
{\textstyle f(x,y)}
on integroituva
S
{\textstyle S}
:ssä ja jos
g
(
u
,
v
)
=
f
(
x
(
u
,
v
)
,
y
(
u
,
v
)
)
{\textstyle g(u,v)=f\left(x(u,v),y(u,v)\right)}
, niin:
∬
S
f
(
x
,
y
)
d
x
d
y
=
∬
D
g
(
u
,
v
)
|
∂
(
x
,
y
)
∂
(
u
,
v
)
|
d
u
d
v
{\displaystyle \iint _{S}f(x,y)\,\mathrm {d} x\,\mathrm {d} y=\iint _{D}g(u,v)\left|{\frac {\partial (x,y)}{\partial (u,v)}}\right|\,\mathrm {d} u\,\mathrm {d} v}
[ 6] [ 7]
Kolminkertaisen integraalin ratkaiseminen muuttujanvaihdolla käy vastaavalla tavalla:
Olkoon
x
=
x
(
u
,
v
,
w
)
{\textstyle x=x(u,v,w)}
,
y
=
y
(
u
,
v
,
w
)
{\textstyle y=y(u,v,w)}
ja
z
=
z
(
u
,
v
,
w
)
{\textstyle z=z(u,v,w)}
injektioita
u
v
w
{\textstyle uvw}
-avaruuden osajoukolta
D
{\textstyle D}
x
y
z
{\textstyle xyz}
-avaruuden osajoukolle
S
{\textstyle S}
. Olkoon funktioilla
x
(
u
,
v
,
w
)
{\textstyle x(u,v,w)}
,
y
(
u
,
v
,
w
)
{\textstyle y(u,v,w)}
ja
z
(
u
,
v
,
w
)
{\textstyle z(u,v,w)}
olemassa joukossa
D
{\textstyle D}
jatkuvat 1. kertaluvun osittaisderivaatat muuttujien
u
{\textstyle u}
,
v
{\textstyle v}
ja
w
{\textstyle w}
suhteen. Jos funktio
f
(
x
,
y
,
z
)
{\textstyle f(x,y,z)}
on integroituva
S
{\textstyle S}
:ssä ja jos
g
(
u
,
v
,
w
)
=
f
(
x
(
u
,
v
,
w
)
,
y
(
u
,
v
,
w
)
,
z
(
u
,
v
,
w
)
)
{\textstyle g(u,v,w)=f\left(x(u,v,w),y(u,v,w),z(u,v,w)\right)}
, niin:
∭
S
f
(
x
,
y
,
z
)
d
x
d
y
d
z
=
∭
D
g
(
u
,
v
,
w
)
|
∂
(
x
,
y
,
z
)
∂
(
u
,
v
,
w
)
|
d
u
d
v
d
w
{\displaystyle \iiint _{S}f(x,y,z)\,\mathrm {d} x\,\mathrm {d} y\,\mathrm {d} z=\iiint _{D}g(u,v,w)\left|{\frac {\partial (x,y,z)}{\partial (u,v,w)}}\right|\,\mathrm {d} u\,\mathrm {d} v\,\mathrm {d} w}
[ 8] [ 9]
Tässä Jacobin determinatti on:
∂
(
x
,
y
,
z
)
∂
(
u
,
v
,
w
)
=
|
∂
x
∂
u
∂
x
∂
v
∂
x
∂
w
∂
y
∂
u
∂
y
∂
v
∂
y
∂
w
∂
z
∂
u
∂
z
∂
v
∂
z
∂
w
|
{\displaystyle {\frac {\partial (x,y,z)}{\partial (u,v,w)}}={\begin{vmatrix}{\frac {\partial x}{\partial u}}&{\frac {\partial x}{\partial v}}&{\frac {\partial x}{\partial w}}\\{\frac {\partial y}{\partial u}}&{\frac {\partial y}{\partial v}}&{\frac {\partial y}{\partial w}}\\{\frac {\partial z}{\partial u}}&{\frac {\partial z}{\partial v}}&{\frac {\partial z}{\partial w}}\end{vmatrix}}}
Muuttujanvaihto yleistyy n :lle muuttujalle:
Olkoon
x
=
x
(
u
)
{\displaystyle \mathbf {x} =\mathbf {x} (\mathbf {u} )}
kuvaus avoimelta joukolta
U
⊂
R
n
{\displaystyle U\subset \mathbb {R} ^{n}}
avoimelle joukolle
S
⊂
R
n
{\displaystyle S\subset \mathbb {R} ^{n}}
siten, että Jacobin determinantti
∂
(
x
1
,
…
,
x
n
)
∂
(
u
1
,
…
,
u
n
)
=
|
∂
x
1
∂
u
1
∂
x
1
∂
u
2
⋯
∂
x
1
∂
u
n
∂
x
2
∂
u
1
∂
x
2
∂
u
2
⋯
∂
x
2
∂
u
n
⋮
⋮
⋱
⋮
∂
x
n
∂
u
1
∂
x
n
∂
u
2
⋯
∂
x
n
∂
u
n
|
≠
0
{\displaystyle {\frac {\partial (x_{1},\dots ,x_{n})}{\partial (u_{1},\dots ,u_{n})}}={\begin{vmatrix}{\frac {\partial x_{1}}{\partial u_{1}}}&{\frac {\partial x_{1}}{\partial u_{2}}}&\dotsb &{\frac {\partial x_{1}}{\partial u_{n}}}\\{\frac {\partial x_{2}}{\partial u_{1}}}&{\frac {\partial x_{2}}{\partial u_{2}}}&\dotsb &{\frac {\partial x_{2}}{\partial u_{n}}}\\\vdots &\vdots &\ddots &\vdots \\{\frac {\partial x_{n}}{\partial u_{1}}}&{\frac {\partial x_{n}}{\partial u_{2}}}&\dotsb &{\frac {\partial x_{n}}{\partial u_{n}}}\end{vmatrix}}\neq 0}
.
Olkoon
D
{\displaystyle D}
rajoitettu
U
{\displaystyle U}
:n osajoukko ja funktio
f
{\displaystyle f}
integroituva kuvajoukossa
x
(
D
)
⊂
S
{\displaystyle \mathbf {x} (D)\subset S}
. Tällöin funktio
f
(
x
(
u
)
)
|
∂
(
x
1
,
…
,
x
n
)
∂
(
u
1
,
…
,
u
n
)
|
{\displaystyle f(\mathbf {x} (\mathbf {u} ))\left|{\frac {\partial (x_{1},\dots ,x_{n})}{\partial (u_{1},\dots ,u_{n})}}\right|}
on integroituva
D
{\displaystyle D}
:ssä ja
∫
…
∫
x
(
D
)
f
(
x
)
d
x
1
⋯
d
x
n
=
∫
…
∫
D
f
(
x
(
u
)
)
|
∂
(
x
1
,
…
,
x
n
)
∂
(
u
1
,
…
,
u
n
)
|
d
u
1
⋯
d
u
n
{\displaystyle \int \dots \int _{\mathbf {x} (D)}f(\mathbf {x} )\,\mathrm {d} x_{1}\dotsb \mathrm {d} x_{n}=\int \dots \int _{D}f(\mathbf {x} (\mathbf {u} ))\left|{\frac {\partial (x_{1},\dots ,x_{n})}{\partial (u_{1},\dots ,u_{n})}}\right|\,\mathrm {d} u_{1}\dotsb \mathrm {d} u_{n}}
.[ 10]
Sama integrointialue kuvattuna vasemmalla karteesisilla koordinaateilla ja oikealla napakoordinaateilla
Eräs käytännöllinen koordinaatistomuunnos kahdessa ulottuvuudessa on vaihtaa karteesiset koordinaatit napakoordinaateiksi :
{
x
(
r
,
θ
)
=
r
cos
θ
y
(
r
,
θ
)
=
r
sin
θ
,
{\displaystyle {\begin{cases}x(r,\theta )=r\cos \theta \\y(r,\theta )=r\sin \theta ,\end{cases}}}
missä
r
≥
0
{\textstyle r\geq 0}
ja
0
≤
θ
<
2
π
{\textstyle 0\leq \theta <2\pi }
(radiaaneina ). Tämän koordinaatistomuunnoksen Jacobin determinantti on
∂
(
x
,
y
)
∂
(
r
,
θ
)
=
|
cos
θ
−
r
sin
θ
sin
θ
r
cos
θ
|
=
r
{\displaystyle {\frac {\partial (x,y)}{\partial (r,\theta )}}={\begin{vmatrix}\cos \theta &-r\sin \theta \\\sin \theta &r\cos \theta \end{vmatrix}}=r}
.
Koordinaatistomuunnos on siis injektio kaikkialla paitsi origossa .
Ratkaistaan integraali
∬
S
(
y
2
/
x
2
)
d
A
{\textstyle \iint _{S}(y^{2}/x^{2})\,\mathrm {d} A}
, missä integrointialue
S
{\textstyle S}
on
x
{\textstyle x}
-akselin ja suoran
y
=
3
⋅
x
{\textstyle y={\sqrt {3}}\cdot x}
väliin jäävä osa ympyrärenkaasta
0
<
a
2
≤
x
2
+
y
2
≤
b
2
{\textstyle 0<a^{2}\leq x^{2}+y^{2}\leq b^{2}}
. Muuttujanvaihdon jälkeen integrandi on:
y
2
x
2
=
r
2
sin
2
θ
r
2
cos
2
θ
=
tan
2
θ
{\displaystyle {\frac {y^{2}}{x^{2}}}={\frac {r^{2}\sin ^{2}\theta }{r^{2}\cos ^{2}\theta }}=\tan ^{2}\theta }
Karteesisissa koordinaateissa integrointialue on osa origokeskistä ympyrärengasta, jonka sisäsäde on
a
>
0
{\textstyle a>0}
ja ulkosäde
b
≥
a
{\textstyle b\geq a}
ja jota raoittavat suorat
y
=
0
{\textstyle y=0}
ja
y
=
3
⋅
x
{\textstyle y={\sqrt {3}}\cdot x}
. Napakoordinaateissa nämä rajat ovat vastaavasti
0
<
a
≤
r
≤
b
{\textstyle 0<a\leq r\leq b}
sekä
θ
=
arctan
(
0
)
=
0
{\textstyle \theta =\arctan(0)=0}
ja
θ
=
arctan
(
3
)
=
π
/
3
{\textstyle \theta =\arctan({\sqrt {3}})=\pi /3}
. Integraali on tällöin:
∬
S
y
2
x
2
d
x
d
y
=
∬
D
tan
2
θ
|
∂
(
x
,
y
)
∂
(
r
,
θ
)
|
d
r
d
θ
=
∬
D
tan
2
(
θ
)
⋅
|
r
|
⏟
=
r
d
r
d
θ
=
∫
θ
=
0
π
/
3
tan
2
θ
d
θ
∫
r
=
a
b
r
d
r
=
1
2
(
b
2
−
a
2
)
∫
θ
=
0
π
/
3
(
sec
2
θ
−
1
)
d
θ
=
1
2
(
b
2
−
a
2
)
|
θ
=
0
π
/
3
(
tan
θ
−
θ
)
=
1
2
(
b
2
−
a
2
)
⋅
(
3
−
π
3
)
,
{\displaystyle {\begin{aligned}\iint _{S}{\frac {y^{2}}{x^{2}}}\,\mathrm {d} x\,\mathrm {d} y&=\iint _{D}\tan ^{2}\theta \left|{\frac {\partial (x,y)}{\partial (r,\theta )}}\right|\,\mathrm {d} r\,\mathrm {d} \theta \\&=\iint _{D}\tan ^{2}(\theta )\cdot \underbrace {|r|} _{=r}\,\mathrm {d} r\,\mathrm {d} \theta \\&=\int _{\theta =0}^{\pi /3}\tan ^{2}\theta \,\mathrm {d} \theta \int _{r=a}^{b}r\,\mathrm {d} r\\&={\frac {1}{2}}(b^{2}-a^{2})\int _{\theta =0}^{\pi /3}(\sec ^{2}\theta -1)\,\mathrm {d} \theta \\&={\frac {1}{2}}(b^{2}-a^{2}){\Bigg |}_{\theta =0}^{\pi /3}(\tan \theta -\theta )\\&={\frac {1}{2}}(b^{2}-a^{2})\cdot \left({\sqrt {3}}-{\frac {\pi }{3}}\right),\end{aligned}}}
missä käytettiin tietoa
sec
2
θ
−
tan
2
θ
=
1
{\textstyle \sec ^{2}\theta -\tan ^{2}\theta =1}
[ 11] .
Määritetään ellipsoidin
E
{\textstyle E}
tilavuus, kun
E
{\textstyle E}
:n määrittelee yhtälö
x
2
a
2
+
y
2
b
2
+
z
2
c
2
≤
1
{\displaystyle {\frac {x^{2}}{a^{2}}}+{\frac {y^{2}}{b^{2}}}+{\frac {z^{2}}{c^{2}}}\leq 1}
,
missä
a
>
0
{\textstyle a>0}
,
b
>
0
{\textstyle b>0}
ja
c
>
0
{\textstyle c>0}
ovat vakioita. Käytetään muuttujanvaihtoa
{
x
(
u
,
v
,
w
)
=
a
u
y
(
u
,
v
,
w
)
=
b
v
z
(
u
,
v
,
w
)
=
c
w
,
{\displaystyle {\begin{cases}x(u,v,w)=au\\y(u,v,w)=bv\\z(u,v,w)=cw,\end{cases}}}
jolloin ellipsoidi
E
{\textstyle E}
kuvautuu yksikköpalloksi
B
{\textstyle B}
, jossa
u
2
+
v
2
+
w
2
≤
1
{\displaystyle u^{2}+v^{2}+w^{2}\leq 1}
.
Jacobin determinantti tälle koordinaatistomuunnokselle on:
∂
(
x
,
y
,
z
)
∂
(
u
,
v
,
w
)
=
|
a
0
0
0
b
0
0
0
c
|
=
a
b
c
{\displaystyle {\frac {\partial (x,y,z)}{\partial (u,v,w)}}={\begin{vmatrix}a&0&0\\0&b&0\\0&0&c\end{vmatrix}}=abc}
Tällöin ellipsoidin
E
{\textstyle E}
tilavuus on (ks. ellipsoidin tilavuus ja pallon tilavuus ):
V
=
∭
E
1
d
V
=
∭
E
1
d
x
d
y
d
z
=
∭
B
|
∂
(
x
,
y
,
z
)
∂
(
u
,
v
,
w
)
|
d
u
d
v
d
w
=
a
b
c
∭
B
d
u
d
v
d
w
=
a
b
c
⋅
(
pallon
B
tilavuus
)
=
4
3
π
a
b
c
{\displaystyle {\begin{aligned}V&=\iiint _{E}1\,\mathrm {d} V=\iiint _{E}1\,\mathrm {d} x\,\mathrm {d} y\,\mathrm {d} z\\&=\iiint _{B}\left|{\frac {\partial (x,y,z)}{\partial (u,v,w)}}\right|\,\mathrm {d} u\,\mathrm {d} v\,\mathrm {d} w\\&=abc\iiint _{B}\mathrm {d} u\,\mathrm {d} v\,\mathrm {d} w\\&=abc\cdot ({\text{pallon }}B{\text{ tilavuus}})\\&={\frac {4}{3}}\pi abc\end{aligned}}}
Tässä esimerkissä sovellettiin sekä symmetrioiden käyttöä, että muuttujanvaihtoa.
Funktion
f
:
D
→
R
{\textstyle f:D\to \mathbb {R} }
, missä
D
⊂
R
2
{\textstyle D\subset \mathbb {R} ^{2}}
on rajoitettu joukko, ja
x
y
{\textstyle xy}
-tason väliin jäävä tilavuus on
V
=
∬
D
f
(
x
,
y
)
d
A
{\displaystyle V=\iint _{D}f(x,y)\,\mathrm {d} A}
Kolmiulotteisen kappaleen, jota rajoittaa avaruus
A
⊂
R
3
{\textstyle A\subset \mathbb {R} ^{3}}
, massa saadaan kolminkertaisella integraalilla
m
=
∭
A
ρ
(
x
,
y
,
z
)
d
V
{\displaystyle m=\iiint _{A}\rho (x,y,z)\,\mathbb {d} V}
,
missä
ρ
(
x
,
y
,
z
)
{\textstyle \rho (x,y,z)}
on kappaleen tiheys pisteessä
(
x
,
y
,
z
)
∈
A
{\textstyle (x,y,z)\in A}
.
Jos
S
{\textstyle {\mathcal {S}}}
on pinta kolmiulotteisessa avaruudessa siten, että
z
=
f
(
x
,
y
)
{\textstyle z=f(x,y)}
on rajoitettu ja jatkuva funktio sekä
(
x
,
y
)
∈
D
{\textstyle (x,y)\in D}
, missä
D
⊂
R
2
{\textstyle D\subset \mathbb {R} ^{2}}
on rajoitettu joukko, niin pinnan
S
{\textstyle {\mathcal {S}}}
pinta-ala on:
A
=
∬
D
1
+
(
∂
z
∂
x
)
2
+
(
∂
z
∂
y
)
2
d
A
{\displaystyle A=\iint _{D}{\sqrt {1+\left({\frac {\partial z}{\partial x}}\right)^{2}+\left({\frac {\partial z}{\partial y}}\right)^{2}}}\,\mathrm {d} A}
[ 12] [ 13]
Kolmiulotteisen kappaleen, jolla on jatkuva tiheysjakauma
ρ
(
x
,
y
,
z
)
{\textstyle \rho (x,y,z)}
alueessa
B
⊂
R
3
{\textstyle B\subset \mathbb {R} ^{3}}
, massakeskipiste
(
x
0
,
y
0
,
z
0
)
{\textstyle (x_{0},y_{0},z_{0})}
määritetään kaavoilla:
{
x
0
=
1
m
∭
B
x
ρ
d
V
y
0
=
1
m
∭
B
y
ρ
d
V
z
0
=
1
m
∭
B
z
ρ
d
V
,
{\displaystyle {\begin{cases}x_{0}=\displaystyle {{\frac {1}{m}}\iiint _{B}x\rho \,\mathrm {d} V}\\y_{0}=\displaystyle {{\frac {1}{m}}\iiint _{B}y\rho \,\mathrm {d} V}\\z_{0}=\displaystyle {{\frac {1}{m}}\iiint _{B}z\rho \,\mathrm {d} V},\end{cases}}}
missä
m
=
∭
B
ρ
(
x
,
y
,
z
)
d
V
{\textstyle m=\iiint _{B}\rho (x,y,z)\,\mathrm {d} V}
on kappaleen kokonaismassa.[ 14] Vektorimuodossa kappaleen massakeskipisteen paikkavektori
r
→
0
=
x
0
i
+
y
0
j
+
z
0
k
{\textstyle {\vec {r}}_{0}=x_{0}\mathbf {i} +y_{0}\mathbf {j} +z_{0}\mathbf {k} }
on
r
→
0
=
1
m
∭
B
ρ
r
→
d
V
{\displaystyle {\vec {r}}_{0}={\frac {1}{m}}\iiint _{B}\rho {\vec {r}}\,\mathrm {d} V}
,
missä
r
→
=
x
i
+
y
j
+
z
k
∈
R
3
{\textstyle {\vec {r}}=x\mathbf {i} +y\mathbf {j} +z\mathbf {k} \in \mathbb {R} ^{3}}
.[ 14]
Kolmiulotteisen kappaleen, jolla on jatkuva tiheysjakauma
ρ
(
x
,
y
,
z
)
{\textstyle \rho (x,y,z)}
alueessa
B
⊂
R
3
{\textstyle B\subset \mathbb {R} ^{3}}
, hitausmomentti pyörimisakselin
L
{\textstyle L}
suhteen on
I
=
∭
B
r
2
d
m
=
∭
B
r
2
ρ
d
V
{\displaystyle I=\iiint _{B}r^{2}\,\mathrm {d} m=\iiint _{B}r^{2}\rho \,\mathrm {d} V}
,
missä
r
{\textstyle r}
on differentiaalisen massaelementin
d
m
=
ρ
d
V
{\textstyle \mathrm {d} m=\rho \,\mathrm {d} V}
etäisyys pyörimisakselista
L
{\textstyle L}
.[ 15]
Kolmiulotteisen kappaleen, jolla on jatkuva tiheysjakauma
ρ
(
x
,
y
,
z
)
{\textstyle \rho (x,y,z)}
alueessa
B
⊂
R
3
{\textstyle B\subset \mathbb {R} ^{3}}
, hitaussäde pyörimisakselin
L
{\textstyle L}
suhteen on
r
g
=
I
m
=
∭
B
r
2
ρ
d
V
∭
B
ρ
d
V
{\displaystyle r_{\text{g}}={\sqrt {\frac {I}{m}}}={\sqrt {\frac {\iiint _{B}r^{2}\rho \,\mathrm {d} V}{\iiint _{B}\rho \,\mathrm {d} V}}}}
.[ 15]
↑ a b c d e f g h i j k l m n o p q r Adams, Robert A. & Essex, Christopher: Calculus: A Complete Course, 8. painos , s. 807−811. Pearson, 2014. ISBN 978-0-321-78107-9 (englanniksi)
↑ Adams & Essex, s. 835
↑ a b c Adams & Essex, s. 855
↑ a b c d e f g h i j k l Ylinen, Kari: Usean muuttujan funktiot II 2010. Turun yliopisto. Viitattu 11.4.2017. [vanhentunut linkki ]
↑ a b c d e f g Adams & Essex, s. 813−815
↑ a b c d e Adams & Essex, s. 829−831
↑ Friedman, Avner: Advanced Calculus , s. 324. Holt, Rinehart and Winston, Inc., 1971. ISBN 0-03-083983-1 (englanniksi)
↑ Adams & Essex, s. 842
↑ Friedman, s. 346
↑ Friedman, s. 277
↑ Spiegel, Murray R. & Lipschutz, Seymour & Liu, John: Schaum's outlines: Mathematical Handbook of Formulas and Tables, 3. painos , s. 44. McGraw-Hill, 2009. ISBN 978-0-07-154855-7 (englanniksi)
↑ Adams & Essex, s. 849
↑ Friedman, s. 331
↑ a b Adams & Essex, s. 851
↑ a b Adams & Essex, s. 853