{ "ode": [
    "dv/dt = (gl*(El-v) - g_na*(m*m*m)*h*(v-ENa) - g_kd*(n*n*n*n)*(v-EK) + I + ih)/Cm",
    "dm/dt = (0.32*(mV**-1)*(13.*mV-v+VT)/(exp((13.*mV-v+VT)/(4.*mV))-1.)/ms*(1-m)-0.28*(mV**-1)*(v-VT-40.*mV)/(exp((v-VT-40.*mV)/(5.*mV))-1.)/ms*m)",
    "dn/dt = 0.032*(mV**-1)*(15.*mV-v+VT)/(exp((15.*mV-v+VT)/(5.*mV))-1.)/ms*(1.-n)-.5*exp((10.*mV-v+VT)/(40.*mV))/ms*n",
    "dh/dt = 0.128*exp((17.*mV-v+VT)/(18.*mV))/ms*(1.-h)-4./(1+exp((40.*mV-v+VT)/(5.*mV)))/ms*h",
    "dr/dt = (rinf-r)/rtau"
  ],
  "init_states":   {
    "v" : "-65 * mV",
    "m" : "0",
    "n" : "0",
    "h" : "0",
    "r" : "0"
  },
  "definitions":   {
    "ih"      : "ghbar*r*(Eh-v)"
  },
  "functions":   {
    "rinf"   : "1. / (1+exp((V_half_h - v/mV) / k_h))",
    "rtau"   : "0.1*(C_base_h + C_amp_h * exp(-(V_max_h-v/mV)**2./sigma_h**2))"
  },
  "parameters":   {
    "area"      : "20000*umetre**2",
    "Cm"        : "(1.0 * (ufarad/cm**2)) * area",
    "gl"        : "5e-5*siemens*cm**-2 * area",
    "El"        : "-65*mV",
    "EK"        : "-90*mV",
    "ENa"       : "50*mV",
    "g_na"      : "100*msiemens*cm**-2 * area",
    "g_kd"      : "30*msiemens*cm**-2 * area",
    "VT"        : "-63*mV",
    "ghbar"     : "40. * nS",
    "Eh"        : "-1.*mV",
    "V_half_h"  : "-95.0",
    "k_h"       : "-5.5",
    "V_max_h"   : "-75.",
    "sigma_h"   : "20.",
    "C_amp_h"   : "50.",
    "C_base_h"  : "10.",
    "I"         : "0.7*nA"
  },
  "bibkey": "none",
  "comments":   {
    "source": "strongly tweaked h-current"
  }
}