  We present the computation, in lattice QCD, of the renormalization constants
and mixing coefficients of operators that measure the first two moments of DIS
Structure Functions. These calculations have been performed using the
Sheikholeslami-Wohlert O(a) improved ``clover'' action, which is known to
reduce the systematic error associated with the finiteness of the lattice
spacing a. Due to the complexities of the computations, we have developed,
using the computer languages Schoonschip and Form, general codes that are able
to automatically carry out all the analytic lattice manipulations.
