Integrazione numerica – Monte Carlo

Esercizio del libro di testo

Data la funzione \displaystyle 3x^2+2x+1 calcolare l’integrale definito della funzione da a=1 a b=5.

Osserva: nell’intervallo [1, 5] la funzione è continua e positiva.
Siccome nell’intervallo è crescente assume il valore massimo nell’estremo a destra, per x=5, MAX=f(5)=86.

Il grafico della funzione è compreso nel rettangolo individuato dai punti

  • A = (a, 0) = (1, 0)
  • B = (b, 0) = (5, 0)
  • C = (b, f(b)) = (5, 86)

e l’area del rettangolo è

  • AREA = (b-a)*f(b) = … = 344

Siano

  • x* ∈ [a, b], ascisse casuali
  • y* ∈ [0, MAX], ordinate casuali
  • PUNTI = numero di punti casuali generati
  • punti = numero di punti casuali con y* <= f(x*)

allora per PUNTI abbastanza grande

\displaystyle \frac{\text{punti}}{\text{PUNTI}}\approx \frac{\text{area}}{\text{AREA}}

e quindi

\displaystyle \text{area}\approx \frac{\text{punti}}{\text{PUNTI}}\cdot \text{AREA}

Valore esatto

Il valore esatto, ricavato con il calcolo integrale, è

area = \displaystyle \int_a^b f(x)\ dx = \displaystyle \int_1^5 (3x^2+2x+1)\ dx = \displaystyle [x^3+x^2+x]_1^5 = \dots = 152

Quindi

\displaystyle \frac{\text{area}}{\text{AREA}} = \displaystyle \frac{\text{152}}{\text{344}} = 0,44186…

Prova!

PUNTI punti \displaystyle \frac{\text{punti}}{\text{PUNTI}} \displaystyle \text{area} \approx \frac{\text{punti}}{\text{PUNTI}}\cdot \text{AREA}
10 ? ? ?
100 ? ? ?
1000 ? ? ?
10000 ? ? ?
100000 ? ? ?
1000000 ? ? ?
0,44186… 152