%% calcul d'une valeur approchee de sqrt (n) par la methode de Newton
/n 2 def
%% la valeur trouvee
/valeur
n %% x0
.00001 %% epsilon
{2 mul} %% f' (x) = 2x
{dup mul n sub} %% f (x) = x^2 - n
newton_solve
def
%% une mesure de l'erreur
/erreur
valeur n sqrt sub abs
def
(valeur =~ ) valeur
(erreur =~ ) erreur
%% une methode brutale pour visualiser le contenu de la pile d'operande
BOUM_!!!
|