C
C There are a total of 12 entries in the algebraic variable array.
C There are a total of 4 entries in each of the rate and state variable arrays.
C There are a total of 20 entries in the constant variable array.
C
C
C VOI is time in component environment (day).
C ALGBRC(7) is rel_LH_E2_P4_RP_LH in component RP_LH (microg_day).
C STATES(1) is RP_LH in component RP_LH (microg).
C ALGBRC(11) is syn_LH_E2_P4 in component RP_LH (microg_day).
C CONSTS(1) is V0_LH in component RP_LH (microg_day).
C CONSTS(2) is V1_LH in component RP_LH (microg_day).
C CONSTS(3) is h in component RP_LH (dimensionless).
C CONSTS(4) is Km_LH in component RP_LH (ng_L).
C CONSTS(5) is Ki_LHP in component RP_LH (nmol_L).
C CONSTS(6) is kLH_rel in component RP_LH (first_order_rate_constant).
C CONSTS(7) is CLH_P in component RP_LH (L_nmol).
C CONSTS(8) is CLH_E in component RP_LH (L_ng).
C ALGBRC(4) is E2 in component E2 (ng_L).
C ALGBRC(5) is E2_dE in component E2_dE (ng_L).
C ALGBRC(6) is P4 in component P4 (nmol_L).
C ALGBRC(9) is P4_dP in component P4_dP (nmol_L).
C STATES(2) is LH in component LH (microg_l).
C CONSTS(9) is v_dis in component LH (litre).
C ALGBRC(1) is clear_LH in component LH (microg_l_day).
C CONSTS(10) is kLH_cl in component LH (first_order_rate_constant).
C ALGBRC(8) is rel_FSH_E2_P4_RP_FSH in component RP_FSH (microg_day).
C STATES(3) is RP_FSH in component RP_FSH (microg).
C ALGBRC(12) is syn_FSH_Ih in component RP_FSH (microg_day).
C CONSTS(11) is V_FSH in component RP_FSH (microg_day).
C CONSTS(12) is Ki_FSH_Ih in component RP_FSH (U_L).
C CONSTS(13) is kFSH_rel in component RP_FSH (first_order_rate_constant).
C CONSTS(14) is CFSH_P in component RP_FSH (L_nmol).
C CONSTS(15) is CFSH_E in component RP_FSH (L_ng2).
C ALGBRC(10) is Ih_dIh in component Ih_dIh (U_L).
C STATES(4) is FSH in component FSH (microg_l).
C CONSTS(16) is v_dis in component FSH (litre).
C ALGBRC(3) is clear_FSH in component FSH (microg_l_day).
C CONSTS(17) is kFSH_cl in component FSH (first_order_rate_constant).
C CONSTS(18) is dE in component E2_dE (day).
C CONSTS(19) is dP in component P4_dP (day).
C ALGBRC(2) is Ih in component Ih (U_L).
C CONSTS(20) is dIh in component Ih_dIh (day).
C RATES(1) is d/dt RP_LH in component RP_LH (microg).
C RATES(2) is d/dt LH in component LH (microg_l).
C RATES(3) is d/dt RP_FSH in component RP_FSH (microg).
C RATES(4) is d/dt FSH in component FSH (microg_l).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = 467.0
      CONSTS(1) = 1400.0
      CONSTS(2) = 95900.0
      CONSTS(3) = 8.0
      CONSTS(4) = 360.0
      CONSTS(5) = 26.0
      CONSTS(6) = 3.0
      CONSTS(7) = 0.024
      CONSTS(8) = 0.008
      STATES(2) = 40.0
      CONSTS(9) = 2.5
      CONSTS(10) = 14.0
      STATES(3) = 0.0
      CONSTS(11) = 4400.0
      CONSTS(12) = 1176.5
      CONSTS(13) = 45.0
      CONSTS(14) = 3.0
      CONSTS(15) = 0.005
      STATES(4) = 150.0
      CONSTS(16) = 2.5
      CONSTS(17) = 8.21
      CONSTS(18) = 0.42
      CONSTS(19) = 2.9
      CONSTS(20) = 2.0
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(4) = (300.000 - ( 240.000*VOI+1.00000 ** 2.00000)/(3.00000+VOI+1.00000 ** 2.00000))+ 90.0000*EXP(- (VOI - 8.00000 ** 2.00000/10.0000))
      ALGBRC(6) =  52.0000*EXP(- (VOI - 7.00000 ** 2.00000/18.0000))
      ALGBRC(7) = ( CONSTS(6)*(1.00000+ CONSTS(7)*ALGBRC(6))*STATES(1))/(1.00000+ CONSTS(8)*ALGBRC(4))
      ALGBRC(1) =  CONSTS(10)*STATES(2)
      RATES(2) = ALGBRC(7)/CONSTS(9) - ALGBRC(1)
      ALGBRC(8) = ( CONSTS(13)*(1.00000+ CONSTS(14)*ALGBRC(6))*STATES(3))/(1.00000+ CONSTS(15)*ALGBRC(4) ** 2.00000)
      ALGBRC(3) =  CONSTS(17)*STATES(4)
      RATES(4) = ALGBRC(8)/CONSTS(16) - ALGBRC(3)
      ALGBRC(5) = (300.000 - ( 240.000*(VOI+1.00000) - CONSTS(18) ** 2.00000)/(3.00000+(VOI+1.00000) - CONSTS(18) ** 2.00000))+ 90.0000*EXP(- (VOI - (CONSTS(18)+8.00000) ** 2.00000/10.0000))
      ALGBRC(9) =  52.0000*EXP(- (VOI - (CONSTS(19)+7.00000) ** 2.00000/18.0000))
      ALGBRC(11) = (CONSTS(1)+( CONSTS(2)*ALGBRC(5) ** CONSTS(3))/(CONSTS(4) ** CONSTS(3)+ALGBRC(5) ** CONSTS(3)))/(1.00000+ALGBRC(9)/CONSTS(5))
      RATES(1) = ALGBRC(11) - ALGBRC(7)
      ALGBRC(10) = 300.000+ 1330.00*EXP(- (VOI - (7.00000+CONSTS(20)) ** 2.00000/19.0000))
      ALGBRC(12) = CONSTS(11)/(1.00000+ALGBRC(10)/CONSTS(12))
      RATES(3) = ALGBRC(12) - ALGBRC(8)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(4) = (300.000 - ( 240.000*VOI+1.00000 ** 2.00000)/(3.00000+VOI+1.00000 ** 2.00000))+ 90.0000*EXP(- (VOI - 8.00000 ** 2.00000/10.0000))
      ALGBRC(6) =  52.0000*EXP(- (VOI - 7.00000 ** 2.00000/18.0000))
      ALGBRC(7) = ( CONSTS(6)*(1.00000+ CONSTS(7)*ALGBRC(6))*STATES(1))/(1.00000+ CONSTS(8)*ALGBRC(4))
      ALGBRC(1) =  CONSTS(10)*STATES(2)
      ALGBRC(8) = ( CONSTS(13)*(1.00000+ CONSTS(14)*ALGBRC(6))*STATES(3))/(1.00000+ CONSTS(15)*ALGBRC(4) ** 2.00000)
      ALGBRC(3) =  CONSTS(17)*STATES(4)
      ALGBRC(5) = (300.000 - ( 240.000*(VOI+1.00000) - CONSTS(18) ** 2.00000)/(3.00000+(VOI+1.00000) - CONSTS(18) ** 2.00000))+ 90.0000*EXP(- (VOI - (CONSTS(18)+8.00000) ** 2.00000/10.0000))
      ALGBRC(9) =  52.0000*EXP(- (VOI - (CONSTS(19)+7.00000) ** 2.00000/18.0000))
      ALGBRC(11) = (CONSTS(1)+( CONSTS(2)*ALGBRC(5) ** CONSTS(3))/(CONSTS(4) ** CONSTS(3)+ALGBRC(5) ** CONSTS(3)))/(1.00000+ALGBRC(9)/CONSTS(5))
      ALGBRC(10) = 300.000+ 1330.00*EXP(- (VOI - (7.00000+CONSTS(20)) ** 2.00000/19.0000))
      ALGBRC(12) = CONSTS(11)/(1.00000+ALGBRC(10)/CONSTS(12))
      ALGBRC(2) = 300.000+ 1330.00*EXP(- (VOI - 7.00000 ** 2.00000/19.0000))
      RETURN
      END