18.01.2015 Views

Appunti di Calcolo Numerico - Esercizi e Dispense - Università degli ...

Appunti di Calcolo Numerico - Esercizi e Dispense - Università degli ...

Appunti di Calcolo Numerico - Esercizi e Dispense - Università degli ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

4. ZERI DI FUNZIONE<br />

Figura 4.6: Il metodo della Regula Falsi applicato alla funzione f (x) = (x/2) 2 − sin(x) con x 0 = 1.3 e x 1 = 1.35<br />

dove p = 1 + 5<br />

2<br />

ɛ n+1 = M 0.618 ɛ 1.618<br />

n<br />

= 1.618 e A è la costante asitontica del metodo <strong>di</strong> Newton-Raphson, da cui<br />

4.8 Confronto tra i meto<strong>di</strong> <strong>di</strong> Newton-Raphson e la Regula Falsi<br />

Sebbene il metodo <strong>di</strong> Newton-Raphson abbia or<strong>di</strong>ne <strong>di</strong> convergenza più elevato della Regula Falsi,<br />

quest’ultimo è computazionalmente più efficiente. Si ha infatti<br />

Metodo p s E<br />

<br />

Newton-Raphson 2 2 2 ≈ 1.414<br />

Regula Falsi 1.618 1 1.618<br />

Ve<strong>di</strong>amo ora come applicare i due meto<strong>di</strong> e le <strong>di</strong>fferenze che si hanno.<br />

Esempio<br />

Esempio 4.8.1 Consideriamo la funzione f (x) = 0 con f (x) = (x/2) 2 − sin(x). La derivata prima è<br />

f ′ (x) = (x/2) − cos(x) Consideriamo come x 0 = 1.3 per entrambi i meto<strong>di</strong> e x 1 = 1.35 per la Regula Falsi.<br />

Come criterio <strong>di</strong> arresto, consideriamo una tolleranza tol = 1.e −8, cioè andremo avanti con le iterazioni<br />

fino a quando troveremo che lo scarto d n = |x n − x n−1 | sarà minore <strong>di</strong> tol . Otteniamo i seguenti risultati<br />

per il metodo <strong>di</strong> Newton-Raphson<br />

n x n f (x n ) f ′ (x n ) d n d n /d<br />

n−1<br />

2<br />

0 1.3 -5.410581854E-01 0.382501171<br />

1 2.714526871831 1.427962127E+00 2.26744846 1.41452687E+00<br />

2 2.084760792766 2.157545986E-01 1.53401376 6.29766079E-01 0.314743565<br />

3 1.944113685369 1.377189572E-02 1.33676314 1.40647107E-01 0.354627390<br />

4 1.933811265085 7.60156095E-05 1.32199993 1.03024203E-02 0.520808008<br />

5 1.933753764621 2.37200355E-09 1.32191743 5.75004640E-05 0.541742396<br />

6 1.933753762827 -1.00668172E-16 1.79436599E-09 0.542710632<br />

50

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!