C
C There are a total of 3 entries in the algebraic variable array.
C There are a total of 10 entries in each of the rate and state variable arrays.
C There are a total of 2 entries in the constant variable array.
C
C
C CONSTS(1) is p in component jung_10_lesions (per_sec).
C CONSTS(2) is c in component jung_10_lesions (per_sec).
C STATES(1) is P0 in component jung_10_lesions (conc).
C STATES(2) is P1 in component jung_10_lesions (conc).
C STATES(3) is P2 in component jung_10_lesions (conc).
C STATES(4) is P3 in component jung_10_lesions (conc).
C STATES(5) is P4 in component jung_10_lesions (conc).
C STATES(6) is P5 in component jung_10_lesions (conc).
C STATES(7) is P6 in component jung_10_lesions (conc).
C STATES(8) is P7 in component jung_10_lesions (conc).
C STATES(9) is P8 in component jung_10_lesions (conc).
C STATES(10) is P9 in component jung_10_lesions (conc).
C VOI is time in component jung_10_lesions (second).
C ALGBRC(1) is A in component jung_10_lesions (conc).
C ALGBRC(2) is D in component jung_10_lesions (conc).
C ALGBRC(3) is S in component jung_10_lesions (conc).
C RATES(1) is d/dt P0 in component jung_10_lesions (conc).
C RATES(2) is d/dt P1 in component jung_10_lesions (conc).
C RATES(3) is d/dt P2 in component jung_10_lesions (conc).
C RATES(4) is d/dt P3 in component jung_10_lesions (conc).
C RATES(5) is d/dt P4 in component jung_10_lesions (conc).
C RATES(6) is d/dt P5 in component jung_10_lesions (conc).
C RATES(7) is d/dt P6 in component jung_10_lesions (conc).
C RATES(8) is d/dt P7 in component jung_10_lesions (conc).
C RATES(9) is d/dt P8 in component jung_10_lesions (conc).
C RATES(10) is d/dt P9 in component jung_10_lesions (conc).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 1
      CONSTS(2) = 1
      STATES(1) = 1
      STATES(2) = 0
      STATES(3) = 0
      STATES(4) = 0
      STATES(5) = 0
      STATES(6) = 0
      STATES(7) = 0
      STATES(8) = 0
      STATES(9) = 0
      STATES(10) = 0
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(1) =  - CONSTS(1)*STATES(1)
      RATES(2) = ( CONSTS(1)*STATES(1) -  CONSTS(1)*STATES(2)) -  CONSTS(2)*STATES(2)
      RATES(3) = ( CONSTS(1)*STATES(2) -  CONSTS(1)*STATES(3)) -  2.00000*CONSTS(2)*STATES(3)
      RATES(4) = ( CONSTS(1)*STATES(3) -  CONSTS(1)*STATES(4)) -  3.00000*CONSTS(2)*STATES(4)
      RATES(5) = ( CONSTS(1)*STATES(4) -  CONSTS(1)*STATES(5)) -  4.00000*CONSTS(2)*STATES(5)
      RATES(6) = ( CONSTS(1)*STATES(5) -  CONSTS(1)*STATES(6)) -  5.00000*CONSTS(2)*STATES(6)
      RATES(7) = ( CONSTS(1)*STATES(6) -  CONSTS(1)*STATES(7)) -  6.00000*CONSTS(2)*STATES(7)
      RATES(8) = ( CONSTS(1)*STATES(7) -  CONSTS(1)*STATES(8)) -  7.00000*CONSTS(2)*STATES(8)
      RATES(9) = ( CONSTS(1)*STATES(8) -  CONSTS(1)*STATES(9)) -  8.00000*CONSTS(2)*STATES(9)
      RATES(10) = ( CONSTS(1)*STATES(9) -  CONSTS(1)*STATES(10)) -  9.00000*CONSTS(2)*STATES(10)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) = STATES(1)+STATES(2)+STATES(3)+STATES(4)+STATES(5)+STATES(6)+STATES(7)+STATES(8)+STATES(9)+STATES(10)
      ALGBRC(2) = 1.00000 - (STATES(1)+STATES(2)+STATES(3)+STATES(4)+STATES(5)+STATES(6)+STATES(7)+STATES(8)+STATES(9)+STATES(10))
      ALGBRC(3) =  1.00000*EXP( (CONSTS(1)/CONSTS(2))*((1.00000 -  CONSTS(2)*VOI) - EXP(- ( CONSTS(2)*VOI))))
      RETURN
      END