  NLO QCD fits using the DGLAP formalism have been made to the high precision
ZEUS $e^+p$ reduced cross-section data and to fixed target structure function
data, in order to detrmine parton distribution functions and the value of
$\alpha_s(M^2_Z)$, taking full account of correlated systematic errors.
Preliminary results are reported.
