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 1 entries in the algebraic variable array. C There are a total of 1 entries in each of the rate and state variable arrays. C There are a total of 6 entries in the constant variable array. C C C VOI is time in component interface (ms). C CONSTS(1) is lambda in component interface (dimensionless). C CONSTS(2) is Ca_max in component interface (mM). C CONSTS(3) is beta in component interface (dimensionless). C CONSTS(4) is c50 in component interface (mM). C CONSTS(5) is h in component interface (dimensionless). C CONSTS(6) is T_ref in component interface (kPa). C ALGBRC(1) is T in component active_contraction (kPa). C STATES(1) is Ca_actn in component interface (mM). C RATES(1) is d/dt Ca_actn in component interface (mM). C SUBROUTINE initConsts(CONSTS, RATES, STATES) REAL CONSTS(*), RATES(*), STATES(*) CONSTS(1) = 1 CONSTS(2) = 1 CONSTS(3) = 1.45 CONSTS(4) = 0.5 CONSTS(5) = 3 CONSTS(6) = 100 STATES(1) = 0 RETURN END SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) RATES(1) = VOI*CONSTS(2) RETURN END SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) ALGBRC(1) = (( CONSTS(2)*STATES(1) ** CONSTS(5))/(STATES(1) ** CONSTS(5)+CONSTS(4) ** CONSTS(5)))*CONSTS(6)*(1.00000+ CONSTS(3)*(CONSTS(1) - 1.00000)) RETURN END