Algoritmo per il calcolo della radice n-esima

La radice  n-esima di di un numero reale positivo A,  ${\displaystyle {\sqrt[{n}]{A}}}$ , e' la soluzione reale positiva  della equazione

${\displaystyle x^{n}=A}$

(per un numero n>1 intero ci sono  n soluzioni nel campo complesso, tuttavia una sola e' reale e positiva se ${\displaystyle A>0}$.

Questa voce - adattata dall' originale in inglese di wikipedia - e' su  un metodo numerico per il calcolo di questa radice  che converge velocemente. I passi dell' algoritmo sono:

1. Si prova a stimare un valore iniziale di partenza  ${\displaystyle x_{0}}$
2. Si pone ${\displaystyle x_{k+1}={\frac {1}{n}}\left[{(n-1)x_{k}+{\frac {A}{x_{k}^{n-1}}}}\right]}$. Questo equivale a  ${\displaystyle \Delta x_{k}={\frac {1}{n}}\left[{\frac {A}{x_{k}^{n-1}}}-x_{k}\right];x_{k+1}=x_{k}+\Delta x_{k}}$.
3. Si ripete lo step 2 fino a che si raggiunge la precisione desiderata, cioe'  ${\displaystyle |\Delta x_{k}|<\epsilon }$ .

Un caso speciale e' il calcolo numerico della radice quadrata ponendo n=2:

${\displaystyle x_{k+1}={\frac {1}{2}}\left(x_{k}+{\frac {A}{x_{k}}}\right)}$

${\displaystyle x_{k+1}={\frac {1}{2}}\left(x_{k}+{\frac {A}{x_{k}}}\right)}$

La derivazione dell' algoritmo si basa sul metodo numerico di Newton-Raphson.

Derivazione  dal metodo di Newton-Raphons

Il metodo delle tangenti o di Newton-Raphson e' un  metodo per trovare numericamente lo zero di una funzione f(x). Lo schema generale e': 

1. Partire da una stima iniziale ${\displaystyle x_{0}}$ 
2. ${\displaystyle x_{k+1}=x_{k}-{\frac {f(x_{k})}{f'(x_{k})}}}$ 
3. Ripetere lo step 2 fino a che si raggiunga la precisione desiderata.

Il calcolo numerico della radice n-esima si puo' concepire come la ricerca di uno zero della funzione:

${\displaystyle f(x)=x^{n}-A}$ 

La cui derivata e':

${\displaystyle f^{\prime }(x)=nx^{n-1}}$ 

In questo modo si costruisce l' iterazione: 

${\displaystyle x_{k+1}=x_{k}-{\frac {f(x_{k})}{f'(x_{k})}}}$ 

${\displaystyle =x_{k}-{\frac {x_{k}^{n}-A}{nx_{k}^{n-1}}}}$ 

${\displaystyle =x_{k}-{\frac {x_{k}}{n}}+{\frac {A}{nx_{k}^{n-1}}}}$ 

${\displaystyle ={\frac {1}{n}}\left[{(n-1)x_{k}+{\frac {A}{x_{k}^{n-1}}}}\right]}$ 

Voci correlate

Metodi per il calcolo della radice quadrata 

Analisi numerica 

Riferimenti bibliografici

• Kendall E. Atkinson, An introduction to numerical analysis, Wiley, 1989..