{ Constante de definition de la precision     }
{ de calcul                                   }

constante EPSILON reel <- 0.00000000001
  
{ Fonction recursive de calcul et retour      }
{ de la racine carree d'un nombre reel        }
{ x : le nombre dont la racine carrée         }
{     est calculée                            }
{ a : la borne minimale de l'intervalle       }
{     de recherche                            }
{ b : la borne minimale de l'intervalle       }
{     de recherche                            }

reel fonction racine(-> reel x,-> reel a,-> reel b)
  reel res
  reel m <- (b+a)/2.0
  si b-a > EPSILON alors
    si m*m > x alors
      res <- racine(x,a,m)
      sinon
      res <- racine(x,m,b)
    fsi
    sinon
    res <- m
  fsi
  retourner res
fin fonction
    
{ Fonction de calcul et retour                }
{ de la racine carree d'un nombre compris     }
{ dans l'intervalle [0.0,1.0]                 }
{ x : le nombre dont la racine carrée         }
{     est calculée                            }

reel fonction racineCarree(-> reel x)
  réel v <- racine(x,0.0,1.0)
  retourner v
fin fonction