# Implementation of Izhikevich neuron model

**URL:** <https://brian.discourse.group/t/implementation-of-izhikevich-neuron-model/578>\
**Category:** Science, Projects, and Showcases\
**Created:** [13 January 2022 22:15 UTC](https://brian.discourse.group/t/implementation-of-izhikevich-neuron-model/578 "2022-01-13T22:15:06Z")\
**Posts on this page:** 1\
**Showing post:** 18

<div class="post-metadata">

**Author:** ![mstimberg](https://yyz2.discourse-cdn.com/free1/user_avatar/brian.discourse.group/mstimberg/32/11_2.png) [@mstimberg](https://brian.discourse.group/u/mstimberg)\
**Post date:** [10 October 2022 09:30 UTC](https://brian.discourse.group/t/implementation-of-izhikevich-neuron-model/578/18 "2022-10-10T09:30:27Z")

</div>

Hi @shehnaz. This is a common reason for confusion, and we mention this issue briefly in our tutorial ([Introduction to Brian part 1: Neurons — Brian 2 2.5.4 documentation](https://brian2.readthedocs.io/en/stable/resources/tutorials/1-intro-to-brian-neurons.html#a-simple-model)). When describing physical systems, it is not meaningful to have a dimensionless right-hand-side for differential equations, even when all variables are dimensionless. An equation like \frac{dx}{dt} = -x does not make sense for a physical system: For a value of x=1, its value is to decrease by 1 , but should it decrease by this amount every second, every millisecond, every year, …? I explained this a bit in a [post asking about the Van-der-Pol oscillator](https://brian.discourse.group/t/coupled-van-der-pol-oscillators/230/2):

> [@Coupled Van der Pol oscillators](https://brian.discourse.group/t/coupled-van-der-pol-oscillators/230/2):
>
> Just to make the link to the mathematic formulation and e.g. the Wikipedia article more clear: By setting `tau = 10*ms`, you are defining your “unit of time” to be 10ms. This means when you simulate for 100ms, you simulate for “10 units of time”. If you want to plot something like [https://en.wikipedia.org/wiki/Van\_der\_Pol\_oscillator#/media/File:Vanderpol\_mu=5.svg](https://en.wikipedia.org/wiki/Van_der_Pol_oscillator#/media/File:Vanderpol_mu=5.svg), you have to divide the time steps provided by Brian by 10ms to get the same time units on the x axis. Finally, note that Brian by default uses a simulation time step of 0.1ms (a reasonable time step for many neuron models), so with your definitions that means a simulation time step of 0.01 time units (often written as “h=0.01” in mathematical formulations).

In Izhikevich’s model, he states that all time constants are measured in ms (and similarly, membrane potentials are measured in mV, etc.), so this is the reason to divide everything by `ms`.

Hope that makes things clearer?

---

_[View the full topic](https://brian.discourse.group/t/implementation-of-izhikevich-neuron-model/578)._