Location: BG_K_ATP @ fe517be0d80e / matlab_parameter_fitting / calc_IGHK.m
- Author:
- Shelley Fong <s.fong@auckland.ac.nz>
- Date:
- 2022-03-29 13:31:02+13:00
- Desc:
- With channel density checked using patch clamp against kinetic model
- Permanent Source URI:
- https://models.physiomeproject.org/workspace/83a/rawfile/fe517be0d80e2e53c6e08b168d460f7e5065d10a/matlab_parameter_fitting/calc_IGHK.m
function I_GHK = calc_IGHK(G_GHK,V,cKi_st,cKo) R = 8.314; T = 310; F = 96485; V_norm = F*V/R/T; GHK_factor = V_norm./(1-exp(-V_norm)); idx_nan = isnan(GHK_factor); GHK_factor(idx_nan) = 1; I_GHK = G_GHK*GHK_factor.*(cKi_st - cKo*exp(-V_norm)); end