C C There are a total of 0 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 3 entries in the constant variable array. C C C VOI is t in component main (dimensionless). C STATES(1) is x in component main (dimensionless). C STATES(2) is y in component main (dimensionless). C STATES(3) is z in component main (dimensionless). C CONSTS(1) is sigma in component main (dimensionless). C CONSTS(2) is rho in component main (dimensionless). C CONSTS(3) is beta in component main (dimensionless). C RATES(1) is d/dt x in component main (dimensionless). C RATES(2) is d/dt y in component main (dimensionless). C RATES(3) is d/dt z in component main (dimensionless). C SUBROUTINE initConsts(CONSTS, RATES, STATES) REAL CONSTS(*), RATES(*), STATES(*) STATES(1) = 1 STATES(2) = 1 STATES(3) = 1 CONSTS(1) = 10 CONSTS(2) = 28 CONSTS(3) = 2.66667 RETURN END SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) RATES(1) = CONSTS(1)*(STATES(2) - STATES(1)) RATES(2) = STATES(1)*(CONSTS(2) - STATES(3)) - STATES(2) RATES(3) = STATES(1)*STATES(2) - CONSTS(3)*STATES(3) RETURN END SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) RETURN END