Next-to-leading order corrections to e^+ e^- \to four jets

I present a general purpose Monte Carlo program for calculating the next-to-leading order corrections to arbitrary four-jet quantities in electron-positron annihilation. In the current version of the program, some subleading in color terms are neglected. As an example, I calculate the four-jet rate in the Durham scheme as well as the Bengtsson-Zerwas angular distribution at $O({\alpha }_s^3)$ and compare the results to existing data.