Generated Code
The following is f77 code generated by the CellML API from this CellML file. (Back to language selection)
The raw code is available.
C
C There are a total of 0 entries in the algebraic variable array.
C There are a total of 5 entries in each of the rate and state variable arrays.
C There are a total of 36 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C STATES(1) is x1 in component x1 (microg_l).
C CONSTS(1) is lambda_1 in component model_parameters (first_order_rate_constant).
C STATES(2) is x3 in component x3 (microg_l).
C CONSTS(2) is a1 in component model_parameters (flux).
C CONSTS(3) is a2 in component model_parameters (flux).
C CONSTS(4) is a3 in component model_parameters (first_order_rate_constant).
C CONSTS(5) is a4 in component model_parameters (second_order_rate_constant).
C CONSTS(6) is a5 in component model_parameters (per_microg_l).
C CONSTS(7) is a6 in component model_parameters (per_microg_l2).
C CONSTS(8) is a7 in component model_parameters (per_microg_l).
C CONSTS(9) is a8 in component model_parameters (per_microg_l2).
C STATES(3) is x2 in component x2 (microg_l).
C CONSTS(10) is lambda_2 in component model_parameters (first_order_rate_constant).
C CONSTS(11) is a9 in component model_parameters (flux).
C CONSTS(12) is a10 in component model_parameters (first_order_rate_constant).
C CONSTS(13) is a11 in component model_parameters (second_order_rate_constant).
C CONSTS(14) is a12 in component model_parameters (per_microg_l).
C CONSTS(15) is a13 in component model_parameters (per_microg_l2).
C CONSTS(16) is a14 in component model_parameters (per_microg_l).
C CONSTS(17) is a15 in component model_parameters (per_microg_l2).
C CONSTS(34) is lambda_3_ in component model_parameters (first_order_rate_constant).
C STATES(4) is x4 in component x4 (microg_l).
C STATES(5) is x5 in component x5 (microg_l).
C CONSTS(18) is a16 in component model_parameters (flux).
C CONSTS(19) is a17 in component model_parameters (first_order_rate_constant).
C CONSTS(20) is a18 in component model_parameters (second_order_rate_constant).
C CONSTS(21) is a19 in component model_parameters (first_order_rate_constant).
C CONSTS(22) is a20 in component model_parameters (second_order_rate_constant).
C CONSTS(23) is a21 in component model_parameters (per_microg_l).
C CONSTS(24) is a22 in component model_parameters (per_microg_l2).
C CONSTS(25) is a23 in component model_parameters (per_microg_l).
C CONSTS(26) is a24 in component model_parameters (per_microg_l2).
C CONSTS(27) is a25 in component model_parameters (first_order_rate_constant).
C CONSTS(28) is a26 in component model_parameters (first_order_rate_constant).
C CONSTS(35) is lambda_4_ in component model_parameters (first_order_rate_constant).
C CONSTS(29) is a27 in component model_parameters (first_order_rate_constant).
C CONSTS(36) is lambda_5_ in component model_parameters (first_order_rate_constant).
C CONSTS(30) is a28 in component model_parameters (first_order_rate_constant).
C CONSTS(31) is lambda_3 in component model_parameters (first_order_rate_constant).
C CONSTS(32) is lambda_4 in component model_parameters (first_order_rate_constant).
C CONSTS(33) is lambda_5 in component model_parameters (first_order_rate_constant).
C RATES(1) is d/dt x1 in component x1 (microg_l).
C RATES(3) is d/dt x2 in component x2 (microg_l).
C RATES(2) is d/dt x3 in component x3 (microg_l).
C RATES(4) is d/dt x4 in component x4 (microg_l).
C RATES(5) is d/dt x5 in component x5 (microg_l).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
STATES(1) = 0.01067
CONSTS(1) = 0.059
STATES(2) = 6.51
CONSTS(2) = 0.000017
CONSTS(3) = 0.0023
CONSTS(4) = 0.6
CONSTS(5) = 45
CONSTS(6) = 36
CONSTS(7) = 216
CONSTS(8) = 0.28
CONSTS(9) = 0.36
STATES(3) = 0.04665
CONSTS(10) = 0.028
CONSTS(11) = 0.0003
CONSTS(12) = 0.18
CONSTS(13) = 150
CONSTS(14) = 18
CONSTS(15) = 460
CONSTS(16) = 0.46
CONSTS(17) = 0.1
STATES(4) = 60.61
STATES(5) = 12.61
CONSTS(18) = 0.04
CONSTS(19) = 150
CONSTS(20) = 3800
CONSTS(21) = 57
CONSTS(22) = 2600
CONSTS(23) = 200
CONSTS(24) = 9400
CONSTS(25) = 10
CONSTS(26) = 320
CONSTS(27) = 0.04
CONSTS(28) = 0.00097
CONSTS(29) = 0.57
CONSTS(30) = 0.0017
CONSTS(31) = 0.0986
CONSTS(32) = 0.024
CONSTS(33) = 3e-5
CONSTS(34) = CONSTS(31)+CONSTS(29)+CONSTS(30)
CONSTS(35) = CONSTS(32)+CONSTS(27)
CONSTS(36) = CONSTS(33)+CONSTS(28)
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
RATES(1) = (CONSTS(2)+(CONSTS(3)+ CONSTS(4)*STATES(1)+ CONSTS(5)*STATES(1) ** 2.00000)/(1.00000+ CONSTS(6)*STATES(1)+ CONSTS(7)*STATES(1) ** 2.00000+ CONSTS(8)*STATES(2)+ CONSTS(9)*STATES(2) ** 2.00000)) - CONSTS(1)*STATES(1)
RATES(3) = (CONSTS(11)+ CONSTS(12)*STATES(1)+ CONSTS(13)*STATES(1) ** 2.00000)/(1.00000+ CONSTS(14)*STATES(1)+ CONSTS(15)*STATES(1) ** 2.00000+ CONSTS(16)*STATES(2)+ CONSTS(17)*STATES(2) ** 2.00000) - CONSTS(10)*STATES(3)
RATES(2) = (CONSTS(18)+( CONSTS(19)*STATES(1)+ CONSTS(20)*STATES(1) ** 2.00000+ CONSTS(21)*STATES(3)+ CONSTS(22)*STATES(3) ** 2.00000)/(1.00000+ CONSTS(23)*STATES(1)+ CONSTS(24)*STATES(1) ** 2.00000+ CONSTS(25)*STATES(3)+ CONSTS(26)*STATES(3) ** 2.00000)+ CONSTS(27)*STATES(4)+ CONSTS(28)*STATES(5)) - CONSTS(34)*STATES(2)
RATES(4) = CONSTS(29)*STATES(2) - CONSTS(35)*STATES(4)
RATES(5) = CONSTS(30)*STATES(2) - CONSTS(36)*STATES(5)
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
RETURN
END
