""" 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()