Quadratic Integrate-and-Fire

$$\tau \frac{d}{dt}u=a\mathrm{\_}0(u-u\mathrm{\_}rest)(u-u\mathrm{\_}c)+RI$$

Default parameters: $\tau =10$ ms, $a_0=0.1m{V}^{-1}$, $u_c=-55$ mV, $u_{rest}=-70$ mV, $R=10M\mathrm{\Omega }$, $I=2.2nA$, firing threshold $u_{thr}=-40mV$.

The source code is Open Source and can be found on GitHub.