  We present a numerical method for calculation of Ruelle-Pollicott resonances
of dynamical systems. It constructs an effective coarse-grained propagator by
considering the correlations of multiple observables over multiple timesteps.
The method is compared to the usual approaches on the example of the perturbed
cat map and is shown to be numerically efficient and robust.
