EERAD3: Event shapes and jet rates in electron-positron annihilation at   order $\alpha_s^3$
  The program EERAD3 computes the parton-level QCD contributions to event shapes and jet rates in electron-positron annihilation through to order $\alpha_s^3$. For three-jet production and related observables, this corresponds to next-to-next-to-leading order corrections, and allows for precision QCD studies. We describe the program and its usage in detail. 
A. Gehrmann-De Ridder, T. Gehrmann, E.W.N. Glover, G. Heinrich
hep-ph