  We have calculated inclusive two-jet cross sections in next-to-leading order
QCD for direct photoproduction in low $Q^2$ $ep$ collisions at HERA. Infrared
and collinear singularities in real and virtual contributions are cancelled
with the phase space slicing method. Analytical formulas for the different
contributions giving the dependence on the slicing parameter are presented.
Various one- and two-jet distributions have been computed demonstrating the
flexibility of the method.
