//@code create soma access soma insert hh print "default el_hh = ", el_hh // set el_hh so that the steady state is exactly -70 mV finitialize(-70) // sets v to -70 and m,h,n to corresponding steady state values fcurrent() // set all assigned variables consistent with states // use current balance: 0 = ina + ik + gl_hh*(v - el_hh) el_hh = (ina + ik + gl_hh*v)/gl_hh print "-70 mV steady state el_hh = ", el_hh fcurrent() // recalculate currents (il_hh)