load_file("nrngui.hoc")


/////// topology ////////
NDEND = 10
create soma, apical, ap[NDEND], oblique[NDEND]
access soma
connect apical(0), soma(1)
connect ap[0](0), apical(1)
connect oblique[0](0), apical(1)

for i=1,NDEND-1 {
connect ap[i](0), ap[i-1](1)
connect oblique[i](0), ap[i-1](1)
}

/////// geometry ////////
forall pt3dclear()
soma {
  pt3dadd(0, 0, 0, 30)
  pt3dadd(30, 0, 0, 30)
}
apical {
  pt3dadd(30, 0, 0, 5)
  pt3dadd(60, 0, 0, 5)
}

for i=0,NDEND-1 {
  ap[i] {
    pt3dadd(60+i*15, 0, 0, 2)
    pt3dadd(60+(i+1)*15, 0, 0, 2)
  } 
  oblique[i] {
    pt3dadd(60+i*15, 0, 0, 1)
    pt3dadd(60+i*15, 15+5*i, 0, 1)
  }
}