Incoherent $pi^0$ photoproduction on a deuteron in the $Delta(1232)$ region is treated in a three-body scattering approach using separable two-body interactions. Results are presented for total and differential cross sections. It turns out that the role of higher order terms beyond the first order in the multiple scattering series is insignificant, and their inclusion cannot explain the existing discrepancy between theory and experiment.