load_file("init.hoc")
load_file(1, "fig4.ses")

forall nseg=50
soma.diam = 400
forall { L = nseg*200}
forall insert hhint
{tstop=6 tstop_changed()}
cvode.atol(1e-6)
cvode_active(1) // for accurate current calculations
run()