Ecrire une fonction PLX qui a pour arguments d’entrée:
x0 : point initial
nmax : nombre d’itérations maximal a ne pas dépasser
err : précision voulue
Ses arguments de sortie sont:
xn : solution approchée après ...
ni : nombre d’itérations de sortie
Cette fonction calcule les termes de la suite suivante
x(k+1)= (1/2)*(x(k)+2/(x(k)))
jusqu’à avoir au moins la précision voulue err dans le cas d'un nombre élevé d’itérations ne pas dépasser nmax
l'erreur entre deux itérations est calculée par erreur=|x(k+1)-x(k)| a comparer avec err
indice: utiliser les commandes suivantes:
[~,ni], [xi,~], format ..., rand, eps, abs, norm, &&
1) Calculer les termes pour différents x0 : 0, 1, 2, rand, sqrt(2), ln(2).
avec nmax=1000; et err=10^(-4) .
2) Quels sont les chiffres significatifs obtenus.