This paper describes the application of the homotopy perturbations method (HPM) in the computation of invariant measures (IMs) of the non-linear dynamical systems which are characterized by the complex, chaotic behavior. The convergence of the HPM is formally investigated and confirmed, and its efficiency is illustrated in several examples of widely used chaotic maps
