This a model from the article:
Which model to use for cortical spiking neurons?
Izhikevich EM. IEEE Trans Neural Netw.2004 Sep;15(5):1063-70. 15484883,
Abstract:
We discuss the biological plausibility and computational efficiency of some of the most useful models of spiking and bursting neurons. We compare their applicability to large-scale simulations of cortical neural networks.

The model is according to the paperWhich Model to Use for Cortical Spiking Neurons? Figure1(I) spike latency has been reproduced by MathSBML. The ODE and the parameters values are taken from the a paper Simple Model of Spiking NeuronsThe original format of the models are encoded in the MATLAB format existed in the ModelDB with Accession number 39948

Figure1 are the simulation results of the same model with different choices of parameters and different stimulus function or events. In this model a=0.02; b=0.2; c=-65; d=6; V=-70; u=b*V=0.2*(-70);

This model originates from BioModels Database: A Database of Annotated Published Models. It is copyright (c) 2005-2010 The BioModels Team.
For more information see the terms of use.
To cite BioModels Database, please use Le Novère N., Bornstein B., Broicher A., Courtot M., Donizelli M., Dharuri H., Li L., Sauro H., Schilstra M., Shapiro B., Snoep J.L., Hucka M. (2006) BioModels Database: A Free, Centralized Database of Curated, Published, Quantitative Kinetic Models of Biochemical and Cellular Systems Nucleic Acids Res., 34: D689-D691.