"""
Sat Apr 13 10:09:10 CEST 2013
"""
from matplotlib.pyplot import figure, show
from numpy import loadtxt
from IF import SimpleIF as IF
# load current resulting from the simulation
current1 = loadtxt('1contact.out')[:,1]
currentburst = loadtxt('bursts.out')[:,1]
# Integrate and fire parameters
mycell = IF(Cm = 4.9, gl = .16, El = -70, Vinit=-70)
voltage1 = mycell.timecourse(-current1)
voltageBurst = mycell.timecourse(-currentburst)
# one single contact
fig1 = figure(1)
ax1 = fig1.add_subplot(111)
ax1.plot(voltage1)
ax1.set_ylim(-80,10)
# ten contacts
fig2 = figure(2)
ax2 = fig2.add_subplot(111)
ax2.plot(voltageBurst)
ax2.set_ylim(-80,10)
show()
"""
Sat Apr 13 00:27:16 CEST 2013
"""
from matplotlib.pyplot import figure, show
from numpy import loadtxt
from IF import SimpleIF as IF
# load current resulting from the simulation
current1 = loadtxt('1contact.out')[:,1]
current10 = loadtxt('10contacts.out')[:,1]
# Integrate and fire parameters
# Cm in ms, gl in uSiemens and El, Vinit in mV
mycell = IF(Cm = 4.9, gl = .16, El = -70, Vinit=-70)
voltage1 = mycell.timecourse(-current1)
voltage10 = mycell.timecourse(-current10)
# one single contact
fig1 = figure(1)
ax1 = fig1.add_subplot(111)
ax1.plot(voltage1)
ax1.set_ylim(-80,10)
# ten contacts
fig2 = figure(2)
ax2 = fig2.add_subplot(111)
ax2.plot(voltage10)
ax2.set_ylim(-80,10)
show()
