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 ////////
soma { L = 30 diam = 30 }
apical { L = 30 diam = 5 }

for i=0,NDEND-1 {
ap[i] { L = 15 diam = 2 }
oblique[i] { L = 15+5*i diam = 1 }
}