We consider position random maps T = {τ 1 (x), τ 2 (x), . . . , τ K (x); p 1 (x), p 2 (x), . . . ,We assume that the random map T posses a density function f * of the unique absolutely continuous invariant measure (acim) µ * . In this paper, first, we present a general numerical algorithm for the approximation of the density function f * . Moreover, we show that Ulam's method is a special case of the general method. Finally, we describe a Monte-Carlo and a Quasi Monte Carlo implementations of Ulam's method for the approximation of f * . The main advantage of these methods is that we do not need to find the inverse images of subsets under the transformations of the random map T .