Studio dell'integrale di h(x) = |Sin(x)/[Graphics:../Images/index_gr_189.gif]|.

L'uso di Mathematica come aiuto in calcoli e ragionamenti non sempre segue regole fisse, come nel caso dello studio di funzioni. Ecco un esempio di come Mathematica può aiutarci a capire le difficoltà nascoste in un calcolo numerico.

Sia h(x)=|Sin(x)/[Graphics:../Images/index_gr_190.gif]|. L'integrale
[Graphics:../Images/index_gr_191.gif]
è difficile da calcolare numericamente. Proviamo a farlo: otterremo una risposta descritta come inaffidabile (con un messaggio che dice, più o meno: "la convergenza dell'integrale è troppo lenta per consentire una stima del attendibile del risultato").

[Graphics:../Images/index_gr_192.gif]
[Graphics:../Images/index_gr_193.gif]
[Graphics:../Images/index_gr_194.gif]
[Graphics:../Images/index_gr_195.gif]
[Graphics:../Images/index_gr_196.gif]
[Graphics:../Images/index_gr_197.gif]

Una risposta più vicina alla realtà ci viene dall'integrale troncato sull'intervallo [π,100]. Quest'ultimo può essere calcolato con grande precisione.

[Graphics:../Images/index_gr_198.gif]
[Graphics:../Images/index_gr_199.gif]

Dobbiamo controllare che la parte non calcolata dell'integrale (su [100,∞]) non modifichi di molto il risultato finale. A tal fine, studiamo l'andamento di h(x).

[Graphics:../Images/index_gr_200.gif]
[Graphics:../Images/index_gr_201.gif]

[Graphics:../Images/index_gr_202.gif]

[Graphics:../Images/index_gr_203.gif]

Il disegno ci mostra una funzione oscillante, di cui il metodo di integrazione fatica a comprendere l'andamento. A seconda di dove cadono i punti scelti come campione per studiare h[x], il valore di h[x] appare più grande o più piccolo. Così l'algoritmo di integrazione  non riesce ad ottenere una stima attendibile di
[Graphics:../Images/index_gr_204.gif]

La funzione h(x) è compresa tra 0 ed (1/[Graphics:../Images/index_gr_205.gif]), dato che 0<=|Sin(x)|<=1. Disegnamo simultaneamente h(x) e (1/[Graphics:../Images/index_gr_206.gif]), per confronto.

[Graphics:../Images/index_gr_207.gif]

[Graphics:../Images/index_gr_208.gif]

[Graphics:../Images/index_gr_209.gif]

Valutando l'integrale di (1/[Graphics:../Images/index_gr_210.gif] una limitazione superiore all'integrale di h(x) su tale intervallo

[Graphics:../Images/index_gr_211.gif]
[Graphics:../Images/index_gr_212.gif]

Ne segue che
[Graphics:../Images/index_gr_213.gif]
è compreso tra:
[Graphics:../Images/index_gr_214.gif]
e
[Graphics:../Images/index_gr_215.gif] + 0.01


Converted by Mathematica      May 26, 2003