Integrale di una funzione dalla serie di potenze

Prima di postare leggi le regole del Forum. Puoi anche leggere le ultime discussioni.
Questa sezione è un contenitore temporaneo per i topic speciali a pagamento (One-Shot).
#101507
avt
Pippo95n
Punto
Ho bisogno di una mano per approssimare l'integrale definito di una funzione esponenziale usando la relativa serie di potenze.

Partendo dalla serie di potenze della funzione f(x) = e^(x^2) calcolare l'integrale

∫_(0)^((1)/(2))f(x)dx

e approssimare il valore dell'integrale a meno di (1)/(10).

Grazie.
#101508
avt
Ifrit
Amministratore
La traccia ci chiede di determinare l'integrale definito

∫_(0)^((1)/(2))e^(x^2)dx

avvalendoci della serie di potenze associata alla funzione esponenziale f(x) = e^(x^2) dopodiché dovremo approssimare il valore dell'integrale a meno di (1)/(10)


Serie di potenze associata a f(x)

Il primo passo prevede di scrivere la serie di potenze associata a f(x) = e^(x^2) e per farlo è sufficiente usare lo sviluppo in serie dell'esponenziale

e^(t) = Σ_(n = 0)^(+∞)(t^n)/(n!) ∀ t∈R

dopodiché rimpiazziamo x^2 a ogni occorrenza di t, ottenendo così la serie di potenze richiesta:

f(x) = e^(x^2) = Σ_(n = 0)^(+∞)((x^2)^(n))/(n!) = Σ_(n = 0)^(+∞)(x^(2n))/(n!) ∀ x∈R


Calcolo dell'integrale

Noto lo sviluppo di f(x), possiamo integrarla tra 0 e (1)/(2) usando il teorema di integrazione per serie, il quale ci permette di scambiare di posizione l'operatore integrale con quello di serie.

∫_(0)^((1)/(2))f(x)dx = ∫_(0)^((1)/(2))e^(x^2)dx = ∫_(0)^((1)/(2))Σ_(n = 0)^(+∞)(x^(2n))/(n!)dx =

Trasportiamo l'integrale all'interno della serie

= Σ_(n = 0)^(+∞)∫_(0)^((1)/(2))(x^(2n))/(n!)dx =

Poiché il fattoriale di n non dipende da x, può tranquillamente uscire dal simbolo di integrazione

= Σ_(n = 0)^(+∞)(1)/(n!)∫_(0)^((1)/(2))x^(2n)dx = Σ_(n = 0)^(+∞)(1)/(n!)[(x^(2n+1))/(2n+1)]_(x = 0)^(x = (1)/(2)) = Σ_(n = 0)^(+∞)(1)/(n!)·(((1)/(2))^(2n+1))/(2n+1) = Σ_(n = 0)^(+∞)(((1)/(2))^(2n+1))/(n!(2n+1)) =

Volendo potremmo migliorare l'estetica del risultato usando le proprietà delle potenze

= Σ_(n = 0)^(+∞)(1)/(2^(2n+1)n!(2n+1))


Approssimazione dell'integrale

Il secondo punto del problema richiede di determinare l'approssimazione dell'integrale a meno di (1)/(10): in questa circostanza procederemo così:

- dall'uguaglianza

∫_(0)^((1)/(2))f(x)dx = Σ_(n = 0)^(+∞)(1)/(2^(2n+1)n!(2n+1))

è possibile scriverne una completamente equivalente: l'unica cosa che cambia è che abbiamo spezzato la serie nella somma tra i primi N+1 addendi (quelli che vanno da 0 e N) e quella che si ottiene addizionando gli infiniti addendi rimanenti (quelli che vanno da N+1 a +∞)

∫_(0)^((1)/(2))f(x)dx = Σ_(n = 0)^(N)(1)/(2^(2n+1)n!(2n+1))+Σ_(n = N+1)^(+∞)(1)/(2^(2n+1)n!(2n+1))

Se trasportiamo al primo membro la prima sommatoria, otteniamo la relazione

∫_(0)^((1)/(2))f(x)dx-Σ_(n = 0)^(N)(1)/(2^(2n+1)n!(2n+1)) = Σ_(n = N+1)^(+∞)(1)/(2^(2n+1)n!(2n+1))

da cui si evince l'espressione dell'errore, dipendente da N, che si commette nel momento in cui sostituiamo il valore dell'integrale con quello della sommatoria:

R_(N) = Σ_(n = N+1)^(+∞)(1)/(2^(2n+1)n!(2n+1))

Nota: per definizione l'errore è il modulo di R_N. In questo caso però stiamo lavorando con quantità positive sicché il valore assoluto non serve.

Il nostro obiettivo è quello di determinare N in modo che R_N sia minore di (1)/(10).

R_(N) < (1)/(10) → Σ_(n = N+1)^(+∞)(1)/(2^(2n+1)n!(2n+1)) < (1)/(10)

Purtroppo un resto in questa forma è del tutto inutile, perché non siamo in grado di ricavare un'espressione analitica che consenta di valutarlo. Non tutto è perduto: possiamo sempre maggiorarlo!

Poiché n! è crescente, così come lo è 2n+1, possiamo scrivere le seguenti disuguaglianze

n! ≥ (N+1)! , 2n+1 ≥ 2(N+1)+1 = 2N+3

valide per ogni n ≥ N+1, pertanto

n!(2n+1) ≥ (N+1)!(2N+3) ∀ n ≥ N+1

Nel momento in cui passiamo ai reciproci, il verso della disuguaglianza si inverte

(1)/(n!(2n+1)) ≤ (1)/((N+1)!(2N+3)) ∀ n ≥ N+1

Infine, se dividiamo i due membri per 2^(2n+1) scopriamo che:

(1)/(2^(2n+1)n!(2n+1)) ≤ (1)/(2^(2n+1)(N+1)!(2N+3)) ∀ n ≥ N+1

Perché abbiamo fatto tutto questo? Perché questa maggiorazione ci permette di maggiorare R_N con qualcosa di più gestibile. Se applichiamo ai due membri la sommatoria che va da N+1 a +∞ ricaviamo

Σ_(n = N+1)^(+∞)(1)/(2^(2n+1)n!(2n+1)) (R_N) ≤ Σ_(n = N+1)^(+∞)(1)/(2^(2n+1)(N+1)!(2N+3)) =

Attenzione! Nella seconda serie, i termini (N+1)! e (2N+3) non dipendono da n, per cui possono abbandonare il simbolo di sommatoria

= (1)/((N+1)! (2N+3))Σ_(n = N+1)^(+∞)(1)/(2^(2n+1)) =

Sempre grazie alle proprietà delle potenze possiamo infine ricondurci a una serie geometrica di ragione (1)/(4).

= (1)/((N+1)! (2N+3))Σ_(n = N+1)^(+∞)(1)/(2^(2n)·2) = (1)/(2(N+1)! (2N+3))Σ_(n = N+1)^(+∞)(1)/((2^2)^(n)) = (1)/(2(N+1)! (2N+3))Σ_(n = N+1)^(+∞)((1)/(4))^n = (1)/(2(N+1)!(2N+3))·(1)/(3·4^(N)) = (1)/(6·4^(N)(N+1)!(2N+3))

In definitiva abbiamo dimostrato la seguente maggiorazione per il resto R_N

R_N ≤ (1)/(6·4^(N)(N+1)!(2N+3)) ∀ N ≥ 0

Abbiamo praticamente concluso: se determiniamo un N che realizza la disuguaglianza

(1)/(6·4^(N)(N+1)!(2N+3)) < (1)/(10)

a maggior ragione questo valore di N farà sì che sussista anche

R_(N) < (1)/(10)

Per trovare il valore di N, basta procedere per tentativi. Non molti a dire il vero, infatti se N = 0 allora

(1)/(6·4^(N)(N+1)!(2N+3)) = (1)/(6·4^(0) (0+1)!(2·0+3)) = (1)/(18) < (1)/(10)

Ottimo! Abbiamo determinato il numero che ci serve e che va sostituito in

Σ_(n = 0)^(N)(1)/(2^(2n+1)n!(2n+1))

ottenendo

Σ_(n = 0)^(0)(1)/(2^(2n+1)n!(2n+1)) = (1)/(2)


Conclusioni

In definitiva l'integrale approssimato a meno di (1)/(10) vale:

∫_(0)^((1)/(2))e^(x^2)dx ≃ (1)/(2) = 0,5

Abbiamo finito!
Ringraziano: Omega
#101510
avt
Pippo95n
Punto
Innanzitutto grazie!

Non ho capito però se c'è un'espressione generale dell'errore, come faccio a ricondurmi alla serie geometrica di ragione 1/4.

Grazie ancora.
#101511
avt
Ifrit
Amministratore
In generale, esiste una formula che consente di scrivere la somma di una serie geometrica di ragione q, il cui indice parte da un valore diverso da zero ed è questa:

Σ_(n = k)^(+∞)q^(n) = (q^k)/(1-q)

Se prendiamo q = (1)/(4) e k = N+1, otteniamo:

Σ_(n = N+1)^(+∞)((1)/(4))^(n) = (((1)/(4))^(N+1))/(1-(1)/(4)) = ((1)/(4^(N+1)))/((3)/(4)) = (1)/(4^(N+1))·(4)/(3) = (1)/(3·4^(N))
  • Pagina:
  • 1