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 2 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 20 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C CONSTS(1) is a_D in component parameters (first_order_rate_constant).
C CONSTS(2) is a_E in component parameters (first_order_rate_constant).
C CONSTS(3) is a_X in component parameters (first_order_rate_constant).
C CONSTS(4) is k in component parameters (dimensionless).
C CONSTS(5) is q_D in component parameters (dimensionless).
C CONSTS(6) is q_E in component parameters (dimensionless).
C CONSTS(7) is q_X in component parameters (dimensionless).
C CONSTS(8) is f in component parameters (first_order_rate_constant).
C CONSTS(9) is g in component parameters (first_order_rate_constant).
C CONSTS(10) is p_S in component parameters (first_order_rate_constant).
C CONSTS(11) is p_D in component parameters (first_order_rate_constant).
C CONSTS(12) is p_E in component parameters (first_order_rate_constant).
C CONSTS(13) is p_X in component parameters (first_order_rate_constant).
C CONSTS(14) is d_D in component parameters (first_order_rate_constant).
C CONSTS(15) is d_E in component parameters (first_order_rate_constant).
C CONSTS(16) is d_X in component parameters (first_order_rate_constant).
C CONSTS(17) is a_f in component parameters (first_order_rate_constant).
C CONSTS(18) is R_T in component parameters (dimensionless).
C CONSTS(19) is theta in component parameters (dimensionless).
C CONSTS(20) is GF in component parameters (dimensionless).
C STATES(1) is D in component D (dimensionless).
C STATES(2) is E in component E (dimensionless).
C STATES(3) is R_S in component R_S (dimensionless).
C STATES(4) is X in component X (dimensionless).
C STATES(5) is R in component R (dimensionless).
C ALGBRC(1) is unpho_RB in component unpho_RB (dimensionless).
C ALGBRC(2) is free_E2F in component free_E2F (dimensionless).
C RATES(1) is d/dt D in component D (dimensionless).
C RATES(2) is d/dt E in component E (dimensionless).
C RATES(5) is d/dt R in component R (dimensionless).
C RATES(3) is d/dt R_S in component R_S (dimensionless).
C RATES(4) is d/dt X in component X (dimensionless).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
CONSTS(1) = 0.4
CONSTS(2) = 0.16
CONSTS(3) = 0.08
CONSTS(4) = 0.05
CONSTS(5) = 0.6
CONSTS(6) = 0.6
CONSTS(7) = 0.8
CONSTS(8) = 0.2
CONSTS(9) = 0.528
CONSTS(10) = 0.6
CONSTS(11) = 0.48
CONSTS(12) = 0.096
CONSTS(13) = 0.48
CONSTS(14) = 0.4
CONSTS(15) = 0.2
CONSTS(16) = 1.04
CONSTS(17) = 0.9
CONSTS(18) = 2.5
CONSTS(19) = 1.5
CONSTS(20) = 6.3
STATES(1) = 0.1
STATES(2) = 0.6
STATES(3) = 1
STATES(4) = 0.7
STATES(5) = 0.5
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
RATES(1) = ( CONSTS(1)*CONSTS(4)*CONSTS(20))/(1.00000+ CONSTS(4)*CONSTS(20)) - CONSTS(14)*STATES(2)*STATES(1)
RATES(2) = CONSTS(2)*(1.00000+ CONSTS(17)*(CONSTS(19) - STATES(3)))*1.00000 - CONSTS(15)*STATES(4)*STATES(2)
RATES(5) = ( CONSTS(13)*STATES(4)*((CONSTS(18) - STATES(3)) - STATES(5)))/(CONSTS(7)+((CONSTS(18) - STATES(3)) - STATES(5))+STATES(4)) - CONSTS(10)*(CONSTS(19) - STATES(3))*STATES(5)
RATES(3) = ( CONSTS(10)*(CONSTS(19) - STATES(3))*STATES(5) - ( CONSTS(11)*STATES(3)*STATES(1))/(CONSTS(5)+STATES(3)+STATES(1))) - ( CONSTS(12)*STATES(3)*STATES(2))/(CONSTS(6)+STATES(3)+STATES(2))
RATES(4) = ( CONSTS(3)*STATES(2)+ CONSTS(8)*(CONSTS(19) - STATES(3))+ CONSTS(9)*STATES(4) ** 2.00000*STATES(2)) - CONSTS(16)*STATES(4)
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(1) = STATES(5)+STATES(3)
ALGBRC(2) = CONSTS(19) - STATES(3)
RETURN
END
