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
|
COMMENT
Test for derivimplicit solver with sympy and Eigen. Based on cacum.mod
ENDCOMMENT
NEURON {
SUFFIX cacum
USEION ca READ ica WRITE cai
RANGE depth, tau, cai0
}
UNITS {
(mM) = (milli/liter)
(mA) = (milliamp)
F = (faraday) (coulombs)
}
PARAMETER {
depth = 1 (nm) : assume volume = area*depth
tau = 10 (ms)
cai0 = 50e-6 (mM) : Requires explicit use in INITIAL
: block for it to take precedence over cai0_ca_ion
: Do not forget to initialize in hoc if different
: from this default.
}
ASSIGNED {
ica (mA/cm2)
}
STATE {
cai (mM)
}
INITIAL {
cai = cai0
extra_solve()
}
BREAKPOINT {
SOLVE integrate METHOD derivimplicit
}
PROCEDURE extra_solve() {
SOLVE integrate METHOD derivimplicit
}
DERIVATIVE integrate {
cai' = -ica/depth/F/2 * (1e7) + (cai0 - cai)/tau
}
|
