- Author:
- WeiweiAi <wai484@aucklanduni.ac.nz>
- Date:
- 2022-03-22 11:47:01+13:00
- Desc:
- Revert Acap to 0.3
- Permanent Source URI:
- https://models.cellml.org/workspace/64f/rawfile/9d4fea2a1b0b60f3c636a2acd2a3c0400c93a8d7/sed-ml/get_init.py
def init_gate(T, V):
import math
phi = 3**((T - 6.3) / 10)
alpha_m = 0.1 * (V + 25) * phi / (math.exp( (V + 25) / 10) -1 )
beta_m = 4 * math.exp( V / 18 ) * phi
alpha_h = 0.07 * math.exp( V / 20 ) * phi
beta_h = phi / (math.exp( (V + 30) /10 ) + 1)
alpha_n = 0.01 * (V + 10) * phi / (math.exp( (V + 10) / 10) -1 )
beta_n = 0.125 * math.exp( V / 80 ) * phi
m_init = alpha_m / (alpha_m + beta_m)
n_init = alpha_n / (alpha_n + beta_n)
h_init = alpha_h / (alpha_h + beta_h)
return m_init, n_init, h_init