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. These models are built over a phenomenological basis and are mainly implemented by the use of iterative two-dimensional maps that can present neuro-computational properties similar to the usual differential models. A piecewise linear two-dimensional map with one fast and one slow coupled variables is used to model spiking