Simulatie van het volledige model van Wilson

Wilson voegde aan zijn basismodel voor neocorticale neuronen van zoogdieren twee termen toe: Een term, het T-kanaal genoemd, is verbonden met een spanningsafhankelijk calciumkanaal en een tweede term, het H-kanaal genoemd, is verbonden met een calcium-afhankelijk kaliumkanaal. Het stelsel van differentiaalvergelijkingen wordt dan:\[\begin{aligned} {C_m} \cdot \frac{{\dd V}}{{\dd t}} &= - {g_{{\rm{Na}}}}(V) \cdot (V - {E_{{\rm{Na}}}}) - {g_{\rm{K}}} \cdot R \cdot (V - {E_{\rm{K}}})\\[0.25cm] &\phantom{=} - {g_{\rm{T}}} \cdot T \cdot (V - {E_{\rm{T}}}) - {g_{\rm{H}}} \cdot H \cdot (V - {E_{\rm{H}}}) + {I_{{\rm{stim}}}} \\[0.25cm] \tau _{\rm{R}} \cdot \frac{{\dd R}}{{\dd t}} &= {R_{\infty}}(V) - R\\[0.25cm] \tau _{\rm{T}} \cdot \frac{{\dd T}}{{\dd t}} &= {T_{\infty}}(V) - T\\[0.25cm] \tau _{\rm{H}} \cdot \frac{{\dd H}}{{\dd t}} &= 3T-H\\[0.25cm] \end{aligned}\] waarbij \(V\) de membraanpotentiaal is, \(E_X\) een referentiewaarde voor een gegeven kanaal \(X\) is, en \[\begin{aligned} {g_{{\rm{Na}}}}(V) &= 17.8 + 0.476\,\,V + 33.8 \cdot {10^{ - 4}}\,{V^2}\\[0.25cm] {R_\infty }(V) &= 1.24 + 0.037\,V + 3.2\, \cdot {10^{ - 4}}\,{V^2}\\[0.25cm] {T_\infty }(V) &= 4.205 + 0.116\,V + 8.1\, \cdot {10^{ - 4}}\,{V^2}\\[0.25cm] \tau_{\rm{T}} &= 14\\[0.25cm] \tau_{\rm{H}}&=45\\[0.25cm] g_{\rm{K}} & =26\end{aligned}\] Het volledige model van Wilson wordt net zo herschaald als het basismodel om numerieke problemen te omzeilen.