TITLE Single channel with random decay constant
:  National Biomedical Simulation Resource, 1-31-93
:  A simple exponential-dacay channel model with a random rate coefficient.
:    The gate variable being set to 1 initially and the channel closes
:    when gate drops below a threshold.
:    $Id:$

PARAMETER
{
  cutoff = 0.3
}

INDEPENDENT
{
  time FROM 0 TO 10 WITH 100
}

PLOT gate,current VS time

STATE
{
  gate FROM 0 TO 1.5 START 1
}

ASSIGNED
{
  decay_rate
  current FROM 0 TO 1.5
}

INITIAL
{
    decay_rate = normrand(0.5, 0.2)
}

KINETIC rates
{
~ gate ->  (decay_rate)
}

BREAKPOINT
{
  SOLVE rates
  IF (gate >= cutoff)
  {
    current = 1
  }
  ELSE
  {
    current = 0
  }
}