  A new next-to-leading order Monte Carlo program for the calculation of fully
differential jet cross sections in photoproduction is described. The
contributions from both resolved and direct components are included. A
comparison between the theoretical predictions and ZEUS data is presented.
