Soluzione degli esercizi
Primo foglio
(1) Calcola
∑ n=1..∞ 1/((n+x)(n+x+1)(n+x+2)) con Maple, ed esprimi il risultato
usando il comando "factor". Il risultato ottenuto vale per ogni intero x? Prova a dimostrare quanto ottenuto
usando la proprietà telescopica.
Con Maple:
> Sum(1/((n+x)*(n+x+1)*(n+x+2)),n=1..infinity) = sum(1/((n+x)*(n+x+1)*(n+x+2)),n=1..infinity);
∑ n=1..∞ 1/((n+x)(n+x+1)(n+x+2)) = 1/(2·(2+3x+x2))
> factor(sum(1/((n+x)*(n+x+1)*(n+x+2)),n=1..infinity));
1/(2(x+2)(x+1))
"A mano":
1/((n+x)(n+x+1)(n+x+2)) = ( 1/((n+x)(n+x+1)) − 1/((n+x+1)(n+x+2)) ) / 2
Quindi si tratta di una serie telescopica, che converge ad a1−L dove an
è
(Nota: ogni serie ∑bi può essere messa in forma telescopica [basta considerare
la serie an+1 = a1−(b1+...+bn), a1 qualunque], ma solo in alcuni casi la cosa
è utile)
Con Maple si può controllare che si ottiene una somma per ogni x intero non negativo, mentre
non si ottiene risultato per x negativo (infatti i temini della somma sono in alcuni casi non definiti). La dimostrazione
precedente vale per ogni x intero non negativo.
(2)
Confronta la convergenza delle serie:
1+1/2+1/4+1/8+1/16+... e 1+1−1/2+1−3/4+1−7/8+1−15/16+...
Attenzione, per quanto 1+1/2 = 1+1−1/2, 1+1/2+1/4 = 1+1−1/2+1−3/4,
le due serie non sono equivalenti in quanto la seconda oscilla tra
1 e 1+1, tra 1+1−1/2 e 1+1−1/2+1, tra 1+1−1/2+1−3/4 e 1+1−1/2+1−3/4+1,
: le somme di posto dispari tendono a 2, quelle di posto pari a 3.
Questo esempio illustra una importante differenza tra le somme finite e quelle infinite:
1+1−1/2+1−3/4+1−7/8+1−15/16+... non equivale a
1+(1−1/2)+(1−3/4)+(1−7/8)+(1−15/16)+...
(3)
Trova con Maple
∑ n=0..∞ xn,
∑ n=0..∞ x2n,
∑ n=0..∞ x2n+1,
∑ n=0..∞ (−x)n,
•
È noto che ∑ n=0..∞ xn = 1/(1−x) converge solo per |x| < 1 (per x=1/2 è la
prima serie dell'es. precedente)
•
∑ n=0..∞ x2n è costituita dai termini della precedente serie
con esponente pari, ma posso ottenerla anche rimpiazando x con x2.
Quindi converge a 1/(1−x2) per |x| < 1.
•
∑ n=0..∞ x2n+1 la ottengo dalla precedente moltpilcandone i termini per x,
quindi converge a x/(1−x2) per |x| < 1.
•
∑ n=0..∞ (−x)n la ottengo da ∑ n=0..∞ xn
rimpiazzando x con −x:
•
∑ n=0..∞ (−x2)n
la ottengo dalla precedente rimpiazzando x con x2: 1/(1+x2) per |x| < 1.
•
∑ n=0..∞ (−1)nx2n+1
la ottengo dalla precedente moltiplicando per x: x/(1+x2) per |x| < 1
(nota:
•
∑ n=0..∞ (4x2)n
la ottengo da ∑ n=0..∞ x2n rimpiazzando x con 2x:
converge a
Tutte le somme (non dove esse valgono) possono essere trovate immediatamente con Maple.
(4) Differenziando e integrando la serie geometrica (e l'espressione che ne esprime il valore)
si ha:
∑ n=0..∞ nxn−1 =
∑ n=0..∞ (−1)nx2n+1/(n+1) =
Per quali x valgono queste eguaglianze?
Si tratta di serie di potenze. Quindi:
∑ n=0..∞ nxn−1 = 1/(1−x)2 se |x| < 1; ho
derivato, rispetto ad x,
∑ n=0..∞ (−1)nx2n+1/(n+1) = log(1+x); vale se |x| < 1,
ma anche se x = 1 −
(5)
Integrando
(attenzione: Maple potrebbe darti una somma "strana").
Per quali x vale questa eguaglianza?
x − x3/3 + x5/5 − x7/7 +
= arctan(x).
La serie sicuramente converge per |x| < 1, dove converge la serie di partenza.
La serie converge anche per |x| = 1:
sum((-1)^n/(2*n+1),n=0..infinity); e
sum(-(-1)^n/(2*n+1),n=0..infinity); convergono a π/4 e −π/4
(le arcotangenti di 1 e −1).
Per |x| > 1 la serie diverge (la somma oscilla tra valori che tendono a ∞ e a −∞).
Verifica: taylor(arctan(x),x,8); dà
Secondo foglio
(1) Usando la formula di MacLaurin del 3° ordine, approssima e0.1+sin(0.1) e maggiora
il valore assoluto dell'errore commesso, e precisa il segno dell'errore.
exp(x) = 1 + x + x2/2 + x3/3! + exp(c1) x4/4!
con c1 tra 0 ed x
sin(x) = x − x3/3! + cos(c2) x5/5!
con c2 tra 0 ed x
exp(x) + sin(x) = 1 + 2x + x2/2 + exp(c1) x4/4! + cos(c2) x5/5!
con c1 e c2 tra 0 ed x
exp(0.1) + sin(0.1) ≈ 1 + 0.2 + 0.12/2 = 1.205
Sia E l'errore. E = exp(c1) 0.14/4! + cos(c2) 0.15/5!
con c1 e c2 tra 0 ed 0.1. Quindi E > 0, ossia 1.205 è una approssimazione per difetto del valore cercato.
E < 3/(104·24) + 1/(105·120) = 1.258333
·10−5
< 0.000013
(2) Determina l'insieme degli x per cui la serie data converge e calcolane la somma:
(a) ∑ n=0..∞ (−1)n x2n
|x| < 1; 1/(1+x2) −
da ∑ n=0..∞ (−1)n xn
rimpiazzando x con x2
(b) ∑ n=0..∞ xn / 3n+1
|x| < 3; 1/(3−x) −
xn/3n+1 = (x/3)n/3
(c) ∑ n=0..∞ n xn
|x| < 1; x/(1−x)2 −
da ∑ n=0..∞ xn = 1/(1−x), derivando e poi moltiplicando per x
(d) ∑ n=0..∞ (−1)n n xn
|x| < 1; −x/(1+x)2 −
dalla precedente rimpiazzando x con −x
(e) ∑ n=0..∞ (−2)n (n+2)/(n+1)·xn
|x| < 1/2; 1/(1+2x) + log(1+2x)/(2x) [vedi, sotto, i commenti]
(f) ∑ n=1..∞ 2n xn / n
−1/2 ≤ x < 1/2; −log(1−2x) −
dallo sviluppo di log(1+x) ricavo quello di log(1-x); usando questo, e riscrivendo
2nxn come (2x)n
(g) ∑ n=0..∞ (−1)n / (2n+1)·(x/2)2n
−2 ≤ x ≤ 2; 2/x·arctan(x/2)
Con Maple si ottiene un'espressione equivalente alla soluzione se si opera in C. Per trovare la
soluzione direttamente occorre svolgere, direttamente, il calcolo:
ricordare che (vedi)
∑ n=0..∞ (−1)n / (2n+1)·x2n+1 = arctan(x),
dedurne che
(h) ∑ n=0..∞ (−1)n x3n/n!
ogni x; exp(−x3)
(i) ∑ n=0..∞ xn / (n+3)!
ogni x; x−3(exp(x)−1−x−x2/2) se x ≠ 0, 0 se x = 0
(usa
simplify(expand(sum(x^n/(n+3)!,n=0..infinity)));
se vuoi riscrivere l'espressione che ti fornisce Maple;
in ogni caso basta consoscere lo sviluppo di exp(x), e osservare che la serie data si ottiene da questo
saltando i primi tre termini e dividendo per x3)
(j) ∑ n=0..∞ (x−1)n / (n+2)!
ogni x; (exp(x−1)−x)/(x−1)2 se x ≠ 1, 1/2 se x = 1
(usa simplify(expand(...)); se vuoi riscrivere l'espressione che ti fornisce Maple)
[traccia: ∑ n=0..∞ xn/n!
→ ∑ n=0..∞ xn/(n+2)!
→ ∑ n=0..∞ (x−1)n/(n+2)!]
(k) ∑ n=0..∞ xn / 2n+1
|x| < 2; 1/(2−x)
(l) ∑ n=0..∞ x2n+1 / (2n+1)!
ogni x; sinh(x) (vedi)
[per l'es. (e) con Maple si ottiene una soluzione strana:
che consente di capire che la serie converge per |x|<1/2, e che è incerto che cosa accade per x = 1/2;
procedendo "a mano", comunque, si ottiene:
∑ n=0..∞ (−2)n (n+2)/(n+1)·xn =
∑ n=0..∞ (−1)n (n+2)/(n+1)·(2x)n
che converge dove converge
(3)
Provare a studiare con Maple le serie considerate negli esempi del
paragrafo
"serie di funzioni".
Es. 1. Sum((-1)^n*x^(3*n)/((n+1)*4^n),n=0..infinity);
Per controllare la risposta posso fare:
Con "simplify" (e la "sum" con la "s" minuscola) ottengo anche la somma
Es. 2. Sum(x^(2*n+1)*(-1)^n/(2*n+1),n=0..infinity);
controllo:
Con
Es. 3. Sum(x^n*(n+2)/n!,n=0..infinity);
Con la "s" minuscola ottengo proprio
Es. 4. In questo caso ci limitiamo
a ritrovare lo stesso valore;
Es. 5.
taylor(tan(x),x,8); fornisce
Es. 6. Per sviluppare (con Maple) in serie di potenze di x-1 mi conviene fare una sostituzione:
con subs(x=u+1, 1/(2+x)); ottengo 1/(u+3)
con taylor(1/(u+3),u); ottengo 1/3-1/9u+1/27u^2-1/81u^3+1/243u^4-1/729u^5+O(u^6)
con subs(u=x-1, convert(taylor(1/(u+3),u),polynom));
(o subs(u=x-1, 1/3- ... -1/729u^5);)
ottengo
Es. 7. Maple di fronte a sum(-(x^n)/(n+1),n=0..infinity);
fornisce subito
Es. 9. Maple di fronte a sum((n+2)*(n+1)/2*x^n,n=0..infinity);
fornisce