C
C There are a total of 18 entries in the algebraic variable array.
C There are a total of 8 entries in each of the rate and state variable arrays.
C There are a total of 12 entries in the constant variable array.
C
C
C VOI is time in component environment (ms).
C STATES(1) is V in component membrane (mV).
C CONSTS(1) is C in component membrane (uF_per_mm2).
C ALGBRC(1) is i_Na in component sodium_current (uA_per_mm2).
C ALGBRC(15) is i_s in component slow_inward_current (uA_per_mm2).
C ALGBRC(16) is i_x1 in component time_dependent_outward_current (uA_per_mm2).
C ALGBRC(17) is i_K1 in component time_independent_outward_current (uA_per_mm2).
C ALGBRC(18) is Istim in component stimulus_protocol (uA_per_mm2).
C CONSTS(2) is g_Na in component sodium_current (mS_per_mm2).
C CONSTS(3) is E_Na in component sodium_current (mV).
C CONSTS(4) is g_Nac in component sodium_current (mS_per_mm2).
C STATES(2) is m in component sodium_current_m_gate (dimensionless).
C STATES(3) is h in component sodium_current_h_gate (dimensionless).
C STATES(4) is j in component sodium_current_j_gate (dimensionless).
C ALGBRC(2) is alpha_m in component sodium_current_m_gate (per_ms).
C ALGBRC(9) is beta_m in component sodium_current_m_gate (per_ms).
C ALGBRC(3) is alpha_h in component sodium_current_h_gate (per_ms).
C ALGBRC(10) is beta_h in component sodium_current_h_gate (per_ms).
C ALGBRC(4) is alpha_j in component sodium_current_j_gate (per_ms).
C ALGBRC(11) is beta_j in component sodium_current_j_gate (per_ms).
C CONSTS(5) is g_s in component slow_inward_current (mS_per_mm2).
C ALGBRC(8) is E_s in component slow_inward_current (mV).
C STATES(5) is Cai in component slow_inward_current (concentration_units).
C STATES(6) is d in component slow_inward_current_d_gate (dimensionless).
C STATES(7) is f in component slow_inward_current_f_gate (dimensionless).
C ALGBRC(5) is alpha_d in component slow_inward_current_d_gate (per_ms).
C ALGBRC(12) is beta_d in component slow_inward_current_d_gate (per_ms).
C ALGBRC(6) is alpha_f in component slow_inward_current_f_gate (per_ms).
C ALGBRC(13) is beta_f in component slow_inward_current_f_gate (per_ms).
C STATES(8) is x1 in component time_dependent_outward_current_x1_gate (dimensionless).
C CONSTS(12) is x1_open in component time_dependent_outward_current (uA_per_mm2).
C CONSTS(6) is x in component time_dependent_outward_current (dimensionless).
C ALGBRC(7) is alpha_x1 in component time_dependent_outward_current_x1_gate (per_ms).
C ALGBRC(14) is beta_x1 in component time_dependent_outward_current_x1_gate (per_ms).
C CONSTS(7) is IstimStart in component stimulus_protocol (ms).
C CONSTS(8) is IstimEnd in component stimulus_protocol (ms).
C CONSTS(9) is IstimAmplitude in component stimulus_protocol (uA_per_mm2).
C CONSTS(10) is IstimPeriod in component stimulus_protocol (ms).
C CONSTS(11) is IstimPulseDuration in component stimulus_protocol (ms).
C RATES(1) is d/dt V in component membrane (mV).
C RATES(2) is d/dt m in component sodium_current_m_gate (dimensionless).
C RATES(3) is d/dt h in component sodium_current_h_gate (dimensionless).
C RATES(4) is d/dt j in component sodium_current_j_gate (dimensionless).
C RATES(5) is d/dt Cai in component slow_inward_current (concentration_units).
C RATES(6) is d/dt d in component slow_inward_current_d_gate (dimensionless).
C RATES(7) is d/dt f in component slow_inward_current_f_gate (dimensionless).
C RATES(8) is d/dt x1 in component time_dependent_outward_current_x1_gate (dimensionless).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = -84.624
      CONSTS(1) = 0.01
      CONSTS(2) = 4e-2
      CONSTS(3) = 50
      CONSTS(4) = 3e-5
      STATES(2) = 0.011
      STATES(3) = 0.988
      STATES(4) = 0.975
      CONSTS(5) = 9e-4
      STATES(5) = 1e-4
      STATES(6) = 0.003
      STATES(7) = 0.994
      STATES(8) = 0.0001
      CONSTS(6) = 0
      CONSTS(7) = 10
      CONSTS(8) = 50000
      CONSTS(9) = 0.5
      CONSTS(10) = 1000
      CONSTS(11) = 1
      CONSTS(12) = SampleUsingPDF(&pdf_0, 0, pdf_roots_0, CONSTANTS, ALGEBRAIC)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(2) = ( - 1.00000*(STATES(1)+47.0000))/(EXP( - 0.100000*(STATES(1)+47.0000)) - 1.00000)
      ALGBRC(9) =  40.0000*EXP( - 0.0560000*(STATES(1)+72.0000))
      RATES(2) =  ALGBRC(2)*(1.00000 - STATES(2)) -  ALGBRC(9)*STATES(2)
      ALGBRC(3) =  0.126000*EXP( - 0.250000*(STATES(1)+77.0000))
      ALGBRC(10) = 1.70000/(EXP( - 0.0820000*(STATES(1)+22.5000))+1.00000)
      RATES(3) =  ALGBRC(3)*(1.00000 - STATES(3)) -  ALGBRC(10)*STATES(3)
      ALGBRC(4) = ( 0.0550000*EXP( - 0.250000*(STATES(1)+78.0000)))/(EXP( - 0.200000*(STATES(1)+78.0000))+1.00000)
      ALGBRC(11) = 0.300000/(EXP( - 0.100000*(STATES(1)+32.0000))+1.00000)
      RATES(4) =  ALGBRC(4)*(1.00000 - STATES(4)) -  ALGBRC(11)*STATES(4)
      ALGBRC(5) = ( 0.0950000*EXP(- (STATES(1) - 5.00000)/100.000))/(1.00000+EXP(- (STATES(1) - 5.00000)/13.8900))
      ALGBRC(12) = ( 0.0700000*EXP(- (STATES(1)+44.0000)/59.0000))/(1.00000+EXP((STATES(1)+44.0000)/20.0000))
      RATES(6) =  ALGBRC(5)*(1.00000 - STATES(6)) -  ALGBRC(12)*STATES(6)
      ALGBRC(6) = ( 0.0120000*EXP(- (STATES(1)+28.0000)/125.000))/(1.00000+EXP((STATES(1)+28.0000)/6.67000))
      ALGBRC(13) = ( 0.00650000*EXP(- (STATES(1)+30.0000)/50.0000))/(1.00000+EXP(- (STATES(1)+30.0000)/5.00000))
      RATES(7) =  ALGBRC(6)*(1.00000 - STATES(7)) -  ALGBRC(13)*STATES(7)
      ALGBRC(7) = ( 0.000500000*EXP((STATES(1)+50.0000)/12.1000))/(1.00000+EXP((STATES(1)+50.0000)/17.5000))
      ALGBRC(14) = ( 0.00130000*EXP(- (STATES(1)+20.0000)/16.6700))/(1.00000+EXP(- (STATES(1)+20.0000)/25.0000))
      RATES(8) =  ALGBRC(7)*(1.00000 - STATES(8)) -  ALGBRC(14)*STATES(8)
      ALGBRC(8) = - 82.3000 -  13.0287*log( STATES(5)*0.00100000)
      ALGBRC(15) =  CONSTS(5)*STATES(6)*STATES(7)*(STATES(1) - ALGBRC(8))
      RATES(5) = ( - 0.0100000*ALGBRC(15))/1.00000+ 0.0700000*(0.000100000 - STATES(5))
      ALGBRC(1) =  ( CONSTS(2)*STATES(2) ** 3.00000*STATES(3)*STATES(4)+CONSTS(4))*(STATES(1) - CONSTS(3))
      ALGBRC(16) = ( STATES(8)*CONSTS(12)*(EXP( 0.0400000*(STATES(1)+77.0000)) - 1.00000))/EXP( 0.0400000*(STATES(1)+35.0000))
      ALGBRC(17) =  0.00350000*(( 4.00000*(EXP( 0.0400000*(STATES(1)+85.0000)) - 1.00000))/(EXP( 0.0800000*(STATES(1)+53.0000))+EXP( 0.0400000*(STATES(1)+53.0000)))+( 0.200000*(STATES(1)+23.0000))/(1.00000 - EXP( - 0.0400000*(STATES(1)+23.0000))))
      ALGBRC(18) = TERNRY(VOI.GE.CONSTS(7).AND.VOI.LE.CONSTS(8).AND.(VOI - CONSTS(7)) -  INT((VOI - CONSTS(7))/CONSTS(10))*CONSTS(10).LE.CONSTS(11), CONSTS(9), 0.00000)
      RATES(1) = (ALGBRC(18) - (ALGBRC(1)+ALGBRC(15)+ALGBRC(16)+ALGBRC(17)))/CONSTS(1)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(2) = ( - 1.00000*(STATES(1)+47.0000))/(EXP( - 0.100000*(STATES(1)+47.0000)) - 1.00000)
      ALGBRC(9) =  40.0000*EXP( - 0.0560000*(STATES(1)+72.0000))
      ALGBRC(3) =  0.126000*EXP( - 0.250000*(STATES(1)+77.0000))
      ALGBRC(10) = 1.70000/(EXP( - 0.0820000*(STATES(1)+22.5000))+1.00000)
      ALGBRC(4) = ( 0.0550000*EXP( - 0.250000*(STATES(1)+78.0000)))/(EXP( - 0.200000*(STATES(1)+78.0000))+1.00000)
      ALGBRC(11) = 0.300000/(EXP( - 0.100000*(STATES(1)+32.0000))+1.00000)
      ALGBRC(5) = ( 0.0950000*EXP(- (STATES(1) - 5.00000)/100.000))/(1.00000+EXP(- (STATES(1) - 5.00000)/13.8900))
      ALGBRC(12) = ( 0.0700000*EXP(- (STATES(1)+44.0000)/59.0000))/(1.00000+EXP((STATES(1)+44.0000)/20.0000))
      ALGBRC(6) = ( 0.0120000*EXP(- (STATES(1)+28.0000)/125.000))/(1.00000+EXP((STATES(1)+28.0000)/6.67000))
      ALGBRC(13) = ( 0.00650000*EXP(- (STATES(1)+30.0000)/50.0000))/(1.00000+EXP(- (STATES(1)+30.0000)/5.00000))
      ALGBRC(7) = ( 0.000500000*EXP((STATES(1)+50.0000)/12.1000))/(1.00000+EXP((STATES(1)+50.0000)/17.5000))
      ALGBRC(14) = ( 0.00130000*EXP(- (STATES(1)+20.0000)/16.6700))/(1.00000+EXP(- (STATES(1)+20.0000)/25.0000))
      ALGBRC(8) = - 82.3000 -  13.0287*log( STATES(5)*0.00100000)
      ALGBRC(15) =  CONSTS(5)*STATES(6)*STATES(7)*(STATES(1) - ALGBRC(8))
      ALGBRC(1) =  ( CONSTS(2)*STATES(2) ** 3.00000*STATES(3)*STATES(4)+CONSTS(4))*(STATES(1) - CONSTS(3))
      ALGBRC(16) = ( STATES(8)*CONSTS(12)*(EXP( 0.0400000*(STATES(1)+77.0000)) - 1.00000))/EXP( 0.0400000*(STATES(1)+35.0000))
      ALGBRC(17) =  0.00350000*(( 4.00000*(EXP( 0.0400000*(STATES(1)+85.0000)) - 1.00000))/(EXP( 0.0800000*(STATES(1)+53.0000))+EXP( 0.0400000*(STATES(1)+53.0000)))+( 0.200000*(STATES(1)+23.0000))/(1.00000 - EXP( - 0.0400000*(STATES(1)+23.0000))))
      ALGBRC(18) = TERNRY(VOI.GE.CONSTS(7).AND.VOI.LE.CONSTS(8).AND.(VOI - CONSTS(7)) -  INT((VOI - CONSTS(7))/CONSTS(10))*CONSTS(10).LE.CONSTS(11), CONSTS(9), 0.00000)
      RETURN
      END
double pdf_0(double bvar, double* CONSTANTS, double* ALGEBRAIC)
{
  return ( (1.00000/ ( 2.00000* 3.14159265358979*1.00000e-06) ** (1.0 / 2))*EXP(- bvar - 0.00800000 ** 2.00000/( 2.00000*1.00000e-06)));
}
double (*pdf_roots_0[])(double bvar, double*, double*) = {};
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END