(* Objectif : calculer la racine carrée de "a" à 1.e−5 près *)
(* méthode de Héron : suite u(n+1) = [ u(n) + a/u(n) ] / 2 *)
PROGRAM heron ;
VAR a, u, u0 : REAL ;
BEGIN
   WriteLn ( "racine carrée de quel nombre ? " ) ;
   ReadLn ( a ) ;
   u := a / 2. ; (* valeur de départ de la suite *)
   REPEAT
      u0 := u ; (* sauvegarde de l'ancienne valeur dans u0 pour le test *)
      u := ( u + (a/u) ) / 2.
   UNTIL ABS(u − u0) < 1.e−5 ;
   WriteLn ( "racine carrée de",a,"=",u );
END

(* Objectif : mettre une formule compliquée en fonction *)
PROGRAM heron ;
FUNCTION fonc ( u : REAL ) : REAL ; (* entrée : u ; sortie : fonc *)
BEGIN
   (* pour une autre suite, il suffit de changer cette fonction *)
   fonc := ( u + (a/u) ) / 2. ;
END
VAR a, u, u0 : REAL ;
BEGIN
   WriteLn ( "racine carrée de quel nombre ? " ) ;
   ReadLn ( a ) ;
   u := a / 2. ; (* valeur de départ de la suite *)
   REPEAT
      u0 := u ; (* sauvegarde de l'ancienne valeur dans u0 pour le test *)
      u := fonc( u )
   UNTIL ABS(u − u0) < 1.e−5 ;
   WriteLn ( "racine carrée de",a,"=",u );
END