! Time course of GABA-B, from Otis, de Koninck & Mody (1993) and proportional
! to that used in Traub et al. 1993 pyramidal cell model, J. Physiol.
                subroutine otis (t,value)

                real*8 t, value

              if (t.le.10.d0) then
                value = 0.d0
              else
            value = (1.d0 - dexp(-(t-10.d0)/38.1d0)) ** 4

       value = value * (10.2d0 * dexp(-(t-10.d0)/122.d0) +
     &    1.1d0 * dexp(-(t-10.d0)/587.d0))
              endif

                 end