void Neuron::derivs(double *y, double *dydt) { dV = (Iapp - G_Na * (V - E_Na) - G_K * (V - E_K) - G_L * (V - E_L)) / Cm; dm = (m_inf(V) - m) / tau_m(V); dh = (h_inf(V) - h) / tau_h(V); dn = (n_inf(V) - n) / tau_n(V); }
void ConnorStevens::derivs(double *y, double *dydt) { dV = (Iapp - input(0) * 1e6 - G_Na * (V - E_Na) - G_K * (V - E_K) - G_L * (V - E_L) - G_A * (V - E_A)) * 1000 / Cm; dm = (m_inf(V) - m) / tau_m(V); dh = (h_inf(V) - h) / tau_h(V); dn = (n_inf(V) - n) / tau_n(V); da = (a_inf(V) - a) / tau_a(V); db = (b_inf(V) - b) / tau_b(V); IKA = G_A * (V - E_A) * 1e-6 ; // A }
void Conductance::derivs(double *y,double *dydt) { dm = (m_inf(Vm)-m)/tau_m(Vm); dh = (h_inf(Vm)-h)/tau_h(Vm); dn = (n_inf(Vm)-n)/tau_n(Vm); }