NEURON {
  SUFFIX kChan
  USEION k READ ek WRITE ik
  RANGE conductance
}

PARAMETER {
  conductance = 0.00000001 (mS)
}


STATE { gates_n_q }

INITIAL {
  LOCAL gates_n_reverseRate_r, gates_n_forwardRate_x, gates_n_forwardRate_r

  gates_n_reverseRate_r = 0.125 * exp(-0.0125 * (65 + v))
  gates_n_forwardRate_x = 0.1 * (55 + v)
  gates_n_forwardRate_r = 0.10000000149011612
  if (gates_n_forwardRate_x != 0) {
    gates_n_forwardRate_r = 0.10000000149011612 * gates_n_forwardRate_x / (1 - exp(-gates_n_forwardRate_x))
  }
  gates_n_q = gates_n_forwardRate_r / (gates_n_forwardRate_r + gates_n_reverseRate_r)
}

DERIVATIVE dstate {
  LOCAL gates_n_reverseRate_r, gates_n_forwardRate_x, gates_n_forwardRate_r

  gates_n_reverseRate_r = 0.125 * exp(-0.0125 * (65 + v))
  gates_n_forwardRate_x = 0.1 * (55 + v)
  gates_n_forwardRate_r = 0.10000000149011612
  if (gates_n_forwardRate_x != 0) {
    gates_n_forwardRate_r = 0.10000000149011612 * gates_n_forwardRate_x / (1 - exp(-gates_n_forwardRate_x))
  }
  gates_n_q' = gates_n_forwardRate_r - gates_n_q * (gates_n_forwardRate_r + gates_n_reverseRate_r)
}

BREAKPOINT {
  SOLVE dstate METHOD cnexp
  LOCAL gates_n_fcond, g

  gates_n_fcond = gates_n_q*gates_n_q
  g = conductance*gates_n_fcond*gates_n_fcond
  ik = g * (v - ek)
}