We use mean field analysis to study bursting in networks of identical, pulse-coupled neurons. The individual neurons are represented using a class of two-dimensional integrate and fire model. The mean field model is a system of switching ordinary differential equations and the transition to bursting involves both standard and nonsmooth bifurcations. The results of the mean field analysis are compared with numerical simulations of large networks. This is joint work with Wilten Nicola.