%% calcul d'une valeur approchee de sqrt (n) par dichotomie

/n 3 def

%% syntaxe : a b epsilon {f} dich_solve

%% la valeur trouvee
/valeur 
   0               %% a
   n               %% b
   .00001          %% epsilon  (pas trop petit !!!)
   {dup mul n sub} %% la fonction x |--> x^2 - n
   dich_solve
def

%% une mesure de l'erreur
/erreur 
   valeur n sqrt sub abs 
def

%% affichage des resultats
/chaine 15 string def
setTimesBold
   (Calcul de racine carrée par dichotomie) 
   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