From da3a335d297c3e2c80a2b2f52a964c9ea404f0f5 Mon Sep 17 00:00:00 2001 From: Jonathan Karr Date: Thu, 17 Feb 2022 01:23:04 -0500 Subject: [PATCH 1/2] Correct y plot variable to `vm` from `v`, which isn't defined --- muscleAP.ode | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/muscleAP.ode b/muscleAP.ode index 548f7ed..f9c907f 100644 --- a/muscleAP.ode +++ b/muscleAP.ode @@ -45,6 +45,6 @@ dVt/dt=((Vm - Vt) / (Rs * Ct)) % Numerical and plotting parameters for xpp @ meth=Euler, dt=0.01, total=25 -@ yp=v, yhi=50, ylo=-110, xlo=0, xhi=25, bounds=5000 +@ yp=vm, yhi=50, ylo=-110, xlo=0, xhi=25, bounds=5000 done From 0f1443f5bf75d60678594538a9ae7dd59d4465df Mon Sep 17 00:00:00 2001 From: Jonathan Karr Date: Thu, 17 Feb 2022 01:25:17 -0500 Subject: [PATCH 2/2] Reduce length of parameter and variable names to pass validation by XPP 8.0 - Replaced `alpha_` with `a_` in parameter names - Replaced `beta_` with `b_` in parameter names --- muscleAP.ode | 26 +++++++++++++------------- 1 file changed, 13 insertions(+), 13 deletions(-) diff --git a/muscleAP.ode b/muscleAP.ode index f9c907f..46763d7 100644 --- a/muscleAP.ode +++ b/muscleAP.ode @@ -13,28 +13,28 @@ iStim= iStim_mag * heav(mod(t,period)-iStim_beg) * heav(iStim_beg+iStim_dur-mod % Values of the model parameters Param gNa_max=1.0, gK_max=0.415, gL_max=0.0024 Param ENa=50.0, EK=-70.0, EL=-95.0, Cm=0.0090, Ct=0.04, Rs=15.0 -Param alpha_m_max=0.208 -Param beta_m_max=2.081 -Param alpha_n_max=0.0229 -Param beta_n_max=0.09616 -Param alpha_h_max=0.0156 -Param beta_h_max=3.382 +Param a_m_max=0.208 +Param b_m_max=2.081 +Param a_n_max=0.0229 +Param b_n_max=0.09616 +Param a_h_max=0.0156 +Param b_h_max=3.382 Param Em=-42.0, En=-40.0, Eh=-41.0 -Param v_alpha_m=10.0, v_alpha_h=14.7, v_alpha_n=7.0 +Param v_alph_m=10.0, v_alph_h=14.7, v_alph_n=7.0 Param v_beta_m=18.0, v_beta_h=7.6, v_beta_n=40.0 % Expressions Ina= (gNa_max * m * m * m * h * (Vm - ENa)) -beta_n= (beta_n_max * exp(((En - Vm) / v_beta_n))) -beta_m= (beta_m_max * exp(((Em - Vm) / v_beta_m))) -beta_h= (beta_h_max / (1.0 + exp(((Eh - Vm) / v_beta_h)))) +beta_n= (b_n_max * exp(((En - Vm) / v_beta_n))) +beta_m= (b_m_max * exp(((Em - Vm) / v_beta_m))) +beta_h= (b_h_max / (1.0 + exp(((Eh - Vm) / v_beta_h)))) IT= ((Vm - Vt) / Rs) IL= (gL_max * (Vm - EL)) IK= (gK_max * n * n * n * n * (Vm - EK)) -alpha_n= (alpha_n_max * (Vm - En) / (1.0 - exp(((En - Vm) / v_alpha_n)))) -alpha_m= (alpha_m_max * (Vm - Em) / (1.0 - exp(((Em - Vm) / v_alpha_m)))) -alpha_h= (alpha_h_max * exp(((Eh - Vm) / v_alpha_h))) +alpha_n= (a_n_max * (Vm - En) / (1.0 - exp(((En - Vm) / v_alph_n)))) +alpha_m= (a_m_max * (Vm - Em) / (1.0 - exp(((Em - Vm) / v_alph_m)))) +alpha_h= (a_h_max * exp(((Eh - Vm) / v_alph_h))) % Differential equations dVm/dt= ((iStim - (INa + IK + IL + IT)) / Cm)