Cos(π/9)

[GioArca67]

Il metodo di Cardano mi fornisce le tre soluzioni dell'equazione di terzo grado x³-3x-1=0, fra di esse c'è 2cos(π/9).
Poiché nel caso particolare 2 sono reali e negative (non ho fatto il grafico qualitativo a caso) possono essere escluse e la terza (positiva e vicina a 2) fornisce il valore desiderato.

[PietroBaima]

La soluzione di Cardano per equazioni in cui il termine di secondo grado è nullo vale:



e nel nostro caso specifico diventa:



che semplificata diventa

non molto utile, no?

[GioArca67]

Concordo. Per la determinazione numerica è del tutto inutile.
Ma è il valore esatto di cos(π/9).
Perché da un punto di vista computazionale le radici continue sono meglio ad es di Newton-Raphson

[PietroBaima]

La soluzione di Cardano restituisce il valore esatto, certo, ma non è utile per la sua determinazione numerica.

Perché da un punto di vista computazionale le radici continue sono meglio ad es di Newton-Raphson

Suppongo sia una domanda?
Le radici continue vanno meglio di NR per il fatto che la approssimazione precedente fornisce un nuovo valore di partenza per il calcolo della radice cubica successiva, per cui posso sfruttarlo per migliorare la approssimazione. In questo modo non devo ricalcolare "tutta" la radice cubica da capo, ma solo le cifre errate (e per capire quali siano faccio il rapporto fra due approssimazioni successive e vedo qual è la posizione della prima cifra diversa da zero).
Facendo questo si può dimostrare che la convergenza è migliore di quella del metodo di NR.

Il problema è che non sempre siamo fortunati e riusciamo a trovare una radice continua che approssimi qualcosa (e talvolta lo approssima peggio di NR, ad esempio), mentre il metodo di NR funziona sempre.

Inoltre il problema di NR è che le derivate di seni e coseni restano sempre seni e coseni, quindi in questo caso non ci aiuta molto.

[GioArca67]

Ok, grazie. Sì era una domanda.
Però in questo caso N-R mi molto semplice visto che abbiamo un banale polinomio (x³-3x-1) e la derivata è altrettanto banale.
Partendo da x0=2 in 3 iterazioni ho 3 cifre decimali corrette

[PietroBaima]

Con la radice in 19 iterazioni hai 10 cifre decimali corrette, partendo da 4.

[GioArca67]

Con la radice in 19 iterazioni hai 10 cifre decimali corrette, partendo da 4.

Perché 4 e non 2?

[PietroBaima]

Puoi assumere quello che vuoi, ma se non specifichi niente Mathematica (che io uso di solito) assume 1 come valore inziale con cui calcolare la prima iterazione.
Ponendo 1 si ha che fa 4

[GioArca67]

Visto che se non vedo non credo, ho fatto una prova (quick&dirty).
Non per dimostrare alcun che.
Newton vs Radici Continue...

[PietroBaima]

Certo, se fai il calcolo a macchina è ovvio...