  We have implemented a systematic procedure for combining parton shower
algorithms with next-to-leading order QCD calculations for the case of jet
production in deep-inelastic electron-proton scattering. Using this method we
have computed inclusive jet cross sections and jet shapes for the case of
single-jet production and compared them to data from the ZEUS collaboration at
HERA. We found good agreement between the data and our calculations, both for
the jet shapes and the inclusive spectra.
