Model Mathematics

Component: environment

Component: Vm_n

ddtimeVm_n=0.001-i_NaT+i_NaP+i_leakNa+i_KDR+i_KA+i_leakK+i_leakf+i_NaKATPase_n+i_appCm i_app=3500iftimet0timet10otherwise

Component: i_NaT

i_NaT=gNaTm3hVm_n-ENa_n J_NaT=gNaTm3hVm_n-ENa_nF

Component: i_NaT_m_gate

alpha_m=0.32-Vm_n-51.9-0.25Vm_n+12.975-1 beta_m=0.28Vm_n+24.890.2Vm_n+4.978-1 ddtimem=alpha_m1-m-beta_mm

Component: i_NaT_h_gate

alpha_h=0.128-0.056Vm_n+2.94 beta_h=4-0.2Vm_n+6+1 ddtimeh=alpha_h1-h-beta_hh

Component: i_NaP

i_NaP=gNaPm2hVm_n-ENa_n J_NaP=gNaPm2hVm_n-ENa_nF

Component: i_NaP_m_gate

alpha_m=1tau_activation1-0.143Vm_n+5.67+1 beta_m=1tau_activation-0.143Vm_n+5.67-0.143Vm_n+5.67+1 ddtimem=alpha_m1-m-beta_mm

Component: i_NaP_h_gate

alpha_h=5.12e-8-0.056Vm_n+2.94 beta_h=1.6e-6-0.2Vm_n+8+1 ddtimeh=alpha_h1-h-beta_hh

Component: i_KDR

i_KDR=gKDRn2Vm_n-EK_n J_KDR=gKDRn2Vm_n-EK_nF

Component: i_KDR_n_gate

alpha_n=0.016-Vm_n-34.9-0.2Vm_n+6.98-1 beta_n=0.25-0.025Vm_n+1.25 ddtimen=alpha_n1-n-beta_nn

Component: i_KA

i_KA=gKAm2hVm_n-EK_n J_KA=gKAm2hVm_n-EK_nF

Component: i_KA_m_gate

alpha_m=0.02-Vm_n-56.9-0.1Vm_n+5.69-1 beta_m=0.0175Vm_n+29.90.1Vm_n+2.99-1 ddtimem=alpha_m1-m-beta_mm

Component: i_KA_h_gate

alpha_h=0.016-0.056Vm_n+4.61 beta_h=0.5-0.2Vm_n+11.98+1 ddtimeh=alpha_h1-h-beta_hh

Component: i_NaKATPase_n

i_NaKATPase_n=J_NaKATPase_nF J_NaKATPase_n=I_NaKATPase_n_maxNan1.5Nan1.5+KmNa1.5KoKo+KmK

Component: i_leakNa

i_leakNa=gleakNaVm_n-ENa_n J_leakNa=gleakNaVm_n-ENa_nF

Component: i_leakK

i_leakK=gleakKVm_n-EK_n J_leakK=gleakKVm_n-EK_nF

Component: i_leakf

i_leakf=gleakfVm_n-Ef_n

Component: Vm_g

Vm_g=gNaENa_g+gKEK_g+gClECl_g+gNBCENBC_ggNa+gK+gCl+gNBC-FJ_NaKATPase_ggNa+gK+gCl+gNBC

Component: J_Na

J_Na=gNaVm_g-ENa_gF

Component: J_K

J_K=gKVm_g-EK_gF

Component: J_NaKATPase_g

J_NaKATPase_g=I_NaKATPase_g_maxNag1.5Nag1.5+KmNa1.5KoKo+KmK

Component: J_NBC

J_NBC=gNBCVm_g-ENBC_gF

Component: J_NKCC1

J_NKCC1=Ko-P_Ko10Ko-P_Ko10+0.0310PsigNKCC1lnKoKgNaoNagCloClg2FifKo>P_Ko0otherwise

Component: N_Nag

ddtimeN_Nag=dN_Nag_dt dN_Nag_dt=0.01-J_Na-3J_NaKATPase_g+J_NKCC1+J_NBC

Component: N_Kg

ddtimeN_Kg=dN_Kg_dt dN_Kg_dt=0.01-J_K+2J_NaKATPase_g+J_NKCC1

Component: wg

ddtimewg=10LpNag+Kg+Clg+HCO3g+Xgwg-Nao+Ko+Clo+HCO3o

Component: wo

wo=P_wg+P_wo-wg

Component: N_Nao

ddtimeN_Nao=0.013J_NaKATPase_n+J_NaT+J_NaP+J_leakNa-100dN_Nag_dt

Component: N_Ko

ddtimeN_Ko=-0.01J_NaT+J_NaP+J_leakNa+3J_NaKATPase_n-dN_Kg_dt

Component: N_HCO3o

ddtimeN_HCO3o=-0.012J_NBC

Component: electric_potentials

ENa_n=PsilnNaoNan EK_n=PsilnKoKn ENa_g=PsilnNaoNag EK_g=PsilnKoKg ECl_g=-1PsilnCloClg ENBC_g=-PsilnNaoNagHCO3oHCO3g2

Component: ion_concentrations

P_Clo=P_Nao+P_Ko-P_HCO3o P_HCO3g=P_HCO3oP_NaoP_NagP_Vm_gPsi P_Clg=P_CloP_Vm_gPsi Kg=N_Kgwg Ko=N_Kowo Kn=P_wnP_Kn+P_wgP_Kg+P_woP_Ko-N_Ko+N_KgP_wn Nag=N_Nagwg Nao=N_Naowo Nan=P_wnP_Nan+P_wgP_Nag+P_woP_Nao-N_Nao+N_NagP_wn Clg=Nag+Kg-HCO3g+rhoXgwg Cln=P_Cln+Nan-P_Nan+Kn-P_Kn Clo=P_wgP_Clg+P_woP_Clo+P_wnP_Cln-wgClg+P_wnClnwo HCO3o=N_HCO3owo HCO3g=P_wgP_HCO3g+P_woP_HCO3o-N_HCO3owg

Component: model_parameters

Psi=1000RTF