import nsrunit; unit conversion on; // unit millivolt predefined unit per_millivolt=1E3 kilogram^(-1)*meter^(-2)*second^3*ampere^1; unit first_order_rate_constant_per_millivolt=1E3 kilogram^(-1)*meter^(-2)*second^2*ampere^1; // unit micromolar predefined unit per_micromolar=1E3 meter^3*mole^(-1); unit per_micromolar2=1E6 meter^6*mole^(-2); // unit picofarad predefined // unit picosiemens predefined // unit picoampere predefined // unit femtoampere predefined unit micromolar_per_coulomb=1E-3 meter^(-3)*second^(-1)*ampere^(-1)*mole^1; unit femtocoulomb_per_zeptomole=1E6 second^1*ampere^1*mole^(-1); unit attojoule_per_zeptomole_kelvin=1E3 kilogram^1*meter^2*second^(-2)*kelvin^(-1)*mole^(-1); unit first_order_rate_constant=1 second^(-1); unit first_order_rate_constant_squared=1 second^(-2); unit flux=1E-3 meter^(-3)*second^(-1)*mole^1; // unit millivolt predefined unit per_millivolt=1E3 kilogram^(-1)*meter^(-2)*second^3*ampere^1; unit per_second_per_millivolt=1E3 kilogram^(-1)*meter^(-2)*second^2*ampere^1; // unit micromolar predefined unit per_micromolar=1E3 meter^3*mole^(-1); unit per_micromolar2=1E6 meter^6*mole^(-2); // unit picofarad predefined // unit picosiemens predefined // unit picoampere predefined // unit femtoampere predefined unit micromolar_per_coulomb=1E-3 meter^(-3)*second^(-1)*ampere^(-1)*mole^1; unit femtocoulomb_per_zeptomole=1E6 second^1*ampere^1*mole^(-1); unit attojoule_per_zeptomole_kelvin=1E3 kilogram^1*meter^2*second^(-2)*kelvin^(-1)*mole^(-1); unit first_order_rate_constant=1 second^(-1); unit per_second_squared=1 second^(-2); unit flux=1E-3 meter^(-3)*second^(-1)*mole^1; math main { realDomain time second; time.min=0; extern time.max; extern time.delta; real I_CaT(time) picoampere; real g_CaT picosiemens; g_CaT=700.0; real ECaT millivolt; ECaT=17.0; real alpha_OT first_order_rate_constant; alpha_OT=800.0; real beta_OT first_order_rate_constant; beta_OT=72; real OCaT(time) dimensionless; when(time=time.min) OCaT=0; real Vm(time) millivolt; when(time=time.min) Vm=-70.0; real dT(time) dimensionless; when(time=time.min) dT=0.010; real fT(time) dimensionless; when(time=time.min) fT=0.001; real infinity_dT(time) dimensionless; real tau_dT1 second; real A_dT1 second; A_dT1=0.0025; real k_dT per_millivolt; k_dT=-0.60; real V_dT millivolt; V_dT=-53; real infinity_fT(time) dimensionless; real tau_fT1(time) second; real A_fT1 second; A_fT1=0.019; real A_fT2 second; A_fT2=6.75; real A_fT3 per_millivolt; A_fT3=2; real A_fT4 millivolt; A_fT4=-40; real k_fT per_millivolt; k_fT=1; real V_fT millivolt; V_fT=-65; real I_CaExt(time) picoampere; real gCaExt femtoampere; gCaExt=100.0; real KCaExt micromolar; KCaExt=1.0; real CCy(time) micromolar; when(time=time.min) CCy=0.12; real I_Kv1_1(time) picoampere; real g_Kv1_1 picosiemens; g_Kv1_1=10.0; real EK millivolt; real dv1_1(time) dimensionless; when(time=time.min) dv1_1=0.000; real fv1_1(time) dimensionless; when(time=time.min) fv1_1=1.0; real alpha_dv1_1(time) first_order_rate_constant; real A_dv1_11(time) first_order_rate_constant; real A_dv1_11a first_order_rate_constant; A_dv1_11a=0.0; real A_dv1_11b first_order_rate_constant; A_dv1_11b=250.0; real A_dv1_11c per_millivolt; A_dv1_11c=-0.077; real A_dv1_11d millivolt; A_dv1_11d=-4.8; real A_dv1_12 dimensionless; A_dv1_12=0.80; real A_dv1_13 per_millivolt; A_dv1_13=-0.13; real A_dv1_14 millivolt; A_dv1_14=-25.0; real beta_dv1_1(time) first_order_rate_constant; real alpha_fv1_1(time) first_order_rate_constant; real A_fv1_11 first_order_rate_constant; A_fv1_11=333.0; real A_fv1_12 dimensionless; A_fv1_12=0.10; real A_fv1_13 per_millivolt; A_fv1_13=0.23; real A_fv1_14 millivolt; A_fv1_14=-44.8; real beta_fv1_1(time) first_order_rate_constant; real I_KERG(time) picoampere; real g_KERG picosiemens; g_KERG=10.0; real dERG(time) dimensionless; when(time=time.min) dERG=0.000; real alpha_dERG(time) first_order_rate_constant; real A_dERG1 first_order_rate_constant; A_dERG1=1000.0; real A_dERG2 dimensionless; A_dERG2=0.70; real A_dERG3 per_millivolt; A_dERG3=-0.56; real A_dERG4 millivolt; A_dERG4=30; real beta_dERG(time) first_order_rate_constant; real I_KB(time) picoampere; real g_KB picosiemens; g_KB=6.0; real I_Na.I_Na(time) picoampere; real g_Na picosiemens; g_Na=2000.0; real ENa_Cy millivolt; real dNa(time) dimensionless; when(time=time.min) dNa=0.000; real fNa(time) dimensionless; when(time=time.min) fNa=0.100; real alpha_dNa(time) first_order_rate_constant; real A_dNa1(time) first_order_rate_constant; real A_dNa1a first_order_rate_constant; A_dNa1a=250.0; real A_dNa1b first_order_rate_constant; A_dNa1b=1000.0; real A_dNa1c per_millivolt; A_dNa1c=-0.2; real A_dNa1d millivolt; A_dNa1d=-15.0; real A_dNa2 dimensionless; A_dNa2=1.0; real A_dNa3 per_millivolt; A_dNa3=-0.147; real A_dNa4 millivolt; A_dNa4=-35.3; real beta_dNa(time) first_order_rate_constant; real alpha_fNa(time) first_order_rate_constant; real A_fNa1(time) first_order_rate_constant; real A_fNa1a first_order_rate_constant; A_fNa1a=10.0; real A_fNa1b first_order_rate_constant; A_fNa1b=50.0; real A_fNa1c per_millivolt; A_fNa1c=-0.036; real A_fNa1d millivolt; A_fNa1d=-38.0; real A_fNa2 dimensionless; A_fNa2=1.00; real A_fNa3 per_millivolt; A_fNa3=0.167; real A_fNa4 millivolt; A_fNa4=-83.0; real beta_fNa(time) first_order_rate_constant; real I_L(time) picoampere; real g_L picosiemens; g_L=0.20; real EL millivolt; EL=0.0; real I_ion_Cy(time) picoampere; real Cm picofarad; Cm=20.0; real I_ion_PU1(time) picoampere; real I_ion_PU2(time) picoampere; real I_ion_PU3(time) picoampere; real I_ion_PU4(time) picoampere; real I_ion_PU5(time) picoampere; real I_ion_PU6(time) picoampere; real I_ion_PU7(time) picoampere; real I_ion_PU8(time) picoampere; real I_ion_PU9(time) picoampere; real I_ion_PU10(time) picoampere; real ECa_Cy(time) millivolt; real Nerst_potentials.T kelvin; Nerst_potentials.T=310.16; real Nerst_potentials.R attojoule_per_zeptomole_kelvin; Nerst_potentials.R=8.314E-3; real Nerst_potentials.F femtocoulomb_per_zeptomole; Nerst_potentials.F=0.09649; real Nerst_potentials.CO micromolar; Nerst_potentials.CO=1.8E3; real Nerst_potentials.KO micromolar; Nerst_potentials.KO=5400.0; real Nerst_potentials.Ki micromolar; Nerst_potentials.Ki=145E3; real Nerst_potentials.NO micromolar; Nerst_potentials.NO=140E3; real Ni micromolar; Ni=10E3; real lambda_S2_Cy dimensionless; lambda_S2_Cy=1.0E-3; real delta_SCy micromolar_per_coulomb; delta_SCy=0.002; real JCy(time) flux; real was_JS2Cy1(time) flux; real was_JS2Cy2(time) flux; real was_JS2Cy3(time) flux; real was_JS2Cy4(time) flux; real was_JS2Cy5(time) flux; real was_JS2Cy6(time) flux; real was_JS2Cy7(time) flux; real was_JS2Cy8(time) flux; real was_JS2Cy9(time) flux; real was_JS2Cy10(time) flux; real mu_S2Cy1 first_order_rate_constant; mu_S2Cy1=0.000; real mu_S2Cy2 first_order_rate_constant; mu_S2Cy2=0.3361; real mu_S2Cy3 first_order_rate_constant; mu_S2Cy3=0.3658; real mu_S2Cy4 first_order_rate_constant; mu_S2Cy4=0.4357; real mu_S2Cy5 first_order_rate_constant; mu_S2Cy5=0.5351; real mu_S2Cy6 first_order_rate_constant; mu_S2Cy6=0.6366; real mu_S2Cy7 first_order_rate_constant; mu_S2Cy7=0.7386; real mu_S2Cy8 first_order_rate_constant; mu_S2Cy8=0.8361; real mu_S2Cy9 first_order_rate_constant; mu_S2Cy9=0.9320; real mu_S2Cy10 first_order_rate_constant; mu_S2Cy10=1.0185; real CS21(time) micromolar; when(time=time.min) CS21=0.023; real CS22(time) micromolar; when(time=time.min) CS22=0.023; real CS23(time) micromolar; when(time=time.min) CS23=0.023; real CS24(time) micromolar; when(time=time.min) CS24=0.023; real CS25(time) micromolar; when(time=time.min) CS25=0.023; real CS26(time) micromolar; when(time=time.min) CS26=0.023; real CS27(time) micromolar; when(time=time.min) CS27=0.023; real CS28(time) micromolar; when(time=time.min) CS28=0.023; real CS29(time) micromolar; when(time=time.min) CS29=0.023; real CS210(time) micromolar; when(time=time.min) CS210=0.023; real was_JS1S2_1(time) flux; real was_JS1S2_2(time) flux; real was_JS1S2_3(time) flux; real was_JS1S2_4(time) flux; real was_JS1S2_5(time) flux; real was_JS1S2_6(time) flux; real was_JS1S2_7(time) flux; real was_JS1S2_8(time) flux; real was_JS1S2_9(time) flux; real was_JS1S2_10(time) flux; real mu_S1S2_1 first_order_rate_constant; real mu_S1S2_2 first_order_rate_constant; real mu_S1S2_3 first_order_rate_constant; real mu_S1S2_4 first_order_rate_constant; real mu_S1S2_5 first_order_rate_constant; real mu_S1S2_6 first_order_rate_constant; real mu_S1S2_7 first_order_rate_constant; real mu_S1S2_8 first_order_rate_constant; real mu_S1S2_9 first_order_rate_constant; real mu_S1S2_10 first_order_rate_constant; real mu_A first_order_rate_constant; mu_A=0.30; real mu_B first_order_rate_constant; mu_B=0.24; real n_PU dimensionless; n_PU=10.0; real CS11(time) micromolar; when(time=time.min) CS11=0.120; real CS12(time) micromolar; when(time=time.min) CS12=0.120; real CS13(time) micromolar; when(time=time.min) CS13=0.120; real CS14(time) micromolar; when(time=time.min) CS14=0.120; real CS15(time) micromolar; when(time=time.min) CS15=0.120; real CS16(time) micromolar; when(time=time.min) CS16=0.120; real CS17(time) micromolar; when(time=time.min) CS17=0.120; real CS18(time) micromolar; when(time=time.min) CS18=0.120; real CS19(time) micromolar; when(time=time.min) CS19=0.120; real CS110(time) micromolar; when(time=time.min) CS110=0.120; real mu_Cy first_order_rate_constant; mu_Cy=1.30; real C_infinity micromolar; C_infinity=0.12; real model_parameters.P micromolar; model_parameters.P=1.0; real PU_membrane.I_iCa(time) picoampere; real PU_membrane.I_iNa(time) picoampere; real PU_membrane.alpha_scale dimensionless; real PU_membrane.I_Ca(time) picoampere; real PU_membrane.I_Na(time) picoampere; real PU_membrane.I_NaP(time) picoampere; real PU_membrane.I_NSCC_Ca(time) picoampere; real PU_membrane.I_PM(time) picoampere; real PU_membrane.I_NSCC_Na(time) picoampere; real I_Ca.gCa_ picosiemens; I_Ca.gCa_=0.074; real I_Ca.kCa per_millivolt; I_Ca.kCa=0.013; real I_Ca.kVCa per_millivolt; I_Ca.kVCa=-0.20; real I_Ca.VhCa millivolt; I_Ca.VhCa=-85.0; real I_Ca.gCa(time) picosiemens; real I_Ca.ECa_PU(time) millivolt; real I_Ca.T kelvin; I_Ca.T=310.16; real I_Ca.R attojoule_per_zeptomole_kelvin; I_Ca.R=8.314E-3; real I_Ca.F femtocoulomb_per_zeptomole; I_Ca.F=0.09649; real I_Ca.CO micromolar; I_Ca.CO=1.8E3; real I_Na_1.gNa picosiemens; I_Na_1.gNa=13.5; real I_Na_1.ENa_PU(time) millivolt; real I_Na_1.NO micromolar; I_Na_1.NO=140.0E3; real I_Na_1.NS1(time) micromolar; when(time=time.min) I_Na_1.NS1=1.01E4; real I_NSCC_Ca.gNSCC_Ca_ picosiemens; I_NSCC_Ca.gNSCC_Ca_=0.10; real I_NSCC_Ca.gNSCC_Ca(time) picosiemens; // Var below replaced by constant in model eqns to satisfy unit correction // real I_NSCC_Ca.hNSCC dimensionless; // I_NSCC_Ca.hNSCC=3.0; real I_NSCC_Ca.ENSCC millivolt; I_NSCC_Ca.ENSCC=0.0; real I_NSCC_Ca.KNSCC micromolar; I_NSCC_Ca.KNSCC=0.12; real I_NSCC_Na.gNSCC_Na_ picosiemens; I_NSCC_Na.gNSCC_Na_=160.0; real I_NSCC_Na.gNSCC_Na(time) picosiemens; real I_PM.gPM femtoampere; I_PM.gPM=675.0; real I_PM.KPM micromolar; I_PM.KPM=1.0; real I_NaP.gNaP picosiemens; I_NaP.gNaP=187.5; real I_NaP.KNaP micromolar; I_NaP.KNaP=1.0E4; // Var below replaced by constant in model eqns to satisfy unit correction // real I_NaP.hNaP dimensionless; // I_NaP.hNaP=4.0; real I_NaP.ENaP millivolt; I_NaP.ENaP=10.0; real was_JSERCA.JSERCA(time) flux; real was_JSERCA.VSERCA first_order_rate_constant; was_JSERCA.VSERCA=1.0E5; real was_JSERCA.A2 dimensionless; was_JSERCA.A2=6E-4; real was_JSERCA.A4 per_micromolar; was_JSERCA.A4=3.57; real was_JSERCA.A5 per_micromolar; was_JSERCA.A5=2.7E-5; real was_JSERCA.A6 per_micromolar2; was_JSERCA.A6=2.31E-5; real was_JSERCA.CER(time) micromolar; when(time=time.min) was_JSERCA.CER=200.0; real JMCU.JMCU(time) flux; real JMCU.VMCU flux; JMCU.VMCU=800.0; real JMCU.KMCU micromolar; JMCU.KMCU=10.0; real JMCU.epsilon_INH(time) dimensionless; real JMCU.KINH micromolar; JMCU.KINH=10.0; // Var below replaced by constant in model eqns to satisfy unit correction // real JMCU.hINH dimensionless; // JMCU.hINH=4.0; real JMCU.CMT(time) micromolar; when(time=time.min) JMCU.CMT=0.200; real JIPR.JIPR(time) flux; real JIPR.kIPR first_order_rate_constant; JIPR.kIPR=2000.0; real JIPR.k_1 flux; JIPR.k_1=6.4; real JIPR.k1 first_order_rate_constant; JIPR.k1=0.0; real JIPR.k2 first_order_rate_constant; JIPR.k2=4.0; real JIPR.r2 first_order_rate_constant; JIPR.r2=250.0; real JIPR.r_2 flux; JIPR.r_2=0.0; real JIPR.r4 first_order_rate_constant; JIPR.r4=750.0; real JIPR.R1 micromolar; JIPR.R1=36.0; real JIPR.R3 micromolar; JIPR.R3=300.0; real JIPR.phi1(time) first_order_rate_constant; real JIPR.phi_1(time) flux; real JIPR.phi2(time) first_order_rate_constant; real JIPR.phi3(time) first_order_rate_constant; real JIPR.H(time) dimensionless; when(time=time.min) JIPR.H=0.200; real JIPR.g_beta first_order_rate_constant; JIPR.g_beta=1500.0; real JIPR.h_beta dimensionless; JIPR.h_beta=4.0; real JIPR.g_alpha first_order_rate_constant; JIPR.g_alpha=0.85; real JIPR.K_beta micromolar; JIPR.K_beta=0.35; real JIPR.alpha_zeta first_order_rate_constant; real JIPR.beta_zeta(time) first_order_rate_constant; real JIPR.zeta(time) dimensionless; when(time=time.min) JIPR.zeta=0.300; real JIPR.K_phi3_act micromolar; JIPR.K_phi3_act=0.1; real JIPR.K_phi3_inh micromolar; JIPR.K_phi3_inh=0.5; real JIPR.h_phi3_act dimensionless; JIPR.h_phi3_act=3.0; real JIPR.h_phi3_inh dimensionless; JIPR.h_phi3_inh=3.0; real JIPR.g_phi3 first_order_rate_constant; JIPR.g_phi3=4.5; real JIPR.P micromolar; JIPR.P=1.0; real JNCX.JNCX(time) flux; real JNCX.VNCX flux; JNCX.VNCX=3.0; real JNCX.KNCX micromolar; JNCX.KNCX=0.3; real CS1.lambda_MT_S1 dimensionless; real CS1.lambda_ER_S1 dimensionless; real CS1.delta_SPU micromolar_per_coulomb; CS1.delta_SPU=18.5; real CS1.gamma_S1 dimensionless; CS1.gamma_S1=100.0; real CS1.gamma_MT dimensionless; CS1.gamma_MT=200.0; real CS1.gamma_ER dimensionless; CS1.gamma_ER=20.0; real CS2.lambda_MT_S2 dimensionless; real CS2.lambda_ER_S2 dimensionless; real CS2.lambda_S1_S2 dimensionless; real CS2.gamma_S2 dimensionless; CS2.gamma_S2=1.0; real CMT.fm(time) dimensionless; real CMT.Km micromolar; CMT.Km=0.01; real CMT.Bm micromolar; CMT.Bm=100.0; real model_parameters_1.KO micromolar; model_parameters_1.KO=5400.0; real model_parameters_1.Ki micromolar; model_parameters_1.Ki=145E3; real model_parameters_1.n_PU_base dimensionless; model_parameters_1.n_PU_base=50.0; real PU_membrane_1.I_iCa(time) picoampere; real PU_membrane_1.I_iNa(time) picoampere; real PU_membrane_1.alpha_scale dimensionless; real PU_membrane_1.I_Ca(time) picoampere; real PU_membrane_1.I_Na(time) picoampere; real PU_membrane_1.I_NaP(time) picoampere; real PU_membrane_1.I_NSCC_Ca(time) picoampere; real PU_membrane_1.I_PM(time) picoampere; real PU_membrane_1.I_NSCC_Na(time) picoampere; real I_Ca_1.gCa_ picosiemens; I_Ca_1.gCa_=0.074; real I_Ca_1.kCa per_millivolt; I_Ca_1.kCa=0.013; real I_Ca_1.kVCa per_millivolt; I_Ca_1.kVCa=-0.20; real I_Ca_1.VhCa millivolt; I_Ca_1.VhCa=-85.0; real I_Ca_1.gCa(time) picosiemens; real I_Ca_1.ECa_PU(time) millivolt; real I_Ca_1.T kelvin; I_Ca_1.T=310.16; real I_Ca_1.R attojoule_per_zeptomole_kelvin; I_Ca_1.R=8.314E-3; real I_Ca_1.F femtocoulomb_per_zeptomole; I_Ca_1.F=0.09649; real I_Ca_1.CO micromolar; I_Ca_1.CO=1.8E3; real I_Na_2.gNa picosiemens; I_Na_2.gNa=13.5; real I_Na_2.ENa_PU(time) millivolt; real I_Na_2.NO micromolar; I_Na_2.NO=140.0E3; real I_Na_2.NS1(time) micromolar; when(time=time.min) I_Na_2.NS1=1.01E4; real I_NSCC_Ca_1.gNSCC_Ca_ picosiemens; I_NSCC_Ca_1.gNSCC_Ca_=0.10; real I_NSCC_Ca_1.gNSCC_Ca(time) picosiemens; // Var below replaced by constant in model eqns to satisfy unit correction // real I_NSCC_Ca_1.hNSCC dimensionless; // I_NSCC_Ca_1.hNSCC=3.0; real I_NSCC_Ca_1.ENSCC millivolt; I_NSCC_Ca_1.ENSCC=0.0; real I_NSCC_Ca_1.KNSCC micromolar; I_NSCC_Ca_1.KNSCC=0.12; real I_NSCC_Na_1.gNSCC_Na_ picosiemens; I_NSCC_Na_1.gNSCC_Na_=160.0; real I_NSCC_Na_1.gNSCC_Na(time) picosiemens; real I_PM_1.gPM femtoampere; I_PM_1.gPM=675.0; real I_PM_1.KPM micromolar; I_PM_1.KPM=1.0; real I_NaP_1.gNaP picosiemens; I_NaP_1.gNaP=187.5; real I_NaP_1.KNaP micromolar; I_NaP_1.KNaP=1.0E4; // Var below replaced by constant in model eqns to satisfy unit correction // real I_NaP_1.hNaP dimensionless; // I_NaP_1.hNaP=4.0; real I_NaP_1.ENaP millivolt; I_NaP_1.ENaP=10.0; real was_JSERCA_1.JSERCA(time) flux; real was_JSERCA_1.VSERCA first_order_rate_constant; was_JSERCA_1.VSERCA=1.0E5; real was_JSERCA_1.A2 dimensionless; was_JSERCA_1.A2=6E-4; real was_JSERCA_1.A4 per_micromolar; was_JSERCA_1.A4=3.57; real was_JSERCA_1.A5 per_micromolar; was_JSERCA_1.A5=2.7E-5; real was_JSERCA_1.A6 per_micromolar2; was_JSERCA_1.A6=2.31E-5; real was_JSERCA_1.CER(time) micromolar; when(time=time.min) was_JSERCA_1.CER=200.0; real JMCU_1.JMCU(time) flux; real JMCU_1.VMCU flux; JMCU_1.VMCU=800.0; real JMCU_1.KMCU micromolar; JMCU_1.KMCU=10.0; real JMCU_1.epsilon_INH(time) dimensionless; real JMCU_1.KINH micromolar; JMCU_1.KINH=10.0; // Var below replaced by constant in model eqns to satisfy unit correction // real JMCU_1.hINH dimensionless; // JMCU_1.hINH=4.0; real JMCU_1.CMT(time) micromolar; when(time=time.min) JMCU_1.CMT=0.200; real JIPR_1.JIPR(time) flux; real JIPR_1.kIPR first_order_rate_constant; JIPR_1.kIPR=2000.0; real JIPR_1.k_1 flux; JIPR_1.k_1=6.4; real JIPR_1.k1 first_order_rate_constant; JIPR_1.k1=0.0; real JIPR_1.k2 first_order_rate_constant; JIPR_1.k2=4.0; real JIPR_1.r2 first_order_rate_constant; JIPR_1.r2=250.0; real JIPR_1.r_2 flux; JIPR_1.r_2=0.0; real JIPR_1.r4 first_order_rate_constant; JIPR_1.r4=750.0; real JIPR_1.R1 micromolar; JIPR_1.R1=36.0; real JIPR_1.R3 micromolar; JIPR_1.R3=300.0; real JIPR_1.phi1(time) first_order_rate_constant; real JIPR_1.phi_1(time) flux; real JIPR_1.phi2(time) first_order_rate_constant; real JIPR_1.phi3(time) first_order_rate_constant; real JIPR_1.H(time) dimensionless; when(time=time.min) JIPR_1.H=0.200; real JIPR_1.g_beta first_order_rate_constant; JIPR_1.g_beta=1500.0; real JIPR_1.h_beta dimensionless; JIPR_1.h_beta=4.0; real JIPR_1.g_alpha first_order_rate_constant; JIPR_1.g_alpha=0.85; real JIPR_1.K_beta micromolar; JIPR_1.K_beta=0.35; real JIPR_1.alpha_zeta first_order_rate_constant; real JIPR_1.beta_zeta(time) first_order_rate_constant; real JIPR_1.zeta(time) dimensionless; when(time=time.min) JIPR_1.zeta=0.300; real JIPR_1.K_phi3_act micromolar; JIPR_1.K_phi3_act=0.1; real JIPR_1.K_phi3_inh micromolar; JIPR_1.K_phi3_inh=0.5; real JIPR_1.h_phi3_act dimensionless; JIPR_1.h_phi3_act=3.0; real JIPR_1.h_phi3_inh dimensionless; JIPR_1.h_phi3_inh=3.0; real JIPR_1.g_phi3 first_order_rate_constant; JIPR_1.g_phi3=4.5; real JIPR_1.P micromolar; JIPR_1.P=1.0; real JNCX_1.JNCX(time) flux; real JNCX_1.VNCX flux; JNCX_1.VNCX=3.0; real JNCX_1.KNCX micromolar; JNCX_1.KNCX=0.3; real CS1_1.lambda_MT_S1 dimensionless; real CS1_1.lambda_ER_S1 dimensionless; real CS1_1.delta_SPU micromolar_per_coulomb; CS1_1.delta_SPU=18.5; real CS1_1.gamma_S1 dimensionless; CS1_1.gamma_S1=100.0; real CS1_1.gamma_MT dimensionless; CS1_1.gamma_MT=200.0; real CS1_1.gamma_ER dimensionless; CS1_1.gamma_ER=20.0; real CS2_1.lambda_MT_S2 dimensionless; real CS2_1.lambda_ER_S2 dimensionless; real CS2_1.lambda_S1_S2 dimensionless; real CS2_1.gamma_S2 dimensionless; CS2_1.gamma_S2=1.0; real CMT_1.fm(time) dimensionless; real CMT_1.Km micromolar; CMT_1.Km=0.01; real CMT_1.Bm micromolar; CMT_1.Bm=100.0; real model_parameters_2.KO micromolar; model_parameters_2.KO=5400.0; real model_parameters_2.Ki micromolar; model_parameters_2.Ki=145E3; real model_parameters_2.n_PU_base dimensionless; model_parameters_2.n_PU_base=50.0; real PU_membrane_2.I_iCa(time) picoampere; real PU_membrane_2.I_iNa(time) picoampere; real PU_membrane_2.alpha_scale dimensionless; real PU_membrane_2.I_Ca(time) picoampere; real PU_membrane_2.I_Na(time) picoampere; real PU_membrane_2.I_NaP(time) picoampere; real PU_membrane_2.I_NSCC_Ca(time) picoampere; real PU_membrane_2.I_PM(time) picoampere; real PU_membrane_2.I_NSCC_Na(time) picoampere; real I_Ca_2.gCa_ picosiemens; I_Ca_2.gCa_=0.074; real I_Ca_2.kCa per_millivolt; I_Ca_2.kCa=0.013; real I_Ca_2.kVCa per_millivolt; I_Ca_2.kVCa=-0.20; real I_Ca_2.VhCa millivolt; I_Ca_2.VhCa=-85.0; real I_Ca_2.gCa(time) picosiemens; real I_Ca_2.ECa_PU(time) millivolt; real I_Ca_2.T kelvin; I_Ca_2.T=310.16; real I_Ca_2.R attojoule_per_zeptomole_kelvin; I_Ca_2.R=8.314E-3; real I_Ca_2.F femtocoulomb_per_zeptomole; I_Ca_2.F=0.09649; real I_Ca_2.CO micromolar; I_Ca_2.CO=1.8E3; real I_Na_3.gNa picosiemens; I_Na_3.gNa=13.5; real I_Na_3.ENa_PU(time) millivolt; real I_Na_3.NO micromolar; I_Na_3.NO=140.0E3; real I_Na_3.NS1(time) micromolar; when(time=time.min) I_Na_3.NS1=1.01E4; real I_NSCC_Ca_2.gNSCC_Ca_ picosiemens; I_NSCC_Ca_2.gNSCC_Ca_=0.10; real I_NSCC_Ca_2.gNSCC_Ca(time) picosiemens; // Var below replaced by constant in model eqns to satisfy unit correction // real I_NSCC_Ca_2.hNSCC dimensionless; // I_NSCC_Ca_2.hNSCC=3.0; real I_NSCC_Ca_2.ENSCC millivolt; I_NSCC_Ca_2.ENSCC=0.0; real I_NSCC_Ca_2.KNSCC micromolar; I_NSCC_Ca_2.KNSCC=0.12; real I_NSCC_Na_2.gNSCC_Na_ picosiemens; I_NSCC_Na_2.gNSCC_Na_=160.0; real I_NSCC_Na_2.gNSCC_Na(time) picosiemens; real I_PM_2.gPM femtoampere; I_PM_2.gPM=675.0; real I_PM_2.KPM micromolar; I_PM_2.KPM=1.0; real I_NaP_2.gNaP picosiemens; I_NaP_2.gNaP=187.5; real I_NaP_2.KNaP micromolar; I_NaP_2.KNaP=1.0E4; // Var below replaced by constant in model eqns to satisfy unit correction // real I_NaP_2.hNaP dimensionless; // I_NaP_2.hNaP=4.0; real I_NaP_2.ENaP millivolt; I_NaP_2.ENaP=10.0; real was_JSERCA_2.JSERCA(time) flux; real was_JSERCA_2.VSERCA first_order_rate_constant; was_JSERCA_2.VSERCA=1.0E5; real was_JSERCA_2.A2 dimensionless; was_JSERCA_2.A2=6E-4; real was_JSERCA_2.A4 per_micromolar; was_JSERCA_2.A4=3.57; real was_JSERCA_2.A5 per_micromolar; was_JSERCA_2.A5=2.7E-5; real was_JSERCA_2.A6 per_micromolar2; was_JSERCA_2.A6=2.31E-5; real was_JSERCA_2.CER(time) micromolar; when(time=time.min) was_JSERCA_2.CER=200.0; real JMCU_2.JMCU(time) flux; real JMCU_2.VMCU flux; JMCU_2.VMCU=800.0; real JMCU_2.KMCU micromolar; JMCU_2.KMCU=10.0; real JMCU_2.epsilon_INH(time) dimensionless; real JMCU_2.KINH micromolar; JMCU_2.KINH=10.0; // Var below replaced by constant in model eqns to satisfy unit correction // real JMCU_2.hINH dimensionless; // JMCU_2.hINH=4.0; real JMCU_2.CMT(time) micromolar; when(time=time.min) JMCU_2.CMT=0.200; real JIPR_2.JIPR(time) flux; real JIPR_2.kIPR first_order_rate_constant; JIPR_2.kIPR=2000.0; real JIPR_2.k_1 flux; JIPR_2.k_1=6.4; real JIPR_2.k1 first_order_rate_constant; JIPR_2.k1=0.0; real JIPR_2.k2 first_order_rate_constant; JIPR_2.k2=4.0; real JIPR_2.r2 first_order_rate_constant; JIPR_2.r2=250.0; real JIPR_2.r_2 flux; JIPR_2.r_2=0.0; real JIPR_2.r4 first_order_rate_constant; JIPR_2.r4=750.0; real JIPR_2.R1 micromolar; JIPR_2.R1=36.0; real JIPR_2.R3 micromolar; JIPR_2.R3=300.0; real JIPR_2.phi1(time) first_order_rate_constant; real JIPR_2.phi_1(time) flux; real JIPR_2.phi2(time) first_order_rate_constant; real JIPR_2.phi3(time) first_order_rate_constant; real JIPR_2.H(time) dimensionless; when(time=time.min) JIPR_2.H=0.200; real JIPR_2.g_beta first_order_rate_constant; JIPR_2.g_beta=1500.0; real JIPR_2.h_beta dimensionless; JIPR_2.h_beta=4.0; real JIPR_2.g_alpha first_order_rate_constant; JIPR_2.g_alpha=0.85; real JIPR_2.K_beta micromolar; JIPR_2.K_beta=0.35; real JIPR_2.alpha_zeta first_order_rate_constant; real JIPR_2.beta_zeta(time) first_order_rate_constant; real JIPR_2.zeta(time) dimensionless; when(time=time.min) JIPR_2.zeta=0.300; real JIPR_2.K_phi3_act micromolar; JIPR_2.K_phi3_act=0.1; real JIPR_2.K_phi3_inh micromolar; JIPR_2.K_phi3_inh=0.5; real JIPR_2.h_phi3_act dimensionless; JIPR_2.h_phi3_act=3.0; real JIPR_2.h_phi3_inh dimensionless; JIPR_2.h_phi3_inh=3.0; real JIPR_2.g_phi3 first_order_rate_constant; JIPR_2.g_phi3=4.5; real JIPR_2.P micromolar; JIPR_2.P=1.0; real JNCX_2.JNCX(time) flux; real JNCX_2.VNCX flux; JNCX_2.VNCX=3.0; real JNCX_2.KNCX micromolar; JNCX_2.KNCX=0.3; real CS1_2.lambda_MT_S1 dimensionless; real CS1_2.lambda_ER_S1 dimensionless; real CS1_2.delta_SPU micromolar_per_coulomb; CS1_2.delta_SPU=18.5; real CS1_2.gamma_S1 dimensionless; CS1_2.gamma_S1=100.0; real CS1_2.gamma_MT dimensionless; CS1_2.gamma_MT=200.0; real CS1_2.gamma_ER dimensionless; CS1_2.gamma_ER=20.0; real CS2_2.lambda_MT_S2 dimensionless; real CS2_2.lambda_ER_S2 dimensionless; real CS2_2.lambda_S1_S2 dimensionless; real CS2_2.gamma_S2 dimensionless; CS2_2.gamma_S2=1.0; real CMT_2.fm(time) dimensionless; real CMT_2.Km micromolar; CMT_2.Km=0.01; real CMT_2.Bm micromolar; CMT_2.Bm=100.0; real model_parameters_3.KO micromolar; model_parameters_3.KO=5400.0; real model_parameters_3.Ki micromolar; model_parameters_3.Ki=145E3; real model_parameters_3.n_PU_base dimensionless; model_parameters_3.n_PU_base=50.0; real PU_membrane_3.I_iCa(time) picoampere; real PU_membrane_3.I_iNa(time) picoampere; real PU_membrane_3.alpha_scale dimensionless; real PU_membrane_3.I_Ca(time) picoampere; real PU_membrane_3.I_Na(time) picoampere; real PU_membrane_3.I_NaP(time) picoampere; real PU_membrane_3.I_NSCC_Ca(time) picoampere; real PU_membrane_3.I_PM(time) picoampere; real PU_membrane_3.I_NSCC_Na(time) picoampere; real I_Ca_3.gCa_ picosiemens; I_Ca_3.gCa_=0.074; real I_Ca_3.kCa per_millivolt; I_Ca_3.kCa=0.013; real I_Ca_3.kVCa per_millivolt; I_Ca_3.kVCa=-0.20; real I_Ca_3.VhCa millivolt; I_Ca_3.VhCa=-85.0; real I_Ca_3.gCa(time) picosiemens; real I_Ca_3.ECa_PU(time) millivolt; real I_Ca_3.T kelvin; I_Ca_3.T=310.16; real I_Ca_3.R attojoule_per_zeptomole_kelvin; I_Ca_3.R=8.314E-3; real I_Ca_3.F femtocoulomb_per_zeptomole; I_Ca_3.F=0.09649; real I_Ca_3.CO micromolar; I_Ca_3.CO=1.8E3; real I_Na_4.gNa picosiemens; I_Na_4.gNa=13.5; real I_Na_4.ENa_PU(time) millivolt; real I_Na_4.NO micromolar; I_Na_4.NO=140.0E3; real I_Na_4.NS1(time) micromolar; when(time=time.min) I_Na_4.NS1=1.01E4; real I_NSCC_Ca_3.gNSCC_Ca_ picosiemens; I_NSCC_Ca_3.gNSCC_Ca_=0.10; real I_NSCC_Ca_3.gNSCC_Ca(time) picosiemens; // Var below replaced by constant in model eqns to satisfy unit correction // real I_NSCC_Ca_3.hNSCC dimensionless; // I_NSCC_Ca_3.hNSCC=3.0; real I_NSCC_Ca_3.ENSCC millivolt; I_NSCC_Ca_3.ENSCC=0.0; real I_NSCC_Ca_3.KNSCC micromolar; I_NSCC_Ca_3.KNSCC=0.12; real I_NSCC_Na_3.gNSCC_Na_ picosiemens; I_NSCC_Na_3.gNSCC_Na_=160.0; real I_NSCC_Na_3.gNSCC_Na(time) picosiemens; real I_PM_3.gPM femtoampere; I_PM_3.gPM=675.0; real I_PM_3.KPM micromolar; I_PM_3.KPM=1.0; real I_NaP_3.gNaP picosiemens; I_NaP_3.gNaP=187.5; real I_NaP_3.KNaP micromolar; I_NaP_3.KNaP=1.0E4; // Var below replaced by constant in model eqns to satisfy unit correction // real I_NaP_3.hNaP dimensionless; // I_NaP_3.hNaP=4.0; real I_NaP_3.ENaP millivolt; I_NaP_3.ENaP=10.0; real was_JSERCA_3.JSERCA(time) flux; real was_JSERCA_3.VSERCA first_order_rate_constant; was_JSERCA_3.VSERCA=1.0E5; real was_JSERCA_3.A2 dimensionless; was_JSERCA_3.A2=6E-4; real was_JSERCA_3.A4 per_micromolar; was_JSERCA_3.A4=3.57; real was_JSERCA_3.A5 per_micromolar; was_JSERCA_3.A5=2.7E-5; real was_JSERCA_3.A6 per_micromolar2; was_JSERCA_3.A6=2.31E-5; real was_JSERCA_3.CER(time) micromolar; when(time=time.min) was_JSERCA_3.CER=200.0; real JMCU_3.JMCU(time) flux; real JMCU_3.VMCU flux; JMCU_3.VMCU=800.0; real JMCU_3.KMCU micromolar; JMCU_3.KMCU=10.0; real JMCU_3.epsilon_INH(time) dimensionless; real JMCU_3.KINH micromolar; JMCU_3.KINH=10.0; // Var below replaced by constant in model eqns to satisfy unit correction // real JMCU_3.hINH dimensionless; // JMCU_3.hINH=4.0; real JMCU_3.CMT(time) micromolar; when(time=time.min) JMCU_3.CMT=0.200; real JIPR_3.JIPR(time) flux; real JIPR_3.kIPR first_order_rate_constant; JIPR_3.kIPR=2000.0; real JIPR_3.k_1 flux; JIPR_3.k_1=6.4; real JIPR_3.k1 first_order_rate_constant; JIPR_3.k1=0.0; real JIPR_3.k2 first_order_rate_constant; JIPR_3.k2=4.0; real JIPR_3.r2 first_order_rate_constant; JIPR_3.r2=250.0; real JIPR_3.r_2 flux; JIPR_3.r_2=0.0; real JIPR_3.r4 first_order_rate_constant; JIPR_3.r4=750.0; real JIPR_3.R1 micromolar; JIPR_3.R1=36.0; real JIPR_3.R3 micromolar; JIPR_3.R3=300.0; real JIPR_3.phi1(time) first_order_rate_constant; real JIPR_3.phi_1(time) flux; real JIPR_3.phi2(time) first_order_rate_constant; real JIPR_3.phi3(time) first_order_rate_constant; real JIPR_3.H(time) dimensionless; when(time=time.min) JIPR_3.H=0.200; real JIPR_3.g_beta first_order_rate_constant; JIPR_3.g_beta=1500.0; real JIPR_3.h_beta dimensionless; JIPR_3.h_beta=4.0; real JIPR_3.g_alpha first_order_rate_constant; JIPR_3.g_alpha=0.85; real JIPR_3.K_beta micromolar; JIPR_3.K_beta=0.35; real JIPR_3.alpha_zeta first_order_rate_constant; real JIPR_3.beta_zeta(time) first_order_rate_constant; real JIPR_3.zeta(time) dimensionless; when(time=time.min) JIPR_3.zeta=0.300; real JIPR_3.K_phi3_act micromolar; JIPR_3.K_phi3_act=0.1; real JIPR_3.K_phi3_inh micromolar; JIPR_3.K_phi3_inh=0.5; real JIPR_3.h_phi3_act dimensionless; JIPR_3.h_phi3_act=3.0; real JIPR_3.h_phi3_inh dimensionless; JIPR_3.h_phi3_inh=3.0; real JIPR_3.g_phi3 first_order_rate_constant; JIPR_3.g_phi3=4.5; real JIPR_3.P micromolar; JIPR_3.P=1.0; real JNCX_3.JNCX(time) flux; real JNCX_3.VNCX flux; JNCX_3.VNCX=3.0; real JNCX_3.KNCX micromolar; JNCX_3.KNCX=0.3; real CS1_3.lambda_MT_S1 dimensionless; real CS1_3.lambda_ER_S1 dimensionless; real CS1_3.delta_SPU micromolar_per_coulomb; CS1_3.delta_SPU=18.5; real CS1_3.gamma_S1 dimensionless; CS1_3.gamma_S1=100.0; real CS1_3.gamma_MT dimensionless; CS1_3.gamma_MT=200.0; real CS1_3.gamma_ER dimensionless; CS1_3.gamma_ER=20.0; real CS2_3.lambda_MT_S2 dimensionless; real CS2_3.lambda_ER_S2 dimensionless; real CS2_3.lambda_S1_S2 dimensionless; real CS2_3.gamma_S2 dimensionless; CS2_3.gamma_S2=1.0; real CMT_3.fm(time) dimensionless; real CMT_3.Km micromolar; CMT_3.Km=0.01; real CMT_3.Bm micromolar; CMT_3.Bm=100.0; real model_parameters_4.KO micromolar; model_parameters_4.KO=5400.0; real model_parameters_4.Ki micromolar; model_parameters_4.Ki=145E3; real model_parameters_4.n_PU_base dimensionless; model_parameters_4.n_PU_base=50.0; real PU_membrane_4.I_iCa(time) picoampere; real PU_membrane_4.I_iNa(time) picoampere; real PU_membrane_4.alpha_scale dimensionless; real PU_membrane_4.I_Ca(time) picoampere; real PU_membrane_4.I_Na(time) picoampere; real PU_membrane_4.I_NaP(time) picoampere; real PU_membrane_4.I_NSCC_Ca(time) picoampere; real PU_membrane_4.I_PM(time) picoampere; real PU_membrane_4.I_NSCC_Na(time) picoampere; real I_Ca_4.gCa_ picosiemens; I_Ca_4.gCa_=0.074; real I_Ca_4.kCa per_millivolt; I_Ca_4.kCa=0.013; real I_Ca_4.kVCa per_millivolt; I_Ca_4.kVCa=-0.20; real I_Ca_4.VhCa millivolt; I_Ca_4.VhCa=-85.0; real I_Ca_4.gCa(time) picosiemens; real I_Ca_4.ECa_PU(time) millivolt; real I_Ca_4.T kelvin; I_Ca_4.T=310.16; real I_Ca_4.R attojoule_per_zeptomole_kelvin; I_Ca_4.R=8.314E-3; real I_Ca_4.F femtocoulomb_per_zeptomole; I_Ca_4.F=0.09649; real I_Ca_4.CO micromolar; I_Ca_4.CO=1.8E3; real I_Na_5.gNa picosiemens; I_Na_5.gNa=13.5; real I_Na_5.ENa_PU(time) millivolt; real I_Na_5.NO micromolar; I_Na_5.NO=140.0E3; real I_Na_5.NS1(time) micromolar; when(time=time.min) I_Na_5.NS1=1.01E4; real I_NSCC_Ca_4.gNSCC_Ca_ picosiemens; I_NSCC_Ca_4.gNSCC_Ca_=0.10; real I_NSCC_Ca_4.gNSCC_Ca(time) picosiemens; // Var below replaced by constant in model eqns to satisfy unit correction // real I_NSCC_Ca_4.hNSCC dimensionless; // I_NSCC_Ca_4.hNSCC=3.0; real I_NSCC_Ca_4.ENSCC millivolt; I_NSCC_Ca_4.ENSCC=0.0; real I_NSCC_Ca_4.KNSCC micromolar; I_NSCC_Ca_4.KNSCC=0.12; real I_NSCC_Na_4.gNSCC_Na_ picosiemens; I_NSCC_Na_4.gNSCC_Na_=160.0; real I_NSCC_Na_4.gNSCC_Na(time) picosiemens; real I_PM_4.gPM femtoampere; I_PM_4.gPM=675.0; real I_PM_4.KPM micromolar; I_PM_4.KPM=1.0; real I_NaP_4.gNaP picosiemens; I_NaP_4.gNaP=187.5; real I_NaP_4.KNaP micromolar; I_NaP_4.KNaP=1.0E4; // Var below replaced by constant in model eqns to satisfy unit correction // real I_NaP_4.hNaP dimensionless; // I_NaP_4.hNaP=4.0; real I_NaP_4.ENaP millivolt; I_NaP_4.ENaP=10.0; real was_JSERCA_4.JSERCA(time) flux; real was_JSERCA_4.VSERCA first_order_rate_constant; was_JSERCA_4.VSERCA=1.0E5; real was_JSERCA_4.A2 dimensionless; was_JSERCA_4.A2=6E-4; real was_JSERCA_4.A4 per_micromolar; was_JSERCA_4.A4=3.57; real was_JSERCA_4.A5 per_micromolar; was_JSERCA_4.A5=2.7E-5; real was_JSERCA_4.A6 per_micromolar2; was_JSERCA_4.A6=2.31E-5; real was_JSERCA_4.CER(time) micromolar; when(time=time.min) was_JSERCA_4.CER=200.0; real JMCU_4.JMCU(time) flux; real JMCU_4.VMCU flux; JMCU_4.VMCU=800.0; real JMCU_4.KMCU micromolar; JMCU_4.KMCU=10.0; real JMCU_4.epsilon_INH(time) dimensionless; real JMCU_4.KINH micromolar; JMCU_4.KINH=10.0; // Var below replaced by constant in model eqns to satisfy unit correction // real JMCU_4.hINH dimensionless; // JMCU_4.hINH=4.0; real JMCU_4.CMT(time) micromolar; when(time=time.min) JMCU_4.CMT=0.200; real JIPR_4.JIPR(time) flux; real JIPR_4.kIPR first_order_rate_constant; JIPR_4.kIPR=2000.0; real JIPR_4.k_1 flux; JIPR_4.k_1=6.4; real JIPR_4.k1 first_order_rate_constant; JIPR_4.k1=0.0; real JIPR_4.k2 first_order_rate_constant; JIPR_4.k2=4.0; real JIPR_4.r2 first_order_rate_constant; JIPR_4.r2=250.0; real JIPR_4.r_2 flux; JIPR_4.r_2=0.0; real JIPR_4.r4 first_order_rate_constant; JIPR_4.r4=750.0; real JIPR_4.R1 micromolar; JIPR_4.R1=36.0; real JIPR_4.R3 micromolar; JIPR_4.R3=300.0; real JIPR_4.phi1(time) first_order_rate_constant; real JIPR_4.phi_1(time) flux; real JIPR_4.phi2(time) first_order_rate_constant; real JIPR_4.phi3(time) first_order_rate_constant; real JIPR_4.H(time) dimensionless; when(time=time.min) JIPR_4.H=0.200; real JIPR_4.g_beta first_order_rate_constant; JIPR_4.g_beta=1500.0; real JIPR_4.h_beta dimensionless; JIPR_4.h_beta=4.0; real JIPR_4.g_alpha first_order_rate_constant; JIPR_4.g_alpha=0.85; real JIPR_4.K_beta micromolar; JIPR_4.K_beta=0.35; real JIPR_4.alpha_zeta first_order_rate_constant; real JIPR_4.beta_zeta(time) first_order_rate_constant; real JIPR_4.zeta(time) dimensionless; when(time=time.min) JIPR_4.zeta=0.300; real JIPR_4.K_phi3_act micromolar; JIPR_4.K_phi3_act=0.1; real JIPR_4.K_phi3_inh micromolar; JIPR_4.K_phi3_inh=0.5; real JIPR_4.h_phi3_act dimensionless; JIPR_4.h_phi3_act=3.0; real JIPR_4.h_phi3_inh dimensionless; JIPR_4.h_phi3_inh=3.0; real JIPR_4.g_phi3 first_order_rate_constant; JIPR_4.g_phi3=4.5; real JIPR_4.P micromolar; JIPR_4.P=1.0; real JNCX_4.JNCX(time) flux; real JNCX_4.VNCX flux; JNCX_4.VNCX=3.0; real JNCX_4.KNCX micromolar; JNCX_4.KNCX=0.3; real CS1_4.lambda_MT_S1 dimensionless; real CS1_4.lambda_ER_S1 dimensionless; real CS1_4.delta_SPU micromolar_per_coulomb; CS1_4.delta_SPU=18.5; real CS1_4.gamma_S1 dimensionless; CS1_4.gamma_S1=100.0; real CS1_4.gamma_MT dimensionless; CS1_4.gamma_MT=200.0; real CS1_4.gamma_ER dimensionless; CS1_4.gamma_ER=20.0; real CS2_4.lambda_MT_S2 dimensionless; real CS2_4.lambda_ER_S2 dimensionless; real CS2_4.lambda_S1_S2 dimensionless; real CS2_4.gamma_S2 dimensionless; CS2_4.gamma_S2=1.0; real CMT_4.fm(time) dimensionless; real CMT_4.Km micromolar; CMT_4.Km=0.01; real CMT_4.Bm micromolar; CMT_4.Bm=100.0; real model_parameters_5.KO micromolar; model_parameters_5.KO=5400.0; real model_parameters_5.Ki micromolar; model_parameters_5.Ki=145E3; real model_parameters_5.n_PU_base dimensionless; model_parameters_5.n_PU_base=50.0; real PU_membrane_5.I_iCa(time) picoampere; real PU_membrane_5.I_iNa(time) picoampere; real PU_membrane_5.alpha_scale dimensionless; real PU_membrane_5.I_Ca(time) picoampere; real PU_membrane_5.I_Na(time) picoampere; real PU_membrane_5.I_NaP(time) picoampere; real PU_membrane_5.I_NSCC_Ca(time) picoampere; real PU_membrane_5.I_PM(time) picoampere; real PU_membrane_5.I_NSCC_Na(time) picoampere; real I_Ca_5.gCa_ picosiemens; I_Ca_5.gCa_=0.074; real I_Ca_5.kCa per_millivolt; I_Ca_5.kCa=0.013; real I_Ca_5.kVCa per_millivolt; I_Ca_5.kVCa=-0.20; real I_Ca_5.VhCa millivolt; I_Ca_5.VhCa=-85.0; real I_Ca_5.gCa(time) picosiemens; real I_Ca_5.ECa_PU(time) millivolt; real I_Ca_5.T kelvin; I_Ca_5.T=310.16; real I_Ca_5.R attojoule_per_zeptomole_kelvin; I_Ca_5.R=8.314E-3; real I_Ca_5.F femtocoulomb_per_zeptomole; I_Ca_5.F=0.09649; real I_Ca_5.CO micromolar; I_Ca_5.CO=1.8E3; real I_Na_6.gNa picosiemens; I_Na_6.gNa=13.5; real I_Na_6.ENa_PU(time) millivolt; real I_Na_6.NO micromolar; I_Na_6.NO=140.0E3; real I_Na_6.NS1(time) micromolar; when(time=time.min) I_Na_6.NS1=1.01E4; real I_NSCC_Ca_5.gNSCC_Ca_ picosiemens; I_NSCC_Ca_5.gNSCC_Ca_=0.10; real I_NSCC_Ca_5.gNSCC_Ca(time) picosiemens; // Var below replaced by constant in model eqns to satisfy unit correction // real I_NSCC_Ca_5.hNSCC dimensionless; // I_NSCC_Ca_5.hNSCC=3.0; real I_NSCC_Ca_5.ENSCC millivolt; I_NSCC_Ca_5.ENSCC=0.0; real I_NSCC_Ca_5.KNSCC micromolar; I_NSCC_Ca_5.KNSCC=0.12; real I_NSCC_Na_5.gNSCC_Na_ picosiemens; I_NSCC_Na_5.gNSCC_Na_=160.0; real I_NSCC_Na_5.gNSCC_Na(time) picosiemens; real I_PM_5.gPM femtoampere; I_PM_5.gPM=675.0; real I_PM_5.KPM micromolar; I_PM_5.KPM=1.0; real I_NaP_5.gNaP picosiemens; I_NaP_5.gNaP=187.5; real I_NaP_5.KNaP micromolar; I_NaP_5.KNaP=1.0E4; // Var below replaced by constant in model eqns to satisfy unit correction // real I_NaP_5.hNaP dimensionless; // I_NaP_5.hNaP=4.0; real I_NaP_5.ENaP millivolt; I_NaP_5.ENaP=10.0; real was_JSERCA_5.JSERCA(time) flux; real was_JSERCA_5.VSERCA first_order_rate_constant; was_JSERCA_5.VSERCA=1.0E5; real was_JSERCA_5.A2 dimensionless; was_JSERCA_5.A2=6E-4; real was_JSERCA_5.A4 per_micromolar; was_JSERCA_5.A4=3.57; real was_JSERCA_5.A5 per_micromolar; was_JSERCA_5.A5=2.7E-5; real was_JSERCA_5.A6 per_micromolar2; was_JSERCA_5.A6=2.31E-5; real was_JSERCA_5.CER(time) micromolar; when(time=time.min) was_JSERCA_5.CER=200.0; real JMCU_5.JMCU(time) flux; real JMCU_5.VMCU flux; JMCU_5.VMCU=800.0; real JMCU_5.KMCU micromolar; JMCU_5.KMCU=10.0; real JMCU_5.epsilon_INH(time) dimensionless; real JMCU_5.KINH micromolar; JMCU_5.KINH=10.0; // Var below replaced by constant in model eqns to satisfy unit correction // real JMCU_5.hINH dimensionless; // JMCU_5.hINH=4.0; real JMCU_5.CMT(time) micromolar; when(time=time.min) JMCU_5.CMT=0.200; real JIPR_5.JIPR(time) flux; real JIPR_5.kIPR first_order_rate_constant; JIPR_5.kIPR=2000.0; real JIPR_5.k_1 flux; JIPR_5.k_1=6.4; real JIPR_5.k1 first_order_rate_constant; JIPR_5.k1=0.0; real JIPR_5.k2 first_order_rate_constant; JIPR_5.k2=4.0; real JIPR_5.r2 first_order_rate_constant; JIPR_5.r2=250.0; real JIPR_5.r_2 flux; JIPR_5.r_2=0.0; real JIPR_5.r4 first_order_rate_constant; JIPR_5.r4=750.0; real JIPR_5.R1 micromolar; JIPR_5.R1=36.0; real JIPR_5.R3 micromolar; JIPR_5.R3=300.0; real JIPR_5.phi1(time) first_order_rate_constant; real JIPR_5.phi_1(time) flux; real JIPR_5.phi2(time) first_order_rate_constant; real JIPR_5.phi3(time) first_order_rate_constant; real JIPR_5.H(time) dimensionless; when(time=time.min) JIPR_5.H=0.200; real JIPR_5.g_beta first_order_rate_constant; JIPR_5.g_beta=1500.0; real JIPR_5.h_beta dimensionless; JIPR_5.h_beta=4.0; real JIPR_5.g_alpha first_order_rate_constant; JIPR_5.g_alpha=0.85; real JIPR_5.K_beta micromolar; JIPR_5.K_beta=0.35; real JIPR_5.alpha_zeta first_order_rate_constant; real JIPR_5.beta_zeta(time) first_order_rate_constant; real JIPR_5.zeta(time) dimensionless; when(time=time.min) JIPR_5.zeta=0.300; real JIPR_5.K_phi3_act micromolar; JIPR_5.K_phi3_act=0.1; real JIPR_5.K_phi3_inh micromolar; JIPR_5.K_phi3_inh=0.5; real JIPR_5.h_phi3_act dimensionless; JIPR_5.h_phi3_act=3.0; real JIPR_5.h_phi3_inh dimensionless; JIPR_5.h_phi3_inh=3.0; real JIPR_5.g_phi3 first_order_rate_constant; JIPR_5.g_phi3=4.5; real JIPR_5.P micromolar; JIPR_5.P=1.0; real JNCX_5.JNCX(time) flux; real JNCX_5.VNCX flux; JNCX_5.VNCX=3.0; real JNCX_5.KNCX micromolar; JNCX_5.KNCX=0.3; real CS1_5.lambda_MT_S1 dimensionless; real CS1_5.lambda_ER_S1 dimensionless; real CS1_5.delta_SPU micromolar_per_coulomb; CS1_5.delta_SPU=18.5; real CS1_5.gamma_S1 dimensionless; CS1_5.gamma_S1=100.0; real CS1_5.gamma_MT dimensionless; CS1_5.gamma_MT=200.0; real CS1_5.gamma_ER dimensionless; CS1_5.gamma_ER=20.0; real CS2_5.lambda_MT_S2 dimensionless; real CS2_5.lambda_ER_S2 dimensionless; real CS2_5.lambda_S1_S2 dimensionless; real CS2_5.gamma_S2 dimensionless; CS2_5.gamma_S2=1.0; real CMT_5.fm(time) dimensionless; real CMT_5.Km micromolar; CMT_5.Km=0.01; real CMT_5.Bm micromolar; CMT_5.Bm=100.0; real model_parameters_6.KO micromolar; model_parameters_6.KO=5400.0; real model_parameters_6.Ki micromolar; model_parameters_6.Ki=145E3; real model_parameters_6.n_PU_base dimensionless; model_parameters_6.n_PU_base=50.0; real PU_membrane_6.I_iCa(time) picoampere; real PU_membrane_6.I_iNa(time) picoampere; real PU_membrane_6.alpha_scale dimensionless; real PU_membrane_6.I_Ca(time) picoampere; real PU_membrane_6.I_Na(time) picoampere; real PU_membrane_6.I_NaP(time) picoampere; real PU_membrane_6.I_NSCC_Ca(time) picoampere; real PU_membrane_6.I_PM(time) picoampere; real PU_membrane_6.I_NSCC_Na(time) picoampere; real I_Ca_6.gCa_ picosiemens; I_Ca_6.gCa_=0.074; real I_Ca_6.kCa per_millivolt; I_Ca_6.kCa=0.013; real I_Ca_6.kVCa per_millivolt; I_Ca_6.kVCa=-0.20; real I_Ca_6.VhCa millivolt; I_Ca_6.VhCa=-85.0; real I_Ca_6.gCa(time) picosiemens; real I_Ca_6.ECa_PU(time) millivolt; real I_Ca_6.T kelvin; I_Ca_6.T=310.16; real I_Ca_6.R attojoule_per_zeptomole_kelvin; I_Ca_6.R=8.314E-3; real I_Ca_6.F femtocoulomb_per_zeptomole; I_Ca_6.F=0.09649; real I_Ca_6.CO micromolar; I_Ca_6.CO=1.8E3; real I_Na_7.gNa picosiemens; I_Na_7.gNa=13.5; real I_Na_7.ENa_PU(time) millivolt; real I_Na_7.NO micromolar; I_Na_7.NO=140.0E3; real I_Na_7.NS1(time) micromolar; when(time=time.min) I_Na_7.NS1=1.01E4; real I_NSCC_Ca_6.gNSCC_Ca_ picosiemens; I_NSCC_Ca_6.gNSCC_Ca_=0.10; real I_NSCC_Ca_6.gNSCC_Ca(time) picosiemens; // Var below replaced by constant in model eqns to satisfy unit correction // real I_NSCC_Ca_6.hNSCC dimensionless; // I_NSCC_Ca_6.hNSCC=3.0; real I_NSCC_Ca_6.ENSCC millivolt; I_NSCC_Ca_6.ENSCC=0.0; real I_NSCC_Ca_6.KNSCC micromolar; I_NSCC_Ca_6.KNSCC=0.12; real I_NSCC_Na_6.gNSCC_Na_ picosiemens; I_NSCC_Na_6.gNSCC_Na_=160.0; real I_NSCC_Na_6.gNSCC_Na(time) picosiemens; real I_PM_6.gPM femtoampere; I_PM_6.gPM=675.0; real I_PM_6.KPM micromolar; I_PM_6.KPM=1.0; real I_NaP_6.gNaP picosiemens; I_NaP_6.gNaP=187.5; real I_NaP_6.KNaP micromolar; I_NaP_6.KNaP=1.0E4; // Var below replaced by constant in model eqns to satisfy unit correction // real I_NaP_6.hNaP dimensionless; // I_NaP_6.hNaP=4.0; real I_NaP_6.ENaP millivolt; I_NaP_6.ENaP=10.0; real was_JSERCA_6.JSERCA(time) flux; real was_JSERCA_6.VSERCA first_order_rate_constant; was_JSERCA_6.VSERCA=1.0E5; real was_JSERCA_6.A2 dimensionless; was_JSERCA_6.A2=6E-4; real was_JSERCA_6.A4 per_micromolar; was_JSERCA_6.A4=3.57; real was_JSERCA_6.A5 per_micromolar; was_JSERCA_6.A5=2.7E-5; real was_JSERCA_6.A6 per_micromolar2; was_JSERCA_6.A6=2.31E-5; real was_JSERCA_6.CER(time) micromolar; when(time=time.min) was_JSERCA_6.CER=200.0; real JMCU_6.JMCU(time) flux; real JMCU_6.VMCU flux; JMCU_6.VMCU=800.0; real JMCU_6.KMCU micromolar; JMCU_6.KMCU=10.0; real JMCU_6.epsilon_INH(time) dimensionless; real JMCU_6.KINH micromolar; JMCU_6.KINH=10.0; // Var below replaced by constant in model eqns to satisfy unit correction // real JMCU_6.hINH dimensionless; // JMCU_6.hINH=4.0; real JMCU_6.CMT(time) micromolar; when(time=time.min) JMCU_6.CMT=0.200; real JIPR_6.JIPR(time) flux; real JIPR_6.kIPR first_order_rate_constant; JIPR_6.kIPR=2000.0; real JIPR_6.k_1 flux; JIPR_6.k_1=6.4; real JIPR_6.k1 first_order_rate_constant; JIPR_6.k1=0.0; real JIPR_6.k2 first_order_rate_constant; JIPR_6.k2=4.0; real JIPR_6.r2 first_order_rate_constant; JIPR_6.r2=250.0; real JIPR_6.r_2 flux; JIPR_6.r_2=0.0; real JIPR_6.r4 first_order_rate_constant; JIPR_6.r4=750.0; real JIPR_6.R1 micromolar; JIPR_6.R1=36.0; real JIPR_6.R3 micromolar; JIPR_6.R3=300.0; real JIPR_6.phi1(time) first_order_rate_constant; real JIPR_6.phi_1(time) flux; real JIPR_6.phi2(time) first_order_rate_constant; real JIPR_6.phi3(time) first_order_rate_constant; real JIPR_6.H(time) dimensionless; when(time=time.min) JIPR_6.H=0.200; real JIPR_6.g_beta first_order_rate_constant; JIPR_6.g_beta=1500.0; real JIPR_6.h_beta dimensionless; JIPR_6.h_beta=4.0; real JIPR_6.g_alpha first_order_rate_constant; JIPR_6.g_alpha=0.85; real JIPR_6.K_beta micromolar; JIPR_6.K_beta=0.35; real JIPR_6.alpha_zeta first_order_rate_constant; real JIPR_6.beta_zeta(time) first_order_rate_constant; real JIPR_6.zeta(time) dimensionless; when(time=time.min) JIPR_6.zeta=0.300; real JIPR_6.K_phi3_act micromolar; JIPR_6.K_phi3_act=0.1; real JIPR_6.K_phi3_inh micromolar; JIPR_6.K_phi3_inh=0.5; real JIPR_6.h_phi3_act dimensionless; JIPR_6.h_phi3_act=3.0; real JIPR_6.h_phi3_inh dimensionless; JIPR_6.h_phi3_inh=3.0; real JIPR_6.g_phi3 first_order_rate_constant; JIPR_6.g_phi3=4.5; real JIPR_6.P micromolar; JIPR_6.P=1.0; real JNCX_6.JNCX(time) flux; real JNCX_6.VNCX flux; JNCX_6.VNCX=3.0; real JNCX_6.KNCX micromolar; JNCX_6.KNCX=0.3; real CS1_6.lambda_MT_S1 dimensionless; real CS1_6.lambda_ER_S1 dimensionless; real CS1_6.delta_SPU micromolar_per_coulomb; CS1_6.delta_SPU=18.5; real CS1_6.gamma_S1 dimensionless; CS1_6.gamma_S1=100.0; real CS1_6.gamma_MT dimensionless; CS1_6.gamma_MT=200.0; real CS1_6.gamma_ER dimensionless; CS1_6.gamma_ER=20.0; real CS2_6.lambda_MT_S2 dimensionless; real CS2_6.lambda_ER_S2 dimensionless; real CS2_6.lambda_S1_S2 dimensionless; real CS2_6.gamma_S2 dimensionless; CS2_6.gamma_S2=1.0; real CMT_6.fm(time) dimensionless; real CMT_6.Km micromolar; CMT_6.Km=0.01; real CMT_6.Bm micromolar; CMT_6.Bm=100.0; real model_parameters_7.KO micromolar; model_parameters_7.KO=5400.0; real model_parameters_7.Ki micromolar; model_parameters_7.Ki=145E3; real model_parameters_7.n_PU_base dimensionless; model_parameters_7.n_PU_base=50.0; real PU_membrane_7.I_iCa(time) picoampere; real PU_membrane_7.I_iNa(time) picoampere; real PU_membrane_7.alpha_scale dimensionless; real PU_membrane_7.I_Ca(time) picoampere; real PU_membrane_7.I_Na(time) picoampere; real PU_membrane_7.I_NaP(time) picoampere; real PU_membrane_7.I_NSCC_Ca(time) picoampere; real PU_membrane_7.I_PM(time) picoampere; real PU_membrane_7.I_NSCC_Na(time) picoampere; real I_Ca_7.gCa_ picosiemens; I_Ca_7.gCa_=0.074; real I_Ca_7.kCa per_millivolt; I_Ca_7.kCa=0.013; real I_Ca_7.kVCa per_millivolt; I_Ca_7.kVCa=-0.20; real I_Ca_7.VhCa millivolt; I_Ca_7.VhCa=-85.0; real I_Ca_7.gCa(time) picosiemens; real I_Ca_7.ECa_PU(time) millivolt; real I_Ca_7.T kelvin; I_Ca_7.T=310.16; real I_Ca_7.R attojoule_per_zeptomole_kelvin; I_Ca_7.R=8.314E-3; real I_Ca_7.F femtocoulomb_per_zeptomole; I_Ca_7.F=0.09649; real I_Ca_7.CO micromolar; I_Ca_7.CO=1.8E3; real I_Na_8.gNa picosiemens; I_Na_8.gNa=13.5; real I_Na_8.ENa_PU(time) millivolt; real I_Na_8.NO micromolar; I_Na_8.NO=140.0E3; real I_Na_8.NS1(time) micromolar; when(time=time.min) I_Na_8.NS1=1.01E4; real I_NSCC_Ca_7.gNSCC_Ca_ picosiemens; I_NSCC_Ca_7.gNSCC_Ca_=0.10; real I_NSCC_Ca_7.gNSCC_Ca(time) picosiemens; // Var below replaced by constant in model eqns to satisfy unit correction // real I_NSCC_Ca_7.hNSCC dimensionless; // I_NSCC_Ca_7.hNSCC=3.0; real I_NSCC_Ca_7.ENSCC millivolt; I_NSCC_Ca_7.ENSCC=0.0; real I_NSCC_Ca_7.KNSCC micromolar; I_NSCC_Ca_7.KNSCC=0.12; real I_NSCC_Na_7.gNSCC_Na_ picosiemens; I_NSCC_Na_7.gNSCC_Na_=160.0; real I_NSCC_Na_7.gNSCC_Na(time) picosiemens; real I_PM_7.gPM femtoampere; I_PM_7.gPM=675.0; real I_PM_7.KPM micromolar; I_PM_7.KPM=1.0; real I_NaP_7.gNaP picosiemens; I_NaP_7.gNaP=187.5; real I_NaP_7.KNaP micromolar; I_NaP_7.KNaP=1.0E4; // Var below replaced by constant in model eqns to satisfy unit correction // real I_NaP_7.hNaP dimensionless; // I_NaP_7.hNaP=4.0; real I_NaP_7.ENaP millivolt; I_NaP_7.ENaP=10.0; real was_JSERCA_7.JSERCA(time) flux; real was_JSERCA_7.VSERCA first_order_rate_constant; was_JSERCA_7.VSERCA=1.0E5; real was_JSERCA_7.A2 dimensionless; was_JSERCA_7.A2=6E-4; real was_JSERCA_7.A4 per_micromolar; was_JSERCA_7.A4=3.57; real was_JSERCA_7.A5 per_micromolar; was_JSERCA_7.A5=2.7E-5; real was_JSERCA_7.A6 per_micromolar2; was_JSERCA_7.A6=2.31E-5; real was_JSERCA_7.CER(time) micromolar; when(time=time.min) was_JSERCA_7.CER=200.0; real JMCU_7.JMCU(time) flux; real JMCU_7.VMCU flux; JMCU_7.VMCU=800.0; real JMCU_7.KMCU micromolar; JMCU_7.KMCU=10.0; real JMCU_7.epsilon_INH(time) dimensionless; real JMCU_7.KINH micromolar; JMCU_7.KINH=10.0; // Var below replaced by constant in model eqns to satisfy unit correction // real JMCU_7.hINH dimensionless; // JMCU_7.hINH=4.0; real JMCU_7.CMT(time) micromolar; when(time=time.min) JMCU_7.CMT=0.200; real JIPR_7.JIPR(time) flux; real JIPR_7.kIPR first_order_rate_constant; JIPR_7.kIPR=2000.0; real JIPR_7.k_1 flux; JIPR_7.k_1=6.4; real JIPR_7.k1 first_order_rate_constant; JIPR_7.k1=0.0; real JIPR_7.k2 first_order_rate_constant; JIPR_7.k2=4.0; real JIPR_7.r2 first_order_rate_constant; JIPR_7.r2=250.0; real JIPR_7.r_2 flux; JIPR_7.r_2=0.0; real JIPR_7.r4 first_order_rate_constant; JIPR_7.r4=750.0; real JIPR_7.R1 micromolar; JIPR_7.R1=36.0; real JIPR_7.R3 micromolar; JIPR_7.R3=300.0; real JIPR_7.phi1(time) first_order_rate_constant; real JIPR_7.phi_1(time) flux; real JIPR_7.phi2(time) first_order_rate_constant; real JIPR_7.phi3(time) first_order_rate_constant; real JIPR_7.H(time) dimensionless; when(time=time.min) JIPR_7.H=0.200; real JIPR_7.g_beta first_order_rate_constant; JIPR_7.g_beta=1500.0; real JIPR_7.h_beta dimensionless; JIPR_7.h_beta=4.0; real JIPR_7.g_alpha first_order_rate_constant; JIPR_7.g_alpha=0.85; real JIPR_7.K_beta micromolar; JIPR_7.K_beta=0.35; real JIPR_7.alpha_zeta first_order_rate_constant; real JIPR_7.beta_zeta(time) first_order_rate_constant; real JIPR_7.zeta(time) dimensionless; when(time=time.min) JIPR_7.zeta=0.300; real JIPR_7.K_phi3_act micromolar; JIPR_7.K_phi3_act=0.1; real JIPR_7.K_phi3_inh micromolar; JIPR_7.K_phi3_inh=0.5; real JIPR_7.h_phi3_act dimensionless; JIPR_7.h_phi3_act=3.0; real JIPR_7.h_phi3_inh dimensionless; JIPR_7.h_phi3_inh=3.0; real JIPR_7.g_phi3 first_order_rate_constant; JIPR_7.g_phi3=4.5; real JIPR_7.P micromolar; JIPR_7.P=1.0; real JNCX_7.JNCX(time) flux; real JNCX_7.VNCX flux; JNCX_7.VNCX=3.0; real JNCX_7.KNCX micromolar; JNCX_7.KNCX=0.3; real CS1_7.lambda_MT_S1 dimensionless; real CS1_7.lambda_ER_S1 dimensionless; real CS1_7.delta_SPU micromolar_per_coulomb; CS1_7.delta_SPU=18.5; real CS1_7.gamma_S1 dimensionless; CS1_7.gamma_S1=100.0; real CS1_7.gamma_MT dimensionless; CS1_7.gamma_MT=200.0; real CS1_7.gamma_ER dimensionless; CS1_7.gamma_ER=20.0; real CS2_7.lambda_MT_S2 dimensionless; real CS2_7.lambda_ER_S2 dimensionless; real CS2_7.lambda_S1_S2 dimensionless; real CS2_7.gamma_S2 dimensionless; CS2_7.gamma_S2=1.0; real CMT_7.fm(time) dimensionless; real CMT_7.Km micromolar; CMT_7.Km=0.01; real CMT_7.Bm micromolar; CMT_7.Bm=100.0; real model_parameters_8.KO micromolar; model_parameters_8.KO=5400.0; real model_parameters_8.Ki micromolar; model_parameters_8.Ki=145E3; real model_parameters_8.n_PU_base dimensionless; model_parameters_8.n_PU_base=50.0; real PU_membrane_8.I_iCa(time) picoampere; real PU_membrane_8.I_iNa(time) picoampere; real PU_membrane_8.alpha_scale dimensionless; real PU_membrane_8.I_Ca(time) picoampere; real PU_membrane_8.I_Na(time) picoampere; real PU_membrane_8.I_NaP(time) picoampere; real PU_membrane_8.I_NSCC_Ca(time) picoampere; real PU_membrane_8.I_PM(time) picoampere; real PU_membrane_8.I_NSCC_Na(time) picoampere; real I_Ca_8.gCa_ picosiemens; I_Ca_8.gCa_=0.074; real I_Ca_8.kCa per_millivolt; I_Ca_8.kCa=0.013; real I_Ca_8.kVCa per_millivolt; I_Ca_8.kVCa=-0.20; real I_Ca_8.VhCa millivolt; I_Ca_8.VhCa=-85.0; real I_Ca_8.gCa(time) picosiemens; real I_Ca_8.ECa_PU(time) millivolt; real I_Ca_8.T kelvin; I_Ca_8.T=310.16; real I_Ca_8.R attojoule_per_zeptomole_kelvin; I_Ca_8.R=8.314E-3; real I_Ca_8.F femtocoulomb_per_zeptomole; I_Ca_8.F=0.09649; real I_Ca_8.CO micromolar; I_Ca_8.CO=1.8E3; real I_Na_9.gNa picosiemens; I_Na_9.gNa=13.5; real I_Na_9.ENa_PU(time) millivolt; real I_Na_9.NO micromolar; I_Na_9.NO=140.0E3; real I_Na_9.NS1(time) micromolar; when(time=time.min) I_Na_9.NS1=1.01E4; real I_NSCC_Ca_8.gNSCC_Ca_ picosiemens; I_NSCC_Ca_8.gNSCC_Ca_=0.10; real I_NSCC_Ca_8.gNSCC_Ca(time) picosiemens; // Var below replaced by constant in model eqns to satisfy unit correction // real I_NSCC_Ca_8.hNSCC dimensionless; // I_NSCC_Ca_8.hNSCC=3.0; real I_NSCC_Ca_8.ENSCC millivolt; I_NSCC_Ca_8.ENSCC=0.0; real I_NSCC_Ca_8.KNSCC micromolar; I_NSCC_Ca_8.KNSCC=0.12; real I_NSCC_Na_8.gNSCC_Na_ picosiemens; I_NSCC_Na_8.gNSCC_Na_=160.0; real I_NSCC_Na_8.gNSCC_Na(time) picosiemens; real I_PM_8.gPM femtoampere; I_PM_8.gPM=675.0; real I_PM_8.KPM micromolar; I_PM_8.KPM=1.0; real I_NaP_8.gNaP picosiemens; I_NaP_8.gNaP=187.5; real I_NaP_8.KNaP micromolar; I_NaP_8.KNaP=1.0E4; // Var below replaced by constant in model eqns to satisfy unit correction // real I_NaP_8.hNaP dimensionless; // I_NaP_8.hNaP=4.0; real I_NaP_8.ENaP millivolt; I_NaP_8.ENaP=10.0; real was_JSERCA_8.JSERCA(time) flux; real was_JSERCA_8.VSERCA first_order_rate_constant; was_JSERCA_8.VSERCA=1.0E5; real was_JSERCA_8.A2 dimensionless; was_JSERCA_8.A2=6E-4; real was_JSERCA_8.A4 per_micromolar; was_JSERCA_8.A4=3.57; real was_JSERCA_8.A5 per_micromolar; was_JSERCA_8.A5=2.7E-5; real was_JSERCA_8.A6 per_micromolar2; was_JSERCA_8.A6=2.31E-5; real was_JSERCA_8.CER(time) micromolar; when(time=time.min) was_JSERCA_8.CER=200.0; real JMCU_8.JMCU(time) flux; real JMCU_8.VMCU flux; JMCU_8.VMCU=800.0; real JMCU_8.KMCU micromolar; JMCU_8.KMCU=10.0; real JMCU_8.epsilon_INH(time) dimensionless; real JMCU_8.KINH micromolar; JMCU_8.KINH=10.0; // Var below replaced by constant in model eqns to satisfy unit correction // real JMCU_8.hINH dimensionless; // JMCU_8.hINH=4.0; real JMCU_8.CMT(time) micromolar; when(time=time.min) JMCU_8.CMT=0.200; real JIPR_8.JIPR(time) flux; real JIPR_8.kIPR first_order_rate_constant; JIPR_8.kIPR=2000.0; real JIPR_8.k_1 flux; JIPR_8.k_1=6.4; real JIPR_8.k1 first_order_rate_constant; JIPR_8.k1=0.0; real JIPR_8.k2 first_order_rate_constant; JIPR_8.k2=4.0; real JIPR_8.r2 first_order_rate_constant; JIPR_8.r2=250.0; real JIPR_8.r_2 flux; JIPR_8.r_2=0.0; real JIPR_8.r4 first_order_rate_constant; JIPR_8.r4=750.0; real JIPR_8.R1 micromolar; JIPR_8.R1=36.0; real JIPR_8.R3 micromolar; JIPR_8.R3=300.0; real JIPR_8.phi1(time) first_order_rate_constant; real JIPR_8.phi_1(time) flux; real JIPR_8.phi2(time) first_order_rate_constant; real JIPR_8.phi3(time) first_order_rate_constant; real JIPR_8.H(time) dimensionless; when(time=time.min) JIPR_8.H=0.200; real JIPR_8.g_beta first_order_rate_constant; JIPR_8.g_beta=1500.0; real JIPR_8.h_beta dimensionless; JIPR_8.h_beta=4.0; real JIPR_8.g_alpha first_order_rate_constant; JIPR_8.g_alpha=0.85; real JIPR_8.K_beta micromolar; JIPR_8.K_beta=0.35; real JIPR_8.alpha_zeta first_order_rate_constant; real JIPR_8.beta_zeta(time) first_order_rate_constant; real JIPR_8.zeta(time) dimensionless; when(time=time.min) JIPR_8.zeta=0.300; real JIPR_8.K_phi3_act micromolar; JIPR_8.K_phi3_act=0.1; real JIPR_8.K_phi3_inh micromolar; JIPR_8.K_phi3_inh=0.5; real JIPR_8.h_phi3_act dimensionless; JIPR_8.h_phi3_act=3.0; real JIPR_8.h_phi3_inh dimensionless; JIPR_8.h_phi3_inh=3.0; real JIPR_8.g_phi3 first_order_rate_constant; JIPR_8.g_phi3=4.5; real JIPR_8.P micromolar; JIPR_8.P=1.0; real JNCX_8.JNCX(time) flux; real JNCX_8.VNCX flux; JNCX_8.VNCX=3.0; real JNCX_8.KNCX micromolar; JNCX_8.KNCX=0.3; real CS1_8.lambda_MT_S1 dimensionless; real CS1_8.lambda_ER_S1 dimensionless; real CS1_8.delta_SPU micromolar_per_coulomb; CS1_8.delta_SPU=18.5; real CS1_8.gamma_S1 dimensionless; CS1_8.gamma_S1=100.0; real CS1_8.gamma_MT dimensionless; CS1_8.gamma_MT=200.0; real CS1_8.gamma_ER dimensionless; CS1_8.gamma_ER=20.0; real CS2_8.lambda_MT_S2 dimensionless; real CS2_8.lambda_ER_S2 dimensionless; real CS2_8.lambda_S1_S2 dimensionless; real CS2_8.gamma_S2 dimensionless; CS2_8.gamma_S2=1.0; real CMT_8.fm(time) dimensionless; real CMT_8.Km micromolar; CMT_8.Km=0.01; real CMT_8.Bm micromolar; CMT_8.Bm=100.0; real model_parameters_9.KO micromolar; model_parameters_9.KO=5400.0; real model_parameters_9.Ki micromolar; model_parameters_9.Ki=145E3; real model_parameters_9.n_PU_base dimensionless; model_parameters_9.n_PU_base=50.0; real PU_membrane_9.I_iCa(time) picoampere; real PU_membrane_9.I_iNa(time) picoampere; real PU_membrane_9.alpha_scale dimensionless; real PU_membrane_9.I_Ca(time) picoampere; real PU_membrane_9.I_Na(time) picoampere; real PU_membrane_9.I_NaP(time) picoampere; real PU_membrane_9.I_NSCC_Ca(time) picoampere; real PU_membrane_9.I_PM(time) picoampere; real PU_membrane_9.I_NSCC_Na(time) picoampere; real I_Ca_9.gCa_ picosiemens; I_Ca_9.gCa_=0.074; real I_Ca_9.kCa per_millivolt; I_Ca_9.kCa=0.013; real I_Ca_9.kVCa per_millivolt; I_Ca_9.kVCa=-0.20; real I_Ca_9.VhCa millivolt; I_Ca_9.VhCa=-85.0; real I_Ca_9.gCa(time) picosiemens; real I_Ca_9.ECa_PU(time) millivolt; real I_Ca_9.T kelvin; I_Ca_9.T=310.16; real I_Ca_9.R attojoule_per_zeptomole_kelvin; I_Ca_9.R=8.314E-3; real I_Ca_9.F femtocoulomb_per_zeptomole; I_Ca_9.F=0.09649; real I_Ca_9.CO micromolar; I_Ca_9.CO=1.8E3; real I_Na_10.gNa picosiemens; I_Na_10.gNa=13.5; real I_Na_10.ENa_PU(time) millivolt; real I_Na_10.NO micromolar; I_Na_10.NO=140.0E3; real I_Na_10.NS1(time) micromolar; when(time=time.min) I_Na_10.NS1=1.01E4; real I_NSCC_Ca_9.gNSCC_Ca_ picosiemens; I_NSCC_Ca_9.gNSCC_Ca_=0.10; real I_NSCC_Ca_9.gNSCC_Ca(time) picosiemens; // Var below replaced by constant in model eqns to satisfy unit correction // real I_NSCC_Ca_9.hNSCC dimensionless; // I_NSCC_Ca_9.hNSCC=3.0; real I_NSCC_Ca_9.ENSCC millivolt; I_NSCC_Ca_9.ENSCC=0.0; real I_NSCC_Ca_9.KNSCC micromolar; I_NSCC_Ca_9.KNSCC=0.12; real I_NSCC_Na_9.gNSCC_Na_ picosiemens; I_NSCC_Na_9.gNSCC_Na_=160.0; real I_NSCC_Na_9.gNSCC_Na(time) picosiemens; real I_PM_9.gPM femtoampere; I_PM_9.gPM=675.0; real I_PM_9.KPM micromolar; I_PM_9.KPM=1.0; real I_NaP_9.gNaP picosiemens; I_NaP_9.gNaP=187.5; real I_NaP_9.KNaP micromolar; I_NaP_9.KNaP=1.0E4; // Var below replaced by constant in model eqns to satisfy unit correction // real I_NaP_9.hNaP dimensionless; // I_NaP_9.hNaP=4.0; real I_NaP_9.ENaP millivolt; I_NaP_9.ENaP=10.0; real was_JSERCA_9.JSERCA(time) flux; real was_JSERCA_9.VSERCA first_order_rate_constant; was_JSERCA_9.VSERCA=1.0E5; real was_JSERCA_9.A2 dimensionless; was_JSERCA_9.A2=6E-4; real was_JSERCA_9.A4 per_micromolar; was_JSERCA_9.A4=3.57; real was_JSERCA_9.A5 per_micromolar; was_JSERCA_9.A5=2.7E-5; real was_JSERCA_9.A6 per_micromolar2; was_JSERCA_9.A6=2.31E-5; real was_JSERCA_9.CER(time) micromolar; when(time=time.min) was_JSERCA_9.CER=200.0; real JMCU_9.JMCU(time) flux; real JMCU_9.VMCU flux; JMCU_9.VMCU=800.0; real JMCU_9.KMCU micromolar; JMCU_9.KMCU=10.0; real JMCU_9.epsilon_INH(time) dimensionless; real JMCU_9.KINH micromolar; JMCU_9.KINH=10.0; // Var below replaced by constant in model eqns to satisfy unit correction // real JMCU_9.hINH dimensionless; // JMCU_9.hINH=4.0; real JMCU_9.CMT(time) micromolar; when(time=time.min) JMCU_9.CMT=0.200; real JIPR_9.JIPR(time) flux; real JIPR_9.kIPR first_order_rate_constant; JIPR_9.kIPR=2000.0; real JIPR_9.k_1 flux; JIPR_9.k_1=6.4; real JIPR_9.k1 first_order_rate_constant; JIPR_9.k1=0.0; real JIPR_9.k2 first_order_rate_constant; JIPR_9.k2=4.0; real JIPR_9.r2 first_order_rate_constant; JIPR_9.r2=250.0; real JIPR_9.r_2 flux; JIPR_9.r_2=0.0; real JIPR_9.r4 first_order_rate_constant; JIPR_9.r4=750.0; real JIPR_9.R1 micromolar; JIPR_9.R1=36.0; real JIPR_9.R3 micromolar; JIPR_9.R3=300.0; real JIPR_9.phi1(time) first_order_rate_constant; real JIPR_9.phi_1(time) flux; real JIPR_9.phi2(time) first_order_rate_constant; real JIPR_9.phi3(time) first_order_rate_constant; real JIPR_9.H(time) dimensionless; when(time=time.min) JIPR_9.H=0.200; real JIPR_9.g_beta first_order_rate_constant; JIPR_9.g_beta=1500.0; real JIPR_9.h_beta dimensionless; JIPR_9.h_beta=4.0; real JIPR_9.g_alpha first_order_rate_constant; JIPR_9.g_alpha=0.85; real JIPR_9.K_beta micromolar; JIPR_9.K_beta=0.35; real JIPR_9.alpha_zeta first_order_rate_constant; real JIPR_9.beta_zeta(time) first_order_rate_constant; real JIPR_9.zeta(time) dimensionless; when(time=time.min) JIPR_9.zeta=0.300; real JIPR_9.K_phi3_act micromolar; JIPR_9.K_phi3_act=0.1; real JIPR_9.K_phi3_inh micromolar; JIPR_9.K_phi3_inh=0.5; real JIPR_9.h_phi3_act dimensionless; JIPR_9.h_phi3_act=3.0; real JIPR_9.h_phi3_inh dimensionless; JIPR_9.h_phi3_inh=3.0; real JIPR_9.g_phi3 first_order_rate_constant; JIPR_9.g_phi3=4.5; real JIPR_9.P micromolar; JIPR_9.P=1.0; real JNCX_9.JNCX(time) flux; real JNCX_9.VNCX flux; JNCX_9.VNCX=3.0; real JNCX_9.KNCX micromolar; JNCX_9.KNCX=0.3; real CS1_9.lambda_MT_S1 dimensionless; real CS1_9.lambda_ER_S1 dimensionless; real CS1_9.delta_SPU micromolar_per_coulomb; CS1_9.delta_SPU=18.5; real CS1_9.gamma_S1 dimensionless; CS1_9.gamma_S1=100.0; real CS1_9.gamma_MT dimensionless; CS1_9.gamma_MT=200.0; real CS1_9.gamma_ER dimensionless; CS1_9.gamma_ER=20.0; real CS2_9.lambda_MT_S2 dimensionless; real CS2_9.lambda_ER_S2 dimensionless; real CS2_9.lambda_S1_S2 dimensionless; real CS2_9.gamma_S2 dimensionless; CS2_9.gamma_S2=1.0; real CMT_9.fm(time) dimensionless; real CMT_9.Km micromolar; CMT_9.Km=0.01; real CMT_9.Bm micromolar; CMT_9.Bm=100.0; real model_parameters_10.KO micromolar; model_parameters_10.KO=5400.0; real model_parameters_10.Ki micromolar; model_parameters_10.Ki=145E3; real model_parameters_10.n_PU_base dimensionless; model_parameters_10.n_PU_base=50.0; // // I_CaT=(g_CaT*OCaT*(Vm-ECaT)); OCaT:time=(alpha_OT*dT*fT*(1-OCaT)-beta_OT*OCaT); // infinity_dT=(1/(1+exp(k_dT*(Vm-V_dT)))); tau_dT1=A_dT1; dT:time=((infinity_dT-dT)/tau_dT1); // infinity_fT=(1/(1+exp(k_fT*(Vm-V_fT)))); tau_fT1=(A_fT1+(A_fT2-A_fT1)/(1+exp(A_fT3*(Vm-A_fT4)))); fT:time=((infinity_fT-fT)/tau_fT1); // I_CaExt=(gCaExt*(CCy/(KCaExt+CCy))); // I_Kv1_1=(g_Kv1_1*dv1_1*fv1_1*(Vm-EK)); // alpha_dv1_1=(A_dv1_11*(A_dv1_12/(1+exp(A_dv1_13*(Vm-A_dv1_14)))+(1-A_dv1_12))); beta_dv1_1=(A_dv1_11*(A_dv1_12-A_dv1_12/(1+exp(A_dv1_13*(Vm-A_dv1_14))))); A_dv1_11=(A_dv1_11a+(A_dv1_11b-A_dv1_11a)/(1+exp(A_dv1_11c*(Vm-A_dv1_11d)))); dv1_1:time=(alpha_dv1_1*(1-dv1_1)-beta_dv1_1*dv1_1); // alpha_fv1_1=(A_fv1_11*(A_fv1_12/(1+exp(A_fv1_13*(Vm-A_fv1_14)))+(1-A_fv1_12))); beta_fv1_1=(A_fv1_11*(A_fv1_12-A_fv1_12/(1+exp(A_fv1_13*(Vm-A_fv1_14))))); fv1_1:time=(alpha_fv1_1*(1-fv1_1)-beta_fv1_1*fv1_1); // I_KERG=(g_KERG*dERG*(Vm-EK)); // alpha_dERG=(A_dERG1*(A_dERG2/(1+exp(A_dERG3*(Vm-A_dERG4)))+(1-A_dERG2))); beta_dERG=(A_dERG1*(A_dERG2-A_dERG2/(1+exp(A_dERG3*(Vm-A_dERG4))))); dERG:time=(alpha_dERG*(1-dERG)-beta_dERG*dERG); // I_KB=(g_KB*(Vm-EK)); // I_Na.I_Na=(g_Na*dNa*fNa*(Vm-ENa_Cy)); // alpha_dNa=(A_dNa1*(A_dNa2/(1+exp(A_dNa3*(Vm-A_dNa4)))+(1-A_dNa2))); beta_dNa=(A_dNa1*(A_dNa2-A_dNa2/(1+exp(A_dNa3*(Vm-A_dNa4))))); A_dNa1=(A_dNa1a+(A_dNa1b-A_dNa1a)/(1+exp(A_dNa1c*(Vm-A_dNa1d)))); dNa:time=(alpha_dNa*(1-dNa)-beta_dNa*dNa); // alpha_fNa=(A_fNa1*(A_fNa2/(1+exp(A_fNa3*(Vm-A_fNa4)))+(1-A_fNa2))); beta_fNa=(A_fNa1*(A_fNa2-A_fNa2/(1+exp(A_fNa3*(Vm-A_fNa4))))); A_fNa1=(A_fNa1a+(A_fNa1b-A_fNa1a)/(1+exp(A_fNa1c*(Vm-A_fNa1d)))); fNa:time=(alpha_fNa*(1-fNa)-beta_fNa*fNa); // I_L=(g_L*(Vm-EL)); // Vm:time=((-1)*(1/Cm)*(I_ion_Cy+I_ion_PU1+I_ion_PU2+I_ion_PU3+I_ion_PU4+I_ion_PU5+I_ion_PU6+I_ion_PU7+I_ion_PU8+I_ion_PU9+I_ion_PU10)); I_ion_Cy=(I_CaT+I_CaExt+I_Kv1_1+I_KERG+I_KB+I_Na.I_Na+I_L); // ECa_Cy=(Nerst_potentials.R*Nerst_potentials.T/(2*Nerst_potentials.F)*ln(Nerst_potentials.CO/CCy)); ENa_Cy=(Nerst_potentials.R*Nerst_potentials.T/Nerst_potentials.F*ln(Nerst_potentials.NO/Ni)); EK=(Nerst_potentials.R*Nerst_potentials.T/Nerst_potentials.F*ln(Nerst_potentials.KO/Nerst_potentials.Ki)); // CCy:time=(JCy-(lambda_S2_Cy*was_JS2Cy1+lambda_S2_Cy*was_JS2Cy2+lambda_S2_Cy*was_JS2Cy3+lambda_S2_Cy*was_JS2Cy4+lambda_S2_Cy*was_JS2Cy5+lambda_S2_Cy*was_JS2Cy6+lambda_S2_Cy*was_JS2Cy7+lambda_S2_Cy*was_JS2Cy8+lambda_S2_Cy*was_JS2Cy9+lambda_S2_Cy*was_JS2Cy10+delta_SCy/2*(I_CaT+I_CaExt))); // was_JS2Cy1=(mu_S2Cy1*(CCy-CS21)); was_JS2Cy2=(mu_S2Cy2*(CCy-CS22)); was_JS2Cy3=(mu_S2Cy3*(CCy-CS23)); was_JS2Cy4=(mu_S2Cy4*(CCy-CS24)); was_JS2Cy5=(mu_S2Cy5*(CCy-CS25)); was_JS2Cy6=(mu_S2Cy6*(CCy-CS26)); was_JS2Cy7=(mu_S2Cy7*(CCy-CS27)); was_JS2Cy8=(mu_S2Cy8*(CCy-CS28)); was_JS2Cy9=(mu_S2Cy9*(CCy-CS29)); was_JS2Cy10=(mu_S2Cy10*(CCy-CS210)); // was_JS1S2_1=(mu_S1S2_1*(CS21-CS11)); was_JS1S2_2=(mu_S1S2_2*(CS22-CS12)); was_JS1S2_3=(mu_S1S2_3*(CS23-CS13)); was_JS1S2_4=(mu_S1S2_4*(CS24-CS14)); was_JS1S2_5=(mu_S1S2_5*(CS25-CS15)); was_JS1S2_6=(mu_S1S2_6*(CS26-CS16)); was_JS1S2_7=(mu_S1S2_7*(CS27-CS17)); was_JS1S2_8=(mu_S1S2_8*(CS28-CS18)); was_JS1S2_9=(mu_S1S2_9*(CS29-CS19)); was_JS1S2_10=(mu_S1S2_10*(CS210-CS110)); mu_S1S2_1=(mu_A+(mu_B-mu_A)*((1-1)/(n_PU-1))); mu_S1S2_2=(mu_A+(mu_B-mu_A)*((2-1)/(n_PU-1))); mu_S1S2_3=(mu_A+(mu_B-mu_A)*((3-1)/(n_PU-1))); mu_S1S2_4=(mu_A+(mu_B-mu_A)*((4-1)/(n_PU-1))); mu_S1S2_5=(mu_A+(mu_B-mu_A)*((5-1)/(n_PU-1))); mu_S1S2_6=(mu_A+(mu_B-mu_A)*((6-1)/(n_PU-1))); mu_S1S2_7=(mu_A+(mu_B-mu_A)*((7-1)/(n_PU-1))); mu_S1S2_8=(mu_A+(mu_B-mu_A)*((8-1)/(n_PU-1))); mu_S1S2_9=(mu_A+(mu_B-mu_A)*((9-1)/(n_PU-1))); mu_S1S2_10=(mu_A+(mu_B-mu_A)*((10-1)/(n_PU-1))); // JCy=(mu_Cy*(C_infinity-CCy)); // // // // // // // // // // // // PU_membrane.I_iCa=(PU_membrane.alpha_scale*(PU_membrane.I_Ca+PU_membrane.I_NSCC_Ca+PU_membrane.I_PM)); PU_membrane.I_iNa=(PU_membrane.alpha_scale*(PU_membrane.I_NSCC_Na+PU_membrane.I_Na+PU_membrane.I_NaP)); I_ion_PU1=(PU_membrane.I_iCa+PU_membrane.I_iNa); // PU_membrane.I_Ca=(I_Ca.gCa*(Vm-I_Ca.ECa_PU)); I_Ca.gCa=(I_Ca.gCa_*exp(I_Ca.kCa*Vm)/(1+exp(I_Ca.kVCa*(Vm-I_Ca.VhCa)))); I_Ca.ECa_PU=(I_Ca.R*I_Ca.T/(2*I_Ca.F)*ln(I_Ca.CO/CS11)); // PU_membrane.I_Na=(I_Na_1.gNa*(Vm-I_Na_1.ENa_PU)); I_Na_1.ENa_PU=(I_Ca.R*I_Ca.T/I_Ca.F*ln(I_Na_1.NO/I_Na_1.NS1)); // I_NSCC_Ca.gNSCC_Ca=(I_NSCC_Ca.gNSCC_Ca_*(I_NSCC_Ca.KNSCC^3/(I_NSCC_Ca.KNSCC^3+CS11^3))); PU_membrane.I_NSCC_Ca=(I_NSCC_Ca.gNSCC_Ca*(Vm-I_NSCC_Ca.ENSCC)); // I_NSCC_Na.gNSCC_Na=(I_NSCC_Na.gNSCC_Na_*(I_NSCC_Ca.KNSCC^3/(I_NSCC_Ca.KNSCC^3+CS11^3))); PU_membrane.I_NSCC_Na=(I_NSCC_Na.gNSCC_Na*(Vm-I_NSCC_Ca.ENSCC)); // PU_membrane.I_PM=(I_PM.gPM*(CS11^2/(I_PM.KPM^2+CS11^2))); // PU_membrane.I_NaP=(I_NaP.gNaP*(I_Na_1.NS1^4/(I_NaP.KNaP^4+I_Na_1.NS1^4))*(I_NaP.ENaP-Vm)); // was_JSERCA.JSERCA=(was_JSERCA.VSERCA*(CS11-was_JSERCA.A2*was_JSERCA.CER)/(1+was_JSERCA.A4*CS11+was_JSERCA.A5*was_JSERCA.CER+was_JSERCA.A6*CS11*was_JSERCA.CER)); // JMCU.JMCU=(JMCU.VMCU*(CS21^2/(JMCU.KMCU^2+CS21^2))*JMCU.epsilon_INH); JMCU.epsilon_INH=(JMCU.KINH^4/(JMCU.KINH^4+JMCU.CMT^4)); // JIPR.H:time=(JIPR.phi3*(1-JIPR.H)-JIPR.P*JIPR.phi1*JIPR.phi2/(JIPR.P*JIPR.phi1+JIPR.phi_1)*JIPR.H); JIPR.JIPR=(JIPR.kIPR*(JIPR.P*JIPR.phi1*JIPR.H/(JIPR.P*JIPR.phi1+JIPR.phi_1))^4*(was_JSERCA.CER-CS21)); JIPR.phi1=((JIPR.k1*JIPR.R1+JIPR.r2*CS21)/(JIPR.R1+CS21)); JIPR.phi_1=((JIPR.k_1+JIPR.r_2)*JIPR.R3/(JIPR.R3+CS21)); JIPR.phi2=((JIPR.k2*JIPR.R3+JIPR.r4*CS21)/(JIPR.R3+CS21)); JIPR.phi3=(JIPR.g_phi3*JIPR.zeta/((1+(JIPR.K_phi3_act/CS21)^JIPR.h_phi3_act)*(1+(CS21/JIPR.K_phi3_inh)^JIPR.h_phi3_inh))); JIPR.alpha_zeta=JIPR.g_alpha; JIPR.beta_zeta=(JIPR.g_beta/(1+(JIPR.K_beta/CS21)^JIPR.h_beta)); JIPR.zeta:time=(JIPR.alpha_zeta*(1-JIPR.zeta)-JIPR.beta_zeta*JIPR.zeta); // JNCX.JNCX=(JNCX.VNCX*(JMCU.CMT/(JMCU.CMT+JNCX.KNCX))); // CS11:time=(was_JS1S2_1+CS1.lambda_MT_S1*JNCX.JNCX-(CS1.delta_SPU/(PU_membrane.alpha_scale*2)*PU_membrane.I_iCa+CS1.lambda_ER_S1*was_JSERCA.JSERCA)); CS1.lambda_MT_S1=(CS1.gamma_MT/CS1.gamma_S1); CS1.lambda_ER_S1=(CS1.gamma_ER/CS1.gamma_S1); // CS21:time=(was_JS2Cy1+CS2.lambda_ER_S2*JIPR.JIPR-(CS2.lambda_S1_S2*was_JS1S2_1+CS2.lambda_MT_S2*JMCU.JMCU)); CS2.lambda_MT_S2=(CS1.gamma_MT/CS2.gamma_S2); CS2.lambda_ER_S2=(CS1.gamma_ER/CS2.gamma_S2); CS2.lambda_S1_S2=(CS1.gamma_S1/CS2.gamma_S2); // was_JSERCA.CER:time=(was_JSERCA.JSERCA-JIPR.JIPR); // JMCU.CMT:time=(CMT.fm*(JMCU.JMCU-JNCX.JNCX)); CMT.fm=(1/(1+CMT.Km*CMT.Bm/(CMT.Km+JMCU.CMT)^2)); // I_Na_1.NS1:time=((-1)*(CS1.delta_SPU/(1*PU_membrane.alpha_scale))*PU_membrane.I_iNa); // PU_membrane.alpha_scale=(model_parameters_1.n_PU_base/n_PU); // PU_membrane_1.I_iCa=(PU_membrane_1.alpha_scale*(PU_membrane_1.I_Ca+PU_membrane_1.I_NSCC_Ca+PU_membrane_1.I_PM)); PU_membrane_1.I_iNa=(PU_membrane_1.alpha_scale*(PU_membrane_1.I_NSCC_Na+PU_membrane_1.I_Na+PU_membrane_1.I_NaP)); I_ion_PU2=(PU_membrane_1.I_iCa+PU_membrane_1.I_iNa); // PU_membrane_1.I_Ca=(I_Ca_1.gCa*(Vm-I_Ca_1.ECa_PU)); I_Ca_1.gCa=(I_Ca_1.gCa_*exp(I_Ca_1.kCa*Vm)/(1+exp(I_Ca_1.kVCa*(Vm-I_Ca_1.VhCa)))); I_Ca_1.ECa_PU=(I_Ca_1.R*I_Ca_1.T/(2*I_Ca_1.F)*ln(I_Ca_1.CO/CS12)); // PU_membrane_1.I_Na=(I_Na_2.gNa*(Vm-I_Na_2.ENa_PU)); I_Na_2.ENa_PU=(I_Ca_1.R*I_Ca_1.T/I_Ca_1.F*ln(I_Na_2.NO/I_Na_2.NS1)); // I_NSCC_Ca_1.gNSCC_Ca=(I_NSCC_Ca_1.gNSCC_Ca_*(I_NSCC_Ca_1.KNSCC^3/(I_NSCC_Ca_1.KNSCC^3+CS12^3))); PU_membrane_1.I_NSCC_Ca=(I_NSCC_Ca_1.gNSCC_Ca*(Vm-I_NSCC_Ca_1.ENSCC)); // I_NSCC_Na_1.gNSCC_Na=(I_NSCC_Na_1.gNSCC_Na_*(I_NSCC_Ca_1.KNSCC^3/(I_NSCC_Ca_1.KNSCC^3+CS12^3))); PU_membrane_1.I_NSCC_Na=(I_NSCC_Na_1.gNSCC_Na*(Vm-I_NSCC_Ca_1.ENSCC)); // PU_membrane_1.I_PM=(I_PM_1.gPM*(CS12^2/(I_PM_1.KPM^2+CS12^2))); // PU_membrane_1.I_NaP=(I_NaP_1.gNaP*(I_Na_2.NS1^4/(I_NaP_1.KNaP^4+I_Na_2.NS1^4))*(I_NaP_1.ENaP-Vm)); // was_JSERCA_1.JSERCA=(was_JSERCA_1.VSERCA*(CS12-was_JSERCA_1.A2*was_JSERCA_1.CER)/(1+was_JSERCA_1.A4*CS12+was_JSERCA_1.A5*was_JSERCA_1.CER+was_JSERCA_1.A6*CS12*was_JSERCA_1.CER)); // JMCU_1.JMCU=(JMCU_1.VMCU*(CS22^2/(JMCU_1.KMCU^2+CS22^2))*JMCU_1.epsilon_INH); JMCU_1.epsilon_INH=(JMCU_1.KINH^4/(JMCU_1.KINH^4+JMCU_1.CMT^4)); // JIPR_1.H:time=(JIPR_1.phi3*(1-JIPR_1.H)-JIPR_1.P*JIPR_1.phi1*JIPR_1.phi2/(JIPR_1.P*JIPR_1.phi1+JIPR_1.phi_1)*JIPR_1.H); JIPR_1.JIPR=(JIPR_1.kIPR*(JIPR_1.P*JIPR_1.phi1*JIPR_1.H/(JIPR_1.P*JIPR_1.phi1+JIPR_1.phi_1))^4*(was_JSERCA_1.CER-CS22)); JIPR_1.phi1=((JIPR_1.k1*JIPR_1.R1+JIPR_1.r2*CS22)/(JIPR_1.R1+CS22)); JIPR_1.phi_1=((JIPR_1.k_1+JIPR_1.r_2)*JIPR_1.R3/(JIPR_1.R3+CS22)); JIPR_1.phi2=((JIPR_1.k2*JIPR_1.R3+JIPR_1.r4*CS22)/(JIPR_1.R3+CS22)); JIPR_1.phi3=(JIPR_1.g_phi3*JIPR_1.zeta/((1+(JIPR_1.K_phi3_act/CS22)^JIPR_1.h_phi3_act)*(1+(CS22/JIPR_1.K_phi3_inh)^JIPR_1.h_phi3_inh))); JIPR_1.alpha_zeta=JIPR_1.g_alpha; JIPR_1.beta_zeta=(JIPR_1.g_beta/(1+(JIPR_1.K_beta/CS22)^JIPR_1.h_beta)); JIPR_1.zeta:time=(JIPR_1.alpha_zeta*(1-JIPR_1.zeta)-JIPR_1.beta_zeta*JIPR_1.zeta); // JNCX_1.JNCX=(JNCX_1.VNCX*(JMCU_1.CMT/(JMCU_1.CMT+JNCX_1.KNCX))); // CS12:time=(was_JS1S2_2+CS1_1.lambda_MT_S1*JNCX_1.JNCX-(CS1_1.delta_SPU/(PU_membrane_1.alpha_scale*2)*PU_membrane_1.I_iCa+CS1_1.lambda_ER_S1*was_JSERCA_1.JSERCA)); CS1_1.lambda_MT_S1=(CS1_1.gamma_MT/CS1_1.gamma_S1); CS1_1.lambda_ER_S1=(CS1_1.gamma_ER/CS1_1.gamma_S1); // CS22:time=(was_JS2Cy2+CS2_1.lambda_ER_S2*JIPR_1.JIPR-(CS2_1.lambda_S1_S2*was_JS1S2_2+CS2_1.lambda_MT_S2*JMCU_1.JMCU)); CS2_1.lambda_MT_S2=(CS1_1.gamma_MT/CS2_1.gamma_S2); CS2_1.lambda_ER_S2=(CS1_1.gamma_ER/CS2_1.gamma_S2); CS2_1.lambda_S1_S2=(CS1_1.gamma_S1/CS2_1.gamma_S2); // was_JSERCA_1.CER:time=(was_JSERCA_1.JSERCA-JIPR_1.JIPR); // JMCU_1.CMT:time=(CMT_1.fm*(JMCU_1.JMCU-JNCX_1.JNCX)); CMT_1.fm=(1/(1+CMT_1.Km*CMT_1.Bm/(CMT_1.Km+JMCU_1.CMT)^2)); // I_Na_2.NS1:time=((-1)*(CS1_1.delta_SPU/(1*PU_membrane_1.alpha_scale))*PU_membrane_1.I_iNa); // PU_membrane_1.alpha_scale=(model_parameters_2.n_PU_base/n_PU); // PU_membrane_2.I_iCa=(PU_membrane_2.alpha_scale*(PU_membrane_2.I_Ca+PU_membrane_2.I_NSCC_Ca+PU_membrane_2.I_PM)); PU_membrane_2.I_iNa=(PU_membrane_2.alpha_scale*(PU_membrane_2.I_NSCC_Na+PU_membrane_2.I_Na+PU_membrane_2.I_NaP)); I_ion_PU3=(PU_membrane_2.I_iCa+PU_membrane_2.I_iNa); // PU_membrane_2.I_Ca=(I_Ca_2.gCa*(Vm-I_Ca_2.ECa_PU)); I_Ca_2.gCa=(I_Ca_2.gCa_*exp(I_Ca_2.kCa*Vm)/(1+exp(I_Ca_2.kVCa*(Vm-I_Ca_2.VhCa)))); I_Ca_2.ECa_PU=(I_Ca_2.R*I_Ca_2.T/(2*I_Ca_2.F)*ln(I_Ca_2.CO/CS13)); // PU_membrane_2.I_Na=(I_Na_3.gNa*(Vm-I_Na_3.ENa_PU)); I_Na_3.ENa_PU=(I_Ca_2.R*I_Ca_2.T/I_Ca_2.F*ln(I_Na_3.NO/I_Na_3.NS1)); // I_NSCC_Ca_2.gNSCC_Ca=(I_NSCC_Ca_2.gNSCC_Ca_*(I_NSCC_Ca_2.KNSCC^3/(I_NSCC_Ca_2.KNSCC^3+CS13^3))); PU_membrane_2.I_NSCC_Ca=(I_NSCC_Ca_2.gNSCC_Ca*(Vm-I_NSCC_Ca_2.ENSCC)); // I_NSCC_Na_2.gNSCC_Na=(I_NSCC_Na_2.gNSCC_Na_*(I_NSCC_Ca_2.KNSCC^3/(I_NSCC_Ca_2.KNSCC^3+CS13^3))); PU_membrane_2.I_NSCC_Na=(I_NSCC_Na_2.gNSCC_Na*(Vm-I_NSCC_Ca_2.ENSCC)); // PU_membrane_2.I_PM=(I_PM_2.gPM*(CS13^2/(I_PM_2.KPM^2+CS13^2))); // PU_membrane_2.I_NaP=(I_NaP_2.gNaP*(I_Na_3.NS1^4/(I_NaP_2.KNaP^4+I_Na_3.NS1^4))*(I_NaP_2.ENaP-Vm)); // was_JSERCA_2.JSERCA=(was_JSERCA_2.VSERCA*(CS13-was_JSERCA_2.A2*was_JSERCA_2.CER)/(1+was_JSERCA_2.A4*CS13+was_JSERCA_2.A5*was_JSERCA_2.CER+was_JSERCA_2.A6*CS13*was_JSERCA_2.CER)); // JMCU_2.JMCU=(JMCU_2.VMCU*(CS23^2/(JMCU_2.KMCU^2+CS23^2))*JMCU_2.epsilon_INH); JMCU_2.epsilon_INH=(JMCU_2.KINH^4/(JMCU_2.KINH^4+JMCU_2.CMT^4)); // JIPR_2.H:time=(JIPR_2.phi3*(1-JIPR_2.H)-JIPR_2.P*JIPR_2.phi1*JIPR_2.phi2/(JIPR_2.P*JIPR_2.phi1+JIPR_2.phi_1)*JIPR_2.H); JIPR_2.JIPR=(JIPR_2.kIPR*(JIPR_2.P*JIPR_2.phi1*JIPR_2.H/(JIPR_2.P*JIPR_2.phi1+JIPR_2.phi_1))^4*(was_JSERCA_2.CER-CS23)); JIPR_2.phi1=((JIPR_2.k1*JIPR_2.R1+JIPR_2.r2*CS23)/(JIPR_2.R1+CS23)); JIPR_2.phi_1=((JIPR_2.k_1+JIPR_2.r_2)*JIPR_2.R3/(JIPR_2.R3+CS23)); JIPR_2.phi2=((JIPR_2.k2*JIPR_2.R3+JIPR_2.r4*CS23)/(JIPR_2.R3+CS23)); JIPR_2.phi3=(JIPR_2.g_phi3*JIPR_2.zeta/((1+(JIPR_2.K_phi3_act/CS23)^JIPR_2.h_phi3_act)*(1+(CS23/JIPR_2.K_phi3_inh)^JIPR_2.h_phi3_inh))); JIPR_2.alpha_zeta=JIPR_2.g_alpha; JIPR_2.beta_zeta=(JIPR_2.g_beta/(1+(JIPR_2.K_beta/CS23)^JIPR_2.h_beta)); JIPR_2.zeta:time=(JIPR_2.alpha_zeta*(1-JIPR_2.zeta)-JIPR_2.beta_zeta*JIPR_2.zeta); // JNCX_2.JNCX=(JNCX_2.VNCX*(JMCU_2.CMT/(JMCU_2.CMT+JNCX_2.KNCX))); // CS13:time=(was_JS1S2_3+CS1_2.lambda_MT_S1*JNCX_2.JNCX-(CS1_2.delta_SPU/(PU_membrane_2.alpha_scale*2)*PU_membrane_2.I_iCa+CS1_2.lambda_ER_S1*was_JSERCA_2.JSERCA)); CS1_2.lambda_MT_S1=(CS1_2.gamma_MT/CS1_2.gamma_S1); CS1_2.lambda_ER_S1=(CS1_2.gamma_ER/CS1_2.gamma_S1); // CS23:time=(was_JS2Cy3+CS2_2.lambda_ER_S2*JIPR_2.JIPR-(CS2_2.lambda_S1_S2*was_JS1S2_3+CS2_2.lambda_MT_S2*JMCU_2.JMCU)); CS2_2.lambda_MT_S2=(CS1_2.gamma_MT/CS2_2.gamma_S2); CS2_2.lambda_ER_S2=(CS1_2.gamma_ER/CS2_2.gamma_S2); CS2_2.lambda_S1_S2=(CS1_2.gamma_S1/CS2_2.gamma_S2); // was_JSERCA_2.CER:time=(was_JSERCA_2.JSERCA-JIPR_2.JIPR); // JMCU_2.CMT:time=(CMT_2.fm*(JMCU_2.JMCU-JNCX_2.JNCX)); CMT_2.fm=(1/(1+CMT_2.Km*CMT_2.Bm/(CMT_2.Km+JMCU_2.CMT)^2)); // I_Na_3.NS1:time=((-1)*(CS1_2.delta_SPU/(1*PU_membrane_2.alpha_scale))*PU_membrane_2.I_iNa); // PU_membrane_2.alpha_scale=(model_parameters_3.n_PU_base/n_PU); // PU_membrane_3.I_iCa=(PU_membrane_3.alpha_scale*(PU_membrane_3.I_Ca+PU_membrane_3.I_NSCC_Ca+PU_membrane_3.I_PM)); PU_membrane_3.I_iNa=(PU_membrane_3.alpha_scale*(PU_membrane_3.I_NSCC_Na+PU_membrane_3.I_Na+PU_membrane_3.I_NaP)); I_ion_PU4=(PU_membrane_3.I_iCa+PU_membrane_3.I_iNa); // PU_membrane_3.I_Ca=(I_Ca_3.gCa*(Vm-I_Ca_3.ECa_PU)); I_Ca_3.gCa=(I_Ca_3.gCa_*exp(I_Ca_3.kCa*Vm)/(1+exp(I_Ca_3.kVCa*(Vm-I_Ca_3.VhCa)))); I_Ca_3.ECa_PU=(I_Ca_3.R*I_Ca_3.T/(2*I_Ca_3.F)*ln(I_Ca_3.CO/CS14)); // PU_membrane_3.I_Na=(I_Na_4.gNa*(Vm-I_Na_4.ENa_PU)); I_Na_4.ENa_PU=(I_Ca_3.R*I_Ca_3.T/I_Ca_3.F*ln(I_Na_4.NO/I_Na_4.NS1)); // I_NSCC_Ca_3.gNSCC_Ca=(I_NSCC_Ca_3.gNSCC_Ca_*(I_NSCC_Ca_3.KNSCC^3/(I_NSCC_Ca_3.KNSCC^3+CS14^3))); PU_membrane_3.I_NSCC_Ca=(I_NSCC_Ca_3.gNSCC_Ca*(Vm-I_NSCC_Ca_3.ENSCC)); // I_NSCC_Na_3.gNSCC_Na=(I_NSCC_Na_3.gNSCC_Na_*(I_NSCC_Ca_3.KNSCC^3/(I_NSCC_Ca_3.KNSCC^3+CS14^3))); PU_membrane_3.I_NSCC_Na=(I_NSCC_Na_3.gNSCC_Na*(Vm-I_NSCC_Ca_3.ENSCC)); // PU_membrane_3.I_PM=(I_PM_3.gPM*(CS14^2/(I_PM_3.KPM^2+CS14^2))); // PU_membrane_3.I_NaP=(I_NaP_3.gNaP*(I_Na_4.NS1^4/(I_NaP_3.KNaP^4+I_Na_4.NS1^4))*(I_NaP_3.ENaP-Vm)); // was_JSERCA_3.JSERCA=(was_JSERCA_3.VSERCA*(CS14-was_JSERCA_3.A2*was_JSERCA_3.CER)/(1+was_JSERCA_3.A4*CS14+was_JSERCA_3.A5*was_JSERCA_3.CER+was_JSERCA_3.A6*CS14*was_JSERCA_3.CER)); // JMCU_3.JMCU=(JMCU_3.VMCU*(CS24^2/(JMCU_3.KMCU^2+CS24^2))*JMCU_3.epsilon_INH); JMCU_3.epsilon_INH=(JMCU_3.KINH^4/(JMCU_3.KINH^4+JMCU_3.CMT^4)); // JIPR_3.H:time=(JIPR_3.phi3*(1-JIPR_3.H)-JIPR_3.P*JIPR_3.phi1*JIPR_3.phi2/(JIPR_3.P*JIPR_3.phi1+JIPR_3.phi_1)*JIPR_3.H); JIPR_3.JIPR=(JIPR_3.kIPR*(JIPR_3.P*JIPR_3.phi1*JIPR_3.H/(JIPR_3.P*JIPR_3.phi1+JIPR_3.phi_1))^4*(was_JSERCA_3.CER-CS24)); JIPR_3.phi1=((JIPR_3.k1*JIPR_3.R1+JIPR_3.r2*CS24)/(JIPR_3.R1+CS24)); JIPR_3.phi_1=((JIPR_3.k_1+JIPR_3.r_2)*JIPR_3.R3/(JIPR_3.R3+CS24)); JIPR_3.phi2=((JIPR_3.k2*JIPR_3.R3+JIPR_3.r4*CS24)/(JIPR_3.R3+CS24)); JIPR_3.phi3=(JIPR_3.g_phi3*JIPR_3.zeta/((1+(JIPR_3.K_phi3_act/CS24)^JIPR_3.h_phi3_act)*(1+(CS24/JIPR_3.K_phi3_inh)^JIPR_3.h_phi3_inh))); JIPR_3.alpha_zeta=JIPR_3.g_alpha; JIPR_3.beta_zeta=(JIPR_3.g_beta/(1+(JIPR_3.K_beta/CS24)^JIPR_3.h_beta)); JIPR_3.zeta:time=(JIPR_3.alpha_zeta*(1-JIPR_3.zeta)-JIPR_3.beta_zeta*JIPR_3.zeta); // JNCX_3.JNCX=(JNCX_3.VNCX*(JMCU_3.CMT/(JMCU_3.CMT+JNCX_3.KNCX))); // CS14:time=(was_JS1S2_4+CS1_3.lambda_MT_S1*JNCX_3.JNCX-(CS1_3.delta_SPU/(PU_membrane_3.alpha_scale*2)*PU_membrane_3.I_iCa+CS1_3.lambda_ER_S1*was_JSERCA_3.JSERCA)); CS1_3.lambda_MT_S1=(CS1_3.gamma_MT/CS1_3.gamma_S1); CS1_3.lambda_ER_S1=(CS1_3.gamma_ER/CS1_3.gamma_S1); // CS24:time=(was_JS2Cy4+CS2_3.lambda_ER_S2*JIPR_3.JIPR-(CS2_3.lambda_S1_S2*was_JS1S2_4+CS2_3.lambda_MT_S2*JMCU_3.JMCU)); CS2_3.lambda_MT_S2=(CS1_3.gamma_MT/CS2_3.gamma_S2); CS2_3.lambda_ER_S2=(CS1_3.gamma_ER/CS2_3.gamma_S2); CS2_3.lambda_S1_S2=(CS1_3.gamma_S1/CS2_3.gamma_S2); // was_JSERCA_3.CER:time=(was_JSERCA_3.JSERCA-JIPR_3.JIPR); // JMCU_3.CMT:time=(CMT_3.fm*(JMCU_3.JMCU-JNCX_3.JNCX)); CMT_3.fm=(1/(1+CMT_3.Km*CMT_3.Bm/(CMT_3.Km+JMCU_3.CMT)^2)); // I_Na_4.NS1:time=((-1)*(CS1_3.delta_SPU/(1*PU_membrane_3.alpha_scale))*PU_membrane_3.I_iNa); // PU_membrane_3.alpha_scale=(model_parameters_4.n_PU_base/n_PU); // PU_membrane_4.I_iCa=(PU_membrane_4.alpha_scale*(PU_membrane_4.I_Ca+PU_membrane_4.I_NSCC_Ca+PU_membrane_4.I_PM)); PU_membrane_4.I_iNa=(PU_membrane_4.alpha_scale*(PU_membrane_4.I_NSCC_Na+PU_membrane_4.I_Na+PU_membrane_4.I_NaP)); I_ion_PU5=(PU_membrane_4.I_iCa+PU_membrane_4.I_iNa); // PU_membrane_4.I_Ca=(I_Ca_4.gCa*(Vm-I_Ca_4.ECa_PU)); I_Ca_4.gCa=(I_Ca_4.gCa_*exp(I_Ca_4.kCa*Vm)/(1+exp(I_Ca_4.kVCa*(Vm-I_Ca_4.VhCa)))); I_Ca_4.ECa_PU=(I_Ca_4.R*I_Ca_4.T/(2*I_Ca_4.F)*ln(I_Ca_4.CO/CS15)); // PU_membrane_4.I_Na=(I_Na_5.gNa*(Vm-I_Na_5.ENa_PU)); I_Na_5.ENa_PU=(I_Ca_4.R*I_Ca_4.T/I_Ca_4.F*ln(I_Na_5.NO/I_Na_5.NS1)); // I_NSCC_Ca_4.gNSCC_Ca=(I_NSCC_Ca_4.gNSCC_Ca_*(I_NSCC_Ca_4.KNSCC^3/(I_NSCC_Ca_4.KNSCC^3+CS15^3))); PU_membrane_4.I_NSCC_Ca=(I_NSCC_Ca_4.gNSCC_Ca*(Vm-I_NSCC_Ca_4.ENSCC)); // I_NSCC_Na_4.gNSCC_Na=(I_NSCC_Na_4.gNSCC_Na_*(I_NSCC_Ca_4.KNSCC^3/(I_NSCC_Ca_4.KNSCC^3+CS15^3))); PU_membrane_4.I_NSCC_Na=(I_NSCC_Na_4.gNSCC_Na*(Vm-I_NSCC_Ca_4.ENSCC)); // PU_membrane_4.I_PM=(I_PM_4.gPM*(CS15^2/(I_PM_4.KPM^2+CS15^2))); // PU_membrane_4.I_NaP=(I_NaP_4.gNaP*(I_Na_5.NS1^4/(I_NaP_4.KNaP^4+I_Na_5.NS1^4))*(I_NaP_4.ENaP-Vm)); // was_JSERCA_4.JSERCA=(was_JSERCA_4.VSERCA*(CS15-was_JSERCA_4.A2*was_JSERCA_4.CER)/(1+was_JSERCA_4.A4*CS15+was_JSERCA_4.A5*was_JSERCA_4.CER+was_JSERCA_4.A6*CS15*was_JSERCA_4.CER)); // JMCU_4.JMCU=(JMCU_4.VMCU*(CS25^2/(JMCU_4.KMCU^2+CS25^2))*JMCU_4.epsilon_INH); JMCU_4.epsilon_INH=(JMCU_4.KINH^4/(JMCU_4.KINH^4+JMCU_4.CMT^4)); // JIPR_4.H:time=(JIPR_4.phi3*(1-JIPR_4.H)-JIPR_4.P*JIPR_4.phi1*JIPR_4.phi2/(JIPR_4.P*JIPR_4.phi1+JIPR_4.phi_1)*JIPR_4.H); JIPR_4.JIPR=(JIPR_4.kIPR*(JIPR_4.P*JIPR_4.phi1*JIPR_4.H/(JIPR_4.P*JIPR_4.phi1+JIPR_4.phi_1))^4*(was_JSERCA_4.CER-CS25)); JIPR_4.phi1=((JIPR_4.k1*JIPR_4.R1+JIPR_4.r2*CS25)/(JIPR_4.R1+CS25)); JIPR_4.phi_1=((JIPR_4.k_1+JIPR_4.r_2)*JIPR_4.R3/(JIPR_4.R3+CS25)); JIPR_4.phi2=((JIPR_4.k2*JIPR_4.R3+JIPR_4.r4*CS25)/(JIPR_4.R3+CS25)); JIPR_4.phi3=(JIPR_4.g_phi3*JIPR_4.zeta/((1+(JIPR_4.K_phi3_act/CS25)^JIPR_4.h_phi3_act)*(1+(CS25/JIPR_4.K_phi3_inh)^JIPR_4.h_phi3_inh))); JIPR_4.alpha_zeta=JIPR_4.g_alpha; JIPR_4.beta_zeta=(JIPR_4.g_beta/(1+(JIPR_4.K_beta/CS25)^JIPR_4.h_beta)); JIPR_4.zeta:time=(JIPR_4.alpha_zeta*(1-JIPR_4.zeta)-JIPR_4.beta_zeta*JIPR_4.zeta); // JNCX_4.JNCX=(JNCX_4.VNCX*(JMCU_4.CMT/(JMCU_4.CMT+JNCX_4.KNCX))); // CS15:time=(was_JS1S2_5+CS1_4.lambda_MT_S1*JNCX_4.JNCX-(CS1_4.delta_SPU/(PU_membrane_4.alpha_scale*2)*PU_membrane_4.I_iCa+CS1_4.lambda_ER_S1*was_JSERCA_4.JSERCA)); CS1_4.lambda_MT_S1=(CS1_4.gamma_MT/CS1_4.gamma_S1); CS1_4.lambda_ER_S1=(CS1_4.gamma_ER/CS1_4.gamma_S1); // CS25:time=(was_JS2Cy5+CS2_4.lambda_ER_S2*JIPR_4.JIPR-(CS2_4.lambda_S1_S2*was_JS1S2_5+CS2_4.lambda_MT_S2*JMCU_4.JMCU)); CS2_4.lambda_MT_S2=(CS1_4.gamma_MT/CS2_4.gamma_S2); CS2_4.lambda_ER_S2=(CS1_4.gamma_ER/CS2_4.gamma_S2); CS2_4.lambda_S1_S2=(CS1_4.gamma_S1/CS2_4.gamma_S2); // was_JSERCA_4.CER:time=(was_JSERCA_4.JSERCA-JIPR_4.JIPR); // JMCU_4.CMT:time=(CMT_4.fm*(JMCU_4.JMCU-JNCX_4.JNCX)); CMT_4.fm=(1/(1+CMT_4.Km*CMT_4.Bm/(CMT_4.Km+JMCU_4.CMT)^2)); // I_Na_5.NS1:time=((-1)*(CS1_4.delta_SPU/(1*PU_membrane_4.alpha_scale))*PU_membrane_4.I_iNa); // PU_membrane_4.alpha_scale=(model_parameters_5.n_PU_base/n_PU); // PU_membrane_5.I_iCa=(PU_membrane_5.alpha_scale*(PU_membrane_5.I_Ca+PU_membrane_5.I_NSCC_Ca+PU_membrane_5.I_PM)); PU_membrane_5.I_iNa=(PU_membrane_5.alpha_scale*(PU_membrane_5.I_NSCC_Na+PU_membrane_5.I_Na+PU_membrane_5.I_NaP)); I_ion_PU6=(PU_membrane_5.I_iCa+PU_membrane_5.I_iNa); // PU_membrane_5.I_Ca=(I_Ca_5.gCa*(Vm-I_Ca_5.ECa_PU)); I_Ca_5.gCa=(I_Ca_5.gCa_*exp(I_Ca_5.kCa*Vm)/(1+exp(I_Ca_5.kVCa*(Vm-I_Ca_5.VhCa)))); I_Ca_5.ECa_PU=(I_Ca_5.R*I_Ca_5.T/(2*I_Ca_5.F)*ln(I_Ca_5.CO/CS16)); // PU_membrane_5.I_Na=(I_Na_6.gNa*(Vm-I_Na_6.ENa_PU)); I_Na_6.ENa_PU=(I_Ca_5.R*I_Ca_5.T/I_Ca_5.F*ln(I_Na_6.NO/I_Na_6.NS1)); // I_NSCC_Ca_5.gNSCC_Ca=(I_NSCC_Ca_5.gNSCC_Ca_*(I_NSCC_Ca_5.KNSCC^3/(I_NSCC_Ca_5.KNSCC^3+CS16^3))); PU_membrane_5.I_NSCC_Ca=(I_NSCC_Ca_5.gNSCC_Ca*(Vm-I_NSCC_Ca_5.ENSCC)); // I_NSCC_Na_5.gNSCC_Na=(I_NSCC_Na_5.gNSCC_Na_*(I_NSCC_Ca_5.KNSCC^3/(I_NSCC_Ca_5.KNSCC^3+CS16^3))); PU_membrane_5.I_NSCC_Na=(I_NSCC_Na_5.gNSCC_Na*(Vm-I_NSCC_Ca_5.ENSCC)); // PU_membrane_5.I_PM=(I_PM_5.gPM*(CS16^2/(I_PM_5.KPM^2+CS16^2))); // PU_membrane_5.I_NaP=(I_NaP_5.gNaP*(I_Na_6.NS1^4/(I_NaP_5.KNaP^4+I_Na_6.NS1^4))*(I_NaP_5.ENaP-Vm)); // was_JSERCA_5.JSERCA=(was_JSERCA_5.VSERCA*(CS16-was_JSERCA_5.A2*was_JSERCA_5.CER)/(1+was_JSERCA_5.A4*CS16+was_JSERCA_5.A5*was_JSERCA_5.CER+was_JSERCA_5.A6*CS16*was_JSERCA_5.CER)); // JMCU_5.JMCU=(JMCU_5.VMCU*(CS26^2/(JMCU_5.KMCU^2+CS26^2))*JMCU_5.epsilon_INH); JMCU_5.epsilon_INH=(JMCU_5.KINH^4/(JMCU_5.KINH^4+JMCU_5.CMT^4)); // JIPR_5.H:time=(JIPR_5.phi3*(1-JIPR_5.H)-JIPR_5.P*JIPR_5.phi1*JIPR_5.phi2/(JIPR_5.P*JIPR_5.phi1+JIPR_5.phi_1)*JIPR_5.H); JIPR_5.JIPR=(JIPR_5.kIPR*(JIPR_5.P*JIPR_5.phi1*JIPR_5.H/(JIPR_5.P*JIPR_5.phi1+JIPR_5.phi_1))^4*(was_JSERCA_5.CER-CS26)); JIPR_5.phi1=((JIPR_5.k1*JIPR_5.R1+JIPR_5.r2*CS26)/(JIPR_5.R1+CS26)); JIPR_5.phi_1=((JIPR_5.k_1+JIPR_5.r_2)*JIPR_5.R3/(JIPR_5.R3+CS26)); JIPR_5.phi2=((JIPR_5.k2*JIPR_5.R3+JIPR_5.r4*CS26)/(JIPR_5.R3+CS26)); JIPR_5.phi3=(JIPR_5.g_phi3*JIPR_5.zeta/((1+(JIPR_5.K_phi3_act/CS26)^JIPR_5.h_phi3_act)*(1+(CS26/JIPR_5.K_phi3_inh)^JIPR_5.h_phi3_inh))); JIPR_5.alpha_zeta=JIPR_5.g_alpha; JIPR_5.beta_zeta=(JIPR_5.g_beta/(1+(JIPR_5.K_beta/CS26)^JIPR_5.h_beta)); JIPR_5.zeta:time=(JIPR_5.alpha_zeta*(1-JIPR_5.zeta)-JIPR_5.beta_zeta*JIPR_5.zeta); // JNCX_5.JNCX=(JNCX_5.VNCX*(JMCU_5.CMT/(JMCU_5.CMT+JNCX_5.KNCX))); // CS16:time=(was_JS1S2_6+CS1_5.lambda_MT_S1*JNCX_5.JNCX-(CS1_5.delta_SPU/(PU_membrane_5.alpha_scale*2)*PU_membrane_5.I_iCa+CS1_5.lambda_ER_S1*was_JSERCA_5.JSERCA)); CS1_5.lambda_MT_S1=(CS1_5.gamma_MT/CS1_5.gamma_S1); CS1_5.lambda_ER_S1=(CS1_5.gamma_ER/CS1_5.gamma_S1); // CS26:time=(was_JS2Cy6+CS2_5.lambda_ER_S2*JIPR_5.JIPR-(CS2_5.lambda_S1_S2*was_JS1S2_6+CS2_5.lambda_MT_S2*JMCU_5.JMCU)); CS2_5.lambda_MT_S2=(CS1_5.gamma_MT/CS2_5.gamma_S2); CS2_5.lambda_ER_S2=(CS1_5.gamma_ER/CS2_5.gamma_S2); CS2_5.lambda_S1_S2=(CS1_5.gamma_S1/CS2_5.gamma_S2); // was_JSERCA_5.CER:time=(was_JSERCA_5.JSERCA-JIPR_5.JIPR); // JMCU_5.CMT:time=(CMT_5.fm*(JMCU_5.JMCU-JNCX_5.JNCX)); CMT_5.fm=(1/(1+CMT_5.Km*CMT_5.Bm/(CMT_5.Km+JMCU_5.CMT)^2)); // I_Na_6.NS1:time=((-1)*(CS1_5.delta_SPU/(1*PU_membrane_5.alpha_scale))*PU_membrane_5.I_iNa); // PU_membrane_5.alpha_scale=(model_parameters_6.n_PU_base/n_PU); // PU_membrane_6.I_iCa=(PU_membrane_6.alpha_scale*(PU_membrane_6.I_Ca+PU_membrane_6.I_NSCC_Ca+PU_membrane_6.I_PM)); PU_membrane_6.I_iNa=(PU_membrane_6.alpha_scale*(PU_membrane_6.I_NSCC_Na+PU_membrane_6.I_Na+PU_membrane_6.I_NaP)); I_ion_PU7=(PU_membrane_6.I_iCa+PU_membrane_6.I_iNa); // PU_membrane_6.I_Ca=(I_Ca_6.gCa*(Vm-I_Ca_6.ECa_PU)); I_Ca_6.gCa=(I_Ca_6.gCa_*exp(I_Ca_6.kCa*Vm)/(1+exp(I_Ca_6.kVCa*(Vm-I_Ca_6.VhCa)))); I_Ca_6.ECa_PU=(I_Ca_6.R*I_Ca_6.T/(2*I_Ca_6.F)*ln(I_Ca_6.CO/CS17)); // PU_membrane_6.I_Na=(I_Na_7.gNa*(Vm-I_Na_7.ENa_PU)); I_Na_7.ENa_PU=(I_Ca_6.R*I_Ca_6.T/I_Ca_6.F*ln(I_Na_7.NO/I_Na_7.NS1)); // I_NSCC_Ca_6.gNSCC_Ca=(I_NSCC_Ca_6.gNSCC_Ca_*(I_NSCC_Ca_6.KNSCC^3/(I_NSCC_Ca_6.KNSCC^3+CS17^3))); PU_membrane_6.I_NSCC_Ca=(I_NSCC_Ca_6.gNSCC_Ca*(Vm-I_NSCC_Ca_6.ENSCC)); // I_NSCC_Na_6.gNSCC_Na=(I_NSCC_Na_6.gNSCC_Na_*(I_NSCC_Ca_6.KNSCC^3/(I_NSCC_Ca_6.KNSCC^3+CS17^3))); PU_membrane_6.I_NSCC_Na=(I_NSCC_Na_6.gNSCC_Na*(Vm-I_NSCC_Ca_6.ENSCC)); // PU_membrane_6.I_PM=(I_PM_6.gPM*(CS17^2/(I_PM_6.KPM^2+CS17^2))); // PU_membrane_6.I_NaP=(I_NaP_6.gNaP*(I_Na_7.NS1^4/(I_NaP_6.KNaP^4+I_Na_7.NS1^4))*(I_NaP_6.ENaP-Vm)); // was_JSERCA_6.JSERCA=(was_JSERCA_6.VSERCA*(CS17-was_JSERCA_6.A2*was_JSERCA_6.CER)/(1+was_JSERCA_6.A4*CS17+was_JSERCA_6.A5*was_JSERCA_6.CER+was_JSERCA_6.A6*CS17*was_JSERCA_6.CER)); // JMCU_6.JMCU=(JMCU_6.VMCU*(CS27^2/(JMCU_6.KMCU^2+CS27^2))*JMCU_6.epsilon_INH); JMCU_6.epsilon_INH=(JMCU_6.KINH^4/(JMCU_6.KINH^4+JMCU_6.CMT^4)); // JIPR_6.H:time=(JIPR_6.phi3*(1-JIPR_6.H)-JIPR_6.P*JIPR_6.phi1*JIPR_6.phi2/(JIPR_6.P*JIPR_6.phi1+JIPR_6.phi_1)*JIPR_6.H); JIPR_6.JIPR=(JIPR_6.kIPR*(JIPR_6.P*JIPR_6.phi1*JIPR_6.H/(JIPR_6.P*JIPR_6.phi1+JIPR_6.phi_1))^4*(was_JSERCA_6.CER-CS27)); JIPR_6.phi1=((JIPR_6.k1*JIPR_6.R1+JIPR_6.r2*CS27)/(JIPR_6.R1+CS27)); JIPR_6.phi_1=((JIPR_6.k_1+JIPR_6.r_2)*JIPR_6.R3/(JIPR_6.R3+CS27)); JIPR_6.phi2=((JIPR_6.k2*JIPR_6.R3+JIPR_6.r4*CS27)/(JIPR_6.R3+CS27)); JIPR_6.phi3=(JIPR_6.g_phi3*JIPR_6.zeta/((1+(JIPR_6.K_phi3_act/CS27)^JIPR_6.h_phi3_act)*(1+(CS27/JIPR_6.K_phi3_inh)^JIPR_6.h_phi3_inh))); JIPR_6.alpha_zeta=JIPR_6.g_alpha; JIPR_6.beta_zeta=(JIPR_6.g_beta/(1+(JIPR_6.K_beta/CS27)^JIPR_6.h_beta)); JIPR_6.zeta:time=(JIPR_6.alpha_zeta*(1-JIPR_6.zeta)-JIPR_6.beta_zeta*JIPR_6.zeta); // JNCX_6.JNCX=(JNCX_6.VNCX*(JMCU_6.CMT/(JMCU_6.CMT+JNCX_6.KNCX))); // CS17:time=(was_JS1S2_7+CS1_6.lambda_MT_S1*JNCX_6.JNCX-(CS1_6.delta_SPU/(PU_membrane_6.alpha_scale*2)*PU_membrane_6.I_iCa+CS1_6.lambda_ER_S1*was_JSERCA_6.JSERCA)); CS1_6.lambda_MT_S1=(CS1_6.gamma_MT/CS1_6.gamma_S1); CS1_6.lambda_ER_S1=(CS1_6.gamma_ER/CS1_6.gamma_S1); // CS27:time=(was_JS2Cy7+CS2_6.lambda_ER_S2*JIPR_6.JIPR-(CS2_6.lambda_S1_S2*was_JS1S2_7+CS2_6.lambda_MT_S2*JMCU_6.JMCU)); CS2_6.lambda_MT_S2=(CS1_6.gamma_MT/CS2_6.gamma_S2); CS2_6.lambda_ER_S2=(CS1_6.gamma_ER/CS2_6.gamma_S2); CS2_6.lambda_S1_S2=(CS1_6.gamma_S1/CS2_6.gamma_S2); // was_JSERCA_6.CER:time=(was_JSERCA_6.JSERCA-JIPR_6.JIPR); // JMCU_6.CMT:time=(CMT_6.fm*(JMCU_6.JMCU-JNCX_6.JNCX)); CMT_6.fm=(1/(1+CMT_6.Km*CMT_6.Bm/(CMT_6.Km+JMCU_6.CMT)^2)); // I_Na_7.NS1:time=((-1)*(CS1_6.delta_SPU/(1*PU_membrane_6.alpha_scale))*PU_membrane_6.I_iNa); // PU_membrane_6.alpha_scale=(model_parameters_7.n_PU_base/n_PU); // PU_membrane_7.I_iCa=(PU_membrane_7.alpha_scale*(PU_membrane_7.I_Ca+PU_membrane_7.I_NSCC_Ca+PU_membrane_7.I_PM)); PU_membrane_7.I_iNa=(PU_membrane_7.alpha_scale*(PU_membrane_7.I_NSCC_Na+PU_membrane_7.I_Na+PU_membrane_7.I_NaP)); I_ion_PU8=(PU_membrane_7.I_iCa+PU_membrane_7.I_iNa); // PU_membrane_7.I_Ca=(I_Ca_7.gCa*(Vm-I_Ca_7.ECa_PU)); I_Ca_7.gCa=(I_Ca_7.gCa_*exp(I_Ca_7.kCa*Vm)/(1+exp(I_Ca_7.kVCa*(Vm-I_Ca_7.VhCa)))); I_Ca_7.ECa_PU=(I_Ca_7.R*I_Ca_7.T/(2*I_Ca_7.F)*ln(I_Ca_7.CO/CS18)); // PU_membrane_7.I_Na=(I_Na_8.gNa*(Vm-I_Na_8.ENa_PU)); I_Na_8.ENa_PU=(I_Ca_7.R*I_Ca_7.T/I_Ca_7.F*ln(I_Na_8.NO/I_Na_8.NS1)); // I_NSCC_Ca_7.gNSCC_Ca=(I_NSCC_Ca_7.gNSCC_Ca_*(I_NSCC_Ca_7.KNSCC^3/(I_NSCC_Ca_7.KNSCC^3+CS18^3))); PU_membrane_7.I_NSCC_Ca=(I_NSCC_Ca_7.gNSCC_Ca*(Vm-I_NSCC_Ca_7.ENSCC)); // I_NSCC_Na_7.gNSCC_Na=(I_NSCC_Na_7.gNSCC_Na_*(I_NSCC_Ca_7.KNSCC^3/(I_NSCC_Ca_7.KNSCC^3+CS18^3))); PU_membrane_7.I_NSCC_Na=(I_NSCC_Na_7.gNSCC_Na*(Vm-I_NSCC_Ca_7.ENSCC)); // PU_membrane_7.I_PM=(I_PM_7.gPM*(CS18^2/(I_PM_7.KPM^2+CS18^2))); // PU_membrane_7.I_NaP=(I_NaP_7.gNaP*(I_Na_8.NS1^4/(I_NaP_7.KNaP^4+I_Na_8.NS1^4))*(I_NaP_7.ENaP-Vm)); // was_JSERCA_7.JSERCA=(was_JSERCA_7.VSERCA*(CS18-was_JSERCA_7.A2*was_JSERCA_7.CER)/(1+was_JSERCA_7.A4*CS18+was_JSERCA_7.A5*was_JSERCA_7.CER+was_JSERCA_7.A6*CS18*was_JSERCA_7.CER)); // JMCU_7.JMCU=(JMCU_7.VMCU*(CS28^2/(JMCU_7.KMCU^2+CS28^2))*JMCU_7.epsilon_INH); JMCU_7.epsilon_INH=(JMCU_7.KINH^4/(JMCU_7.KINH^4+JMCU_7.CMT^4)); // JIPR_7.H:time=(JIPR_7.phi3*(1-JIPR_7.H)-JIPR_7.P*JIPR_7.phi1*JIPR_7.phi2/(JIPR_7.P*JIPR_7.phi1+JIPR_7.phi_1)*JIPR_7.H); JIPR_7.JIPR=(JIPR_7.kIPR*(JIPR_7.P*JIPR_7.phi1*JIPR_7.H/(JIPR_7.P*JIPR_7.phi1+JIPR_7.phi_1))^4*(was_JSERCA_7.CER-CS28)); JIPR_7.phi1=((JIPR_7.k1*JIPR_7.R1+JIPR_7.r2*CS28)/(JIPR_7.R1+CS28)); JIPR_7.phi_1=((JIPR_7.k_1+JIPR_7.r_2)*JIPR_7.R3/(JIPR_7.R3+CS28)); JIPR_7.phi2=((JIPR_7.k2*JIPR_7.R3+JIPR_7.r4*CS28)/(JIPR_7.R3+CS28)); JIPR_7.phi3=(JIPR_7.g_phi3*JIPR_7.zeta/((1+(JIPR_7.K_phi3_act/CS28)^JIPR_7.h_phi3_act)*(1+(CS28/JIPR_7.K_phi3_inh)^JIPR_7.h_phi3_inh))); JIPR_7.alpha_zeta=JIPR_7.g_alpha; JIPR_7.beta_zeta=(JIPR_7.g_beta/(1+(JIPR_7.K_beta/CS28)^JIPR_7.h_beta)); JIPR_7.zeta:time=(JIPR_7.alpha_zeta*(1-JIPR_7.zeta)-JIPR_7.beta_zeta*JIPR_7.zeta); // JNCX_7.JNCX=(JNCX_7.VNCX*(JMCU_7.CMT/(JMCU_7.CMT+JNCX_7.KNCX))); // CS18:time=(was_JS1S2_8+CS1_7.lambda_MT_S1*JNCX_7.JNCX-(CS1_7.delta_SPU/(PU_membrane_7.alpha_scale*2)*PU_membrane_7.I_iCa+CS1_7.lambda_ER_S1*was_JSERCA_7.JSERCA)); CS1_7.lambda_MT_S1=(CS1_7.gamma_MT/CS1_7.gamma_S1); CS1_7.lambda_ER_S1=(CS1_7.gamma_ER/CS1_7.gamma_S1); // CS28:time=(was_JS2Cy8+CS2_7.lambda_ER_S2*JIPR_7.JIPR-(CS2_7.lambda_S1_S2*was_JS1S2_8+CS2_7.lambda_MT_S2*JMCU_7.JMCU)); CS2_7.lambda_MT_S2=(CS1_7.gamma_MT/CS2_7.gamma_S2); CS2_7.lambda_ER_S2=(CS1_7.gamma_ER/CS2_7.gamma_S2); CS2_7.lambda_S1_S2=(CS1_7.gamma_S1/CS2_7.gamma_S2); // was_JSERCA_7.CER:time=(was_JSERCA_7.JSERCA-JIPR_7.JIPR); // JMCU_7.CMT:time=(CMT_7.fm*(JMCU_7.JMCU-JNCX_7.JNCX)); CMT_7.fm=(1/(1+CMT_7.Km*CMT_7.Bm/(CMT_7.Km+JMCU_7.CMT)^2)); // I_Na_8.NS1:time=((-1)*(CS1_7.delta_SPU/(1*PU_membrane_7.alpha_scale))*PU_membrane_7.I_iNa); // PU_membrane_7.alpha_scale=(model_parameters_8.n_PU_base/n_PU); // PU_membrane_8.I_iCa=(PU_membrane_8.alpha_scale*(PU_membrane_8.I_Ca+PU_membrane_8.I_NSCC_Ca+PU_membrane_8.I_PM)); PU_membrane_8.I_iNa=(PU_membrane_8.alpha_scale*(PU_membrane_8.I_NSCC_Na+PU_membrane_8.I_Na+PU_membrane_8.I_NaP)); I_ion_PU9=(PU_membrane_8.I_iCa+PU_membrane_8.I_iNa); // PU_membrane_8.I_Ca=(I_Ca_8.gCa*(Vm-I_Ca_8.ECa_PU)); I_Ca_8.gCa=(I_Ca_8.gCa_*exp(I_Ca_8.kCa*Vm)/(1+exp(I_Ca_8.kVCa*(Vm-I_Ca_8.VhCa)))); I_Ca_8.ECa_PU=(I_Ca_8.R*I_Ca_8.T/(2*I_Ca_8.F)*ln(I_Ca_8.CO/CS19)); // PU_membrane_8.I_Na=(I_Na_9.gNa*(Vm-I_Na_9.ENa_PU)); I_Na_9.ENa_PU=(I_Ca_8.R*I_Ca_8.T/I_Ca_8.F*ln(I_Na_9.NO/I_Na_9.NS1)); // I_NSCC_Ca_8.gNSCC_Ca=(I_NSCC_Ca_8.gNSCC_Ca_*(I_NSCC_Ca_8.KNSCC^3/(I_NSCC_Ca_8.KNSCC^3+CS19^3))); PU_membrane_8.I_NSCC_Ca=(I_NSCC_Ca_8.gNSCC_Ca*(Vm-I_NSCC_Ca_8.ENSCC)); // I_NSCC_Na_8.gNSCC_Na=(I_NSCC_Na_8.gNSCC_Na_*(I_NSCC_Ca_8.KNSCC^3/(I_NSCC_Ca_8.KNSCC^3+CS19^3))); PU_membrane_8.I_NSCC_Na=(I_NSCC_Na_8.gNSCC_Na*(Vm-I_NSCC_Ca_8.ENSCC)); // PU_membrane_8.I_PM=(I_PM_8.gPM*(CS19^2/(I_PM_8.KPM^2+CS19^2))); // PU_membrane_8.I_NaP=(I_NaP_8.gNaP*(I_Na_9.NS1^4/(I_NaP_8.KNaP^4+I_Na_9.NS1^4))*(I_NaP_8.ENaP-Vm)); // was_JSERCA_8.JSERCA=(was_JSERCA_8.VSERCA*(CS19-was_JSERCA_8.A2*was_JSERCA_8.CER)/(1+was_JSERCA_8.A4*CS19+was_JSERCA_8.A5*was_JSERCA_8.CER+was_JSERCA_8.A6*CS19*was_JSERCA_8.CER)); // JMCU_8.JMCU=(JMCU_8.VMCU*(CS29^2/(JMCU_8.KMCU^2+CS29^2))*JMCU_8.epsilon_INH); JMCU_8.epsilon_INH=(JMCU_8.KINH^4/(JMCU_8.KINH^4+JMCU_8.CMT^4)); // JIPR_8.H:time=(JIPR_8.phi3*(1-JIPR_8.H)-JIPR_8.P*JIPR_8.phi1*JIPR_8.phi2/(JIPR_8.P*JIPR_8.phi1+JIPR_8.phi_1)*JIPR_8.H); JIPR_8.JIPR=(JIPR_8.kIPR*(JIPR_8.P*JIPR_8.phi1*JIPR_8.H/(JIPR_8.P*JIPR_8.phi1+JIPR_8.phi_1))^4*(was_JSERCA_8.CER-CS29)); JIPR_8.phi1=((JIPR_8.k1*JIPR_8.R1+JIPR_8.r2*CS29)/(JIPR_8.R1+CS29)); JIPR_8.phi_1=((JIPR_8.k_1+JIPR_8.r_2)*JIPR_8.R3/(JIPR_8.R3+CS29)); JIPR_8.phi2=((JIPR_8.k2*JIPR_8.R3+JIPR_8.r4*CS29)/(JIPR_8.R3+CS29)); JIPR_8.phi3=(JIPR_8.g_phi3*JIPR_8.zeta/((1+(JIPR_8.K_phi3_act/CS29)^JIPR_8.h_phi3_act)*(1+(CS29/JIPR_8.K_phi3_inh)^JIPR_8.h_phi3_inh))); JIPR_8.alpha_zeta=JIPR_8.g_alpha; JIPR_8.beta_zeta=(JIPR_8.g_beta/(1+(JIPR_8.K_beta/CS29)^JIPR_8.h_beta)); JIPR_8.zeta:time=(JIPR_8.alpha_zeta*(1-JIPR_8.zeta)-JIPR_8.beta_zeta*JIPR_8.zeta); // JNCX_8.JNCX=(JNCX_8.VNCX*(JMCU_8.CMT/(JMCU_8.CMT+JNCX_8.KNCX))); // CS19:time=(was_JS1S2_9+CS1_8.lambda_MT_S1*JNCX_8.JNCX-(CS1_8.delta_SPU/(PU_membrane_8.alpha_scale*2)*PU_membrane_8.I_iCa+CS1_8.lambda_ER_S1*was_JSERCA_8.JSERCA)); CS1_8.lambda_MT_S1=(CS1_8.gamma_MT/CS1_8.gamma_S1); CS1_8.lambda_ER_S1=(CS1_8.gamma_ER/CS1_8.gamma_S1); // CS29:time=(was_JS2Cy9+CS2_8.lambda_ER_S2*JIPR_8.JIPR-(CS2_8.lambda_S1_S2*was_JS1S2_9+CS2_8.lambda_MT_S2*JMCU_8.JMCU)); CS2_8.lambda_MT_S2=(CS1_8.gamma_MT/CS2_8.gamma_S2); CS2_8.lambda_ER_S2=(CS1_8.gamma_ER/CS2_8.gamma_S2); CS2_8.lambda_S1_S2=(CS1_8.gamma_S1/CS2_8.gamma_S2); // was_JSERCA_8.CER:time=(was_JSERCA_8.JSERCA-JIPR_8.JIPR); // JMCU_8.CMT:time=(CMT_8.fm*(JMCU_8.JMCU-JNCX_8.JNCX)); CMT_8.fm=(1/(1+CMT_8.Km*CMT_8.Bm/(CMT_8.Km+JMCU_8.CMT)^2)); // I_Na_9.NS1:time=((-1)*(CS1_8.delta_SPU/(1*PU_membrane_8.alpha_scale))*PU_membrane_8.I_iNa); // PU_membrane_8.alpha_scale=(model_parameters_9.n_PU_base/n_PU); // PU_membrane_9.I_iCa=(PU_membrane_9.alpha_scale*(PU_membrane_9.I_Ca+PU_membrane_9.I_NSCC_Ca+PU_membrane_9.I_PM)); PU_membrane_9.I_iNa=(PU_membrane_9.alpha_scale*(PU_membrane_9.I_NSCC_Na+PU_membrane_9.I_Na+PU_membrane_9.I_NaP)); I_ion_PU10=(PU_membrane_9.I_iCa+PU_membrane_9.I_iNa); // PU_membrane_9.I_Ca=(I_Ca_9.gCa*(Vm-I_Ca_9.ECa_PU)); I_Ca_9.gCa=(I_Ca_9.gCa_*exp(I_Ca_9.kCa*Vm)/(1+exp(I_Ca_9.kVCa*(Vm-I_Ca_9.VhCa)))); I_Ca_9.ECa_PU=(I_Ca_9.R*I_Ca_9.T/(2*I_Ca_9.F)*ln(I_Ca_9.CO/CS110)); // PU_membrane_9.I_Na=(I_Na_10.gNa*(Vm-I_Na_10.ENa_PU)); I_Na_10.ENa_PU=(I_Ca_9.R*I_Ca_9.T/I_Ca_9.F*ln(I_Na_10.NO/I_Na_10.NS1)); // I_NSCC_Ca_9.gNSCC_Ca=(I_NSCC_Ca_9.gNSCC_Ca_*(I_NSCC_Ca_9.KNSCC^3/(I_NSCC_Ca_9.KNSCC^3+CS110^3))); PU_membrane_9.I_NSCC_Ca=(I_NSCC_Ca_9.gNSCC_Ca*(Vm-I_NSCC_Ca_9.ENSCC)); // I_NSCC_Na_9.gNSCC_Na=(I_NSCC_Na_9.gNSCC_Na_*(I_NSCC_Ca_9.KNSCC^3/(I_NSCC_Ca_9.KNSCC^3+CS110^3))); PU_membrane_9.I_NSCC_Na=(I_NSCC_Na_9.gNSCC_Na*(Vm-I_NSCC_Ca_9.ENSCC)); // PU_membrane_9.I_PM=(I_PM_9.gPM*(CS110^2/(I_PM_9.KPM^2+CS110^2))); // PU_membrane_9.I_NaP=(I_NaP_9.gNaP*(I_Na_10.NS1^4/(I_NaP_9.KNaP^4+I_Na_10.NS1^4))*(I_NaP_9.ENaP-Vm)); // was_JSERCA_9.JSERCA=(was_JSERCA_9.VSERCA*(CS110-was_JSERCA_9.A2*was_JSERCA_9.CER)/(1+was_JSERCA_9.A4*CS110+was_JSERCA_9.A5*was_JSERCA_9.CER+was_JSERCA_9.A6*CS110*was_JSERCA_9.CER)); // JMCU_9.JMCU=(JMCU_9.VMCU*(CS210^2/(JMCU_9.KMCU^2+CS210^2))*JMCU_9.epsilon_INH); JMCU_9.epsilon_INH=(JMCU_9.KINH^4/(JMCU_9.KINH^4+JMCU_9.CMT^4)); // JIPR_9.H:time=(JIPR_9.phi3*(1-JIPR_9.H)-JIPR_9.P*JIPR_9.phi1*JIPR_9.phi2/(JIPR_9.P*JIPR_9.phi1+JIPR_9.phi_1)*JIPR_9.H); JIPR_9.JIPR=(JIPR_9.kIPR*(JIPR_9.P*JIPR_9.phi1*JIPR_9.H/(JIPR_9.P*JIPR_9.phi1+JIPR_9.phi_1))^4*(was_JSERCA_9.CER-CS210)); JIPR_9.phi1=((JIPR_9.k1*JIPR_9.R1+JIPR_9.r2*CS210)/(JIPR_9.R1+CS210)); JIPR_9.phi_1=((JIPR_9.k_1+JIPR_9.r_2)*JIPR_9.R3/(JIPR_9.R3+CS210)); JIPR_9.phi2=((JIPR_9.k2*JIPR_9.R3+JIPR_9.r4*CS210)/(JIPR_9.R3+CS210)); JIPR_9.phi3=(JIPR_9.g_phi3*JIPR_9.zeta/((1+(JIPR_9.K_phi3_act/CS210)^JIPR_9.h_phi3_act)*(1+(CS210/JIPR_9.K_phi3_inh)^JIPR_9.h_phi3_inh))); JIPR_9.alpha_zeta=JIPR_9.g_alpha; JIPR_9.beta_zeta=(JIPR_9.g_beta/(1+(JIPR_9.K_beta/CS210)^JIPR_9.h_beta)); JIPR_9.zeta:time=(JIPR_9.alpha_zeta*(1-JIPR_9.zeta)-JIPR_9.beta_zeta*JIPR_9.zeta); // JNCX_9.JNCX=(JNCX_9.VNCX*(JMCU_9.CMT/(JMCU_9.CMT+JNCX_9.KNCX))); // CS110:time=(was_JS1S2_10+CS1_9.lambda_MT_S1*JNCX_9.JNCX-(CS1_9.delta_SPU/(PU_membrane_9.alpha_scale*2)*PU_membrane_9.I_iCa+CS1_9.lambda_ER_S1*was_JSERCA_9.JSERCA)); CS1_9.lambda_MT_S1=(CS1_9.gamma_MT/CS1_9.gamma_S1); CS1_9.lambda_ER_S1=(CS1_9.gamma_ER/CS1_9.gamma_S1); // CS210:time=(was_JS2Cy10+CS2_9.lambda_ER_S2*JIPR_9.JIPR-(CS2_9.lambda_S1_S2*was_JS1S2_10+CS2_9.lambda_MT_S2*JMCU_9.JMCU)); CS2_9.lambda_MT_S2=(CS1_9.gamma_MT/CS2_9.gamma_S2); CS2_9.lambda_ER_S2=(CS1_9.gamma_ER/CS2_9.gamma_S2); CS2_9.lambda_S1_S2=(CS1_9.gamma_S1/CS2_9.gamma_S2); // was_JSERCA_9.CER:time=(was_JSERCA_9.JSERCA-JIPR_9.JIPR); // JMCU_9.CMT:time=(CMT_9.fm*(JMCU_9.JMCU-JNCX_9.JNCX)); CMT_9.fm=(1/(1+CMT_9.Km*CMT_9.Bm/(CMT_9.Km+JMCU_9.CMT)^2)); // I_Na_10.NS1:time=((-1)*(CS1_9.delta_SPU/(1*PU_membrane_9.alpha_scale))*PU_membrane_9.I_iNa); // PU_membrane_9.alpha_scale=(model_parameters_10.n_PU_base/n_PU); }