Photoproduction of multi-jet events at HERA: a Monte Carlo simulation