Salta al contenuto
Note per Studenti Formula di Taylor con resto di Peano

Formula di Taylor con resto di Peano

In questa pagina 7

Capitolo 4, "Calcolo differenziale" (lezione 31). Serve a calcolare limiti difficili, a trovare l'ordine di infinitesimi e infiniti e a calcolare derivate alte senza derivare. Tabella degli sviluppi già pronti: Sviluppi di Mac-Laurin notevoliTabella degli sviluppi di Mac-Laurin da sapere a memoria (e^x, sin, cos, log(1+x), (1+x)^alpha, arctan, sinh, cosh, tan) e regole per combinarli: algebra degli o piccoli, prodotti, funzioni composte, quanti termini tenere.Sviluppi di Mac-Laurin notevoli →. Versione con l'errore scritto esplicitamente: Formula di Taylor con resto di LagrangeSe f ha n+1 derivate continue su [a,b], l'errore tra f e il suo polinomio di Taylor si scrive esattamente come f^(n+1)(c)(x-x0)^(n+1)/(n+1)! per un c opportuno: serve a stimare l'errore su tutto un intervallo e a scegliere il grado del polinomio.Formula di Taylor con resto di Lagrange →. Esercizi: Esercizio 102 · sviluppo di Mac-Laurin di log(cos x), Esercizio 103 · sviluppo di Mac-Laurin di (1+2x)sinh x, Esercizio 104 · derivate quinta e sesta in 0 senza derivare.

Cosa serve sapere prima. Le derivate (anche successive: f(2)=f′′f^{(2)} = f'', f(3)=f′′′f^{(3)} = f''', …), il simbolo o(⋅)o(\cdot) ("o piccolo")Per x→x0x \to x_0, g(x)=o(h(x))g(x) = o(h(x)) significa che g(x)/h(x)→0g(x)/h(x) \to 0: gg diventa trascurabile rispetto a hh. Per esempio x3=o(x2)x^3 = o(x^2) per x→0x \to 0., il simbolo ∼\sim (asintotico)g∼hg \sim h per x→x0x \to x_0 significa g(x)/h(x)→1g(x)/h(x) \to 1: le due funzioni, vicino a x0x_0, sono "quasi uguali". e il teorema di de l'HôpitalSe f/gf/g è una forma 0/00/0 (o ∞/∞\infty/\infty) e il limite di f′/g′f'/g' esiste, allora il limite di f/gf/g è lo stesso. Nei passaggi della prof è segnato con una H cerchiata..

La domanda

Data ff derivabile nn volte in x0=0x_0 = 0, si può trovare un polinomio Tn(x)T_n(x) di grado ≤n\le n tale che

f(x)=Tn(x)+o(xn)per x→0 ?f(x) = T_n(x) + o(x^n) \qquad \text{per } x \to 0 \ ?

Cioè: vicino a 00, ff è uguale a un polinomio a meno di un errore che va a zero più in fretta di xnx^n. I polinomi sono le funzioni più facili da trattare (si sommano, si moltiplicano, si derivano a mano), quindi sostituire ff con un polinomio semplifica moltissimo i conti.

Il caso n=1n = 1 lo conosciamo già. Che ff sia derivabile in 00 equivale a dire che

f(x)=f(0)+f′(0) x+o(x)per x→0,f(x) = f(0) + f'(0)\,x + o(x) \qquad \text{per } x \to 0,

cioè che la retta tangente T1(x)=f(0)+f′(0) xT_1(x) = f(0) + f'(0)\,x approssima ff con un errore o(x)o(x) (è la caratterizzazione della derivabilità con l'o piccolo). La formula di Taylor è lo stesso discorso con polinomi di grado più alto, che approssimano meglio.

L'esempio che fa capire: f(x)=sin⁡xf(x) = \sin x

n=1n = 1. sin⁡0=0\sin 0 = 0 e cos⁡0=1\cos 0 = 1, quindi sin⁡x=x+o(x)\sin x = x + o(x): T1(x)=xT_1(x) = x.

n=2n = 2. Vogliamo capire quanto vale l'errore sin⁡x−T1(x)=sin⁡x−x\sin x - T_1(x) = \sin x - x confrontato con x2x^2. Calcoliamo il limite del rapporto; è una forma 00\frac00, applichiamo de l'Hôpital (deriviamo sopra e sotto):

lim⁡x→0sin⁡x−xx2=Hlim⁡x→0cos⁡x−12x=Hlim⁡x→0−sin⁡x2=0\lim_{x \to 0} \frac{\sin x - x}{x^2} \overset{H}{=} \lim_{x \to 0} \frac{\cos x - 1}{2x} \overset{H}{=} \lim_{x \to 0} \frac{-\sin x}{2} = 0

Il rapporto tende a 00: per definizione di o piccolo, sin⁡x−x=o(x2)\sin x - x = o(x^2). Quindi sin⁡x=x+o(x2)\sin x = x + o(x^2) e T2(x)=T1(x)=xT_2(x) = T_1(x) = x: il termine di grado 2 non c'è (ha coefficiente 00). Osserviamo che 0=12f′′(0)0 = \frac12 f''(0), perché f′′(x)=−sin⁡xf''(x) = -\sin x vale 00 in 00.

n=3n = 3. Ora confrontiamo l'errore con x3x^3:

lim⁡x→0sin⁡x−xx3=Hlim⁡x→0cos⁡x−13x2=Hlim⁡x→0−sin⁡x3⋅2 x=Hlim⁡x→0−cos⁡x3!=−13!\lim_{x \to 0} \frac{\sin x - x}{x^3} \overset{H}{=} \lim_{x \to 0} \frac{\cos x - 1}{3x^2} \overset{H}{=} \lim_{x \to 0} \frac{-\sin x}{3 \cdot 2\, x} \overset{H}{=} \lim_{x \to 0} \frac{-\cos x}{3!} = -\frac{1}{3!}

Il limite è un numero diverso da zero: vuol dire che sin⁡x−x∼−13!x3\sin x - x \sim -\frac{1}{3!}x^3, cioè

sin⁡x−x=−13!x3+o(x3)⟹sin⁡x=x−13!x3⏟T3(x)+o(x3)\sin x - x = -\frac{1}{3!}x^3 + o(x^3) \quad\Longrightarrow\quad \sin x = \underbrace{x - \frac{1}{3!}x^3}_{T_3(x)} + o(x^3)

E di nuovo il coefficiente trovato è −13!=13!f′′′(0)-\frac{1}{3!} = \frac{1}{3!}f'''(0), perché f′′′(x)=−cos⁡xf'''(x) = -\cos x vale −1-1 in 00.

Lo schema che emerge. Ogni volta il nuovo coefficiente di xkx^k è 1k!f(k)(0)\frac{1}{k!}f^{(k)}(0). Andando avanti (confrontando con x4x^4, x5x^5, …) si continua allo stesso modo. Questo è il contenuto del teorema.

Teorema (formula di Taylor in x0=0x_0 = 0, o sviluppo di Mac-Laurin)

Se ff è derivabile nn volte in x0=0x_0 = 0, esiste uno e un solo polinomio Tn(x)T_n(x) di grado ≤n\le n tale che f(x)=Tn(x)+o(xn)per x→0,f(x) = T_n(x) + o(x^n) \qquad \text{per } x \to 0, ed è Tn(x)=f(0)+f′(0) x+12f′′(0) x2+13!f′′′(0) x3+⋯+1n!f(n)(0) xn=∑k=0nf(k)(0)k! xkT_n(x) = f(0) + f'(0)\,x + \frac{1}{2}f''(0)\,x^2 + \frac{1}{3!}f'''(0)\,x^3 + \dots + \frac{1}{n!}f^{(n)}(0)\,x^n = \sum_{k=0}^{n} \frac{f^{(k)}(0)}{k!}\,x^k (con la convenzione f(0)=ff^{(0)} = f).

Nomi.

  • Tn(x)T_n(x) è il polinomio di Taylor di ordine nn di ff (in 00; detto anche polinomio di Mac-Laurin).
  • Rn(x)=f(x)−Tn(x)=o(xn)R_n(x) = f(x) - T_n(x) = o(x^n) è il resto nn-esimo di Peano: è l'errore che si commette sostituendo ff con TnT_n, e la formula dice solo che è trascurabile rispetto a xnx^n, senza dire quanto vale.
  • TnT_n è anche l'unico polinomio di grado ≤n\le n che ha le stesse derivate di ff in 00 fino all'ordine nn: Tn(k)(0)=f(k)(0)T_n^{(k)}(0) = f^{(k)}(0) per k=0,1,…,nk = 0, 1, \dots, n.

Come si legge la formula. Il coefficiente di xkx^k è "derivata kk-esima in 00 divisa per k!k!Il fattoriale: k!=1⋅2⋯kk! = 1 \cdot 2 \cdots k, con 0!=10! = 1. Per esempio 3!=63! = 6, 4!=244! = 24, 5!=1205! = 120.Fattoriale e coefficienti binomiali →". Il fattoriale cresce molto in fretta, quindi i coefficienti di solito diventano piccoli.

Dimostrazione (per induzione su nn)

Il ragionamento è quello fatto sopra con sin⁡x\sin x, ripetuto per una ff qualsiasi.

n=1n = 1. È la caratterizzazione della derivabilità: f(x)=f(0)+f′(0) x+o(x)f(x) = f(0) + f'(0)\,x + o(x), quindi T1(x)=f(0)+f′(0) xT_1(x) = f(0) + f'(0)\,x.

n=2n = 2. Si calcola

lim⁡x→0f(x)−T1(x)x2=lim⁡x→0f(x)−f(0)−f′(0) xx2=Hlim⁡x→0f′(x)−f′(0)2x=12f′′(0)\lim_{x \to 0} \frac{f(x) - T_1(x)}{x^2} = \lim_{x \to 0} \frac{f(x) - f(0) - f'(0)\,x}{x^2} \overset{H}{=} \lim_{x \to 0} \frac{f'(x) - f'(0)}{2x} = \frac12 f''(0)

Perché si può usare de l'Hôpital: siccome esiste f′′(0)f''(0), le funzioni ff e f′f' sono derivabili (quindi continue) in 00; allora sopra si ha f(x)−f(0)−f′(0)x→0f(x) - f(0) - f'(0)x \to 0 e sotto x2→0x^2 \to 0, forma 00\frac00. L'ultimo limite è 12\frac12 per il rapporto incrementale di f′f'f′(x)−f′(0)x−0\frac{f'(x) - f'(0)}{x - 0}: il suo limite per x→0x \to 0 è, per definizione, la derivata di f′f' in 00, cioè f′′(0)f''(0)., che tende a f′′(0)f''(0) per definizione di derivata seconda.

Quindi f(x)−T1(x)∼12f′′(0) x2f(x) - T_1(x) \sim \frac12 f''(0)\,x^2 (se f′′(0)≠0f''(0) \neq 0; se è 00 il limite dice comunque che la differenza è o(x2)o(x^2)). In entrambi i casi:

f(x)−T1(x)=12f′′(0) x2+o(x2) ⟹ f(x)=f(0)+f′(0) x+12f′′(0) x2+o(x2)f(x) - T_1(x) = \frac12 f''(0)\,x^2 + o(x^2) \ \Longrightarrow\ f(x) = f(0) + f'(0)\,x + \frac12 f''(0)\,x^2 + o(x^2)

n=3n = 3. Stessa cosa con un passaggio di de l'Hôpital in più (esiste f′′′(0)f'''(0), quindi f,f′,f′′f, f', f'' sono continue in 00 e ogni volta si ha 00\frac00):

lim⁡x→0f(x)−f(0)−f′(0) x−12f′′(0) x2x3=Hlim⁡x→0f′(x)−f′(0)−f′′(0) x3x2=Hlim⁡x→0f′′(x)−f′′(0)3⋅2 x=13!f′′′(0)\lim_{x \to 0} \frac{f(x) - f(0) - f'(0)\,x - \frac12 f''(0)\,x^2}{x^3} \overset{H}{=} \lim_{x \to 0} \frac{f'(x) - f'(0) - f''(0)\,x}{3x^2} \overset{H}{=} \lim_{x \to 0} \frac{f''(x) - f''(0)}{3 \cdot 2\,x} = \frac{1}{3!}f'''(0)

e quindi f(x)=T2(x)+13!f′′′(0) x3+o(x3)f(x) = T_2(x) + \frac{1}{3!}f'''(0)\,x^3 + o(x^3).

Passo induttivo. Si suppone che Tn−1(x)=f(0)+f′(0) x+⋯+1(n−1)!f(n−1)(0) xn−1T_{n-1}(x) = f(0) + f'(0)\,x + \dots + \frac{1}{(n-1)!}f^{(n-1)}(0)\,x^{n-1} e si calcola, applicando de l'Hôpital n−1n - 1 volte,

lim⁡x→0f(x)−Tn−1(x)xn=⋯=1n!f(n)(0)\lim_{x \to 0} \frac{f(x) - T_{n-1}(x)}{x^n} = \dots = \frac{1}{n!}f^{(n)}(0)

da cui f(x)=Tn−1(x)+1n!f(n)(0) xn+o(xn)=Tn(x)+o(xn)f(x) = T_{n-1}(x) + \frac{1}{n!}f^{(n)}(0)\,x^n + o(x^n) = T_n(x) + o(x^n) ∎.

Perché il polinomio è unico (quello con le stesse derivate)

Prendiamo un polinomio qualsiasi di grado ≤n\le n, Pn(x)=a0+a1x+a2x2+⋯+anxnP_n(x) = a_0 + a_1 x + a_2 x^2 + \dots + a_n x^n, e imponiamo che abbia in 00 le stesse derivate di ff: Pn(k)(0)=f(k)(0)P_n^{(k)}(0) = f^{(k)}(0) per k=0,…,nk = 0, \dots, n. Vediamo che cosa obbliga a fare.

  • k=0k = 0: Pn(0)=a0P_n(0) = a_0 (tutti gli altri termini contengono xx e si annullano). Quindi a0=f(0)a_0 = f(0).
  • k=1k = 1: Pn′(x)=a1+2a2x+3a3x2+⋯+n anxn−1P_n'(x) = a_1 + 2a_2 x + 3a_3 x^2 + \dots + n\,a_n x^{n-1}, quindi Pn′(0)=a1P_n'(0) = a_1. Imponendo a1=f′(0)a_1 = f'(0).
  • k=2k = 2: Pn′′(x)=2a2+3⋅2 a3x+…P_n''(x) = 2a_2 + 3 \cdot 2\, a_3 x + \dots, quindi Pn′′(0)=2a2=f′′(0)P_n''(0) = 2a_2 = f''(0), cioè a2=12f′′(0)a_2 = \frac12 f''(0).
  • k=3k = 3: Pn′′′(0)=3⋅2⋅a3=3! a3=f′′′(0)P_n'''(0) = 3 \cdot 2 \cdot a_3 = 3!\,a_3 = f'''(0), cioè a3=13!f′′′(0)a_3 = \frac{1}{3!}f'''(0).
  • In generale (per induzione): derivando kk volte akxka_k x^k resta k! akk!\,a_k, i termini di grado più basso spariscono (derivati troppe volte) e quelli di grado più alto contengono ancora xx (valgono 00 in 00). Quindi ak=1k!f(k)(0)a_k = \frac{1}{k!}f^{(k)}(0).

I coefficienti sono obbligati: Pn=TnP_n = T_n ∎.

Conseguenza pratica fondamentale. Se in qualunque modo (anche senza derivare!) si riesce a scrivere

f(x)=a0+a1x+⋯+anxn+o(xn)per x→0,f(x) = a_0 + a_1 x + \dots + a_n x^n + o(x^n) \qquad \text{per } x \to 0,

allora quello è il polinomio di Taylor, e quindi ak=f(k)(0)k!a_k = \frac{f^{(k)}(0)}{k!}. È il motivo per cui gli sviluppi si costruiscono componendo quelli noti invece di derivare: Esercizio 102 · sviluppo di Mac-Laurin di log(cos x), e per cui si possono ricavare le derivate dai coefficienti: f(k)(0)=k! akf^{(k)}(0) = k!\,a_k (Esercizio 104 · derivate quinta e sesta in 0 senza derivare).

Primi esempi

f(x)=exf(x) = e^x

Tutte le derivate di exe^x sono exe^x, e e0=1e^0 = 1: f(k)(0)=1f^{(k)}(0) = 1 per ogni kk. Quindi ogni coefficiente è 1k!\frac{1}{k!}:

ex=1+x+x22+x33!+⋯+xnn!+o(xn)per x→0e^x = 1 + x + \frac{x^2}{2} + \frac{x^3}{3!} + \dots + \frac{x^n}{n!} + o(x^n) \qquad \text{per } x \to 0

Verifica numerica con x=0,1x = 0{,}1: T3(0,1)=1+0,1+0,005+0,000167=1,105167T_3(0{,}1) = 1 + 0{,}1 + 0{,}005 + 0{,}000167 = 1{,}105167, mentre e0,1=1,105171…e^{0{,}1} = 1{,}105171\ldots: errore circa 4⋅10−64 \cdot 10^{-6}, molto più piccolo di x3=0,001x^3 = 0{,}001.

f(x)=sin⁡xf(x) = \sin x

Le derivate si ripetono ogni quattro:

kk f(k)(x)f^{(k)}(x) f(k)(0)f^{(k)}(0)
00 sin⁡x\sin x 00
11 cos⁡x\cos x 11
22 −sin⁡x-\sin x 00
33 −cos⁡x-\cos x −1-1
44 sin⁡x\sin x 00
  • Derivate pari (k=2hk = 2h): valgono ±sin⁡0=0\pm\sin 0 = 0, quindi le potenze pari non compaiono.
  • Derivate dispari (k=2h+1k = 2h + 1): valgono alternativamente +1+1 (h=0h = 0), −1-1 (h=1h = 1), +1+1, …, cioè (−1)h(-1)^h.

sin⁡x=x−x33!+x55!−⋯+(−1)nx2n+1(2n+1)!+o(x2n+1)\sin x = x - \frac{x^3}{3!} + \frac{x^5}{5!} - \dots + (-1)^n \frac{x^{2n+1}}{(2n+1)!} + o(x^{2n+1})

Finezza utile. Il termine successivo, di grado 2n+22n + 2, ha coefficiente 00 (è una derivata pari). Quindi si può scrivere anche + 0⋅x2n+2+o(x2n+2)+\,0 \cdot x^{2n+2} + o(x^{2n+2}): lo stesso polinomio approssima con errore o(x2n+2)o(x^{2n+2}), un ordine meglio. Per esempio sin⁡x=x−x36+o(x4)\sin x = x - \frac{x^3}{6} + o(x^4). Vale per ogni funzione dispari (sviluppo con sole potenze dispari); per una funzione pari (come cos⁡x\cos x) ci sono solo potenze pari.

Gli altri sviluppi da sapere (cos⁡x\cos x, log⁡(1+x)\log(1+x), (1+x)α(1+x)^\alpha, arctan⁡x\arctan x, sinh⁡x\sinh x, cosh⁡x\cosh x) sono raccolti in Sviluppi di Mac-Laurin notevoliTabella degli sviluppi di Mac-Laurin da sapere a memoria (e^x, sin, cos, log(1+x), (1+x)^alpha, arctan, sinh, cosh, tan) e regole per combinarli: algebra degli o piccoli, prodotti, funzioni composte, quanti termini tenere.Sviluppi di Mac-Laurin notevoli →.

Teorema (formula di Taylor in un punto x0∈Rx_0 \in \mathbb{R} qualsiasi)

Lo stesso si può fare attorno a qualunque punto, usando le potenze di (x−x0)(x - x_0) al posto di quelle di xx.

Se ff è derivabile nn volte in x0∈Rx_0 \in \mathbb{R}, esiste uno e un solo polinomio Tn,x0(x)T_{n,x_0}(x) di grado ≤n\le n tale che f(x)=Tn,x0(x)+o((x−x0)n)f(x) = T_{n,x_0}(x) + o\big((x - x_0)^n\big) per x→x0x \to x_0, ed è Tn,x0(x)=f(x0)+f′(x0)(x−x0)+12f′′(x0)(x−x0)2+⋯+1n!f(n)(x0)(x−x0)nT_{n,x_0}(x) = f(x_0) + f'(x_0)(x - x_0) + \frac12 f''(x_0)(x - x_0)^2 + \dots + \frac{1}{n!}f^{(n)}(x_0)(x - x_0)^n

Tn,x0T_{n,x_0} è il polinomio di Taylor di ordine nn di ff in x0x_0, Rn,x0(x)=o((x−x0)n)R_{n,x_0}(x) = o\big((x - x_0)^n\big) il resto di Peano; Tn,x0T_{n,x_0} è l'unico polinomio con Tn,x0(k)(x0)=f(k)(x0)T_{n,x_0}^{(k)}(x_0) = f^{(k)}(x_0) per k=0,…,nk = 0, \dots, n.

Attenzione. I polinomi di Taylor in due punti diversi x0≠x1x_0 \neq x_1 sono indipendenti: non si può ricavare l'uno dall'altro in nessun modo (in particolare, non "spostando" quello in 00).

Esempio: sin⁡x\sin x in x0=π4x_0 = \frac{\pi}{4}, ordine 3

Servono le derivate fino alla terza, calcolate in π4\frac\pi4 (dove seno e coseno valgono 22\frac{\sqrt2}{2}, vedi Funzioni trigonometriche e loro inverseSeno, coseno e tangente definiti sulla circonferenza goniometrica (angoli in radianti), le loro proprietà e le inverse arcoseno, arcocoseno e arcotangente.Funzioni trigonometriche e loro inverse →):

  • f(π4)=sin⁡π4=22f(\frac\pi4) = \sin\frac\pi4 = \frac{\sqrt2}{2};
  • f′(x)=cos⁡xf'(x) = \cos x, f′(π4)=22f'(\frac\pi4) = \frac{\sqrt2}{2};
  • f′′(x)=−sin⁡xf''(x) = -\sin x, f′′(π4)=−22f''(\frac\pi4) = -\frac{\sqrt2}{2};
  • f′′′(x)=−cos⁡xf'''(x) = -\cos x, f′′′(π4)=−22f'''(\frac\pi4) = -\frac{\sqrt2}{2}.

Sostituendo (con 12\frac12 e 13!=16\frac{1}{3!} = \frac16 davanti al secondo e terzo termine) e raccogliendo 22\frac{\sqrt2}{2}:

sin⁡x=22(1+(x−π4)−12(x−π4)2−16(x−π4)3)+o ⁣((x−π4)3)per x→π4\sin x = \frac{\sqrt2}{2}\left(1 + \left(x - \tfrac\pi4\right) - \frac12\left(x - \tfrac\pi4\right)^2 - \frac16\left(x - \tfrac\pi4\right)^3\right) + o\!\left(\left(x - \tfrac\pi4\right)^3\right) \quad \text{per } x \to \tfrac\pi4

Errore tipico. Partire da sin⁡x=x−x36+o(x3)\sin x = x - \frac{x^3}{6} + o(x^3) e sostituire xx con x−π4x - \frac\pi4 dà lo sviluppo di sin⁡(x−π4)\sin\left(x - \frac\pi4\right), non quello di sin⁡x\sin x attorno a π4\frac\pi4: è un'altra funzione (il grafico di sin⁡\sin traslato a destra di π4\frac\pi4).

Grafico interattivo: Errore sin x − T₃(x) del polinomio in π/4 ≈ 0,785: è praticamente zero vicino a π/4 e cresce allontanandosi

A che cosa serve

  1. Calcolare limiti. Sostituendo ogni funzione con il suo polinomio di Taylor (fino all'ordine giusto) le forme indeterminate diventano confronti tra potenze. Esempio: lim⁡x→0sin⁡x−xx3=lim⁡x→0−x36+o(x3)x3=−16\displaystyle\lim_{x \to 0}\frac{\sin x - x}{x^3} = \lim_{x \to 0}\frac{-\frac{x^3}{6} + o(x^3)}{x^3} = -\frac16, senza de l'Hôpital.
  2. Trovare l'ordine di infinitesimo di una funzione (il primo termine non nullo dello sviluppo dice a quale potenza è asintotica): Esercizio 106 · ordine di infinito al variare di α.
  3. Calcolare derivate alte in un punto leggendo i coefficienti: Esercizio 104 · derivate quinta e sesta in 0 senza derivare.

È un'approssimazione locale. Il resto di Peano dice che l'errore è piccolo solo per xx vicino a x0x_0 e non dà nessuna stima numerica. Per sapere quanto vale l'errore su un intervallo intero serve la Formula di Taylor con resto di LagrangeSe f ha n+1 derivate continue su [a,b], l'errore tra f e il suo polinomio di Taylor si scrive esattamente come f^(n+1)(c)(x-x0)^(n+1)/(n+1)! per un c opportuno: serve a stimare l'errore su tutto un intervallo e a scegliere il grado del polinomio.Formula di Taylor con resto di Lagrange →.

Errori comuni

Esercizi su questo argomento

Lezioni in cui compare

Teoria collegata