76879/wc.mod at master · ModelDBRepository/76879 · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
: $Id: wc.mod,v 1.4 2006/07/28 17:19:41 billl Exp $
TITLE wc.mod Wilson-Cowan equations

COMMENT
Wilson-Cowan equations
ENDCOMMENT

NEURON {
  POINT_PROCESS WC
  RANGE ii0,ii1,ie0,ie1,w,aee,aie,aei,aii,ze,zi,tau,i_e,i_i
}

PARAMETER {
  aee = 10
  aie = 8
  aei = 12
  aii = 3
  ze = 0.2
  zi = 4
  tau = 1
  ie0 = 0
  ie1 = 0
  ii0 = 0
  ii1 = 0
  w = 0.25
}

STATE { uu vv }

ASSIGNED {
  i_e
  i_i
}

BREAKPOINT {
  SOLVE states METHOD derivimplicit
}

INITIAL {
  uu = 0.1
  vv = 0.05
}

DERIVATIVE states {
  i_e = ie0+ie1*sin(w*t)
  i_i = ii0+ii1*sin(w*t)
  uu' = -uu+f(aee*uu-aie*vv-ze+i_e)
  vv' = (-vv+f(aei*uu-aii*vv-zi+i_i))/tau
}

FUNCTION f(x) {
  f = 1/(1+exp(-x))
}