TITLE: A parameterization method for the computation of invariant tori and their whiskers in quasi-periodic maps: numerical implementation and examples. AUTHORS: Alex Haro (1) and Rafael de la Llave (2) (1) Departament de Matematica Aplicada i Analisi, Universitat de Barcelona, Gran Via 585, 08007 Barcelona (Spain) E-mail: haro@mat.ub.es (2) Department of Mathematics, University of Texas at Austin, Austin, TX 78712 (USA). E-mails: llave@math.utexas.edu ABSTRACT: In this paper we describe the implementation of the numerical algorithms for the computation of invariant manifolds (both tori and their whiskers) in quasi-periodically forced systems based on the parameterization method introduced in two companion papers. We apply the implemented algorithms to some examples considered already in the literature and report on efficiency, accuracy, storage requirements, running times, etc. The methods allow us to continue invariant objects close to the breakdown of their hyperbolicity properties. We find that some of the systems loose hyperbolicity because the stable and unstable bundles approach each other but the Lyapunov multipliers remain away from $1$. Computing several measures of hyperbolicity (the distance between the invariant bundles, and the Lyapunov multipliers) we find power laws with universal exponents. We also observe that, even if the rigorous justifications in the companion papers are developed only for hyperbolic tori, the algorithms work also for elliptic tori in Hamiltonian systems. We can continue these tori and also compute some bifurcations at resonance which may lead to the existence of hyperbolic tori with non-orientable bundles, which we compute.