Introduzione

Se prendiamo una funzione non costante, la regione delimitata dal suo grafico ed il
piano xy non è un parallelepipedo e dobbiamo definire il suo volume. Cominciamo
con il caso in cui il dominio della funzione sia comunque un rettangolo. Si definisce
cilindroide la regione sopra il rettangolo R e sotto il grafico di f :
𝑉 = (𝑥, 𝑦, 𝑧 ∈ 𝑅3 : 𝑎 ≤ 𝑥 ≤ 𝑏, 𝑐 ≤ 𝑦 ≤ 𝑑, 0 ≤ 𝑧 ≤ 𝑓(𝑥, 𝑦)
 Somme di Riemann
Analogamente al caso del calcolo di un’area sottesa ad una curva, possiamo calcolare il
volume di un cilindroide come il limite dei volumi approssimanti: le somme di Riemann.
Consideriamo una funzione f definita e limitata su un rettangolo chiuso e limitato.
Supponiamo (solo inizialmente) che la funzione sia positiva e chiamiamo V il solido
compreso tra R e il grafico di f:
𝑉 = (𝑥, 𝑦, 𝑧 ∈ 𝑅3 : 𝑎 ≤ 𝑥 ≤ 𝑏, 𝑐 ≤ 𝑦 ≤ 𝑑, 0 ≤ 𝑧 ≤ 𝑓(𝑥, 𝑦)
Il primo passo consiste nel dividere R in sottorettangoli, dividendo l’intervallo
[a, b] in n sottointervalli [𝑥𝑖−1 , 𝑥𝑖 ] con n + 1 punti 𝑎 = 𝑥0 < 𝑥1 < 𝑥2 < ⋯ < 𝑥𝑛 = 𝑏
Analogamente dividiamo l’intervallo [c, d] (sull’asse delle y) in m sottointervalli
[𝑦𝑗−1 , 𝑦𝑗 ] con m+1 punti 𝑐 = 𝑦0 < 𝑦1 < 𝑦2 < ⋯ < 𝑦𝑚 = 𝑑
Tracciando le n+1 rette parallele all’asse y di equazione 𝑥 = 𝑥𝑖 e le m+1 rette parallele
all’asse x di equazione 𝑦 = 𝑦𝑗 , otteniamo una partizione di R in 𝑛 ∙ 𝑚 sottorettangoli di
R,
𝑅𝑖𝑗 = 𝑥𝑖−1 , 𝑥𝑖 × 𝑦𝑗−1 , 𝑦𝑗 = 𝑥, 𝑦 : 𝑥𝑖−1 ≤ 𝑥 ≤ 𝑥𝑖 ; 𝑦𝑗−1 ≤ 𝑦 ≤ 𝑦𝑗
ciascuno di area 𝐴𝑖𝑗 = 𝑥𝑖 − 𝑥𝑖−1 𝑦𝑗 − 𝑦𝑗−1 .
 Somme di Riemann
Se f ≥ 0 una sua somma di Riemann rappresenta il volume della regione formata da
parallelepipedi affiancati (pluri-parallelepipedo).

Definizione:
chiamiamo ampiezza della partizione la lunghezza della più grande diagonale di un
sottorettangolo della suddivisione:
2 2
𝛿𝑃 = 𝑚𝑎𝑥 𝑥𝑖 − 𝑥𝑖−1 + 𝑦𝑗 − 𝑦𝑗−1 : 𝑖 = 1, … , 𝑛, 𝑗 = 1, … , 𝑚
 Somme di Riemann
∗ ∗
Scegliamo ora un punto base (𝑥𝑖𝑗 , 𝑦𝑖𝑗 ) in ognuno dei rettangoli 𝑅𝑖𝑗 , e consideriamo un
∗ ∗
parallelepipedo di base 𝑅𝑖𝑗 e altezza f(𝑥𝑖𝑗 , 𝑦𝑖𝑗 ). Il volume di questo parallelepipedo è
∗ ∗
pari all’area del rettangolo di base per l’altezza, vale a dire 𝑓(𝑥𝑖𝑗 , 𝑦𝑖𝑗 ) ∙ 𝐴𝑖𝑗 . Ora
compiamo questa operazione per ciascuno dei rettangoli 𝑅𝑖𝑗 della suddivisione e
sommiamo i volumi dei parallelepipedi corrispondenti: otteniamo la quantità
𝑛 𝑚
∗ ∗
𝑓(𝑥𝑖𝑗 , 𝑦𝑖𝑗 ) ∙ 𝐴𝑖𝑗
𝑖=1 𝑗=1

che possiamo interpretare come un’approssimazione del volume totale di V.

Definizione
Chiamiamo somma di Riemann della funzione f , corrispondente alla suddivisione di R
∗ ∗
nei sottorettangoli 𝑅𝑖𝑗 e alla scelta di punti base (𝑥𝑖𝑗 , 𝑦𝑖𝑗 ) ∈ 𝑅𝑖𝑗 , la somma
𝑛 𝑚
∗ ∗
𝑓(𝑥𝑖𝑗 , 𝑦𝑖𝑗 ) ∙ 𝐴𝑖𝑗
𝑖=1 𝑗=1
 Somme di Riemann
Oltre alle somme di Riemann possiamo definire la somma inferiore e la
somma superiore, relative alla partizione P:
𝑛 𝑚

𝑠 𝑃 = 𝑚𝑖𝑗 ∙ 𝐴𝑖𝑗 𝑚𝑖𝑗 = inf 𝑓

𝑅𝑖𝑗
𝑖=1 𝑗=1
𝑛 𝑚

𝑆 𝑃 = 𝑀𝑖𝑗 ∙ 𝐴𝑖𝑗 𝑀𝑖𝑗 = sup 𝑓

𝑅𝑖𝑗
𝑖=1 𝑗=1
 Convergenza delle somme di Riemann
Esse rappresentano l’approssimazione per difetto e per eccesso del volume V: si ha
infatti, per ogni scelta dei punti base,
𝑛 𝑚
∗ ∗
𝑠(𝑃) ≤ 𝑓 𝑥𝑖𝑗 , 𝑦𝑖𝑗 ∙ 𝐴𝑖𝑗 ≤ 𝑆(𝑃)
𝑖=1 𝑗=1

Infatti, per definizione

∗ ∗
𝑚𝑖𝑗 ≤ 𝑓 𝑥𝑖𝑗 , 𝑦𝑖𝑗 ≤ 𝑀𝑖𝑗

Diciamo che le somme di Riemann della funzione f convergono al limite S quando

l’ampiezza della suddivisione tende a zero se, fissato un qualunque ε > 0, possiamo
trovare δ > 0 per cui, comunque si prenda una suddivisione di ampiezza 𝛿𝑃 < 𝛿 e
∗ ∗
comunque si scelgano i punti base 𝑥𝑖𝑗 , 𝑦𝑖𝑗 , la somma di Riemann corrispondente
differisca da S per una quantità minore di ε:
𝑛 𝑚
∗ ∗
𝛿𝑃 < 𝛿 ⟹ 𝑓 𝑥𝑖𝑗 , 𝑦𝑖𝑗 ∙ 𝐴𝑖𝑗 − 𝑆 < 𝜀
𝑖=1 𝑗=1
Quando tale limite esiste, il suo valore definisce l’integrale doppio della funzione f sul
rettangolo R.
 Definizione di integrale doppio
Definizione
L’integrale doppio di f sul rettangolo R è
𝑛 𝑚
∗ ∗
𝑓 𝑥, 𝑦 𝑑𝑥𝑑𝑦 = lim 𝑓 𝑥𝑖𝑗 , 𝑦𝑖𝑗 ∙ 𝐴𝑖𝑗
𝛿𝑃
𝑖=1 𝑗=1
quando il limite delle somme di Riemann esiste. In tal caso diciamo che la funzione è
integrabile sul rettangolo.

Chiamiamo volume (con segno) della regione solida V compresa tra il grafico di z = f (x,
y) e il rettangolo R il valore del limite. Nel caso di funzioni positive l’integrale definisce
il volume del solido V:

𝑓≥0 ⟹ 𝑉𝑜𝑙𝑢𝑚𝑒 𝑉 = 𝑓 𝑥, 𝑦 𝑑𝑥𝑑𝑦

Teorema. Una funzione limitata è integrabile sul rettangolo se e soltanto se l’estremo

superiore delle somme inferiori, fra tutte partizioni del rettangolo, è uguale all’estremo
inferiore delle somme superiori fra le partizioni. Tali estremi sono a loro volta uguali
all’integrale doppio.
L’esistenza dell’integrale, dunque, significa che, al tendere a zero della
massima ampiezza dei rettangoli della suddivisione di R, i volumi (con
segno) di tutti i possibili pluri-parallelepipedi tendono ad un limite
comune.
 Il calcolo degli integrali doppi su
rettangoli
Per calcolare gli integrali doppi, così come capita per quelli semplici, non si applica
praticamente mai la definizione di somma di Riemann, ma ci si riduce al calcolo di
due integrali semplici. Per fare questa riduzione occorre la nozione di integrazione
parziale, operazione che corrisponde alla derivazione parziale.
Data la funzione
𝑓: 𝑎, 𝑏 × [𝑐, 𝑑] → 𝑅
integrare parzialmente rispetto a x significa integrare rispetto a 𝑦 ∈ [𝑐, 𝑑] la famiglia
delle tracce di f ad 𝑥 ∈ [𝑎, 𝑏] fissato.
 Esempio
Se 𝑓 𝑥, 𝑦 = 𝑥 2 𝑦 3 e 𝑎, 𝑏 × 𝑐, 𝑑 = 1,3 × 0,1 otteniamo, per ogni valore fissato
di 𝑥 ∈ 1,3
1
𝑦 4 𝑦=1 𝑥2
2 3 2
𝑥 𝑦 𝑑𝑦 = 𝑥 =
0 4 𝑦=0
4

Integrando parzialmente la funzione 𝑧 = 𝑓(𝑥, 𝑦) rispetto alla variabile y otteniamo

una nuova funzione A della sola variabile x:
𝑑
𝐴 𝑥 = 𝑓 𝑥, 𝑦 𝑑𝑦, 𝑥 ∈ [𝑎, 𝑏]
𝑐
Questo integrale viene chiamato integrale iterato e il suo valore è un numero, che ha
un preciso significato geometrico. Infatti, se prendiamo per buono il Principio di
Cavalieri secondo il quale se due regioni solide hanno sezioni trasversali
corrispondenti di uguale area, allora hanno lo stesso volume, possiamo interpretare
l’integrale iterato come il volume del solido V:
𝑏
𝑉𝑜𝑙𝑢𝑚𝑒 𝑉 = 𝐴 𝑥 𝑑𝑥
𝑎
 Teorema (formule di riduzione)
Teorema

Se f è continua sul rettangolo 𝑅 = 𝑎, 𝑏 × [𝑐, 𝑑] allora è integrabile e l’integrale

doppio è uguale all’integrale iterato:

𝑏 𝑑 𝑑 𝑏
𝑓 𝑥, 𝑦 𝑑𝑥𝑑𝑦 = 𝑓 𝑥, 𝑦 𝑑𝑦 𝑑𝑥 = 𝑓 𝑥, 𝑦 𝑑𝑥 𝑑𝑦
𝑎 𝑐 𝑐 𝑎

Esempio 1 Esempio 2
Su R = [1, 2] × [0, π], si calcoli l’integrale Si calcoli il volume del solido V che giace
della funzione f (x, y) = x sin y sotto l grafico di 𝑧 = 𝑦 2 − 𝑥 2 e sopra il
quadrato 𝑅 = −1,1 × 1,3
2 𝜋
𝑥 𝑠𝑖𝑛𝑦𝑑𝑦 𝑑𝑥
1 0
 Integrazione su domini generici

La definizione e la teoria che abbiamo sviluppato finora, e in particolare il

Teorema di riduzione, ci permettono di integrare le funzioni continue sui
rettangoli. D’altra parte, nelle applicazioni, capita di frequente di voler
integrare delle funzioni su domini che non siano rettangoli.

Teorema. Un sottoinsieme limitato D del piano è misurabile secondo Peano–Jordan se e

soltanto se l’estremo superiore delle aree dei poligoni contenuti in D eguaglia l’estremo
inferiore delle aree di quelli contenenti D e i due valori sono uguali all’area.
Proprietà dell’integrale doppio
Proprietà di
linearità

Proprietà di
monotonia
 Integrazione su regioni semplici

Definizione
Una regione 𝐷 ⊂ 𝑅2 è detta y–semplice se è compresa tra i grafici di due funzioni della
variabile x, cioè se è del tipo
𝐷 = 𝑥, 𝑦 ∈ 𝑅2 : 𝑎 ≤ 𝑥 ≤ 𝑏, 𝑔1 (𝑥) ≤ 𝑦 ≤ 𝑔2 (𝑥)
Quindi l’area di una regione semplice è data dalla formula:
𝑏
𝐴𝑟𝑒𝑎 𝐷 = [𝑔2 𝑥 − 𝑔1 (𝑥)]𝑑𝑥
𝑎
Analogamente si possono definire le regioni x–semplici come quelle del tipo
𝐷 = 𝑥, 𝑦 ∈ 𝑅2 : 𝑐 ≤ 𝑦 ≤ 𝑑 ℎ1 (𝑦) ≤ 𝑥 ≤ ℎ2 (𝑦)
 Teorema (formule di riduzione)
Ogni funzione f continua su un insieme semplice D è integrabile su D. inoltre valgono le
formule di riduzione:
𝑏 𝑔2 (𝑥)
𝐷 𝑦 − 𝑠𝑒𝑚𝑝𝑙𝑖𝑐𝑒 ⟹ 𝑓 𝑥, 𝑦 𝑑𝑥𝑑𝑦 = 𝑓 𝑥, 𝑦 𝑑𝑦 𝑑𝑥
𝑎 𝑔1 (𝑥)

𝑑 ℎ2 (𝑥)
𝐷 𝑥 − 𝑠𝑒𝑚𝑝𝑙𝑖𝑐𝑒 ⟹ 𝑓 𝑥, 𝑦 𝑑𝑥𝑑𝑦 = 𝑓 𝑥, 𝑦 𝑑𝑥 𝑑𝑦
𝑐 ℎ1 (𝑥)
Esempi
Calcoliamo

𝑥 + 𝑦 2 𝑑𝑥𝑑𝑦

nel dominio D rappresentato in figura

Calcoliamo
2
𝑒 𝑦 𝑑𝑥𝑑𝑦

nel dominio D rappresentato in figura

 Integrazione su domini non semplici
Teorema. Supponiamo che gli insiemi semplici 𝐷1 , 𝐷2 , … , 𝐷𝑘 non abbiano, a due a due,
punti in comune oltre ad una parte della frontiera. Allora ogni funzione continua
sull’unione 𝐷 = 𝐷1 ∪ 𝐷2 ∪ ⋯ ∪ 𝐷𝑘 è integrabile e si ha:

𝑓 𝑥, 𝑦 𝑑𝑥𝑑𝑦 = 𝑓 𝑥, 𝑦 𝑑𝑥𝑑𝑦 + ⋯ + 𝑓 𝑥, 𝑦 𝑑𝑥𝑑𝑦

𝐷 𝐷1 𝐷2

Possiamo usare il Teorema di additività rispetto al dominio per estendere le nostre

capacità di calcolo a tutte le regioni che non sono necessariamente semplici, ma che
possono essere suddivise in un numero finito di sottoregioni semplici .Tali domini
verranno detti semplicemente decomponibili.
 Esempio

Il dominio D può essere visto come unione di tre

domini semplici, così definiti:

𝐷1 = 𝑥, 𝑦 : −1 ≤ 𝑥 ≤ 0, −1 ≤ 𝑦 ≤ 1 + 𝑥 2 ,
𝐷2 = 𝑥, 𝑦 : 0 ≤ 𝑥 ≤ 1, −1 ≤ 𝑦 ≤ −√𝑥
𝐷3 = 𝑥, 𝑦 : 0 ≤ 𝑥 ≤ 1, 𝑥 ≤ 𝑦 ≤ 1 + 𝑥 2

10/3

𝐴𝑟𝑒𝑎 𝐷 𝑑𝑥𝑑𝑦 = 𝑑𝑥𝑑𝑦 + 𝑑𝑥𝑑𝑦 + 𝑑𝑥𝑑𝑦

𝐷 𝐷1 𝐷2 𝐷3
 Applicazioni dell’integrale doppio: la
densità
Il procedimento di passaggio al limite utilizzato per definire l’integrale suggerisce una
metodologia generale per calcolare quantità di interesse in molte applicazioni, dalla
geometria, alla probabilità e alla fisica, attraverso il passaggio al limite da un
problema discreto (cioè con un numero finito di parti) ad un problema continuo.

Consideriamo una lamina che occupa una regione D del piano. Dire che la densità
(superficiale) della lamina è costantemente uguale a 𝜌0 , significa dire che rapporto
fra la massa contenuta in un rettangolo R ⊂ D e l’area di R è costante:
𝑚𝑎𝑠𝑠𝑎 (𝑅)
= 𝜌0
𝐴𝑟𝑒𝑎 (𝑅)

Se invece tale rapporto non è costante, possiamo considerare il rapporto

massa(R)/Area(R) come una densità media e definire la densità nel punto (x, y) come
il limite (se esiste)
𝑚𝑎𝑠𝑠𝑎 𝑅𝛿 (𝑥, 𝑦)
𝜌 𝑥, 𝑦 = lim
𝛿→0 𝐴𝑟𝑒𝑎 𝑅𝛿 (𝑥, 𝑦)

al tendere a zero della diagonale δ del rettangolo 𝑅𝛿 (𝑥, 𝑦) contenente il punto (x, y),
 Densità

Ora, rovesciando il punto di vista, supponiamo di conoscere la funzione densità

puntuale e di volerne dedurre la massa totale della lamina, o di una sua porzione:
possiamo utilizzare il procedimento di approssimazione tramite le somme di Riemann.
Suddividiamo infatti un rettangolo contenente D in sottorettangoli 𝑅𝑖𝑗 e poniamo ρ
uguale a zero fuori da D.
Se il rettangolino 𝑅𝑖𝑗 è molto piccolo, e se la densità è una funzione continua,
possiamo pensare che ρ sia approssimativamente costante in 𝑅𝑖𝑗 .
La massa contenuta nel rettangolino sarà quindi circa
𝑚𝑎𝑠𝑠𝑎 𝑅𝑖𝑗 ~𝜌 𝑥𝑖𝑗 , 𝑦𝑖𝑗 ∙ 𝐴𝑟𝑒𝑎(𝑅𝑖𝑗 )
e la massa totale sarà approssimativamente una somma di Riemann:
𝑛 𝑚 𝑛 𝑚

𝑚𝑎𝑠𝑠𝑎 𝐷 = 𝑚𝑎𝑠𝑠𝑎 𝑅𝑖𝑗 ~ 𝜌 𝑥𝑖𝑗 , 𝑦𝑖𝑗 𝐴𝑟𝑒𝑎 𝑅𝑖𝑗

𝑖=1 𝑗=1 𝑖=1 𝑗=1
 Dal caso discreto al caso continuo
Se la funzione ρ è integrabile su D, le somme di Riemann convergono, al tendere a zero
dell’ampiezza della partizione, verso l’integrale doppio di f su D, che è quindi la massa
della lamina:

𝑚𝑎𝑠𝑠𝑎 𝐷 = 𝜌 𝑥, 𝑦 𝑑𝑥𝑑𝑦
𝐷

Il ragionamento esposto qui sopra si può ripetere, oltre al caso della massa, anche ad
altre grandezze fisiche. Per esempio, se una carica elettrica è distribuita su una regione
piana D, si parlerà della funzione densità di carica (prima media e poi, passando al limite,
puntuale) σ(x, y) (in unità di carica per unità di superficie). Nota σ, si può calcolare la
carica totale attraverso l’integrale:

𝑄= 𝜎 𝑥, 𝑦 𝑑𝑥𝑑𝑦
𝐷
 Baricentro
Per un sistema di n punti nel piano, di masse 𝑚𝑖 e posizioni 𝑃𝑖 𝑥𝑖 , 𝑦𝑖 , le coordinate del
baricentro sono:
𝑛 𝑛
1 1
𝑥= 𝑚𝑖 𝑥𝑖 , 𝑦= 𝑚𝑖 𝑦𝑖
𝑀 𝑀
𝑖=1 𝑖=1
Consideriamo una lamina piana, avente densità superficiale di massa variabile (e
continua) ρ(x, y),che occupa una regione D: determiniamone il baricentro.
A tale fine, ragioniamo approssimando la lamina con un sistema finito di punti materiali
(il cui numero tenderà poi all’infinito);
Consideriamo una partizione di D in rettangolini 𝑅𝑖𝑗 . Il sistema di punti approssimante è
dato dai punti 𝑃𝑖𝑗 (𝑥𝑖𝑗 , 𝑦𝑖𝑗 ) aventi masse 𝑚𝑖𝑗 = 𝜌 𝑥𝑖𝑗 , 𝑦𝑖𝑗 ∙ 𝐴𝑟𝑒𝑎(𝑅𝑖𝑗 ), dove i punti
𝑃𝑖𝑗 (𝑥𝑖𝑗 , 𝑦𝑖𝑗 ) sono punti base arbitrariamente scelti nei rettangolini 𝑅𝑖𝑗 . La coordinata x
del baricentro di tale sistema di punti è:
𝑛 𝑚 𝑛 𝑚
1 1
𝑥= 𝑚𝑖𝑗 ∙ 𝑥𝑖𝑗 = 𝑥 = 𝜌 𝑥𝑖𝑗 , 𝑦𝑖𝑗 ∙ 𝐴𝑟𝑒𝑎(𝑅𝑖𝑗 ) ∙ 𝑥𝑖𝑗
𝑚𝑎𝑠𝑠𝑎 (𝐷) 𝑚𝑎𝑠𝑠𝑎 (𝐷)
𝑖=1 𝑗=1 𝑖=1 𝑗=1
Nel secondo membro si può riconoscere una somma di Riemann della funzione
𝑓 𝑥, 𝑦 = 𝑥 ∙ 𝜌(𝑥, 𝑦)
 Dal caso discreto al caso continuo
Se la funzione 𝑓 𝑥, 𝑦 = 𝑥 ∙ 𝜌(𝑥, 𝑦) è integrabile, possiamo passare al limite al
tendere a zero dell’ampiezza della suddivisione, ottenendo
1
𝑥= 𝜌(𝑥, 𝑦) ∙ 𝑥 ∙ 𝑑𝑥𝑑𝑦
𝑚𝑎𝑠𝑠𝑎(𝐷) 𝐷
1
𝑦= 𝜌(𝑥, 𝑦) ∙ 𝑦 ∙ 𝑑𝑥𝑑𝑦
𝑚𝑎𝑠𝑠𝑎(𝐷) 𝐷
Ricordando l’espressione della massa D, otteniamo:

𝐷
𝜌(𝑥, 𝑦) ∙ 𝑥 ∙ 𝑑𝑥𝑑𝑦 𝐷
𝜌(𝑥, 𝑦) ∙ 𝑦 ∙ 𝑑𝑥𝑑𝑦
𝑥= 𝑦=
𝐷
𝜌 𝑥, 𝑦 𝑑𝑥𝑑𝑦 𝐷
𝜌 𝑥, 𝑦 𝑑𝑥𝑑𝑦

In particolare, se la densità è costante, allora si elide dalla frazione e si ottengono

le coordinate del baricentro geometrico o centroide di D:

𝐷
𝑥 ∙ 𝑑𝑥𝑑𝑦 𝐷
𝑦 ∙ 𝑑𝑥𝑑𝑦
𝑥= 𝑦=
𝐷
𝑑𝑥𝑑𝑦 𝐷
𝑑𝑥𝑑𝑦
 Momenti d’inerzia
Con un ragionamento del tutto analogo si perviene alla definizione dei momenti di
inerzia di una distribuzione piana di massa. Il momento di inerzia di una particella di
massa m rispetto a un asse è definito come 𝑚𝑟 2 dove r è la distanza della particella
dall’asse; ad esempio il momento di inerzia rispetto all’asse y è 𝑚𝑥 2 . Fisicamente
esso misura la resistenza del punto a mettersi in rotazione attorno all’asse in
considerazione. Se consideriamo una lamina piana di densità ρ(x, y) che occupa la
regione D, i momenti di inerzia rispetto agli assi x e y della lamina sono:

𝐼𝑥 = 𝜌(𝑥, 𝑦) ∙ 𝑦 2 𝑑𝑥𝑑𝑦 , 𝐼𝑦 = 𝜌(𝑥, 𝑦) ∙ 𝑥 2 𝑑𝑥𝑑𝑦

𝐷 𝐷

Il momento di inerzia polare, relativo all’asse ortogonale al piano e passante per

l’origine:

𝐼𝑂 = 𝜌(𝑥, 𝑦) ∙ 𝑥 2 + 𝑦 2 2 𝑑𝑥𝑑𝑦 ,
𝐷
 Cambiamento di variabili negli integrali
doppi
Spesso conviene esprimere un Dominio D rispetto a cui calcoliamo un integrale doppio
delle variabili x e y, attraverso un nuovo dominio T, sfruttando eventuali simmetrie sia
del dominio di integrazione che della funzione integranda.
Un cambiamento di variabili è una funzione

𝑇: 𝑢, 𝑣 ∈ 𝑅2 ⟼ (𝑥 𝑢, 𝑣 , 𝑦(𝑢, 𝑣)) ∈ 𝑅2

Consideriamo una funzione T, lineare dal piano uv al piano xy:

𝑢 𝑎𝑢 + 𝑏𝑣 𝑢 𝑎 𝑏
𝑇 = =𝐴 , 𝐴=
𝑣 𝑐𝑢 + 𝑑𝑣 𝑣 𝑐 𝑑

Per via della linearità, T trasforma parallelogrammi del piano uv in altri parallelogrammi
del piano xy.
Vi è un rapporto costante tra l’area di un parallelogramma S e quella della sua
immagine T(S), il cui fattore di proporzionalità è il modulo del determinante della
matrice associata A:
𝐴𝑟𝑒𝑎 𝑇 𝑆 = 𝑑𝑒𝑡𝐴 ∙ 𝐴𝑟𝑒𝑎(𝑆)
Cosa succede se S non è un parallelogramma? Dato che una funzione lineare manda
parallelogrammi in parallelogrammi, allora manda poligoni in poligoni, mantenendo il
medesimo fattore di proporzionalità fra l’area del poligono e quello della sua
immagine. In particolare, manda poligoni inclusi in S in poligoni inclusi in T(S) e
ugualmente per quelli contenenti S, che hanno immagine contenente T(S). Perciò
possiamo concludere in base alla definizione di misura secondo Peano–Jordan che la
relazione di proporzionalità fra l’area di un insieme e quella della sua immagine vale
per tutti gli insiemi dotati di area.
Definizione
SeT è una funzione lineare, allora, per ogni insieme S, misurabile secondo Peano
Jordan, si ha
𝐴𝑟𝑒𝑎 𝑇 𝑆 = 𝑑𝑒𝑡 𝐽𝑇 ∙ 𝐴𝑟𝑒𝑎(𝑆)

dove 𝐽𝑇 = 𝐴 è la matrice associata, nonché la matrice jacobiana di T.

La matrice jacobiana della funzione

𝑇: 𝑢, 𝑣 ∈ 𝑅 2 ⟼ (𝑥 𝑢, 𝑣 , 𝑦(𝑢, 𝑣)) ∈ 𝑅2
è:
𝜕𝑥 𝜕𝑥
𝜕𝑢 𝜕𝑣
𝐽𝑇 𝑢, 𝑣 = 𝜕𝑦 𝜕𝑦
𝜕𝑢 𝜕𝑣
 Teorema del cambiamento di variabili
Se 𝑇: 𝐴 ⊂ 𝑅2 → 𝑅2 è una trasformazione di classe 𝐶 1 , 𝑇 𝑢, 𝑣 = (𝑥 𝑢, 𝑣 , 𝑦(𝑢, 𝑣) che
sia:
 biunivoca da un aperto misurabile S del piano uv in un aperto misurabile 𝑇 = 𝑇(𝑆) del
piano xy;
 Con jacobiano sempre non nullo e limitato
(una trasformazione che gode di queste proprietà si chiama anche diffeomorfismo)
Allora, ∀𝑓: 𝑇 → 𝑅 integrabile su T, vale la formula:
𝜕(𝑥, 𝑦)
𝑓 𝑥, 𝑦 𝑑𝑥𝑑𝑦 = 𝑓(𝑥 𝑢, 𝑣 , 𝑦 𝑢, 𝑣 𝑑𝑢𝑑𝑣
𝑆 𝜕(𝑢, 𝑣)
𝑇
 Integrali doppi in coordinate polari
Quando il dominio T è un disco, una corona, un settore circolare, per risolvere
l’integrale doppio , conviene utilizzare le coordinate polari. In questo modo le regioni
corrispondono, attraverso il cambiamento di variabili 𝑇 𝜌, 𝜃 = 𝜌𝑐𝑜𝑠𝜃, 𝜌𝑠𝑖𝑛𝜃 a
rettangoli del piano 𝜌𝜃 detti anche rettangoli polari.
La trasformazione in coordinata polari ha jacobiano:

𝜕(𝑥, 𝑦) 𝑐𝑜𝑠𝜃 −𝜌𝑠𝑖𝑛𝜃

= =𝜌
𝜕(𝜌, 𝜃) 𝑠𝑖𝑛𝜃 𝜌𝑐𝑜𝑠𝜃

Lo jacobiano è dunque una funzione limitata

sui limitati del piano 𝜌𝜃 ed è diverso da zero
tranne nei punti con 𝜌 = 0. La trasformazione è
biunivoca, tranne nei punti del piano 𝜌𝜃 con
𝜌 = 0, che vengono tutti mandati nell’origine
del piano xy. Dunque sono soddisfatte le
ipotesi del Teorema relativo al cambiamento di
variabili nell’integrale.
 Teorema del passaggio a coordinate
polari
Sia 𝑇 𝜌, 𝜃 = 𝜌𝑐𝑜𝑠𝜃, 𝜌𝑠𝑖𝑛𝜃 . Sia 𝑆 ⊂ 0, +∞ × 0, 2𝜋 un aperto misurabile del
piano 𝜌𝜃 e sia 𝑇 = 𝑇(𝑆).
Allora ∀𝐹: 𝑇 → 𝑅 integrabile su T, vale la relazione

𝑓 𝑥, 𝑦 𝑑𝑥𝑑𝑦 = 𝑓 𝜌𝑐𝑜𝑠𝜃, 𝜌𝑠𝑖𝑛𝜃 𝜌𝑑𝜌𝑑𝜃

𝑇
𝑆
 Relazione tra integrali curvilinei e integrali
doppi (formula di Gauss – Green nel piano)
Definizione
Se D è un dominio limitato in 𝑅2 , la cui frontiera 𝜕𝐷 sia una curva di Jordan regolare a
tratti, si dice che 𝜕𝐷 è orientata positivamente se è orientata in senso antiorario e si
denoterà col simbolo +𝜕𝐷

Lemma
Sia D un dominio semplice (o normale) rispetto agli assi e 𝐹 𝑥, 𝑦 = 𝑃 𝑥, 𝑦 +
𝑄 𝑥, 𝑦 un campo vettoriale definito dalla chiusura di D di classe 𝐶 1 :
Se 𝐷 = 𝑥, 𝑦 ∈ 𝑅2 : 𝑎 < 𝑥 < 𝑏, 𝜑1 𝑥 < 𝑦 < 𝜑2 (𝑥) con 𝜑1 𝑥 𝑒 𝜑2 (𝑥) regolari a
tratti, allora

𝑃𝑦 𝑑𝑥𝑑𝑦 = − 𝑃𝑑𝑥
𝐷 +𝜕𝐷
Se 𝐷 = 𝑥, 𝑦 ∈ 𝑅2 : 𝑐 < 𝑦 < 𝑑, 𝜓1 𝑦 < 𝑥 < 𝜓2 (𝑦) con 𝜓1 𝑦 𝑒 𝜓2 (𝑦) regolari a
tratti, allora

𝑄𝑥 𝑑𝑥𝑑𝑦 = 𝑄𝑑𝑦
𝐷 +𝜕𝐷
 Dimostrazione

Ricordando il teorema di riduzione per integrali doppi:

𝑏 𝜑2 (𝑥) 𝑏

𝑃𝑦 𝑑𝑥𝑑𝑦 = 𝑑𝑥 𝑃𝑦 𝑑𝑦 = 𝑃 𝑥, 𝜑2 (𝑥 − 𝑃 𝑥, 𝜑1 (𝑥 𝑑𝑥
𝐷 𝑎 𝜑1 (𝑥) 𝑎
Essendo 𝑑𝑥 = 0 sui tratti rettilinei del bordo:
𝑏 𝑎 𝑏

𝑃𝑑𝑥 = 𝑃 𝑥, 𝜑1 (𝑥) 𝑑𝑥 + 𝑃 𝑥, 𝜑2 𝑥 𝑑𝑥 = 𝑃 𝑥, 𝜑1 (𝑥) − 𝑃 𝑥, 𝜑2 (𝑥)

+𝜕𝐷 𝑎 𝑏 𝑎
Dalle due relazioni soprascritte, si ricava

𝑃𝑦 𝑑𝑥𝑑𝑦 = − 𝑃𝑑𝑥
𝐷 +𝜕𝐷
Analogamente si dimostra la seconda delle due formule.
 Integrale triplo
Il caso più semplice è quello di una funzione 𝑓(𝑥, 𝑦, 𝑧) definita su un parallelepipedo R
= [a, b] × [c, d] × [r, s] .
Per definire l’integrale triplo iniziamo col sezionare il parallelepipedo R mediante piani
paralleli a ciascuno dei piani coordinati. A tale fine consideriamo una partizione
dell’intervallo [a, b], 𝑎 = 𝑥0 < 𝑥1 < 𝑥2 < ⋯ < 𝑥𝑛 = 𝑏, una dell’intervallo [c,d],
𝑐 = 𝑦0 < 𝑦1 < 𝑦2 < · · · < 𝑦𝑚 = 𝑑 e infine una dell’intervallo [r, s],
r = 𝑧0 < 𝑧1 < 𝑧2 <· · ·< 𝑧𝑙 =s. Si ottengono in questo modo n × m × l parallelepipedi
𝑅𝑖𝑗𝑘 = 𝑥𝑖−1 , 𝑥𝑖 × 𝑦𝑗−1 , 𝑦𝑗 × 𝑧𝑘−1 , 𝑧𝑘 ciascuno di volume
𝑉𝑖𝑗𝑘 = 𝑥𝑖−1 − 𝑥𝑖 𝑦𝑗−1 − 𝑦𝑗 𝑧𝑘−1 − 𝑧𝑘 e diagonale
2 2 2
𝛿𝑖𝑗𝑘 = 𝑥𝑖−1 − 𝑥𝑖 + 𝑦𝑗−1 − 𝑦𝑗 + 𝑧𝑘−1 − 𝑧𝑘
 Integrale triplo
Scegliamo ora un punto base 𝑥𝑖𝑗 , 𝑦𝑖𝑗 , 𝑧𝑖𝑗 in ognuno dei parallelepipedi 𝑅𝑖𝑗 :
possiamo così definire la tripla somma di Riemann:
𝑛 𝑚 𝑙

𝑓 𝑥𝑖𝑗 , 𝑦𝑖𝑗 , 𝑧𝑖𝑗 𝑉𝑖𝑗𝑘

𝑖=1 𝑗=1 𝑘=1

Definizione
L’integrale triplo di f sul parallelepipedo R è il limite delle somme di Riemann al
tendere a zero dell’ampiezza della partizione:

𝑛 𝑚 𝑙

𝑓 𝑥, 𝑦, 𝑧 𝑑𝑥𝑑𝑦𝑑𝑧 = lim 𝑓 𝑥𝑖𝑗 , 𝑦𝑖𝑗 , 𝑧𝑖𝑗 𝑉𝑖𝑗𝑘

𝛿𝑃 →0
𝑅 𝑖=1 𝑗=1 𝑘=1

quando tale limite esiste.

 Integrali tripli
Come nel caso delle funzioni di due variabili, l’integrale si può definire attraverso il
seguente teorema:
Una funzione limitata è integrabile sul rettangolo se e soltanto se l’estremo superiore
delle somme inferiori, fra tutte partizioni del rettangolo, è uguale all’estremo inferiore
delle somme superiori fra le partizioni. Tali estremi sono a loro volta uguali all’integrale
doppio.
Si può asserire che l’integrale esiste sicuramente quando f è continua sul parallelepipedo
e che valgono le formule di riduzione analoghe al quelle per gli integrali doppi.
Tutte le proprietà degli integrali doppi hanno un analogo per gli integrali tripli, così come
le applicazioni in campo fisico.
Come per gli integrali doppi, definiamo l’integrale di una funzione limitata f : E → R su
una regione limitata𝐸 ⊂ 𝑅3 come l’integrale su un parallelepipedo 𝑅 ⊃ 𝐸 della funzione
f che coincide con f su E e vale zero al di fuori di E. Il volume di Peano–Jordan
tridimensionale di un insieme E come l’integrale su E della funzione costante uno:

𝑉𝑜𝑙𝑢𝑚𝑒 𝐸 = 𝑑𝑥𝑑𝑦𝑑𝑧
𝐸
sempre che tale integrale esista (in tale caso diciamo che E è misurabile secondo Peano–
Jordan). Valgono nel caso tridimensionale le considerazioni teoriche del tutto analoghe a
quelle già fatte nel caso della misura piana,
 Formule di riduzione per gli integrali tripli
Per il calcolo concreto dell’integrale in tre variabili, abbiamo diversi modi di ordinare le
variabili: a ciascuna di esse corrisponde una formula di riduzione a tre integrazioni
semplici successive. Inoltre possiamo interpretare tali formule di riduzione pensando
di eseguire dapprima una integrazione doppia e successivamente una semplice (o
viceversa). Per gli integrali tripli abbiamo pertanto tre tipi diversi di formule di
riduzione:

Teorema
Se f è continua sul parallelepipedo R = [a, b]×[c, d]×[r, s] allora:
𝑠 𝑏

𝑓 𝑥, 𝑦, 𝑧 𝑑𝑥𝑑𝑦𝑑𝑧 = 𝑓 𝑥, 𝑦, 𝑧 𝑑𝑧 𝑑𝑥𝑑𝑦 = 𝑓 𝑥, 𝑦, 𝑧 𝑑𝑦𝑑𝑧 𝑑𝑥

𝑎,𝑏 ×[𝑐,𝑑]
𝑅 𝑟 𝑎 𝑐,𝑑 ×[𝑟,𝑠]
𝑏 𝑑 𝑠

= 𝑓 𝑥, 𝑦, 𝑧 𝑑𝑧𝑑𝑦𝑑𝑥
𝑎 𝑐 𝑟
 Integrazione per fili e per strati

Il significato del primo integrale iterato è facile da comprendere: si esegue prima

l’integrazione più interna, rispetto alla variabile z, tenendo fissate x e y. Il risultato è
una funzione delle due variabili x e y,
𝑠

𝐿 𝑥, 𝑦 = 𝑓 𝑥, 𝑦, 𝑧 𝑑𝑧
𝑟
che verrà poi integrata in un secondo tempo.
Poiché la prima integrazione avviene lungo il segmento in cui x e y sono fissati e z
varia fra r e s tale metodo di calcolo prende il nome di integrazione per fili.

Nell’integrazione per strati invece, si esegue prima un integrale doppio sul rettangolo
del piano in cui x è fissata, {x} × [c, d] × [r, s], ottenendo così una funzione della sola
variabile x, che viene successivamente integrata sull’intervallo [a, b].
In 𝑅3 ci sono tre tipi di regioni semplici, in corrispondenza
ai tre assi coordinati.

Definizione.
Una regione solida E è detta z–semplice se è compresa tra i
grafici di due funzione continue di x e y, cioè è di tipo

𝐸= 𝑥, 𝑦, 𝑧 ∈ 𝑅3 : 𝑥, 𝑦 ∈ 𝐷, 𝑢1 (𝑥, 𝑦) ≤ 𝑧 ≤ 𝑢2 (𝑥, 𝑦)

dove 𝐷 ⊂ 𝑅2 , proiezione di E sul piano z = 0, è un insieme

piano semplice oppure semplicemente decomponibile.
 Cambiamento di variabile negli
integrali tripli
Anche le considerazioni relative al cambiamento di variabili negli integrali doppi si
possono estendere agli integrali tripli (e a integrali di ordine più elevato). Consideriamo
una trasformazione T di classe 𝐶 1 dallo spazio 𝑅3 in sé, descritta dalle relazioni
𝑥 𝑢 𝑥(𝑢, 𝑣, 𝑤)
𝑦 = 𝑇 𝑣 = 𝑦(𝑢, 𝑣, 𝑤)
𝑧 𝑤 𝑧(𝑢, 𝑣, 𝑤)
Definiamo lo jacobiano di T come il determinante della matrice jacobiana
𝜕𝑥 𝜕𝑥 𝜕𝑥
𝜕𝑢 𝜕𝑣 𝜕𝑤
𝜕(𝑥, 𝑦, 𝑧) 𝜕𝑦 𝜕𝑦 𝜕𝑦
= 𝑑𝑒𝑡𝐽𝑇 = 𝑑𝑒𝑡
𝜕(𝑢, 𝑣, 𝑤) 𝜕𝑢 𝜕𝑣 𝜕𝑤
𝜕𝑧 𝜕𝑧 𝜕𝑧
𝜕𝑢 𝜕𝑣 𝜕𝑤

Per le trasformazioni del piano, lo jacobiano rappresenta il fattore (puntuale) di

trasformazione delle aree secondo T; analogamente, in dimensione tre lo jacobiano
rappresenta il fattore infinitesimale di trasformazione dei volumi.
 Teorema del cambiamento di variabili per
gli integrali tripli
Se T sé biunivoca da un aperto misurabile S in un aperto misurabile T = T(S), (la
misurabilità è intesa nel senso della teoria di Peano–Jordan in dimensione tre),
con determinante jacobiano sempre non nullo e limitato. Allora, per ogni f : S → T
continua, vale la formula

𝑓 𝑥, 𝑦, 𝑧 𝑑𝑥𝑑𝑦𝑑𝑧 = 𝑓(𝑥 𝑢, 𝑣, 𝑤 , 𝑦 𝑢, 𝑣, 𝑤 , 𝑧(𝑢, 𝑣, 𝑤) 𝐽 𝑑𝑢𝑑𝑣𝑑𝑤

𝑅 𝑆
In 𝑅3 le coordinate cilindriche (r, 𝜃, z) sono legate alle coordinate cartesiane dalle
relazioni
𝑥 = 𝑟𝑐𝑜𝑠𝜃
𝑦 = 𝑟𝑠𝑖𝑛𝜃
𝑧=𝑧
dove 𝑟 ≥ 0, 𝜃 ∈ 0,2𝜋 , 𝑧 ∈ 𝑅
Lo Jacobiano associato alla trasformazione è:

𝜕(𝑥, 𝑦, 𝑧) 𝑐𝑜𝑠𝜃 −𝑟𝑠𝑖𝑛𝜃 0

= 𝑑𝑒𝑡 𝑠𝑖𝑛𝜃 𝑟𝑐𝑜𝑠𝜃 0 = 𝑟
𝜕(𝑟, 𝜃, 𝑧)
0 0 1
L’uso delle coordinate cilindriche è particolarmente utile per trattare domini dello spazio
con simmetria assiale rispetto all’asse z.
 Coordinate sferiche
Un altro sistema di coordinate in 𝑅3 è quello delle
coordinate sferiche, a in cui ogni punto P è
rappresentato dalla terna (r, φ, ϑ):

𝑥 = 𝑟𝑠𝑖𝑛φ𝑐𝑜𝑠𝜃
𝑦 = 𝑟𝑠𝑖𝑛φ𝑠𝑖𝑛𝜃
𝑧 = 𝑟𝑐𝑜𝑠φ

Il valore di r è la distanza dall’origine di P, φ (la colatitudine) è l’angolo che la retta

passante per O e P forma col semiasse positivo delle z; infine ϑ (la longitudine) è la
coordinata angolare della proiezione del segmento OP sul piano z = 0. ϑ ha quindi lo
stesso significato che nelle coordinate cilindriche. Si ha r ≥ 0, φ ∈ [0, π] e ϑ ∈ [0, 2π).