Generated Code

C
C There are a total of 4 entries in the algebraic variable array.
C There are a total of 3 entries in each of the rate and state variable arrays.
C There are a total of 15 entries in the constant variable array.
C
C
C VOI is time in component environment (hour).
C ALGBRC(1) is Qv in component v (per_second).
C STATES(1) is Vv in component v (mV).
C CONSTS(1) is tau_v in component v (second).
C CONSTS(2) is v_vm in component v (mV_second).
C ALGBRC(2) is Qm in component m (per_second).
C CONSTS(3) is Qmax in component model_parameters (per_second).
C ALGBRC(4) is D in component D (mV).
C CONSTS(4) is theta in component model_parameters (mV).
C CONSTS(5) is sigma in component model_parameters (mV).
C STATES(2) is Vm in component m (mV).
C CONSTS(6) is tau_m in component m (second).
C CONSTS(7) is v_mv in component m (mV_second).
C CONSTS(8) is A in component m (mV).
C STATES(3) is H in component H (nM).
C CONSTS(9) is chi in component H (hour).
C CONSTS(10) is mu in component H (nM_second).
C ALGBRC(3) is C in component D (dimensionless).
C CONSTS(11) is c0 in component D (dimensionless).
C CONSTS(15) is omega in component D (per_hour).
C CONSTS(12) is alpha in component D (hour).
C CONSTS(13) is v_vc in component D (mV).
C CONSTS(14) is v_vh in component D (mV_per_nM).
C RATES(1) is d/dt Vv in component v (mV).
C RATES(2) is d/dt Vm in component m (mV).
C RATES(3) is d/dt H in component H (nM).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = -6.0
      CONSTS(1) = 10.0
      CONSTS(2) = -1.9
      CONSTS(3) = 100.0
      CONSTS(4) = 10.0
      CONSTS(5) = 3.0
      STATES(2) = 0.0
      CONSTS(6) = 10.0
      CONSTS(7) = -1.9
      CONSTS(8) = 1.0
      STATES(3) = 15.0
      CONSTS(9) = 11.0
      CONSTS(10) = 3.6
      CONSTS(11) = 1.0
      CONSTS(12) = 16.0
      CONSTS(13) = -3.2
      CONSTS(14) = 0.19
      CONSTS(15) = ( 2.00000* 3.14159265358979)/24.0000
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) = CONSTS(3)/(1.00000+EXP(- (STATES(1) - CONSTS(4))/CONSTS(5)))
      RATES(2) = ((CONSTS(8)+ CONSTS(7)*ALGBRC(1)) - STATES(2))/(CONSTS(6)/3600.00)
      ALGBRC(2) = CONSTS(3)/(1.00000+EXP(- (STATES(2) - CONSTS(4))/CONSTS(5)))
      RATES(3) = ( CONSTS(10)*ALGBRC(2) - STATES(3))/CONSTS(9)
      ALGBRC(3) = CONSTS(11)+cos( CONSTS(15)*(VOI - CONSTS(12)))
      ALGBRC(4) =  CONSTS(13)*ALGBRC(3)+ CONSTS(14)*STATES(3)
      RATES(1) = (( CONSTS(2)*ALGBRC(2)+ALGBRC(4)) - STATES(1))/(CONSTS(1)/3600.00)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) = CONSTS(3)/(1.00000+EXP(- (STATES(1) - CONSTS(4))/CONSTS(5)))
      ALGBRC(2) = CONSTS(3)/(1.00000+EXP(- (STATES(2) - CONSTS(4))/CONSTS(5)))
      ALGBRC(3) = CONSTS(11)+cos( CONSTS(15)*(VOI - CONSTS(12)))
      ALGBRC(4) =  CONSTS(13)*ALGBRC(3)+ CONSTS(14)*STATES(3)
      RETURN
      END