AnQCD: A Mathematica package for calculations in general analytic QCD models

We provide a Mathematica package that evaluates the QCD analytic couplings (in the complex domain) A_ν(Q₂), which are analytic analogs of the powers a(Q²)ν of the underlying perturbative QCD (pQCD) coupling a(Q₂)≡α^s(Q₂)π, in three analytic QCD models (anQCD): Fractional Analytic Perturbation Theory (FAPT), Two-delta analytic QCD (2δanQCD), and Massive Perturbation Theory (MPT). The analytic (holomorphic) running couplings A_ν(Q₂), in contrast to the corresponding pQCD expressions a(Q²)^ν, reflect correctly the analytic properties of the spacelike observables D(Q₂) in the complex Q₂ plane as dictated by the general principles of quantum field theory. They are thus more suited for evaluations of such physical quantities, especially at low momenta |Q₂|∼1GeV².