.. currentmodule:: brian

.. index::
   pair: example usage; raster_plot
   pair: example usage; NeuronGroup
   pair: example usage; run
   pair: example usage; SpikeMonitor

.. _example-misc_non_reliability:

Example: non_reliability (misc)
===============================

Reliability of spike timing.
See e.g. Mainen & Sejnowski (1995) for experimental results in vitro.

Here: a constant current is injected in all trials.

R. Brette

::

    from brian import *
    
    N = 25
    tau = 20 * ms
    sigma = .015
    eqs_neurons = '''
    dx/dt=(1.1-x)/tau+sigma*(2./tau)**.5*xi:1
    '''
    neurons = NeuronGroup(N, model=eqs_neurons, threshold=1, reset=0, refractory=5 * ms)
    spikes = SpikeMonitor(neurons)
    
    run(500 * ms)
    raster_plot(spikes)
    show()