  We have implemented in the RAPGAP program a previously derived subtraction
method for next-to-leading order (NLO) corrections in Monte-Carlo event
generators, and we show results for jet production in deep inelastic
scattering. At small x_{bj}, NLO corrections are comparable to the LO cross
section, because of the large gluon density. We devise a jet observable that is
particularly sensitive to the treatment of parton kinematics, which our method
is intended to treat correctly. We compare the results of the calculation with
LO calculations and with the previously existing treatment of NLO processes in
RAPGAP. Substantial corrections to the event distribution are found.
