Pulse broadening in multimode optical fibres with large Δn/n: numerical results