FacetsPythonCourse2008: interactive_session_part_1_1.py

from brian import *

tau = 10*ms
Vt = -55*mV
Vr = -70*mV
El = -54*mV

eqs = '''
    dV/dt = -(V-El)/tau : volt
    '''
    
model = Model(equations=eqs,
        threshold=Vt, reset=Vr)

neuron = NeuronGroup(1, model)

M = StateMonitor(neuron,'V', record=True)

run(200*ms)

plot(M.times, M[0])
show()