Models of neurons based on iterative maps allows the simulation of big networks of coupled neurons without loss of biophysical properties such as spiking, bursting or tonic bursting and with an affordable computational effort. A piecewise linear two dimensional map with one fast and one slow variable is used to model spiking–bursting neural behavior. This map shows oscillations similar to other phenomenological models based on maps that require a much higher computational effort. The dynamics of coupled neurons is studied for different coupling strengths.