  We report on a high statistics quenched lattice QCD calculation of the
deep-inelastic structure functions $F_1$, $F_2$, $g_1$ and $g_2$ of the proton
and neutron. The theoretical basis for the calculation is the operator product
expansion. We consider the moments of the leading twist operators up to spin
four. Using Wilson fermions the calculation is done for three values of
$\kappa$, and we perform the extrapolation to the chiral limit. The
renormalization constants, which lead us from lattice to continuum operators,
are calculated in perturbation theory to one loop order.
