The mechanism of incoherent π 0 and η photoproduction from complex nuclei is investigated from 4 to 12 GeV with an extended version of the multicollisional Monte Carlo (MCMC) intranuclear cascade model. The calculations take into account the elementary photoproduction amplitudes via a Regge model and the nuclear effects of photon shadowing, Pauli blocking, and meson-nucleus final-state interactions. The results for π 0 photoproduction reproduced for the first time the magnitude and energy dependence of the measured rations σ γ A /σ γ N for several nuclei (Be, C, Al, Cu, Ag, and Pb) from a Cornell experiment. The results for η photoproduction fitted the inelastic background in Cornell's yields remarkably well, which is clearly not isotropic as previously considered in Cornell's analysis. With this constraint for the background, the η → γ γ decay width was extracted using the Primakoff method, combining Be and Cu data [ η→γ γ = 0.476 (62)