69
and
vω = 2vmid - v^-1∖ (3.19)
where, for current injection, vmid is the solution to the linear system
Lωvmid = rω - Uτ(G^.(v - Eβ)), (3.20)
where
LÛ) = 211 At — H ÷ Rdiag(Φ(wθ))e)Z
r0) = 2v(>-1)∕∆f + RΦ(w^)Ei.
In the case of synaptic conductance, vmjd is the solution to the linear system
(1ω + Uτdiag(Gw)U)vmid = rω + Uτ(Gω.E.s). (3.21)
3.2.3 Numerical Solution of the Reduced System
We have written a MATLAB software suite that loads morphology and channel
kinetics and distributions, simulates the full system and obtains POD and DEIM
snapshots, and then simulates the reduced system and displays the response of these
models to the desired inputs. Gating variable evaluation was performed via a look-
up table similar to (Hines, 1984). These codes are available from the authors upon