Sim1.hoc

load_file("nrngui.hoc")

/* axon morphology */

create s, a[10], b[10]

forall {
	Ra = 110
	cm = 1
}

access s
    nseg = 1
	diam = 10
	L    = 10

for i = 0, 9 a[i] {
	nseg = 100
	diam = 0.2
	L    = 100
}

for i = 0, 9 b[i] {
	nseg = 4
	diam = 4
	L    = 4
}

       connect s(1), a[0](0)

for i = 0, 9 {connect b[i](0), a[i](1)		
}
		
for i = 0, 8 {connect a[i+1](0), b[i](1)		
}

/* channel */

forall {insert hhmfb
        insert KIn
        gnabar_hhmfb = 0.05 /* S/cm^2 */
		gkbar_hhmfb = 0 /* S/cm^2 */
		gl_hhmfb = 0.0001 /* S/cm^2 */
		gkbar_KIn = 0.036 /* S/cm^2 */
		ena = 50    /* mV */
        ek = -85    /* mV */
        celsius = 25
}

forsec "s" gnabar_hhmfb = 0.01 /* S/cm^2 */ 

// forsec "b[7]" gkbar_KIn = 0    /* S/cm^2*/
// forsec "b[8]" gkbar_KIn = 0    /* S/cm^2*/
// forsec "b[9]" gkbar_KIn = 0    /* S/cm^2*/

// forsec "a[8]" gkbar_KIn = 0    /* S/cm^2*/
// forsec "a[9]" gkbar_KIn = 0    /* S/cm^2*/

/* single shock */

objectvar stim
s stim = new IClamp(0.5)
stim.del = 50
stim.dur = 0.2
stim.amp = 1

load_file("rec1.ses")