I will discuss the generalized unitarity algorithm used in the public C++ code NJet which
provides a framework for numerical evaluation of one-loop QCD matrix elements. I will then
discuss phenomenological applications at the LHC by combining NJet with the Sherpa
Monte-Carlo program to study observables in 3 and 4 jet production.