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

%% affichage des resulats
/chaine 15 string def
setTimesBold
   (Calcul de racine carrée par Newton) 
   xmin xmax add 2 div
   .7 ymax mul
   cctext

setTimes
   (racine carrée de ) xmin 0 brtext
   n chaine cvs show
   ( => ) show
   valeur chaine cvs show

   (erreur => ) xmin -1 brtext
   erreur chaine cvs show