def __init__ (self, type='RS', host=None, cellid=-1):
self.type=type
if host is None: # need to set up a sec for this
self.sec=h.Section(name='izhi2007'+type+str(cellid))
self.sec.L, self.sec.diam, self.sec.cm = 10, 10, 31.831 # empirically tuned
self.izh = h.Izhi2007b(0.5, sec=self.sec)
self.vinit = -60
else:
self.sec = host
self.izh = h.Izhi2007a(0.5, sec=host) # Create a new u,V 2007 neuron at location 0.5 (doesn't matter where)
self.izh.C,self.izh.k,self.izh.vr,self.izh.vt,self.izh.vpeak,self.izh.a,self.izh.b,self.izh.c,self.izh.d,self.izh.celltype = type2007[type]
self.izh.cellid = cellid # Cell ID for keeping track which cell this is
def iz03b_init():
sec03b(0.5).v, iz03b.u = -65, -65 * iz03b.b
fih.append(h.FInitializeHandler(iz03b_init))
# make a 2007a (NMODL) cell
iz07a = h.Izhi2007a(0.5, sec=dummy)
iz07a.Iin = 70
# make a 2007b (section) cell
sec07b = h.Section()
sec07b.L, sec07b.diam = 5, 6.3
iz07b = h.Izhi2007b(0.5, sec=sec07b)
iz07b.Iin = 70
def iz07b_init():
sec07b.v = -60
fih.append(h.FInitializeHandler(iz07b_init))
# make a 2007b (section) cell using the Wrapper
iz07bw = izh07.IzhiCell() # defaults to RS
iz07bw.izh.Iin = 70
fih.append(h.FInitializeHandler(iz07bw.init))
# vectors and plot