  Many natural phenomena exhibit a stochastic nature that one attempts at
modeling by using stochastic processes of different types. In this context,
often one is interested in investigating the memory properties of the natural
phenomenon at hand. This is usually accomplished by computing the
autocorrelation function of the numerical series describing the considered
phenomenon. Often, especially when considering real world data, the
autocorrelation function must be computed starting from a single numerical
series: i.e. with a time-average approach. Hereafter, we will propose a novel
way of evaluating the time-average autocorrelation function, based on the
preliminary evaluation of the quantity N({\tau},g{\mu},g{\nu}), that, apart
from normalization factors, represents a numerical estimate, based on a single
realization of the process, of the 2-point joint probability density function
P(x2,{\tau};x1,0). The main advantage of the proposed method is that it allows
to quantitatively assess what is the error that one makes when numerically
evaluating the autocorrelation function due to the fact that any simulated time
series is necessarily bounded. In fact, we show that, for a wide class of
stochastic processes admitting a nonlinear Langevin equation with white noise
and that can be described by using a Fokker-Planck equation, the way the
numerical estimate of the autocorrelation function converges to its theoretical
prediction depends on the pdf tails. Moreover, the knowledge of
N({\tau},g{\mu},g{\nu}) allows to easily compute the process histogram and to
characterize processes with multiple timescales. We will show the effectiveness
of our new methodology by considering three stochastic processes whose
autocorrelation function and two-point probability density function are both
known in an analytical or numerical form, thus allowing direct comparisons.
