🚀 go-pugleaf

Some time ago I wrote some procedures for rounding real numbers:

PROCEDURE round (x:LONGREAL; n:INTEGER): LONGREAL;
(* Round x to n digits. *)

PROCEDURE roundp (x:LONGREAL; n:INTEGER): LONGREAL;
(* Round x to n digits behind the decimal point. *)

PROCEDURE RealToStr (x:LONGREAL; n:INTEGER; VAR s:ARRAY OF CHAR);
(* Convert x to s, rounding to a total of n characters.
   n>0: floating point, else exponential.
*)

I am not so happy with my solutions, because I use logarithms and 
exponential functions to generate power-of-10-factors for left- and 
right-shifting and ENTIER for truncation. Sometimes arithmetic errors in 
those functions give me a result of 9.99 when I want to round 10. 

Does anybody know procedures or algorithms or references about that topic? 


Thanks in advance.

-----------------------------------------------------------
Fundacao CERTI			Tel.:	+55 / 48 / 234-3000
Dr. Frank Hrebabetzky		Fax:	+55 / 48 / 234-1514
Campus da UFSC, CP 5053		email:	fh@mbox1.ufsc.br
Florianopolis, SC		        fh@certi.ufsc.br
88040-970
Brasil