  In this talk we review the GOLEM approach to one-loop calculations and
present an automated implementation of this technique. This method is based on
Feynman diagrams and an advanced reduction of one-loop tensor integrals which
avoids numerical instabilities. We have extended our one-loop integral library
golem95 with an automated one-loop matrix element generator to compute the
virtual corrections of the process $q\bar{q}\to b\bar{b}b\bar{b}$. The
implementation of the virtual matrix element has been interfaced with
tree-level Monte Carlo programs to provide the full result for the above
process.
