The Basic Izhikevich Spiking-Neuron Model.

gistfile1.py

for i in range(args.n_timesteps):
    spike_train.append(v)
    
    if v > 0.035:
        v = c
        u = u + d
        
        spikes.append(i)
    
    d_v = (1 / 10.0) * (k * v**2 + 5 * v + 140 - u + 1)
    d_u = (a / 10.0) * (b * v - u)
    
    v += d_v
    u += d_u