An improvement on precision of recursive function simulation in IEEE floating point standard is presented. It is shown that the average of rounding towards negative infinite and rounding towards positive infinite yields a better result than the usual standard rounding to the nearest in the simulation of recursive functions. In general, the method improves one digit of precision and it has also been useful to avoid divergence from a correct stationary regime in the logistic map. Numerical studies are presented to illustrate the method.