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

constante EPSILON : reel <- 0.00000000001
  
{ Methode recursive de calcul                 }
{ de la racine carree d'un nombre reel        }

reel fonction racine(x,a,b)
  Données
    x : reel
    a : reel
    b : reel
  Locales
    res : reel
    m : 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 action
    
{ Calcul d'une racine carree                  }

reel fonction racineCarree(x)
  Données
    x : reel
  Locales
    v : réel
  v <- racine(x,0.0,1.0)
  retourner v
fin action