FacetsPythonCourse2008: interactive_session_part_1_2.py

from brian import *

tau = 10*ms
tau2 = 100*ms
Vt = -55*mV
Vr = -70*mV
El = -54*mV
weight = 5*mV

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

neuron = NeuronGroup(1, model)
neuron2 = NeuronGroup(1, model2)

C = Connection(neuron, neuron2, 'V')
C.connect_full(neuron, neuron2, weight=weight)

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

run(1*second)

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