Calcolare ∬Dydxdy, dove D è la regione delimitata dalla curva γ(t)=(tcost,sint), t∈[0,2π], e dall'asse y.
Soluzione
Geometria. Per t=0 si ha γ=(0,0); per t=2π si ha γ=(0,1). In mezzo x=tcost>0 e y=sint cresce da 0 a 1. Dunque D sta fra l'asse y (a sinistra) e la curva (a destra), per 0≤y≤1.
Metodo 1: riduzione (integrazione per fili orizzontali)
Poiché y=sint è invertibile, t=arcsiny e lungo la curva x=arcsin(y)1−y2 (da cost=1−y2). Quindi
D={0≤y≤1,0≤x≤arcsiny1−y2},∬Dydxdy=∫01yarcsiny1−y2dy.
Con y=sint, dy=costdt: l'integrale diventa ∫0π/2tsintcos2tdt. Per parti, con u=t e v′=sintcos2t (quindi v=−31cos3t): il termine di bordo [−3tcos3t]0π/2 è nullo e resta 31∫0π/2cos3tdt=31[sint−3sin3t]0π/2=31⋅32=92.
Metodo 2: Gauss-Green
Si sceglie F1=−2y2, F2=0, così che ∂xF2−∂yF1=y e
∬Dydxdy=∮∂+D(−2y2)dx.
Il bordo positivo ha D a sinistra: si percorre la curva da (0,0) a (0,1) (la regione sta a sinistra, verso l'asse y) e poi l'asse y da (0,1) a (0,0). Sull'asse y, x=0 e dx=0: nessun contributo. Sulla curva dx=(cost−tsint)dt e y=sint:
−21∫0π/2sin2t(cost−tsint)dt=−21∫0π/2sin2tcostdt+21∫0π/2tsin3tdt.
Il primo: −21[3sin3t]0π/2=−61.
Il secondo: da sin3t=43sint−sin3t e dalle integrazioni per parti ∫0π/2tsintdt=1, ∫0π/2tsin3tdt=−91 si ottiene ∫0π/2tsin3tdt=43+1/9=97, quindi il secondo pezzo vale 187.
Somma: 187−61=187−3=92.
Risultato.∬Dydxdy=92, con entrambi i metodi (e sympy).
Versione ripasso
Testo.∬Dydxdy, D delimitata da (tcost,sint), t∈[0,2π], e dall'asse y.
Metodo 1.D={0≤y≤1,0≤x≤arcsiny1−y2}; con y=sint: ∫0π/2tsintcos2tdt=31∫0π/2cos3tdt=92 (per parti).