Il seguente algoritmo PhP calcola l'integrale definito di una funzione con
la formula dei trapezi.
La funzione, inserita nella classe $funzione ha tre parametri di ingresso: x1 e x2
(estremi di integrazione), numero (numero di trapezi); ritorna in uscita l'approssimazione x.
function Trapezi($x1, $x2, $numero){
$h = ($x2 - $x1)/$numero;
$x = $x2;
$area = $funzione->Valore($x2);
$x = $x1;
$area += $funzione->Valore($x1);
$area = $area/2;
for ($conta = 1; $conta < $numero; $conta++){
$x += $h;
$area += $funzione->Valore($x);
}
return $area * $h;
}
Nel riquadro sopra a destra è possibile provare l'algoritmo; è consigliabile partire con un basso numero di intervalli e quindi aumentarlo gradualmente osservando come il risultato tende a stabilizzarsi.
Il grafico a sinistra visualizza i trapezi e la curva. Il numero massimo di trapezi
visualizzabili è 100; i calcoli vengono comunque eseguiti anche per numeri di trapezi superiori.
|