def beta_h(self, V): """Channel gating kinetics. Functions of membrane voltage""" return 1.0/(1.0 + sp.exp(-(V+35.0) / 10.0))