Introduzione all'uso di Mathematica

Marco Bramanti
Questi appunti sono scritti per chi vuole imparare velocemente a usare Mathematica per fare le pi comuni operazioni che hanno a che fare con il contenuto dei corsi di Analisi Matematica 1 e 2. Scorrendo il Sommario qui sotto saprete subito cosa potrete imparare e cosa no, leggendo questi appunti. Facendo le operazioni che sono spiegate qui, vi capiter di chiedervi come potreste farne molte altre, che qui non sono spiegate. Le risposte a queste domande si trovano nei manuali ufficiali di Mathematica, che sono i vari libri di Stephen Wolfram, il creatore del programma: c' un libro per ogni versione di Mathematica (attualmente, ottobre 2000, sono in circolazione le versioni 2, 3 e 4). Ciascuno di questi libri ha dalle 900 pagine in su, il che potrebbe scoraggiare qualcuno a consultarli: in realt sono libri ben scritti, in cui non difficile orientarsi. Comunque, questi appunti dovrebbero servire proprio per evitare di sfogliare i manuali, a chi desideri imparare in fretta l'essenziale.

Sommario
1. Nozioni di base su Mathematica
(Conviene leggere integralmente questa prima parte) 1.1. 1.2. 1.3. 1.4. 1.5. 1.6. Primo impatto Calcoli numerici Funzioni elementari Costanti notevoli Definire una funzione o un numero Matematica discreta 1.6.1. Successioni, tabulazioni 1.6.2. Sommatorie

2. Argomenti specifici di Analisi Matematica

(Sono abbastanza indipendenti gli uni dagli altri, si pu leggere solo quello che interessa; richiedono per le Nozioni di base)
2.1. 2.2. 2.3. 2.4. 2.5. 2.6. 2.7. Serie numeriche Numeri complessi 2.2.1. Radici n-esime in C Risoluzione di equazioni, esatta o approssimata Grafici di funzioni Funzioni definite a tratti Uso dei comandi Sum e Table con funzioni Calcolo infinitesimale per funzioni di una variabile 2.7.1. Calcolo dei limiti 2.7.2. Calcolo delle derivate

2.8.

2.7.3. Formula di Taylor 2.7.4. Calcolo delle primitive 2.7.5. Calcolo degli integrali definiti, esatto o approssimato 2.7.6. Integrali generalizzati 2.7.7. Funzioni integrali 2.7.8. Equazioni differenziali 2.7.9. Serie di Fourier Funzioni di pi variabili (questa parte non ancora stata scritta) 2.8.1. Grafici di funzioni reali di due variabili 2.8.2. Grafici di linee in forma parametrica,nel piano o nello spazio 2.8.3. Grafici di superfici in forma parametrica 2.8.4. Calcolo di integrali doppi,esatto o approssimato

1. Nozioni di base su Mathematica

1.1. Primo impatto
All'avvio di Mathematica compare una finestra principale bianca (Untitled), pi eventualmente altre finestre con raccolte di simboli e comandi. Non ostante l'aspetto grafico sia quello di un comune programma windows (con men ecc.), per fare qualcosa di significativo non basta cliccare sui vari men con il mouse: occorre usare la tastiera e digitare istruzioni secondo una certa sintassi. Quanta sintassi sia necessario imparare dipende quindi da quale e quanta matematica si vuole fare con Mathematica. Pi precisamente, per fare qualsiasi operazione, si digita dalla tastiera una opportuna stringa di testo (input) nella finestra principale, e quindi si preme shift+invio: Mathematica a questo punto interpreta come input quello che si appena scritto (pi precisamente, quello che scritto nel blocco di righe in cui sta lampeggiando il cursore), esegue quanto richiesto e fornisce il risultato (output): un numero, un'espressione, un grafico..., a seconda di cosa si chiesto. L'input va digitato secondo una precisa sintassi: in particolare, sono importanti le maiuscole/minuscole, gli spazi, l'uso corretto delle parentesi graffe, quadre e tonde. Ricordiamo che sulla tastiera italiana: le parentesi quadre [ ] si trovano sulla tastiera, ma sono il terzo carattere di un tasto, perci si ottengono premendo quel tasto insieme ad AltGr; le parentesi graffe {} non ci sono sulla tastiera, ma si ottengono digitando i numeri 123 e 125 col tastierino numerico, mentre si tiene premuto il tasto Alt. Nella versione 2 di Mathematica, gli input si digitano esclusivamente da tastiera. Nelle versioni 3 o 4 invece, ci sono varie barre di strumenti cliccando sui quali si pu immettere almeno parte dell'input, usando un po' meno sintassi. In questi appunti si insegna a usare solo la tastiera; l'uso del mouse e dei toolbar, del resto, un'alternativa intuitiva che non richiede molte spiegazioni. Nota tecnica: alcune delle operazioni che pu capitare di fare con Mathematica sono piuttosto impegnative per il PC; se usate una vecchia macchina (con processore "lento", tipo 486 o poca RAM, tipo 8Mb), potrebbe bloccarsi o metterci parecchio. Nel caso si abbia questo dubbio, meglio salvare prima di premere shift+invio.

1.2. Calcoli numerici

Esempio. Se scriviamo: 2+3*4^2

e quindi premiamo shift+invio, succede questo:

In[1]:= 2 + 3 * 4 ^ 2 Out[1]= 50

Come si vede, Mathematica scrive In[n]= davanti all'input (n il numero progressivo degli input digitati nella sessione), e fornisce un output. (Quando un file Mathematica viene salvato, chiuso e riaperto, le scritte In[n]= e Out[n]= non ricompaiono, finch eventualmente l'input non viene valutato di nuovo; questo il motivo per cui in questo file non tutti gli input e output sono numerati). Regole di base: le 4 operazioni si indicano con
+- *

Il prodotto si pu indicare anche con uno spazio, ad esempio:

2x

lo stesso che
2*x

ma diverso da 2x (senza spazio tra 2 e x), che denoterebbe una nuova variabile di nome "2x", e non 2 volte la variabile x. L'elevamento a potenza "a alla b" si scrive cos:
a^ b

Si possono usare tutte le parentesi che si vuole, sempre tonde, le une dentro le altre. Anzi, nel dubbio su quale sia la precedenza con cui Mathematica interpreta le operazioni, meglio abbondare, ad esempio:
In[2]:= H2 ^ H3 + xLL H5 x ^ 2 + 1L Out[2]=

23+x 1 + 5 x2

Le parentesi quadre e graffe hanno un significato sintattico particolare che si vedr in seguito, e non vanno usate quindi come parentesi "normali". Esempio.
In[3]:= H2 45L ^ 36 Out[3]=

68719476736 327831360408733226531187361612285036244429647922515869140625

Se il risultato esatto di un calcolo un numero razionale, Mathematica lo scrive come frazione esatta! Se vogliamo un valore numerico approssimato, dobbiamo scrivere:
N@operazione che voglio fareD

Ad esempio:
In[4]:= N@H2 45L ^ 36D Out[4]= 2.09618 10-49

Trucco. Col simbolo % si pu indicare, in un input, l'ultimo output precedente. Ad esempio:

In[5]:= % * 1000 Out[5]= 2.09618 10-46

Ad esempio, per ottenere un valore numerico approssimato dell'ultimo output, basta scrivere N[%]. Invece, per indicare l'output numero n si scrive %n. Ad esempio,
In[6]:= 1 %3 Out[6]=

327831360408733226531187361612285036244429647922515869140625 68719476736

1.3. Funzioni elementari:

In generale una funzione elementare si scrive F[x], cio ha un nome che comincia sempre per maiuscola, e il suo argomento scritto tra parentesi quadre. Funzioni comuni:
Exp@xD

(esponenziale di base e)
Log@xD

(logaritmo in base e)
Log@a, xD

(logaritmo in base a)
Sin@xD, Cos@xD, Tan@xD, ArcSin@xD, ArcCos@xD, ArcTan@xD

(funzioni trigonometriche e loro inverse)

Sinh@xD, Cosh@xD, Tanh@xD

(funzioni iperboliche)
Sqrt@xD

(radice quadrata di x; ricordare che si pu comunque scrivere x^a per avere potenze a esponente reale qualunque; Mathematica, per, richiede sempre che la base sia positiva. Ad esempio, se si chiede a Mathematica di disegnare il grafico di x^(1/3), lo disegna solo per x>0. Per averlo tutto occorre usare un trucco che si spiegher poi).
Abs@xD

(valore assoluto di x)

1.4. Costanti notevoli:

pi greco il numero e di Nepero l'unit immaginaria pi infinito meno infinito si scrive si scrive si scrive si scrive si scrive Pi E I Infinity -Infinity

1.5. Definire una funzione o un numero

Se si vogliono fare pi operazioni su una funzione non elementare, conviene darle un nome, in modo da poterla poi richiamare con quel nome. Ad esempio:
In[7]:= F@x_D = HSin@2 xD + 2 Cos@3 xDL H1 + x ^ 2L Out[7]=

2 Cos@3 xD + Sin@2 xD 1 + x2

Se ora scriviamo, ad esempio:

F@2 Pi 3D
3 2 - 2 4 p2 1 + 9 !!!!

Mathematica calcola il valore di F in x=2 Pi/3. La sintassi dunque :

Nome@x_ D = Qualcosa che contiene@xD

Il nome della funzione pu essere una lettera (maiuscola) o una parola (con l'iniziale maiuscola). Alcune lettere sono proibite perch hanno gi un altro significato. Ad esempio, si possono usare:
F@x_D, G@x_ D, H@x_ D, M @x_ D

ma non si possono usare:

N@x_ D, O@x_ D, D@x_ D.

(Motivo: N indica la valutazione numerica, O il simbolo di "o grande", D la derivazione). Se si prova a farlo, Mathematica ci avverte che non lecito. Naturalemente il nome della variabile indipendente arbitrario. Ad esempio potremmo digitare:
In[8]:= F@t_D = HSin@2 tD + 2 Cos@3 tDL H1 + t ^ 2L Out[8]=

2 Cos@3 tD + Sin@2 tD 1 + t2

e sarebbe stata la stessa funzione dell'esempio precedente (In[7]) (t o x sono "variabili mute"). Analogamente, ma pi semplicemente, si pu dare un nome a un numero che in seguito comparir in parecchie operazioni. Si pone ad esempio:

6
a = 2 Pi Sqrt@3D 2 !!!! 3 p

A questo punto possiamo ad esempio calcolare:

Exp@aD 2
!!!! 3 p

N@%D 53252.3

Avvertenza. Quando si definisce una funzione o un numero dandogli dei nomi, per tutta la sessione di lavoro Mathematica dar a quel nome quel significato. Dopo un po' questo potrebbe diventare fastidioso (non si sa pi che lettere usare perch sono gi tutte "compromesse"). Per dire a Mathematica "ora dimentica che x vuol dire la tal cosa" il comando :
Clear@xD

Lo stesso per i nomi di funzione: con

Clear@FD

Mathematica fa s che F non abbia pi il significato di funzione che gli era stato dato.

1.6. Matematica discreta

1.6.1. Successioni, tabulazioni
Definire una successione esattamente analogo a definire una funzione:
A @k_ D = 1 Hk ^ 2L 1 k2

Fin qui, anzi, Mathematica non sa se noi intendiamo che k sia intero o reale. Potremmo chiedergli di calcolare:
A @PiD 1 p2

Fattoriale e coefficiente binomiale si indicano, naturalmente, con n! e con Binomial[n,k]

Binomial@10, 4D 210 10 ! 3628800

Se vogliamo tabulare i primi (o un po' di) valori della successione, il comando :

Table@ A @kD, 8k, n1, n2<D

che produce i valori di A[k] quando k intero tra n1 e n2. Questo comando, quindi, forza Mathematica a intepretare k come intero.
Table@ A @kD, 8k, 1, 10<D 1 1 1 1 1 1 1 1 1 91, , , , , , , , , = 4 9 16 25 36 49 64 81 100

Al solito, se volessimo valori numerici in forma decimale, dovremmo chiedere:

Table@ N@ A @kDD, 8k, 1, 10<D 81., 0.25, 0.111111, 0.0625, 0.04, 0.0277778, 0.0204082, 0.015625, 0.0123457, 0.01<

Ossia: "tabula i valori numerici approssimati di A[k] per k da 1 a 10" Se vogliamo avere un'idea dell'andamento di A[k] per k grande, pu essere utile tabulare i valori di A[k] fino a un valore grande, ma incrementando k ogni volta di un passo maggiore di 1 (per non visualizzare troppi numeri). Ad esempio:
Table@ N@ A @kDD, 8k, 1, 100, 10<D 81., 0.00826446, 0.00226757, 0.00104058, 0.000594884, 0.000384468, 0.000268745, 0.000198373, 0.000152416, 0.000120758<

produce i valori di A[k] per k che si incrementa di 10 alla volta, da 1 a 100. Per rappresentare graficamente i valori di (un po' di termini di) una successione, il comando :
ListPlot@Table@ A @kD, 8k, n1, n2<DD

che produce il grafico dei punti di coordinate (k,A[k]) per k da n1 a n2:

ListPlot@Table@ A @kD, 8k, 5, 10<DD
0.04 0.035 0.03 0.025 0.02 0.015 2 3 4 5 6

Graphics

Se volete punti pi grossi, fate cos (e non chiedetemi perch):

8
ListPlot@Table@ A @kD, 8k, 5, 10<D, Prolog -> AbsolutePointSize@5DD
0.04 0.035 0.03 0.025 0.02 0.015 2 3 4 5 6

Graphics

Pi in generale, ListPlot[{lista di numeri}], produce un grafico dei punti di ascisse 1,2,3,... e ordinate questi numeri
ListPlot@83, 4, 1, 9<D

1.5

2.5

3.5

Graphics

1.6.2. Sommatorie
L'operazione di sommatoria per k da n1 a n2 di A[k] si scrive naturalmente cos:
Sum @A @kD, 8k, n1, n2<D Sum @1 k, 8k, 1, 20<D 55835135 15519504

N@%D

3.59774

L'ultimo numero il valore numerico approssimato del precedente. Come il comando Table, cos il comando Sum forza Mathematica a interpretare la variabile k (o come l'abbiamo chiamata) come intero. Con ci abbiamo introdotto la sintassi principale di Mathematica, che si utilizza in qualsiasi contesto. Ora passiamo agli argomenti specifici di calcolo infinitesimale.

2. Argomenti specifici di Analisi Matematica

2.1. Serie numeriche
Si pu definire una serie con lo stesso operatore Sum di sommatoria, facendo variare l'indice fino a infinito; se ne capace, Mathematica fornisce la somma corretta della serie...
Sum @1 Hk ^ 2L, 8k, 1, Infinity<D p2 6

Sum @1 k, 8k, 1, Infinity<D

Sum::div : Sum does not converge.

k=1

1 k

...oppure dice che la serie non converge!

Sum @1 Hk ^ 2 + 2 ^ kL, 8k, 1, Infinity<D 1 2 + 2k k k=1

...oppure ricopia l'espressione perch non ne sa calcolare la somma, ma non ha neanche la prova che la serie diverga. Per indagare numericamente il comportamento di una serie pu essere utile definire la successione delle somme parziali, e tabularne e/o rappresentarne graficamente i valori. Consideriamo ad esempio la serie il cui termine generale :
In[9]:= B@k_D = H2 ^ kL Hk + k !L Out[9]=

In[10]:= S@n_ D = Sum @B@kD, 8k, 1, n<D Out[10]= B@kD

n k=1

2k k + k!

In[11]:= Table@ N@S@kDD, 8k, 1, 10<D

Out[11]= 81., 2., 2.88889, 3.46032, 3.71632, 3.80447, 3.82983, 3.83618, 3.83759, 3.83787<

...dopo i primi 10 termini appaiono gi stabilizzate le prime 3 cifre decimali della somma...

10
In[12]:= ListPlot@%11D

3.5 3 2.5 2 1.5

4 Out[12]= Graphics

10

2.2. Numeri complessi

I numeri complessi si scrivono a+I b. Si possono eseguire le operazioni algebriche:
H3 + 2 IL H2 - IL ^ 3 28 - 29

Inoltre, si possono calcolare parte reale, parte immaginaria, coniugato, modulo, argomento di un numero. I comandi corrsipondenti sono, rispettivamente:
Re@zD, Im @zD, Conjugate@zD, Abs@zD, Arg @zD

Esempi:
Re@H1 + 3 IL ^ 5D 316 Arg @1 + ID p 4

Arg @1 - 2 ID - ArcTan@2D Abs@1 + ID !!!! 2 Conjugate@H3 + 2 IL ^ 2D 5 - 12

Esempio di calcolo della somma di una serie nel campo complesso (v. "Serie numeriche"):
Sum @1 H1 + IL ^ k, 8k, 0, Infinity<D 1-

11

2.2.1. Radici ennesime in C.

Se si calcola la radice quadrata o ennesima di un numero complesso direttamente, Mathematica ne calcola una sola (anzich n). Questo perch il simbolo Sqrt[ ] oppure ( )^n denota per Mathematica una funzione (univoca). Esempi:
Sqrt@2 ID 1+ I ^ H1 3L H- 1L16

Il risultato non soddisfacente. Per trovare davvero le radici ennesime di z, bisogna chiedere esplicitamente a Mathematica di applicare la formula di De Moivre, cio dare una sequenza di istruzioni del tipo:
x=1+ I 1+ n=6 6 r = Abs@xD !!!! 2 t = Arg @xD p 4

z@k_D = Hr ^ H1 nLL HCos@Ht + 2 k PiL nD + I Sin@Ht + 2 k PiL nDL 1 p 1 p 2112 JCosA I + 2 k p ME + SinA I + 2 k p MEN 6 4 6 4 Table@ N@z@kDD, 8k, 0, n - 1<D

81.0504 + 0.138288 , 0.405439 + 0.978816 , - 0.64496 + 0.840529 , - 1.0504 - 0.138288 , - 0.405439 - 0.978816 , 0.64496 - 0.840529 <

Plottare gli n punti nel campo complesso richiede invece un'istruzione un po'... complessa, si rimanda al Notebook specifico su questo.

2.3. Risoluzione di equazioni, esatta o approssimata

La risoluzione di un'equazione pu essere esatta o approssimata. 1. Risolvere un'equazione esattamente significa, con passaggi algebrici opportuni, trovare il valore esatto di tutte le soluzioni (come quando risolviamo un'equazione di secondo grado applicando la formula relativa). 2. Risolvere un'equazione con approssimazione significa determinare, coi metodi dell'Analisi Numerica (ossia con qualche tipo di algoritmo iterativo) una soluzione approssimata (o tutte le soluzioni, approssimate) dell'equazione stessa.

12

Siamo noi che dobbiamo dire a Mathematica quale delle 2 cose fare: se gli chiediamo di risolvere esattamente un'equazione quando ci non possibile, non far nulla (in particolare, non ci dar soluzioni approssimate). Se chiediamo soluzioni approssimate, queste saranno "solo" approssimate, e non esatte, anche qualora fosse possibile trovarle esatte. Il comando per trovare le soluzioni esatte Solve; il comando per trovare le soluzioni approssimate FindRoot. Sintassi di Solve:
Solve@ A @xD B@xD, xD

Esempio:
Solve@x ^ 2 + 3 x + 1 0, xD 1 1 !!!! !!!! 99x I- 3 - 5 M=, 9x I- 3 + 5 M== 2 2

Con questo comando Mathematica trova le radici reali o complesse dell'equazione. Esempi:
Solve@x ^ 3 + x ^ 2 + 1 0, xD 1 1 i 2 y z 99x - - j j !z 3 3 k 29 - 3 !!!!!! 93 {
13

!!!! 1 1 1 1- 3 !!!! !!!!!! ! 13 9x - + I1 + 3 M J I29 - 3 93 MN + =, !!!!!! ! 13 3 6 2 3 223 I29 - 3 93 M !!!! 1 1 1 1+ 3 !!!! !!!!!! ! 13 9x - + I1 - 3 M J I29 - 3 93 MN + == !!!!!! ! 13 3 6 2 3 223 I29 - 3 93 M

1 1 !!!!!! ! 13 - J I29 - 3 93 MN =, 3 2

L'equazione di terzo grado precedente ha una radice reale e due radici complesse coniugate.
Solve@Exp@2 xD - 3 Exp@xD + 1 0, xD
Solve::ifun : Inverse functions are being used by Solve, so some solutions may not be found.

1 1 !!!! !!!! 99x LogA I3 - 5 ME=, 9x LogA I3 + 5 ME== 2 2

L'equazione precedente un'equazione di secondo grado in Exp[x], perci pu essere risolta esattamente. Si possono risolvere anche sistemi, lineari e (talvolta!) non lineari:
Solve@8x + 2 y == 1, 3 x - y == 0<, 8x, y<D 1 3 99x , y == 7 7

Solve@8x H y + 1L 0, 3 x - y == 0<, 8x, y<D 1 99x - , y - 1=, 8x 0, y 0<= 3

Sintassi di FindRoot:
FindRoot@ A @xD 0, 8x, x0<D

13

Questo comando trova una soluzione approssimata dell'equazione in un intorno del punto x0 (pi precisamente, x0 viene assunto come primo valore dell'algoritmo iterativo che definisce la soluzione approssimata) Ad esempio, l'equazione Tan[x]x ha infinite soluzioni (una zero). Possiamo trovarne qualcuna dando valori diversi a x0:
FindRoot@Tan@xD x, 8x, 2<D 8x 4.49341< 8x 7.72525< 8x 10.9041< FindRoot@Tan@xD x, 8x, 3<D FindRoot@Tan@xD x, 8x, 8<D

Per cercare le soluzioni con criterio, anzich andando a casaccio, pu essere utile prima tracciare con Mathematica il grafico di F[x]=Tan[x]-x, e vedere approssimativamente dove cadono i suoi zeri.Vedremo nel prossimo paragrafo come si fa.

2.4. Grafici di funzioni

La sintassi di base per tracciare il grafico di una funzione di una variabile :
Plot@F@xD, 8x, a, b<D

Cos facendo, Mathematica traccia il grafico di F[x] per x nell'intervallo [a,b]. Se in quell'intervallo la funzione non sempre definita, Mathematica protesta, ma fa il possibile.
Plot@ Abs@Log @Abs@xDDD, 8x, - 3, 3<D
6 5 4 3 2 1 -3 -2 -1 1 2 3

Graphics

Notare che noi abbiamo detto ben poco a Mathematica su come tracciare il grafico, e lui ha fatto una serie di scelte (quale scala mettere su ciascun asse, ad esempio). Talvolta non saremo d'accordo con le sue scelte, e vorremo chiedere qualcosa di pi preciso. Questo si fa con vari tipi di opzione. Esempi:

14
Plot@1 + Sin@3 xD ^ 2, 8x, - Pi, Pi<D
2 1.8 1.6 1.4 1.2

-3

-2

-1

Graphics

Mathematica ha tracciato il grafico scegliendo come asse orizzontale la retta y=1 (guardare i numeri scritti sull'asse!). Questo pu essere fuorviante. Se vogliamo che il grafico mostri le y tra 0 e 2, ad esempio (in particolare, che l'asse orizzontale sia la retta y=0), possiamo usare l'opzione PlotRange->{a,b} dove a,b sono i valori estremi delle y che vogliamo compaiano nel grafico:
Plot@1 + Sin@3 xD ^ 2, 8x, - Pi, Pi<, PlotRange 80, 2<D
2 1.75 1.5 1.25 1 0.75 0.5 0.25 -3 -2 -1 1 2 3

Graphics

Altro problema:
Plot@x ^ 2, 8x, - 3, 3<D

8 6 4 2

-3

-2

-1

Graphics

15

Nel grafico precedente la scala non la stessa sugli assi. L'opzione AspectRatioa (a numero positivo, maggiore o minore di 1) permette di dilatare un asse rispetto all'altro. Per tentativi, si pu trovare la scala che si vuole:
Plot@x ^ 2, 8x, - 3, 3<, AspectRatio 1D

-3

-2

-1

Graphics Plot@x ^ 2, 8x, - 3, 3<, AspectRatio 0.5D

8 6 4 2

-3

-2

-1

Graphics

16
Plot@x ^ 2, 8x, - 3, 3<, AspectRatio 2D

-3

-2

-1

Graphics

Vi sono molte altre opzioni possibili per il comando Plot, ma non ve le spiego per non riscrivere qui il manuale di Mathematica... Si osservi il prossimo trucco per tracciare il grafico di x^(2/3) anche per x<0:
Plot@Hx ^ 2L ^ H1 3L, 8x, - 3, 3<D
2

1.5

0.5

-3

-2

-1

Graphics

2.5. Funzioni definite a tratti

Si possono ottenere con un comando del tipo "se allora". Ad esempio la funzione che vale x se x<0 e x^2 se x>0 si pu ottenere cos:
In[13]:= If@x < 0, x, x ^ 2D Out[13]= If@x < 0, x, x2 D

La sintassi : If[Condizione, Allora cos, Altrimenti cos]

17
In[14]:= Plot@If@x < 0, x, x ^ 2D, 8x, - 1, 1<D 1

0.5

-1

-0.5 -0.5

0.5

-1 Out[14]= Graphics

Si possono annidare tanti If[...] uno dentro l'altro, e ottenere funzioni definite in 3 o pi pezzi diversi, ad esempio:
In[15]:= Plot@If@x < 0, x, If@x < 1, x ^ 2, 2 - xDD, 8x, - 1, 2<D 1

0.5

-1

-0.5 -0.5

0.5

1.5

-1 Out[15]= Graphics

L'ultimo input si pu leggere cos: "Disegna, per x in [-1,2], la funzione che se x<0 vale x, altrimenti (se x<1 vale x^2, altrimenti vale 2-x)". Ad esempio, per ottenere il grafico della radice cubica di x anche per x<0 si pu fare cos:
In[16]:= Plot@If@x > 0, x ^ H1 3L, - H- xL ^ H1 3LD, 8x, - 2, 2<D

1 0.5

-2

-1 -0.5 -1

Out[16]= Graphics

Si possono tracciare i grafici di pi funzioni su un solo grafico. Questo si pu fare in due modi:

18
Plot@8F@xD, G@xD, H@xD<, 8x, a, b<D

disegna sullo stesso grafico le funzioni F, G, H. Se invece abbiamo gi prodotto separatamente i grafici di F, G e H, e ci piacerebbe ora vederli tutti insieme, usiamo il comando Show[%n,%m,%r], dove n, m, r sono i numeri progressivi degli input in cui abbiamo fatto i grafici di F, G, H. Ad esempio:
In[17]:= Show@%15, %16D

1 0.5

-2

-1 -0.5 -1

Out[17]= Graphics

Si noti che, col comando Show, Mathematica non esegue nuovi calcoli, ma si limita a mostrare sullo stesso sistema d'assi due o pi grafici gi ottenuti; se le due funzioni erano state disegnate su intervalli diversi, col comando show vedremo solo i tratti di grafico che erano gi stati disegnati (nell'ultimo esempio, si noti che la retta x "non esiste" per x<-1, perch cos era nel grafico originale). Per ovviare a questo problema, dovrei usare il primo metodo, col quale Mathematica ricalcola i due grafici:
In[18]:= Plot@8If@x < 0, x, If@x < 1, x ^ 2, 2 - xDD,

If@x > 0, x ^ H1 3L, - H- xL ^ H1 3LD<, 8x, - 2, 2<D

1 0.5

-2

-1 -0.5 -1 -1.5 -2

Out[18]= Graphics

2.6. Uso dei comandi Sum e Table con funzioni

Col comando Sum si possono definire anche funzioni che sono somme parziali di serie di funzioni, o serie di funzioni:
F@x_, n_ D = Sum @Hx ^ kL k !, 8k, 0, n<D x H1 + nL Gamma@1 + n, xD Gamma@2 + nD

19

Quella ottenuta cos una funzione di 2 variabili: x e n. Fissando un valore di n, una funzione di x:
Plot@8Exp@xD, F@x, 4D<, 8x, - 3, 3<D
10 8 6 4 2

-3

-2

-1

Graphics

Il grafico precedente mostra Exp[x] insieme al suo "polinomio di Taylor di grado 4": vicino a x=0 c' buona approssimazione tra i due. Si possono produrre con una sola istruzione gruppi di grafici col comando Table. Ad esempio, vogliamo confrontare Exp[x] e il suo polinomio di Taylor F[x,n] sull'intervallo [-3,3] per valori crescenti di n:
Table@Plot@8Exp@xD, F@x, nD<, 8x, - 3, 3<D, 8n, 1, 5<D

8 6 4 2

-3

-2

-1 -2 8

-3

-2

-1

20

20

15

10

-3

-2

-1 10 8 6 4 2

-3

-2

-1 20

15

10

-3

8 Graphics , Graphics , Graphics , Graphics , Graphics <

-2

-1

(Un solo input ha prodotto 5 grafici). Al crescere di n l'approssimazione sempre migliore, su un intervallo fissato. Oppure, vogliamo confrontare Exp[x] e il suo polinomio di Taylor F[x,5] sull'intervallo[-n,n] per valori crescenti di n:
Table@Plot@8Exp@xD, F@x, 5D<, 8x, - n, n<D, 8n, 1, 5<D

21

2.5 2 1.5 1 0.5

-1

-0.5 7 6 5 4 3 2 1

0.5

-2

-1 20

15

10

-3

-2

-1

30

20

10

-4

-2

22

40 30 20 10

-4

-2 -10

8 Graphics , Graphics , Graphics , Graphics , Graphics <

Aumentando l'ampiezza dell'intervallo, l'approssimazione sempre peggiore.

2.7. Calcolo infinitesimale per funzioni di una variabile

2.7.1. Calcolo dei limiti
I prossimi esempi mostrano la sintassi:
Limit@Sin@xD x, x 0D 1 Limit@Hx ^ 2 + 2 ^ xL H1 + Log @xD ^ 3L, x InfinityD Limit@Hx ^ 2 + 2 ^ xL H1 + x ^ 3L, x - InfinityD 0

(La freccia si ottiene digitando - seguito da >). Se il limite destro e sinistro sono diversi, Mathematica per default calcola il limite destro (questa una caratteristica non molto bella del programma! Se si trova un risultato non si sa se il limite c' o c' solo da destra):
Limit@Exp@1 xD, x 0D

Per chiedere esplicitamente a Mathematica di calcolare un limite destro o sinistro usiamo l'opzione Direction1 per il limite sinistro, Direction-1 per il limite destro:
Limit@x Abs@xD, x 0, Direction 1D -1 Limit@x Abs@xD, x 0, Direction - 1D 1

23
Limit@Exp@1 xD, x 0, Direction 1D 0

Se si prova a calcolare un limite che non esiste, cosa succede?

Limit@Sin@1 xD, x 0D Interval@8- 1, 1<D

L'output significa che il limite non esiste e la "classe limite" l'intervallo [-1,1].
Limit@x Sin@xD, x InfinityD Interval@8-, <D

L'output significa che il limite non esiste e la "classe limite" l'intervallo (-Infinito,Infinito).

2.7.2. Calcolo delle derivate

La sintassi per il calcolo della derivata prima di una funzione :
D@F@xD, xD

Esempio:
D@Sin@2 xD, xD 2 Cos@2 xD

A questo punto per calcolare la derivata seconda, cio la derivata prima di quanto appena trovato, si pu fare cos:
D@%, xD - 4 Sin@2 xD

Oppure:
D@D@Sin@2 xD, xD, xD - 4 Sin@2 xD

2.7.3. Formula di Taylor

Per scrivere la formula di Taylor di F[x] all'ordine n, centrata in x0, il comando :
Series@F@xD, 8x, x0, n<D

Notare che, non ostante il comando si chiami Series, serve a scrivere la formula di Taylor, non la serie di Taylor. Inoltre, notare che il resto (secondo Peano) scritto in modo insolito: al posto del simbolo di o piccolo viene usato il simbolo di O grande. In pratica, ci significa che se Mathematica scrive O[x]^4 bisogna interpretare o[x]^3, e cos via.

24
Series@Sin@xD, 8x, 0, 9<D

x3 x5 x7 x9 x - + - + + O@xD10 6 120 5040 362880

Qui O[x]^10=o[x]^9
Series@Exp@x ^ 2D - x Sin@xD, 8x, 0, 5<D 2 x4 1 + + O@xD6 3

Se si vuole solo il polinomio di Taylor, anzich la formula di Taylor, il comando

Normal@Series@F@xD, 8x, x0, n<DD Normal@Series@Exp@xD, 8x, 1, 4<DD

1 1 1 + H- 1 + xL + H- 1 + xL2 + H- 1 + xL3 + H- 1 + xL4 2 6 24

2.7.4. Calcolo delle primitive

La sintassi per il calcolo di una primitiva di una funzione :
Integrate@F@xD, xD Integrate@H1 + xL H1 + x ^ 2 + xL, xD
1+2 x ArcTanA !!!! E 1 3 + Log@1 + x + x2 D !!!! 2 3

(Notare che Mathematica non aggiunge la costante arbitraria). Naturalmente se la primitiva non una funzione elementare, Mathematica non trova nulla:
Integrate@1 ArcTan@xD, xD 1 x ArcTan@xD

che come dire che Mathematica non ha fatto nulla. Talvolta la primitiva di una funzione una funzione non elementare che per in Analisi ha un nome ed studiata; se Mathematica la conosce, fornisce questo risultato. Ad esempio:
Integrate@1 Log @xD, xD LogIntegral@xD

1 La funzione LogIntegral@xD definita proprio come primitiva di . Mathematica Log@xD conosce questa funzione ed in grado di calcolarne valori. Ad esempio :
N@LogIntegral@2DD 1.04516

25
Plot@LogIntegral@xD, 8x, 0, 1<D
0.2 -2 -4 -6 -8 -10 0.4 0.6 0.8 1

Graphics

2.7.5. Calcolo degli integrali definiti, esatto ed approssimato

La sintassi per calcolare l'integrale definito di F sull'intervallo [a,b] :
Integrate@F@xD, 8x, a, b<D Integrate@1 H1 + x ^ 2L, 8x, 0, 1<D p 4

Integrate@H1 + xL H1 + x ^ 2 + xL, 8x, 0, 1<D

1 !!!! !!!! I3 Log@3D + 3 LogA3 I1 - 3 ME 6 !!!! !!!! !!!! !!!! !!!! !!!! 3 LogA3 - 3 E - 3 LogA3 I1 + 3 ME + 3 LogA3 + 3 EM

Talvolta, come nell'ultimo caso, l'output infelice: qui ad esempio espresso attraverso numeri complessi (anche se in realt reale). Comunque possiamo sempre chiedergli quanto fa:
N@%D 0.851606 + 0.

cio un numero reale, circa 0.851... Il comando Integrate cerca di calcolare un integrale definito in modo esatto, passando attraverso la primitiva (o con altri metodi esatti): se si prova a calcolare l'integrale definito di una funzione che non ha primitiva conosciuta, Mathematica non calcola nulla:
Integrate@1 H1 + x + Exp@xDL, 8x, 0, 1<D
1 1 x x 0 1+ +x

C' una strada diversa: calcolare un valore approssimato dell'integrale definito, coi metodi dell'analisi numerica. Questo si fa con un comando diverso:
NIntegrate@F@xD, 8x, a, b<D

26
NIntegrate@1 H1 + x + Exp@xDL, 8x, 0, 1<D 0.329971

Si confronti il risultato trovato con Integrate e NIntegrate in un caso in cui entrambi i metodi sono applicabili:
Integrate@1 H1 + x ^ 2L, 8x, 0, 1<D p 4

NIntegrate@1 H1 + x ^ 2L, 8x, 0, 1<D 0.785398

Si osservi che:
N@Pi 4D 0.785398

cio il valore "approssimato" approssimato davvero bene.

2.7.6. Integrali generalizzati

Si possono calcolare integrali definiti generalizzati:
Integrate@1 Sqrt@xD, 8x, 0, 1<D 2

L'integrale precedente converge ed calcolabile esattamente.

Integrate@1 x, 8x, 0, 1<D
1 Integrate::idiv : Integral of does not converge on 80, 1<. x
1 1 x 0 x

L'integrale precedente diverge.

Integrate@Exp@- x ^ 2D, 8x, - Infinity, Infinity<D !!!! p

L'integrale precedente converge ed calcolabile esattamente.

Integrate@1 Sqrt@Sin@xDD, 8x, 0, 1<D 1 1 5 !!!!!!!!!!!!!!!!! 2 Hypergeometric2F1A , , , Sin@1D2 E Sin@1D 4 2 4

Questa una tipica risposta "da Mathematica": ha calcolato il valore esatto dell'integrale, esprimendolo con complicate funzioni. Quanto vale?

27
N@%D 2.03481 Integrate@1 Sqrt@x + Sin@xDD, 8x, 0, 1<D
1 1 x !!!!!!!!!!!!!!!!!!!!!!!!! 0 x + Sin@xD

Questa volta non ce l'ha fatta. Se per gli chiediamo di calcolare l'integrale numericamente:
NIntegrate@1 Sqrt@x + Sin@xDD, 8x, 0, 1<D 1.42608

si capisce che l'integrale converge. Se si provasse a calcolare numericamente un integrale divergente:

NIntegrate@1 Hx + Sin@xDL, 8x, 0, 1<D
NIntegrate::slwcon : Numerical integration converging too slowly; suspect one of the following: singularity, value of the integration being 0, oscillatory integrand, or insufficient WorkingPrecision. If your integrand is oscillatory try using the option Method->Oscillatory in NIntegrate. NIntegrate::ncvb : NIntegrate failed to converge to prescribed accuracy after 7 recursive bisections in x near x = 4.369993747903707`*^-57.

11976.6

Il messaggio di Mathematica mette in guardia sul fatto che la convergenza dell'algoritmo di approssimazione lenta, il che deve far sospettare che l'integrale possa divergere. L'output comunque un valore numerico (in questo caso inattendibile).

2.7.7. Funzioni integrali

Una funzione integrale di F[x] una funzione del tipo:
G@x_D = Integrate@F@tD, 8t, 0, x<D

(Attenzione all'uso corretto della variabile muta t, di integrazione). Proviamo a definire, ad esempio:
G@x_D = Integrate@1 H1 + t ^ 2L, 8t, 0, x<D ArcTan@xD

La primitiva qui calcolabile esattamente e la funzione integrale calcolata esplicitamente.

H@x_D = Integrate@1 HExp@tD + t ^ 2L, 8t, 0, x<D
x 1 t t + t2 0

Ovviamente, non riesce a fare nulla in questo caso. Se proviamo a valutare:

28
H@10D
10 0

1 t t + t2

non troviamo niente. In questi casi (quando cio la primitiva non calcolabile elementarmente, e non nemmeno una funzione speciale, conosciuta da Mathematica), conviene definire la funzione integrale come integrale calcolato numericamente , cio porre:
H@x_D = NIntegrate@1 HExp@tD + t ^ 2L, 8t, 0, x<D
NIntegrate::nlim : t = x is not a valid limit of integration. NIntegrate::nlim : t = x is not a valid limit of integration.

1 NIntegrateA , 8t, 0, x<E t + t2

Mathematica ha protestato, perch gli abbiamo chiesto di calcolare numericamente l'integrale definito tra due estremi di cui il secondo non un numero, ma x! Ovviamente non pu farlo; per ora Mathematica sa chi H[x] e pu calcolarla in un punto, o anche tracciarne il grafico:
H@3D 0.780494 Plot@H@xD, 8x, - 5, 5<D

0.5

-4

-2 -0.5 -1 -1.5

Graphics

Nota: il grafico di una funzione integrale un esempio di operazione che pu risultare "impegnativa" per un vecchio PC, come accennato all'inizio. Nel dubbio, meglio salvare prima di premere shift+invio.

2.7.8. Equazioni differenziali

Mathematica pu risolvere le equazioni differenziali (ed anche i problemi di Cauchy) che noi sapremmo risolvere a mano, con carta e penna. Ma ci mette molto meno. La sintassi con cui si scrive un'equazione differenziale illustrata dall'esempio seguente:
y "@xD+3 y'@xD-2 y@xD==Sin@xD

29

Notare che la funzione incognita va sempre scritta nella forma y[x] (cio indicando la variabile indipendente, che potrebbe anche non chiamarsi x), le derivate si indicano con y'[x], y''[x], y'''[x],..., e tra i due membri dell'equazione ci vogliono DUE segni di =. Per la derivata seconda, usare due segni ' di fila, e NON il segno di virgolette ". La sintassi con cui si chiede a Mathematica di risolvere un'equazione differenziale :
DSolve@equazione, y@xD, xD

Esempi. Equazione del second'ordine, lineare, a coefficienti costanti, anche non omogenea:
DSolve@ y ''@xD + 3 y '@xD - 2 y@xD == Sin@xD, y@xD, xD
2 I-399y@xD 1

!!!!!!! I-3+ 17 M x !!!!!!! 1 2 i 2 I-3+ 17 M - Sin@xD y 4 Cos@xD 2 I-3+ 17 M x j - z j- z !!!!!!! !!!!!!! !!!!!!! !!!!!!! IH-3-2 L+ 17 M IH-3+2 L+ 17 M IH-3-2 L+ 17 M IH-3+2 L+ 17 M { k !!!!!! ! 17
1 I-3+!!!!!!! - 17 M x 2 1 !!!!!!! 1 I3+!!!!!!! !!!!!!! !!!!!!! 17 M x 1 2 I3+ 17 M x Sin@xD y i 2 I3+ 17 M 2 4 Cos@xD 2 I-3- 17 M x j + z j- z !!!!!!! !!!!!!! !!!!!!! !!!!!!! IH 3 2 L + 17 M IH 3 + 2 L + 17 M IH 3 2 L + 17 M IH3+2 L+ 17 M { k == !!!!!! ! 17 1 !!!!!!!

!!!!!!! 17 M x

2 I-3+ C@1D + 1

!!!!!!! 17 M x

C@2D +

Come si vede, Mathematica d l'integrale generale, chiamando C[1] e C[2] le costanti arbitrarie di integrazione. Equazione del prim'ordine, lineare:
DSolve@ y '@xD + x y@xD x, y@xD, xD
2 99y@xD 1 + - C@1D==
x2

Equazione del prim'ordine a variabili separabili:

DSolve@ y@xD y '@xD x ^ 2, y@xD, xD !!!!!!!!!!!!!!!!!!!!!!!!! ! !!!!!!!!!!!!!!!!!!!!!!!!! ! 2 x3 + C@1D 2 x3 + C@1D 99y@xD - =, 9y@xD == !!!! !!!! 3 3

Problema di Cauchy: La sintassi :

DSolve@8equazione, condizioni<, y@xD, xD

dove le condizioni sono scritte con la stessa sintassi dell'equazione:

DSolve@8 y ''@xD + 3 y '@xD - 2 y@xD 0, y@0D 1, y '@0D 2<, y@xD, xD
!!!!!!! 1 1 7 7 i1 y i1 y z z 99y@xD j - j z 2 I-3- 17 M x + j j + z2 !!!!!! ! !!!!!! ! k2 k2 2 17 { 2 17 { !!!!!!! I-3+ 17 M x

==

30

2.7.9. Serie di Fourier

Non esistono comandi appositi per scrivere la serie di Fourie di una funzione: si tratta di combinare le varie tecniche viste fin qui: si definiscono le due successioni dei coefficienti, come opportuni integrali (esatti o approssimati? Dipende!), poi si definisce la successione delle somme parziali della serie di Fourier, quindi si pu plottare questa per certi n... Si veda il notebook apposito.