Neuronal Network

rete.hoc

create soma[3]

soma[0] {pt3dclear()
        pt3dadd(-2.5, 0, 0, 5)
        pt3dadd(2.5, 0, 0, 5)
}

soma[1] {pt3dclear()
  pt3dadd(-17.5, 0, 0, 5)
  pt3dadd(-12.5, 0, 0, 5)
}

soma[2] {pt3dclear()
  pt3dadd(12.5, 0, 0, 5)
  pt3dadd(17.5, 0, 0, 5)
}

SimulationNetwork.hoc

load_file("rete.hoc")

celsius=35
objref stim[3], ns, nc[4]
access soma[0]
ns = new NetStim()
ns.start=10
ns.number=1
ns.interval=50
ns.noise=1
k=0

for k=0, 2 {
  soma[k] {
    stim[k] = new Exp2Syn(.5)
    stim[k].tau1=1
    stim[k].tau2=10
    stim[k].e=0
  }
}

nc[3]= new NetCon(ns, stim, 0, 0, 2.5e-5)

for k=0, 1 {
  soma[k] nc[k]= new NetCon(&v(.5), stim[k+1], 0, 0, 2.5e-5)
}

soma[2] nc[2]= new NetCon(&v(.5), stim[0], 0, 0, 2.5e-5)

forall { insert pas cm=1 g_pas=1/30000 e_pas=-65 }
forall { insert nax gbar_nax=0.01 }
forall { insert kdr gkdrbar_kdr=0.01 }