$d V d time = − (i_K_dr + i_K_Ca + i_K_ATP + i_fast + i_Ca + i_NS + i_NaL) / Cm$

$i_fast = g_fast m_infinity 3.0 h (V − V_fast)$

$m_infinity = 1.0 / (1.0 + e (V − Vm) / Sm)$

$d h d time = (h_infinity − h) / tau_h h_infinity = 1.0 / (1.0 + e (V − Vh) / Sh) tau_h = 1.0 / lamda_h (e (Vh − V) / 2.0 Sh + e (V − Vh) / 2.0 Sh)$

$i_Ca = P_Ca d f_infinity 2.0 F V / R T (Ca_o − Ca_i e 2.0 F V / R T) / (1.0 − e 2.0 F V / R T)$

$d d d time = (d_infinity − d) / tau_d d_infinity = 1.0 / (1.0 + e (V − Vd) / Sd) tau_d = 1.0 / lamda_d (e (Vd − V) / 2.0 Sd + e (V − Vd) / 2.0 Sd)$

$f_infinity = K_Ca / (K_Ca + Ca_i)$

$i_NS = g_NS K_NS 2.0 / (K_NS 2.0 + Ca_lum 2.0) ((V − VNS) / (1.0 − e 0.1 (VNS − V)) − 10.0)$

$i_K_dr = g_K_dr n (V − V_K)$

$d n d time = (n_infinity − n) / tau_n n_infinity = 1.0 / (1.0 + e (V − Vn) / Sn) tau_n = 1.0 / lamda_n (e (Vn − V) / 2.0 Sn + e (V − Vn) / 2.0 Sn)$

$i_K_Ca = g_K_Ca Ca_i 3.0 / (K_Ca 3.0 + Ca_i 3.0) (V − V_K)$

$i_K_ATP = g_K_ATP (V − V_K)$

$i_NaL = g_NaL (V − V_Na)$

$d Ca_i d time = k_rel (Ca_lum − Ca_i) − omega i_Ca + k_Ca Ca_i + k_pump Ca_i d Ca_lum d time = − (k_rel (Ca_lum − Ca_i)) + k_pump Ca_i$