TAUDECAY is a library of helicity amplitudes to simulate polarized tau decays, constructed in the FEYN-RULES and MADGRAPH5 framework. Together with the leptonic mode, the decay library includes the main hadronic modes, τ → ν τ + π, 2π, and 3π, which are introduced as effective vertices by using FEYNRULES. The model file allows us to simulate tau decays when the on-shell tau production is kinematically forbidden. We also demonstrate that all possible correlations among the decay products of pairproduced taus through a Z boson and a scalar/pseudoscalar Higgs boson are produced automatically. The program has been tested carefully by making use of the standard tau decay library TAUOLA.