%% 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_!!!