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 160 entries in the algebraic variable array.
C There are a total of 50 entries in each of the rate and state variable arrays.
C There are a total of 370 entries in the constant variable array.
C
C
C ALGBRC(1) is sum_g_syn_e_s_1 in component synaptic_coupling (nanoS).
C ALGBRC(6) is sum_g_syn_e_s_2 in component synaptic_coupling (nanoS).
C ALGBRC(11) is sum_g_syn_e_s_3 in component synaptic_coupling (nanoS).
C ALGBRC(16) is sum_g_syn_e_s_4 in component synaptic_coupling (nanoS).
C ALGBRC(21) is sum_g_syn_e_s_5 in component synaptic_coupling (nanoS).
C ALGBRC(26) is sum_g_syn_e_s_6 in component synaptic_coupling (nanoS).
C ALGBRC(31) is sum_g_syn_e_s_7 in component synaptic_coupling (nanoS).
C ALGBRC(36) is sum_g_syn_e_s_8 in component synaptic_coupling (nanoS).
C ALGBRC(41) is sum_g_syn_e_s_9 in component synaptic_coupling (nanoS).
C ALGBRC(46) is sum_g_syn_e_s_10 in component synaptic_coupling (nanoS).
C CONSTS(1) is g_syn_e_1_2 in component synaptic_coupling (nanoS).
C CONSTS(2) is g_syn_e_1_3 in component synaptic_coupling (nanoS).
C CONSTS(3) is g_syn_e_1_4 in component synaptic_coupling (nanoS).
C CONSTS(4) is g_syn_e_1_5 in component synaptic_coupling (nanoS).
C CONSTS(5) is g_syn_e_1_6 in component synaptic_coupling (nanoS).
C CONSTS(6) is g_syn_e_1_7 in component synaptic_coupling (nanoS).
C CONSTS(7) is g_syn_e_1_8 in component synaptic_coupling (nanoS).
C CONSTS(8) is g_syn_e_1_9 in component synaptic_coupling (nanoS).
C CONSTS(9) is g_syn_e_1_10 in component synaptic_coupling (nanoS).
C CONSTS(10) is g_syn_e_2_1 in component synaptic_coupling (nanoS).
C CONSTS(11) is g_syn_e_2_3 in component synaptic_coupling (nanoS).
C CONSTS(12) is g_syn_e_2_4 in component synaptic_coupling (nanoS).
C CONSTS(13) is g_syn_e_2_5 in component synaptic_coupling (nanoS).
C CONSTS(14) is g_syn_e_2_6 in component synaptic_coupling (nanoS).
C CONSTS(15) is g_syn_e_2_7 in component synaptic_coupling (nanoS).
C CONSTS(16) is g_syn_e_2_8 in component synaptic_coupling (nanoS).
C CONSTS(17) is g_syn_e_2_9 in component synaptic_coupling (nanoS).
C CONSTS(18) is g_syn_e_2_10 in component synaptic_coupling (nanoS).
C CONSTS(19) is g_syn_e_3_1 in component synaptic_coupling (nanoS).
C CONSTS(20) is g_syn_e_3_2 in component synaptic_coupling (nanoS).
C CONSTS(21) is g_syn_e_3_4 in component synaptic_coupling (nanoS).
C CONSTS(22) is g_syn_e_3_5 in component synaptic_coupling (nanoS).
C CONSTS(23) is g_syn_e_3_6 in component synaptic_coupling (nanoS).
C CONSTS(24) is g_syn_e_3_7 in component synaptic_coupling (nanoS).
C CONSTS(25) is g_syn_e_3_8 in component synaptic_coupling (nanoS).
C CONSTS(26) is g_syn_e_3_9 in component synaptic_coupling (nanoS).
C CONSTS(27) is g_syn_e_3_10 in component synaptic_coupling (nanoS).
C CONSTS(28) is g_syn_e_4_1 in component synaptic_coupling (nanoS).
C CONSTS(29) is g_syn_e_4_2 in component synaptic_coupling (nanoS).
C CONSTS(30) is g_syn_e_4_3 in component synaptic_coupling (nanoS).
C CONSTS(31) is g_syn_e_4_5 in component synaptic_coupling (nanoS).
C CONSTS(32) is g_syn_e_4_6 in component synaptic_coupling (nanoS).
C CONSTS(33) is g_syn_e_4_7 in component synaptic_coupling (nanoS).
C CONSTS(34) is g_syn_e_4_8 in component synaptic_coupling (nanoS).
C CONSTS(35) is g_syn_e_4_9 in component synaptic_coupling (nanoS).
C CONSTS(36) is g_syn_e_4_10 in component synaptic_coupling (nanoS).
C CONSTS(37) is g_syn_e_5_1 in component synaptic_coupling (nanoS).
C CONSTS(38) is g_syn_e_5_2 in component synaptic_coupling (nanoS).
C CONSTS(39) is g_syn_e_5_3 in component synaptic_coupling (nanoS).
C CONSTS(40) is g_syn_e_5_4 in component synaptic_coupling (nanoS).
C CONSTS(41) is g_syn_e_5_6 in component synaptic_coupling (nanoS).
C CONSTS(42) is g_syn_e_5_7 in component synaptic_coupling (nanoS).
C CONSTS(43) is g_syn_e_5_8 in component synaptic_coupling (nanoS).
C CONSTS(44) is g_syn_e_5_9 in component synaptic_coupling (nanoS).
C CONSTS(45) is g_syn_e_5_10 in component synaptic_coupling (nanoS).
C CONSTS(46) is g_syn_e_6_1 in component synaptic_coupling (nanoS).
C CONSTS(47) is g_syn_e_6_2 in component synaptic_coupling (nanoS).
C CONSTS(48) is g_syn_e_6_3 in component synaptic_coupling (nanoS).
C CONSTS(49) is g_syn_e_6_4 in component synaptic_coupling (nanoS).
C CONSTS(50) is g_syn_e_6_5 in component synaptic_coupling (nanoS).
C CONSTS(51) is g_syn_e_6_7 in component synaptic_coupling (nanoS).
C CONSTS(52) is g_syn_e_6_8 in component synaptic_coupling (nanoS).
C CONSTS(53) is g_syn_e_6_9 in component synaptic_coupling (nanoS).
C CONSTS(54) is g_syn_e_6_10 in component synaptic_coupling (nanoS).
C CONSTS(55) is g_syn_e_7_1 in component synaptic_coupling (nanoS).
C CONSTS(56) is g_syn_e_7_2 in component synaptic_coupling (nanoS).
C CONSTS(57) is g_syn_e_7_3 in component synaptic_coupling (nanoS).
C CONSTS(58) is g_syn_e_7_4 in component synaptic_coupling (nanoS).
C CONSTS(59) is g_syn_e_7_5 in component synaptic_coupling (nanoS).
C CONSTS(60) is g_syn_e_7_6 in component synaptic_coupling (nanoS).
C CONSTS(61) is g_syn_e_7_8 in component synaptic_coupling (nanoS).
C CONSTS(62) is g_syn_e_7_9 in component synaptic_coupling (nanoS).
C CONSTS(63) is g_syn_e_7_10 in component synaptic_coupling (nanoS).
C CONSTS(64) is g_syn_e_8_1 in component synaptic_coupling (nanoS).
C CONSTS(65) is g_syn_e_8_2 in component synaptic_coupling (nanoS).
C CONSTS(66) is g_syn_e_8_3 in component synaptic_coupling (nanoS).
C CONSTS(67) is g_syn_e_8_4 in component synaptic_coupling (nanoS).
C CONSTS(68) is g_syn_e_8_5 in component synaptic_coupling (nanoS).
C CONSTS(69) is g_syn_e_8_6 in component synaptic_coupling (nanoS).
C CONSTS(70) is g_syn_e_8_7 in component synaptic_coupling (nanoS).
C CONSTS(71) is g_syn_e_8_9 in component synaptic_coupling (nanoS).
C CONSTS(72) is g_syn_e_8_10 in component synaptic_coupling (nanoS).
C CONSTS(73) is g_syn_e_9_1 in component synaptic_coupling (nanoS).
C CONSTS(74) is g_syn_e_9_2 in component synaptic_coupling (nanoS).
C CONSTS(75) is g_syn_e_9_3 in component synaptic_coupling (nanoS).
C CONSTS(76) is g_syn_e_9_4 in component synaptic_coupling (nanoS).
C CONSTS(77) is g_syn_e_9_5 in component synaptic_coupling (nanoS).
C CONSTS(78) is g_syn_e_9_6 in component synaptic_coupling (nanoS).
C CONSTS(79) is g_syn_e_9_7 in component synaptic_coupling (nanoS).
C CONSTS(80) is g_syn_e_9_8 in component synaptic_coupling (nanoS).
C CONSTS(81) is g_syn_e_9_10 in component synaptic_coupling (nanoS).
C CONSTS(82) is g_syn_e_10_1 in component synaptic_coupling (nanoS).
C CONSTS(83) is g_syn_e_10_2 in component synaptic_coupling (nanoS).
C CONSTS(84) is g_syn_e_10_3 in component synaptic_coupling (nanoS).
C CONSTS(85) is g_syn_e_10_4 in component synaptic_coupling (nanoS).
C CONSTS(86) is g_syn_e_10_5 in component synaptic_coupling (nanoS).
C CONSTS(87) is g_syn_e_10_6 in component synaptic_coupling (nanoS).
C CONSTS(88) is g_syn_e_10_7 in component synaptic_coupling (nanoS).
C CONSTS(89) is g_syn_e_10_8 in component synaptic_coupling (nanoS).
C CONSTS(90) is g_syn_e_10_9 in component synaptic_coupling (nanoS).
C STATES(1) is s in component synaptic_input (dimensionless).
C STATES(2) is s in component synaptic_input (dimensionless).
C STATES(3) is s in component synaptic_input (dimensionless).
C STATES(4) is s in component synaptic_input (dimensionless).
C STATES(5) is s in component synaptic_input (dimensionless).
C STATES(6) is s in component synaptic_input (dimensionless).
C STATES(7) is s in component synaptic_input (dimensionless).
C STATES(8) is s in component synaptic_input (dimensionless).
C STATES(9) is s in component synaptic_input (dimensionless).
C STATES(10) is s in component synaptic_input (dimensionless).
C VOI is time in component environment (millisecond).
C STATES(11) is V in component membrane (millivolt).
C CONSTS(91) is C in component membrane (picoF).
C CONSTS(92) is i_app in component membrane (picoA).
C ALGBRC(121) is i_NaP in component persistent_sodium_current (picoA).
C ALGBRC(91) is i_Na in component fast_sodium_current (picoA).
C ALGBRC(101) is i_K in component potassium_current (picoA).
C ALGBRC(131) is i_L in component leakage_current (picoA).
C ALGBRC(141) is i_tonic_e in component tonic_current (picoA).
C ALGBRC(151) is i_syn_e in component synaptic_input (picoA).
C CONSTS(93) is E_Na in component fast_sodium_current (millivolt).
C CONSTS(94) is g_Na in component fast_sodium_current (nanoS).
C ALGBRC(51) is m_infinity in component fast_sodium_current_m_gate (dimensionless).
C STATES(12) is n in component fast_sodium_current_n_gate (dimensionless).
C CONSTS(95) is theta_m in component fast_sodium_current_m_gate (millivolt).
C CONSTS(96) is sigma_m in component fast_sodium_current_m_gate (millivolt).
C ALGBRC(2) is n_infinity in component fast_sodium_current_n_gate (dimensionless).
C ALGBRC(52) is tau_n in component fast_sodium_current_n_gate (millisecond).
C CONSTS(97) is tau_n_max in component fast_sodium_current_n_gate (millisecond).
C CONSTS(98) is theta_n in component fast_sodium_current_n_gate (millivolt).
C CONSTS(99) is sigma_n in component fast_sodium_current_n_gate (millivolt).
C CONSTS(100) is g_K in component potassium_current (nanoS).
C CONSTS(101) is E_K in component potassium_current (millivolt).
C STATES(13) is n in component potassium_current_n_gate (dimensionless).
C ALGBRC(3) is n_infinity in component potassium_current_n_gate (dimensionless).
C ALGBRC(53) is tau_n in component potassium_current_n_gate (millisecond).
C CONSTS(102) is tau_n_max in component potassium_current_n_gate (millisecond).
C CONSTS(103) is theta_n in component potassium_current_n_gate (millivolt).
C CONSTS(104) is sigma_n in component potassium_current_n_gate (millivolt).
C CONSTS(105) is g_NaP in component persistent_sodium_current (nanoS).
C ALGBRC(111) is m_infinity in component persistent_sodium_current_m_gate (dimensionless).
C STATES(14) is h in component persistent_sodium_current_h_gate (dimensionless).
C CONSTS(106) is theta_m in component persistent_sodium_current_m_gate (millivolt).
C CONSTS(107) is sigma_m in component persistent_sodium_current_m_gate (millivolt).
C ALGBRC(4) is h_infinity in component persistent_sodium_current_h_gate (dimensionless).
C ALGBRC(54) is tau_h in component persistent_sodium_current_h_gate (millisecond).
C CONSTS(108) is tau_h_max in component persistent_sodium_current_h_gate (millisecond).
C CONSTS(109) is theta_h in component persistent_sodium_current_h_gate (millivolt).
C CONSTS(110) is sigma_h in component persistent_sodium_current_h_gate (millivolt).
C CONSTS(111) is g_L in component leakage_current (nanoS).
C CONSTS(112) is E_L in component leakage_current (millivolt).
C CONSTS(113) is E_syn_e in component tonic_current (millivolt).
C CONSTS(114) is g_tonic_e in component tonic_current (nanoS).
C ALGBRC(5) is s_infinity in component synaptic_input (dimensionless).
C CONSTS(115) is kr in component synaptic_input (dimensionless).
C CONSTS(116) is tau_s in component synaptic_input (millisecond).
C CONSTS(117) is sigma_s in component synaptic_input (millivolt).
C CONSTS(118) is theta_s in component synaptic_input (millivolt).
C VOI is time in component environment (millisecond).
C STATES(15) is V in component membrane (millivolt).
C CONSTS(119) is C in component membrane (picoF).
C CONSTS(120) is i_app in component membrane (picoA).
C ALGBRC(122) is i_NaP in component persistent_sodium_current (picoA).
C ALGBRC(92) is i_Na in component fast_sodium_current (picoA).
C ALGBRC(102) is i_K in component potassium_current (picoA).
C ALGBRC(132) is i_L in component leakage_current (picoA).
C ALGBRC(142) is i_tonic_e in component tonic_current (picoA).
C ALGBRC(152) is i_syn_e in component synaptic_input (picoA).
C CONSTS(121) is E_Na in component fast_sodium_current (millivolt).
C CONSTS(122) is g_Na in component fast_sodium_current (nanoS).
C ALGBRC(55) is m_infinity in component fast_sodium_current_m_gate (dimensionless).
C STATES(16) is n in component fast_sodium_current_n_gate (dimensionless).
C CONSTS(123) is theta_m in component fast_sodium_current_m_gate (millivolt).
C CONSTS(124) is sigma_m in component fast_sodium_current_m_gate (millivolt).
C ALGBRC(7) is n_infinity in component fast_sodium_current_n_gate (dimensionless).
C ALGBRC(56) is tau_n in component fast_sodium_current_n_gate (millisecond).
C CONSTS(125) is tau_n_max in component fast_sodium_current_n_gate (millisecond).
C CONSTS(126) is theta_n in component fast_sodium_current_n_gate (millivolt).
C CONSTS(127) is sigma_n in component fast_sodium_current_n_gate (millivolt).
C CONSTS(128) is g_K in component potassium_current (nanoS).
C CONSTS(129) is E_K in component potassium_current (millivolt).
C STATES(17) is n in component potassium_current_n_gate (dimensionless).
C ALGBRC(8) is n_infinity in component potassium_current_n_gate (dimensionless).
C ALGBRC(57) is tau_n in component potassium_current_n_gate (millisecond).
C CONSTS(130) is tau_n_max in component potassium_current_n_gate (millisecond).
C CONSTS(131) is theta_n in component potassium_current_n_gate (millivolt).
C CONSTS(132) is sigma_n in component potassium_current_n_gate (millivolt).
C CONSTS(133) is g_NaP in component persistent_sodium_current (nanoS).
C ALGBRC(112) is m_infinity in component persistent_sodium_current_m_gate (dimensionless).
C STATES(18) is h in component persistent_sodium_current_h_gate (dimensionless).
C CONSTS(134) is theta_m in component persistent_sodium_current_m_gate (millivolt).
C CONSTS(135) is sigma_m in component persistent_sodium_current_m_gate (millivolt).
C ALGBRC(9) is h_infinity in component persistent_sodium_current_h_gate (dimensionless).
C ALGBRC(58) is tau_h in component persistent_sodium_current_h_gate (millisecond).
C CONSTS(136) is tau_h_max in component persistent_sodium_current_h_gate (millisecond).
C CONSTS(137) is theta_h in component persistent_sodium_current_h_gate (millivolt).
C CONSTS(138) is sigma_h in component persistent_sodium_current_h_gate (millivolt).
C CONSTS(139) is g_L in component leakage_current (nanoS).
C CONSTS(140) is E_L in component leakage_current (millivolt).
C CONSTS(141) is E_syn_e in component tonic_current (millivolt).
C CONSTS(142) is g_tonic_e in component tonic_current (nanoS).
C ALGBRC(10) is s_infinity in component synaptic_input (dimensionless).
C CONSTS(143) is kr in component synaptic_input (dimensionless).
C CONSTS(144) is tau_s in component synaptic_input (millisecond).
C CONSTS(145) is sigma_s in component synaptic_input (millivolt).
C CONSTS(146) is theta_s in component synaptic_input (millivolt).
C VOI is time in component environment (millisecond).
C STATES(19) is V in component membrane (millivolt).
C CONSTS(147) is C in component membrane (picoF).
C CONSTS(148) is i_app in component membrane (picoA).
C ALGBRC(123) is i_NaP in component persistent_sodium_current (picoA).
C ALGBRC(93) is i_Na in component fast_sodium_current (picoA).
C ALGBRC(103) is i_K in component potassium_current (picoA).
C ALGBRC(133) is i_L in component leakage_current (picoA).
C ALGBRC(143) is i_tonic_e in component tonic_current (picoA).
C ALGBRC(153) is i_syn_e in component synaptic_input (picoA).
C CONSTS(149) is E_Na in component fast_sodium_current (millivolt).
C CONSTS(150) is g_Na in component fast_sodium_current (nanoS).
C ALGBRC(59) is m_infinity in component fast_sodium_current_m_gate (dimensionless).
C STATES(20) is n in component fast_sodium_current_n_gate (dimensionless).
C CONSTS(151) is theta_m in component fast_sodium_current_m_gate (millivolt).
C CONSTS(152) is sigma_m in component fast_sodium_current_m_gate (millivolt).
C ALGBRC(12) is n_infinity in component fast_sodium_current_n_gate (dimensionless).
C ALGBRC(60) is tau_n in component fast_sodium_current_n_gate (millisecond).
C CONSTS(153) is tau_n_max in component fast_sodium_current_n_gate (millisecond).
C CONSTS(154) is theta_n in component fast_sodium_current_n_gate (millivolt).
C CONSTS(155) is sigma_n in component fast_sodium_current_n_gate (millivolt).
C CONSTS(156) is g_K in component potassium_current (nanoS).
C CONSTS(157) is E_K in component potassium_current (millivolt).
C STATES(21) is n in component potassium_current_n_gate (dimensionless).
C ALGBRC(13) is n_infinity in component potassium_current_n_gate (dimensionless).
C ALGBRC(61) is tau_n in component potassium_current_n_gate (millisecond).
C CONSTS(158) is tau_n_max in component potassium_current_n_gate (millisecond).
C CONSTS(159) is theta_n in component potassium_current_n_gate (millivolt).
C CONSTS(160) is sigma_n in component potassium_current_n_gate (millivolt).
C CONSTS(161) is g_NaP in component persistent_sodium_current (nanoS).
C ALGBRC(113) is m_infinity in component persistent_sodium_current_m_gate (dimensionless).
C STATES(22) is h in component persistent_sodium_current_h_gate (dimensionless).
C CONSTS(162) is theta_m in component persistent_sodium_current_m_gate (millivolt).
C CONSTS(163) is sigma_m in component persistent_sodium_current_m_gate (millivolt).
C ALGBRC(14) is h_infinity in component persistent_sodium_current_h_gate (dimensionless).
C ALGBRC(62) is tau_h in component persistent_sodium_current_h_gate (millisecond).
C CONSTS(164) is tau_h_max in component persistent_sodium_current_h_gate (millisecond).
C CONSTS(165) is theta_h in component persistent_sodium_current_h_gate (millivolt).
C CONSTS(166) is sigma_h in component persistent_sodium_current_h_gate (millivolt).
C CONSTS(167) is g_L in component leakage_current (nanoS).
C CONSTS(168) is E_L in component leakage_current (millivolt).
C CONSTS(169) is E_syn_e in component tonic_current (millivolt).
C CONSTS(170) is g_tonic_e in component tonic_current (nanoS).
C ALGBRC(15) is s_infinity in component synaptic_input (dimensionless).
C CONSTS(171) is kr in component synaptic_input (dimensionless).
C CONSTS(172) is tau_s in component synaptic_input (millisecond).
C CONSTS(173) is sigma_s in component synaptic_input (millivolt).
C CONSTS(174) is theta_s in component synaptic_input (millivolt).
C VOI is time in component environment (millisecond).
C STATES(23) is V in component membrane (millivolt).
C CONSTS(175) is C in component membrane (picoF).
C CONSTS(176) is i_app in component membrane (picoA).
C ALGBRC(124) is i_NaP in component persistent_sodium_current (picoA).
C ALGBRC(94) is i_Na in component fast_sodium_current (picoA).
C ALGBRC(104) is i_K in component potassium_current (picoA).
C ALGBRC(134) is i_L in component leakage_current (picoA).
C ALGBRC(144) is i_tonic_e in component tonic_current (picoA).
C ALGBRC(154) is i_syn_e in component synaptic_input (picoA).
C CONSTS(177) is E_Na in component fast_sodium_current (millivolt).
C CONSTS(178) is g_Na in component fast_sodium_current (nanoS).
C ALGBRC(63) is m_infinity in component fast_sodium_current_m_gate (dimensionless).
C STATES(24) is n in component fast_sodium_current_n_gate (dimensionless).
C CONSTS(179) is theta_m in component fast_sodium_current_m_gate (millivolt).
C CONSTS(180) is sigma_m in component fast_sodium_current_m_gate (millivolt).
C ALGBRC(17) is n_infinity in component fast_sodium_current_n_gate (dimensionless).
C ALGBRC(64) is tau_n in component fast_sodium_current_n_gate (millisecond).
C CONSTS(181) is tau_n_max in component fast_sodium_current_n_gate (millisecond).
C CONSTS(182) is theta_n in component fast_sodium_current_n_gate (millivolt).
C CONSTS(183) is sigma_n in component fast_sodium_current_n_gate (millivolt).
C CONSTS(184) is g_K in component potassium_current (nanoS).
C CONSTS(185) is E_K in component potassium_current (millivolt).
C STATES(25) is n in component potassium_current_n_gate (dimensionless).
C ALGBRC(18) is n_infinity in component potassium_current_n_gate (dimensionless).
C ALGBRC(65) is tau_n in component potassium_current_n_gate (millisecond).
C CONSTS(186) is tau_n_max in component potassium_current_n_gate (millisecond).
C CONSTS(187) is theta_n in component potassium_current_n_gate (millivolt).
C CONSTS(188) is sigma_n in component potassium_current_n_gate (millivolt).
C CONSTS(189) is g_NaP in component persistent_sodium_current (nanoS).
C ALGBRC(114) is m_infinity in component persistent_sodium_current_m_gate (dimensionless).
C STATES(26) is h in component persistent_sodium_current_h_gate (dimensionless).
C CONSTS(190) is theta_m in component persistent_sodium_current_m_gate (millivolt).
C CONSTS(191) is sigma_m in component persistent_sodium_current_m_gate (millivolt).
C ALGBRC(19) is h_infinity in component persistent_sodium_current_h_gate (dimensionless).
C ALGBRC(66) is tau_h in component persistent_sodium_current_h_gate (millisecond).
C CONSTS(192) is tau_h_max in component persistent_sodium_current_h_gate (millisecond).
C CONSTS(193) is theta_h in component persistent_sodium_current_h_gate (millivolt).
C CONSTS(194) is sigma_h in component persistent_sodium_current_h_gate (millivolt).
C CONSTS(195) is g_L in component leakage_current (nanoS).
C CONSTS(196) is E_L in component leakage_current (millivolt).
C CONSTS(197) is E_syn_e in component tonic_current (millivolt).
C CONSTS(198) is g_tonic_e in component tonic_current (nanoS).
C ALGBRC(20) is s_infinity in component synaptic_input (dimensionless).
C CONSTS(199) is kr in component synaptic_input (dimensionless).
C CONSTS(200) is tau_s in component synaptic_input (millisecond).
C CONSTS(201) is sigma_s in component synaptic_input (millivolt).
C CONSTS(202) is theta_s in component synaptic_input (millivolt).
C VOI is time in component environment (millisecond).
C STATES(27) is V in component membrane (millivolt).
C CONSTS(203) is C in component membrane (picoF).
C CONSTS(204) is i_app in component membrane (picoA).
C ALGBRC(125) is i_NaP in component persistent_sodium_current (picoA).
C ALGBRC(95) is i_Na in component fast_sodium_current (picoA).
C ALGBRC(105) is i_K in component potassium_current (picoA).
C ALGBRC(135) is i_L in component leakage_current (picoA).
C ALGBRC(145) is i_tonic_e in component tonic_current (picoA).
C ALGBRC(155) is i_syn_e in component synaptic_input (picoA).
C CONSTS(205) is E_Na in component fast_sodium_current (millivolt).
C CONSTS(206) is g_Na in component fast_sodium_current (nanoS).
C ALGBRC(67) is m_infinity in component fast_sodium_current_m_gate (dimensionless).
C STATES(28) is n in component fast_sodium_current_n_gate (dimensionless).
C CONSTS(207) is theta_m in component fast_sodium_current_m_gate (millivolt).
C CONSTS(208) is sigma_m in component fast_sodium_current_m_gate (millivolt).
C ALGBRC(22) is n_infinity in component fast_sodium_current_n_gate (dimensionless).
C ALGBRC(68) is tau_n in component fast_sodium_current_n_gate (millisecond).
C CONSTS(209) is tau_n_max in component fast_sodium_current_n_gate (millisecond).
C CONSTS(210) is theta_n in component fast_sodium_current_n_gate (millivolt).
C CONSTS(211) is sigma_n in component fast_sodium_current_n_gate (millivolt).
C CONSTS(212) is g_K in component potassium_current (nanoS).
C CONSTS(213) is E_K in component potassium_current (millivolt).
C STATES(29) is n in component potassium_current_n_gate (dimensionless).
C ALGBRC(23) is n_infinity in component potassium_current_n_gate (dimensionless).
C ALGBRC(69) is tau_n in component potassium_current_n_gate (millisecond).
C CONSTS(214) is tau_n_max in component potassium_current_n_gate (millisecond).
C CONSTS(215) is theta_n in component potassium_current_n_gate (millivolt).
C CONSTS(216) is sigma_n in component potassium_current_n_gate (millivolt).
C CONSTS(217) is g_NaP in component persistent_sodium_current (nanoS).
C ALGBRC(115) is m_infinity in component persistent_sodium_current_m_gate (dimensionless).
C STATES(30) is h in component persistent_sodium_current_h_gate (dimensionless).
C CONSTS(218) is theta_m in component persistent_sodium_current_m_gate (millivolt).
C CONSTS(219) is sigma_m in component persistent_sodium_current_m_gate (millivolt).
C ALGBRC(24) is h_infinity in component persistent_sodium_current_h_gate (dimensionless).
C ALGBRC(70) is tau_h in component persistent_sodium_current_h_gate (millisecond).
C CONSTS(220) is tau_h_max in component persistent_sodium_current_h_gate (millisecond).
C CONSTS(221) is theta_h in component persistent_sodium_current_h_gate (millivolt).
C CONSTS(222) is sigma_h in component persistent_sodium_current_h_gate (millivolt).
C CONSTS(223) is g_L in component leakage_current (nanoS).
C CONSTS(224) is E_L in component leakage_current (millivolt).
C CONSTS(225) is E_syn_e in component tonic_current (millivolt).
C CONSTS(226) is g_tonic_e in component tonic_current (nanoS).
C ALGBRC(25) is s_infinity in component synaptic_input (dimensionless).
C CONSTS(227) is kr in component synaptic_input (dimensionless).
C CONSTS(228) is tau_s in component synaptic_input (millisecond).
C CONSTS(229) is sigma_s in component synaptic_input (millivolt).
C CONSTS(230) is theta_s in component synaptic_input (millivolt).
C VOI is time in component environment (millisecond).
C STATES(31) is V in component membrane (millivolt).
C CONSTS(231) is C in component membrane (picoF).
C CONSTS(232) is i_app in component membrane (picoA).
C ALGBRC(126) is i_NaP in component persistent_sodium_current (picoA).
C ALGBRC(96) is i_Na in component fast_sodium_current (picoA).
C ALGBRC(106) is i_K in component potassium_current (picoA).
C ALGBRC(136) is i_L in component leakage_current (picoA).
C ALGBRC(146) is i_tonic_e in component tonic_current (picoA).
C ALGBRC(156) is i_syn_e in component synaptic_input (picoA).
C CONSTS(233) is E_Na in component fast_sodium_current (millivolt).
C CONSTS(234) is g_Na in component fast_sodium_current (nanoS).
C ALGBRC(71) is m_infinity in component fast_sodium_current_m_gate (dimensionless).
C STATES(32) is n in component fast_sodium_current_n_gate (dimensionless).
C CONSTS(235) is theta_m in component fast_sodium_current_m_gate (millivolt).
C CONSTS(236) is sigma_m in component fast_sodium_current_m_gate (millivolt).
C ALGBRC(27) is n_infinity in component fast_sodium_current_n_gate (dimensionless).
C ALGBRC(72) is tau_n in component fast_sodium_current_n_gate (millisecond).
C CONSTS(237) is tau_n_max in component fast_sodium_current_n_gate (millisecond).
C CONSTS(238) is theta_n in component fast_sodium_current_n_gate (millivolt).
C CONSTS(239) is sigma_n in component fast_sodium_current_n_gate (millivolt).
C CONSTS(240) is g_K in component potassium_current (nanoS).
C CONSTS(241) is E_K in component potassium_current (millivolt).
C STATES(33) is n in component potassium_current_n_gate (dimensionless).
C ALGBRC(28) is n_infinity in component potassium_current_n_gate (dimensionless).
C ALGBRC(73) is tau_n in component potassium_current_n_gate (millisecond).
C CONSTS(242) is tau_n_max in component potassium_current_n_gate (millisecond).
C CONSTS(243) is theta_n in component potassium_current_n_gate (millivolt).
C CONSTS(244) is sigma_n in component potassium_current_n_gate (millivolt).
C CONSTS(245) is g_NaP in component persistent_sodium_current (nanoS).
C ALGBRC(116) is m_infinity in component persistent_sodium_current_m_gate (dimensionless).
C STATES(34) is h in component persistent_sodium_current_h_gate (dimensionless).
C CONSTS(246) is theta_m in component persistent_sodium_current_m_gate (millivolt).
C CONSTS(247) is sigma_m in component persistent_sodium_current_m_gate (millivolt).
C ALGBRC(29) is h_infinity in component persistent_sodium_current_h_gate (dimensionless).
C ALGBRC(74) is tau_h in component persistent_sodium_current_h_gate (millisecond).
C CONSTS(248) is tau_h_max in component persistent_sodium_current_h_gate (millisecond).
C CONSTS(249) is theta_h in component persistent_sodium_current_h_gate (millivolt).
C CONSTS(250) is sigma_h in component persistent_sodium_current_h_gate (millivolt).
C CONSTS(251) is g_L in component leakage_current (nanoS).
C CONSTS(252) is E_L in component leakage_current (millivolt).
C CONSTS(253) is E_syn_e in component tonic_current (millivolt).
C CONSTS(254) is g_tonic_e in component tonic_current (nanoS).
C ALGBRC(30) is s_infinity in component synaptic_input (dimensionless).
C CONSTS(255) is kr in component synaptic_input (dimensionless).
C CONSTS(256) is tau_s in component synaptic_input (millisecond).
C CONSTS(257) is sigma_s in component synaptic_input (millivolt).
C CONSTS(258) is theta_s in component synaptic_input (millivolt).
C VOI is time in component environment (millisecond).
C STATES(35) is V in component membrane (millivolt).
C CONSTS(259) is C in component membrane (picoF).
C CONSTS(260) is i_app in component membrane (picoA).
C ALGBRC(127) is i_NaP in component persistent_sodium_current (picoA).
C ALGBRC(97) is i_Na in component fast_sodium_current (picoA).
C ALGBRC(107) is i_K in component potassium_current (picoA).
C ALGBRC(137) is i_L in component leakage_current (picoA).
C ALGBRC(147) is i_tonic_e in component tonic_current (picoA).
C ALGBRC(157) is i_syn_e in component synaptic_input (picoA).
C CONSTS(261) is E_Na in component fast_sodium_current (millivolt).
C CONSTS(262) is g_Na in component fast_sodium_current (nanoS).
C ALGBRC(75) is m_infinity in component fast_sodium_current_m_gate (dimensionless).
C STATES(36) is n in component fast_sodium_current_n_gate (dimensionless).
C CONSTS(263) is theta_m in component fast_sodium_current_m_gate (millivolt).
C CONSTS(264) is sigma_m in component fast_sodium_current_m_gate (millivolt).
C ALGBRC(32) is n_infinity in component fast_sodium_current_n_gate (dimensionless).
C ALGBRC(76) is tau_n in component fast_sodium_current_n_gate (millisecond).
C CONSTS(265) is tau_n_max in component fast_sodium_current_n_gate (millisecond).
C CONSTS(266) is theta_n in component fast_sodium_current_n_gate (millivolt).
C CONSTS(267) is sigma_n in component fast_sodium_current_n_gate (millivolt).
C CONSTS(268) is g_K in component potassium_current (nanoS).
C CONSTS(269) is E_K in component potassium_current (millivolt).
C STATES(37) is n in component potassium_current_n_gate (dimensionless).
C ALGBRC(33) is n_infinity in component potassium_current_n_gate (dimensionless).
C ALGBRC(77) is tau_n in component potassium_current_n_gate (millisecond).
C CONSTS(270) is tau_n_max in component potassium_current_n_gate (millisecond).
C CONSTS(271) is theta_n in component potassium_current_n_gate (millivolt).
C CONSTS(272) is sigma_n in component potassium_current_n_gate (millivolt).
C CONSTS(273) is g_NaP in component persistent_sodium_current (nanoS).
C ALGBRC(117) is m_infinity in component persistent_sodium_current_m_gate (dimensionless).
C STATES(38) is h in component persistent_sodium_current_h_gate (dimensionless).
C CONSTS(274) is theta_m in component persistent_sodium_current_m_gate (millivolt).
C CONSTS(275) is sigma_m in component persistent_sodium_current_m_gate (millivolt).
C ALGBRC(34) is h_infinity in component persistent_sodium_current_h_gate (dimensionless).
C ALGBRC(78) is tau_h in component persistent_sodium_current_h_gate (millisecond).
C CONSTS(276) is tau_h_max in component persistent_sodium_current_h_gate (millisecond).
C CONSTS(277) is theta_h in component persistent_sodium_current_h_gate (millivolt).
C CONSTS(278) is sigma_h in component persistent_sodium_current_h_gate (millivolt).
C CONSTS(279) is g_L in component leakage_current (nanoS).
C CONSTS(280) is E_L in component leakage_current (millivolt).
C CONSTS(281) is E_syn_e in component tonic_current (millivolt).
C CONSTS(282) is g_tonic_e in component tonic_current (nanoS).
C ALGBRC(35) is s_infinity in component synaptic_input (dimensionless).
C CONSTS(283) is kr in component synaptic_input (dimensionless).
C CONSTS(284) is tau_s in component synaptic_input (millisecond).
C CONSTS(285) is sigma_s in component synaptic_input (millivolt).
C CONSTS(286) is theta_s in component synaptic_input (millivolt).
C VOI is time in component environment (millisecond).
C STATES(39) is V in component membrane (millivolt).
C CONSTS(287) is C in component membrane (picoF).
C CONSTS(288) is i_app in component membrane (picoA).
C ALGBRC(128) is i_NaP in component persistent_sodium_current (picoA).
C ALGBRC(98) is i_Na in component fast_sodium_current (picoA).
C ALGBRC(108) is i_K in component potassium_current (picoA).
C ALGBRC(138) is i_L in component leakage_current (picoA).
C ALGBRC(148) is i_tonic_e in component tonic_current (picoA).
C ALGBRC(158) is i_syn_e in component synaptic_input (picoA).
C CONSTS(289) is E_Na in component fast_sodium_current (millivolt).
C CONSTS(290) is g_Na in component fast_sodium_current (nanoS).
C ALGBRC(79) is m_infinity in component fast_sodium_current_m_gate (dimensionless).
C STATES(40) is n in component fast_sodium_current_n_gate (dimensionless).
C CONSTS(291) is theta_m in component fast_sodium_current_m_gate (millivolt).
C CONSTS(292) is sigma_m in component fast_sodium_current_m_gate (millivolt).
C ALGBRC(37) is n_infinity in component fast_sodium_current_n_gate (dimensionless).
C ALGBRC(80) is tau_n in component fast_sodium_current_n_gate (millisecond).
C CONSTS(293) is tau_n_max in component fast_sodium_current_n_gate (millisecond).
C CONSTS(294) is theta_n in component fast_sodium_current_n_gate (millivolt).
C CONSTS(295) is sigma_n in component fast_sodium_current_n_gate (millivolt).
C CONSTS(296) is g_K in component potassium_current (nanoS).
C CONSTS(297) is E_K in component potassium_current (millivolt).
C STATES(41) is n in component potassium_current_n_gate (dimensionless).
C ALGBRC(38) is n_infinity in component potassium_current_n_gate (dimensionless).
C ALGBRC(81) is tau_n in component potassium_current_n_gate (millisecond).
C CONSTS(298) is tau_n_max in component potassium_current_n_gate (millisecond).
C CONSTS(299) is theta_n in component potassium_current_n_gate (millivolt).
C CONSTS(300) is sigma_n in component potassium_current_n_gate (millivolt).
C CONSTS(301) is g_NaP in component persistent_sodium_current (nanoS).
C ALGBRC(118) is m_infinity in component persistent_sodium_current_m_gate (dimensionless).
C STATES(42) is h in component persistent_sodium_current_h_gate (dimensionless).
C CONSTS(302) is theta_m in component persistent_sodium_current_m_gate (millivolt).
C CONSTS(303) is sigma_m in component persistent_sodium_current_m_gate (millivolt).
C ALGBRC(39) is h_infinity in component persistent_sodium_current_h_gate (dimensionless).
C ALGBRC(82) is tau_h in component persistent_sodium_current_h_gate (millisecond).
C CONSTS(304) is tau_h_max in component persistent_sodium_current_h_gate (millisecond).
C CONSTS(305) is theta_h in component persistent_sodium_current_h_gate (millivolt).
C CONSTS(306) is sigma_h in component persistent_sodium_current_h_gate (millivolt).
C CONSTS(307) is g_L in component leakage_current (nanoS).
C CONSTS(308) is E_L in component leakage_current (millivolt).
C CONSTS(309) is E_syn_e in component tonic_current (millivolt).
C CONSTS(310) is g_tonic_e in component tonic_current (nanoS).
C ALGBRC(40) is s_infinity in component synaptic_input (dimensionless).
C CONSTS(311) is kr in component synaptic_input (dimensionless).
C CONSTS(312) is tau_s in component synaptic_input (millisecond).
C CONSTS(313) is sigma_s in component synaptic_input (millivolt).
C CONSTS(314) is theta_s in component synaptic_input (millivolt).
C VOI is time in component environment (millisecond).
C STATES(43) is V in component membrane (millivolt).
C CONSTS(315) is C in component membrane (picoF).
C CONSTS(316) is i_app in component membrane (picoA).
C ALGBRC(129) is i_NaP in component persistent_sodium_current (picoA).
C ALGBRC(99) is i_Na in component fast_sodium_current (picoA).
C ALGBRC(109) is i_K in component potassium_current (picoA).
C ALGBRC(139) is i_L in component leakage_current (picoA).
C ALGBRC(149) is i_tonic_e in component tonic_current (picoA).
C ALGBRC(159) is i_syn_e in component synaptic_input (picoA).
C CONSTS(317) is E_Na in component fast_sodium_current (millivolt).
C CONSTS(318) is g_Na in component fast_sodium_current (nanoS).
C ALGBRC(83) is m_infinity in component fast_sodium_current_m_gate (dimensionless).
C STATES(44) is n in component fast_sodium_current_n_gate (dimensionless).
C CONSTS(319) is theta_m in component fast_sodium_current_m_gate (millivolt).
C CONSTS(320) is sigma_m in component fast_sodium_current_m_gate (millivolt).
C ALGBRC(42) is n_infinity in component fast_sodium_current_n_gate (dimensionless).
C ALGBRC(84) is tau_n in component fast_sodium_current_n_gate (millisecond).
C CONSTS(321) is tau_n_max in component fast_sodium_current_n_gate (millisecond).
C CONSTS(322) is theta_n in component fast_sodium_current_n_gate (millivolt).
C CONSTS(323) is sigma_n in component fast_sodium_current_n_gate (millivolt).
C CONSTS(324) is g_K in component potassium_current (nanoS).
C CONSTS(325) is E_K in component potassium_current (millivolt).
C STATES(45) is n in component potassium_current_n_gate (dimensionless).
C ALGBRC(43) is n_infinity in component potassium_current_n_gate (dimensionless).
C ALGBRC(85) is tau_n in component potassium_current_n_gate (millisecond).
C CONSTS(326) is tau_n_max in component potassium_current_n_gate (millisecond).
C CONSTS(327) is theta_n in component potassium_current_n_gate (millivolt).
C CONSTS(328) is sigma_n in component potassium_current_n_gate (millivolt).
C CONSTS(329) is g_NaP in component persistent_sodium_current (nanoS).
C ALGBRC(119) is m_infinity in component persistent_sodium_current_m_gate (dimensionless).
C STATES(46) is h in component persistent_sodium_current_h_gate (dimensionless).
C CONSTS(330) is theta_m in component persistent_sodium_current_m_gate (millivolt).
C CONSTS(331) is sigma_m in component persistent_sodium_current_m_gate (millivolt).
C ALGBRC(44) is h_infinity in component persistent_sodium_current_h_gate (dimensionless).
C ALGBRC(86) is tau_h in component persistent_sodium_current_h_gate (millisecond).
C CONSTS(332) is tau_h_max in component persistent_sodium_current_h_gate (millisecond).
C CONSTS(333) is theta_h in component persistent_sodium_current_h_gate (millivolt).
C CONSTS(334) is sigma_h in component persistent_sodium_current_h_gate (millivolt).
C CONSTS(335) is g_L in component leakage_current (nanoS).
C CONSTS(336) is E_L in component leakage_current (millivolt).
C CONSTS(337) is E_syn_e in component tonic_current (millivolt).
C CONSTS(338) is g_tonic_e in component tonic_current (nanoS).
C ALGBRC(45) is s_infinity in component synaptic_input (dimensionless).
C CONSTS(339) is kr in component synaptic_input (dimensionless).
C CONSTS(340) is tau_s in component synaptic_input (millisecond).
C CONSTS(341) is sigma_s in component synaptic_input (millivolt).
C CONSTS(342) is theta_s in component synaptic_input (millivolt).
C VOI is time in component environment (millisecond).
C STATES(47) is V in component membrane (millivolt).
C CONSTS(343) is C in component membrane (picoF).
C CONSTS(344) is i_app in component membrane (picoA).
C ALGBRC(130) is i_NaP in component persistent_sodium_current (picoA).
C ALGBRC(100) is i_Na in component fast_sodium_current (picoA).
C ALGBRC(110) is i_K in component potassium_current (picoA).
C ALGBRC(140) is i_L in component leakage_current (picoA).
C ALGBRC(150) is i_tonic_e in component tonic_current (picoA).
C ALGBRC(160) is i_syn_e in component synaptic_input (picoA).
C CONSTS(345) is E_Na in component fast_sodium_current (millivolt).
C CONSTS(346) is g_Na in component fast_sodium_current (nanoS).
C ALGBRC(87) is m_infinity in component fast_sodium_current_m_gate (dimensionless).
C STATES(48) is n in component fast_sodium_current_n_gate (dimensionless).
C CONSTS(347) is theta_m in component fast_sodium_current_m_gate (millivolt).
C CONSTS(348) is sigma_m in component fast_sodium_current_m_gate (millivolt).
C ALGBRC(47) is n_infinity in component fast_sodium_current_n_gate (dimensionless).
C ALGBRC(88) is tau_n in component fast_sodium_current_n_gate (millisecond).
C CONSTS(349) is tau_n_max in component fast_sodium_current_n_gate (millisecond).
C CONSTS(350) is theta_n in component fast_sodium_current_n_gate (millivolt).
C CONSTS(351) is sigma_n in component fast_sodium_current_n_gate (millivolt).
C CONSTS(352) is g_K in component potassium_current (nanoS).
C CONSTS(353) is E_K in component potassium_current (millivolt).
C STATES(49) is n in component potassium_current_n_gate (dimensionless).
C ALGBRC(48) is n_infinity in component potassium_current_n_gate (dimensionless).
C ALGBRC(89) is tau_n in component potassium_current_n_gate (millisecond).
C CONSTS(354) is tau_n_max in component potassium_current_n_gate (millisecond).
C CONSTS(355) is theta_n in component potassium_current_n_gate (millivolt).
C CONSTS(356) is sigma_n in component potassium_current_n_gate (millivolt).
C CONSTS(357) is g_NaP in component persistent_sodium_current (nanoS).
C ALGBRC(120) is m_infinity in component persistent_sodium_current_m_gate (dimensionless).
C STATES(50) is h in component persistent_sodium_current_h_gate (dimensionless).
C CONSTS(358) is theta_m in component persistent_sodium_current_m_gate (millivolt).
C CONSTS(359) is sigma_m in component persistent_sodium_current_m_gate (millivolt).
C ALGBRC(49) is h_infinity in component persistent_sodium_current_h_gate (dimensionless).
C ALGBRC(90) is tau_h in component persistent_sodium_current_h_gate (millisecond).
C CONSTS(360) is tau_h_max in component persistent_sodium_current_h_gate (millisecond).
C CONSTS(361) is theta_h in component persistent_sodium_current_h_gate (millivolt).
C CONSTS(362) is sigma_h in component persistent_sodium_current_h_gate (millivolt).
C CONSTS(363) is g_L in component leakage_current (nanoS).
C CONSTS(364) is E_L in component leakage_current (millivolt).
C CONSTS(365) is E_syn_e in component tonic_current (millivolt).
C CONSTS(366) is g_tonic_e in component tonic_current (nanoS).
C ALGBRC(50) is s_infinity in component synaptic_input (dimensionless).
C CONSTS(367) is kr in component synaptic_input (dimensionless).
C CONSTS(368) is tau_s in component synaptic_input (millisecond).
C CONSTS(369) is sigma_s in component synaptic_input (millivolt).
C CONSTS(370) is theta_s in component synaptic_input (millivolt).
C RATES(11) is d/dt V in component membrane (millivolt).
C RATES(12) is d/dt n in component fast_sodium_current_n_gate (dimensionless).
C RATES(13) is d/dt n in component potassium_current_n_gate (dimensionless).
C RATES(14) is d/dt h in component persistent_sodium_current_h_gate (dimensionless).
C RATES(1) is d/dt s in component synaptic_input (dimensionless).
C RATES(15) is d/dt V in component membrane (millivolt).
C RATES(16) is d/dt n in component fast_sodium_current_n_gate (dimensionless).
C RATES(17) is d/dt n in component potassium_current_n_gate (dimensionless).
C RATES(18) is d/dt h in component persistent_sodium_current_h_gate (dimensionless).
C RATES(2) is d/dt s in component synaptic_input (dimensionless).
C RATES(19) is d/dt V in component membrane (millivolt).
C RATES(20) is d/dt n in component fast_sodium_current_n_gate (dimensionless).
C RATES(21) is d/dt n in component potassium_current_n_gate (dimensionless).
C RATES(22) is d/dt h in component persistent_sodium_current_h_gate (dimensionless).
C RATES(3) is d/dt s in component synaptic_input (dimensionless).
C RATES(23) is d/dt V in component membrane (millivolt).
C RATES(24) is d/dt n in component fast_sodium_current_n_gate (dimensionless).
C RATES(25) is d/dt n in component potassium_current_n_gate (dimensionless).
C RATES(26) is d/dt h in component persistent_sodium_current_h_gate (dimensionless).
C RATES(4) is d/dt s in component synaptic_input (dimensionless).
C RATES(27) is d/dt V in component membrane (millivolt).
C RATES(28) is d/dt n in component fast_sodium_current_n_gate (dimensionless).
C RATES(29) is d/dt n in component potassium_current_n_gate (dimensionless).
C RATES(30) is d/dt h in component persistent_sodium_current_h_gate (dimensionless).
C RATES(5) is d/dt s in component synaptic_input (dimensionless).
C RATES(31) is d/dt V in component membrane (millivolt).
C RATES(32) is d/dt n in component fast_sodium_current_n_gate (dimensionless).
C RATES(33) is d/dt n in component potassium_current_n_gate (dimensionless).
C RATES(34) is d/dt h in component persistent_sodium_current_h_gate (dimensionless).
C RATES(6) is d/dt s in component synaptic_input (dimensionless).
C RATES(35) is d/dt V in component membrane (millivolt).
C RATES(36) is d/dt n in component fast_sodium_current_n_gate (dimensionless).
C RATES(37) is d/dt n in component potassium_current_n_gate (dimensionless).
C RATES(38) is d/dt h in component persistent_sodium_current_h_gate (dimensionless).
C RATES(7) is d/dt s in component synaptic_input (dimensionless).
C RATES(39) is d/dt V in component membrane (millivolt).
C RATES(40) is d/dt n in component fast_sodium_current_n_gate (dimensionless).
C RATES(41) is d/dt n in component potassium_current_n_gate (dimensionless).
C RATES(42) is d/dt h in component persistent_sodium_current_h_gate (dimensionless).
C RATES(8) is d/dt s in component synaptic_input (dimensionless).
C RATES(43) is d/dt V in component membrane (millivolt).
C RATES(44) is d/dt n in component fast_sodium_current_n_gate (dimensionless).
C RATES(45) is d/dt n in component potassium_current_n_gate (dimensionless).
C RATES(46) is d/dt h in component persistent_sodium_current_h_gate (dimensionless).
C RATES(9) is d/dt s in component synaptic_input (dimensionless).
C RATES(47) is d/dt V in component membrane (millivolt).
C RATES(48) is d/dt n in component fast_sodium_current_n_gate (dimensionless).
C RATES(49) is d/dt n in component potassium_current_n_gate (dimensionless).
C RATES(50) is d/dt h in component persistent_sodium_current_h_gate (dimensionless).
C RATES(10) is d/dt s in component synaptic_input (dimensionless).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
CONSTS(1) = 0.10
CONSTS(2) = 0.10
CONSTS(3) = 0.10
CONSTS(4) = 0.10
CONSTS(5) = 0.10
CONSTS(6) = 0.10
CONSTS(7) = 0.10
CONSTS(8) = 0.10
CONSTS(9) = 0.10
CONSTS(10) = 0.10
CONSTS(11) = 0.10
CONSTS(12) = 0.10
CONSTS(13) = 0.10
CONSTS(14) = 0.10
CONSTS(15) = 0.10
CONSTS(16) = 0.10
CONSTS(17) = 0.10
CONSTS(18) = 0.10
CONSTS(19) = 0.10
CONSTS(20) = 0.10
CONSTS(21) = 0.10
CONSTS(22) = 0.10
CONSTS(23) = 0.10
CONSTS(24) = 0.10
CONSTS(25) = 0.10
CONSTS(26) = 0.10
CONSTS(27) = 0.10
CONSTS(28) = 0.10
CONSTS(29) = 0.10
CONSTS(30) = 0.10
CONSTS(31) = 0.10
CONSTS(32) = 0.10
CONSTS(33) = 0.10
CONSTS(34) = 0.10
CONSTS(35) = 0.10
CONSTS(36) = 0.10
CONSTS(37) = 0.10
CONSTS(38) = 0.10
CONSTS(39) = 0.10
CONSTS(40) = 0.10
CONSTS(41) = 0.10
CONSTS(42) = 0.10
CONSTS(43) = 0.10
CONSTS(44) = 0.10
CONSTS(45) = 0.10
CONSTS(46) = 0.10
CONSTS(47) = 0.10
CONSTS(48) = 0.10
CONSTS(49) = 0.10
CONSTS(50) = 0.10
CONSTS(51) = 0.10
CONSTS(52) = 0.10
CONSTS(53) = 0.10
CONSTS(54) = 0.10
CONSTS(55) = 0.10
CONSTS(56) = 0.10
CONSTS(57) = 0.10
CONSTS(58) = 0.10
CONSTS(59) = 0.10
CONSTS(60) = 0.10
CONSTS(61) = 0.10
CONSTS(62) = 0.10
CONSTS(63) = 0.10
CONSTS(64) = 0.10
CONSTS(65) = 0.10
CONSTS(66) = 0.10
CONSTS(67) = 0.10
CONSTS(68) = 0.10
CONSTS(69) = 0.10
CONSTS(70) = 0.10
CONSTS(71) = 0.10
CONSTS(72) = 0.10
CONSTS(73) = 0.10
CONSTS(74) = 0.10
CONSTS(75) = 0.10
CONSTS(76) = 0.10
CONSTS(77) = 0.10
CONSTS(78) = 0.10
CONSTS(79) = 0.10
CONSTS(80) = 0.10
CONSTS(81) = 0.10
CONSTS(82) = 0.10
CONSTS(83) = 0.10
CONSTS(84) = 0.10
CONSTS(85) = 0.10
CONSTS(86) = 0.10
CONSTS(87) = 0.10
CONSTS(88) = 0.10
CONSTS(89) = 0.10
CONSTS(90) = 0.10
STATES(1) = 1
STATES(2) = 1
STATES(3) = 1
STATES(4) = 1
STATES(5) = 1
STATES(6) = 1
STATES(7) = 1
STATES(8) = 1
STATES(9) = 1
STATES(10) = 1
STATES(11) = -50.0
CONSTS(91) = 21.0
CONSTS(92) = 0.0
CONSTS(93) = 50.0
CONSTS(94) = 28.0
STATES(12) = 0.01
CONSTS(95) = -34.0
CONSTS(96) = -5.0
CONSTS(97) = 10.0
CONSTS(98) = -29.0
CONSTS(99) = -4.0
CONSTS(100) = 11.2
CONSTS(101) = -85.0
STATES(13) = 0.01
CONSTS(102) = 10.0
CONSTS(103) = -29.0
CONSTS(104) = -4.0
CONSTS(105) = 2.8
STATES(14) = 0.46
CONSTS(106) = -40.0
CONSTS(107) = -6.0
CONSTS(108) = 10000.0
CONSTS(109) = -48.0
CONSTS(110) = 6.0
CONSTS(111) = 2.8
CONSTS(112) = -57.5
CONSTS(113) = 0.0
CONSTS(114) = 0.0
CONSTS(115) = 1.0
CONSTS(116) = 5.0
CONSTS(117) = -5.0
CONSTS(118) = -10.0
STATES(15) = -50.0
CONSTS(119) = 21.0
CONSTS(120) = 0.0
CONSTS(121) = 50.0
CONSTS(122) = 28.0
STATES(16) = 0.01
CONSTS(123) = -34.0
CONSTS(124) = -5.0
CONSTS(125) = 10.0
CONSTS(126) = -29.0
CONSTS(127) = -4.0
CONSTS(128) = 11.2
CONSTS(129) = -85.0
STATES(17) = 0.01
CONSTS(130) = 10.0
CONSTS(131) = -29.0
CONSTS(132) = -4.0
CONSTS(133) = 2.8
STATES(18) = 0.46
CONSTS(134) = -40.0
CONSTS(135) = -6.0
CONSTS(136) = 10000.0
CONSTS(137) = -48.0
CONSTS(138) = 6.0
CONSTS(139) = 2.8
CONSTS(140) = -57.5
CONSTS(141) = 0.0
CONSTS(142) = 0.0
CONSTS(143) = 1.0
CONSTS(144) = 5.0
CONSTS(145) = -5.0
CONSTS(146) = -10.0
STATES(19) = -50.0
CONSTS(147) = 21.0
CONSTS(148) = 0.0
CONSTS(149) = 50.0
CONSTS(150) = 28.0
STATES(20) = 0.01
CONSTS(151) = -34.0
CONSTS(152) = -5.0
CONSTS(153) = 10.0
CONSTS(154) = -29.0
CONSTS(155) = -4.0
CONSTS(156) = 11.2
CONSTS(157) = -85.0
STATES(21) = 0.01
CONSTS(158) = 10.0
CONSTS(159) = -29.0
CONSTS(160) = -4.0
CONSTS(161) = 2.8
STATES(22) = 0.46
CONSTS(162) = -40.0
CONSTS(163) = -6.0
CONSTS(164) = 10000.0
CONSTS(165) = -48.0
CONSTS(166) = 6.0
CONSTS(167) = 2.8
CONSTS(168) = -57.5
CONSTS(169) = 0.0
CONSTS(170) = 0.0
CONSTS(171) = 1.0
CONSTS(172) = 5.0
CONSTS(173) = -5.0
CONSTS(174) = -10.0
STATES(23) = -50.0
CONSTS(175) = 21.0
CONSTS(176) = 0.0
CONSTS(177) = 50.0
CONSTS(178) = 28.0
STATES(24) = 0.01
CONSTS(179) = -34.0
CONSTS(180) = -5.0
CONSTS(181) = 10.0
CONSTS(182) = -29.0
CONSTS(183) = -4.0
CONSTS(184) = 11.2
CONSTS(185) = -85.0
STATES(25) = 0.01
CONSTS(186) = 10.0
CONSTS(187) = -29.0
CONSTS(188) = -4.0
CONSTS(189) = 2.8
STATES(26) = 0.46
CONSTS(190) = -40.0
CONSTS(191) = -6.0
CONSTS(192) = 10000.0
CONSTS(193) = -48.0
CONSTS(194) = 6.0
CONSTS(195) = 2.8
CONSTS(196) = -57.5
CONSTS(197) = 0.0
CONSTS(198) = 0.0
CONSTS(199) = 1.0
CONSTS(200) = 5.0
CONSTS(201) = -5.0
CONSTS(202) = -10.0
STATES(27) = -50.0
CONSTS(203) = 21.0
CONSTS(204) = 0.0
CONSTS(205) = 50.0
CONSTS(206) = 28.0
STATES(28) = 0.01
CONSTS(207) = -34.0
CONSTS(208) = -5.0
CONSTS(209) = 10.0
CONSTS(210) = -29.0
CONSTS(211) = -4.0
CONSTS(212) = 11.2
CONSTS(213) = -85.0
STATES(29) = 0.01
CONSTS(214) = 10.0
CONSTS(215) = -29.0
CONSTS(216) = -4.0
CONSTS(217) = 2.8
STATES(30) = 0.46
CONSTS(218) = -40.0
CONSTS(219) = -6.0
CONSTS(220) = 10000.0
CONSTS(221) = -48.0
CONSTS(222) = 6.0
CONSTS(223) = 2.8
CONSTS(224) = -57.5
CONSTS(225) = 0.0
CONSTS(226) = 0.0
CONSTS(227) = 1.0
CONSTS(228) = 5.0
CONSTS(229) = -5.0
CONSTS(230) = -10.0
STATES(31) = -50.0
CONSTS(231) = 21.0
CONSTS(232) = 0.0
CONSTS(233) = 50.0
CONSTS(234) = 28.0
STATES(32) = 0.01
CONSTS(235) = -34.0
CONSTS(236) = -5.0
CONSTS(237) = 10.0
CONSTS(238) = -29.0
CONSTS(239) = -4.0
CONSTS(240) = 11.2
CONSTS(241) = -85.0
STATES(33) = 0.01
CONSTS(242) = 10.0
CONSTS(243) = -29.0
CONSTS(244) = -4.0
CONSTS(245) = 2.8
STATES(34) = 0.46
CONSTS(246) = -40.0
CONSTS(247) = -6.0
CONSTS(248) = 10000.0
CONSTS(249) = -48.0
CONSTS(250) = 6.0
CONSTS(251) = 2.8
CONSTS(252) = -57.5
CONSTS(253) = 0.0
CONSTS(254) = 0.0
CONSTS(255) = 1.0
CONSTS(256) = 5.0
CONSTS(257) = -5.0
CONSTS(258) = -10.0
STATES(35) = -50.0
CONSTS(259) = 21.0
CONSTS(260) = 0.0
CONSTS(261) = 50.0
CONSTS(262) = 28.0
STATES(36) = 0.01
CONSTS(263) = -34.0
CONSTS(264) = -5.0
CONSTS(265) = 10.0
CONSTS(266) = -29.0
CONSTS(267) = -4.0
CONSTS(268) = 11.2
CONSTS(269) = -85.0
STATES(37) = 0.01
CONSTS(270) = 10.0
CONSTS(271) = -29.0
CONSTS(272) = -4.0
CONSTS(273) = 2.8
STATES(38) = 0.46
CONSTS(274) = -40.0
CONSTS(275) = -6.0
CONSTS(276) = 10000.0
CONSTS(277) = -48.0
CONSTS(278) = 6.0
CONSTS(279) = 2.8
CONSTS(280) = -57.5
CONSTS(281) = 0.0
CONSTS(282) = 0.0
CONSTS(283) = 1.0
CONSTS(284) = 5.0
CONSTS(285) = -5.0
CONSTS(286) = -10.0
STATES(39) = -50.0
CONSTS(287) = 21.0
CONSTS(288) = 0.0
CONSTS(289) = 50.0
CONSTS(290) = 28.0
STATES(40) = 0.01
CONSTS(291) = -34.0
CONSTS(292) = -5.0
CONSTS(293) = 10.0
CONSTS(294) = -29.0
CONSTS(295) = -4.0
CONSTS(296) = 11.2
CONSTS(297) = -85.0
STATES(41) = 0.01
CONSTS(298) = 10.0
CONSTS(299) = -29.0
CONSTS(300) = -4.0
CONSTS(301) = 2.8
STATES(42) = 0.46
CONSTS(302) = -40.0
CONSTS(303) = -6.0
CONSTS(304) = 10000.0
CONSTS(305) = -48.0
CONSTS(306) = 6.0
CONSTS(307) = 2.8
CONSTS(308) = -57.5
CONSTS(309) = 0.0
CONSTS(310) = 0.0
CONSTS(311) = 1.0
CONSTS(312) = 5.0
CONSTS(313) = -5.0
CONSTS(314) = -10.0
STATES(43) = -50.0
CONSTS(315) = 21.0
CONSTS(316) = 0.0
CONSTS(317) = 50.0
CONSTS(318) = 28.0
STATES(44) = 0.01
CONSTS(319) = -34.0
CONSTS(320) = -5.0
CONSTS(321) = 10.0
CONSTS(322) = -29.0
CONSTS(323) = -4.0
CONSTS(324) = 11.2
CONSTS(325) = -85.0
STATES(45) = 0.01
CONSTS(326) = 10.0
CONSTS(327) = -29.0
CONSTS(328) = -4.0
CONSTS(329) = 2.8
STATES(46) = 0.46
CONSTS(330) = -40.0
CONSTS(331) = -6.0
CONSTS(332) = 10000.0
CONSTS(333) = -48.0
CONSTS(334) = 6.0
CONSTS(335) = 2.8
CONSTS(336) = -57.5
CONSTS(337) = 0.0
CONSTS(338) = 0.0
CONSTS(339) = 1.0
CONSTS(340) = 5.0
CONSTS(341) = -5.0
CONSTS(342) = -10.0
STATES(47) = -50.0
CONSTS(343) = 21.0
CONSTS(344) = 0.0
CONSTS(345) = 50.0
CONSTS(346) = 28.0
STATES(48) = 0.01
CONSTS(347) = -34.0
CONSTS(348) = -5.0
CONSTS(349) = 10.0
CONSTS(350) = -29.0
CONSTS(351) = -4.0
CONSTS(352) = 11.2
CONSTS(353) = -85.0
STATES(49) = 0.01
CONSTS(354) = 10.0
CONSTS(355) = -29.0
CONSTS(356) = -4.0
CONSTS(357) = 2.8
STATES(50) = 0.46
CONSTS(358) = -40.0
CONSTS(359) = -6.0
CONSTS(360) = 10000.0
CONSTS(361) = -48.0
CONSTS(362) = 6.0
CONSTS(363) = 2.8
CONSTS(364) = -57.5
CONSTS(365) = 0.0
CONSTS(366) = 0.0
CONSTS(367) = 1.0
CONSTS(368) = 5.0
CONSTS(369) = -5.0
CONSTS(370) = -10.0
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(5) = 1.00000/(1.00000+EXP((STATES(11) - CONSTS(118))/CONSTS(117)))
RATES(1) = ( (1.00000 - STATES(1))*ALGBRC(5) - - ( CONSTS(115)*STATES(1)))/CONSTS(116)
ALGBRC(10) = 1.00000/(1.00000+EXP((STATES(15) - CONSTS(146))/CONSTS(145)))
RATES(2) = ( (1.00000 - STATES(2))*ALGBRC(10) - - ( CONSTS(143)*STATES(2)))/CONSTS(144)
ALGBRC(15) = 1.00000/(1.00000+EXP((STATES(19) - CONSTS(174))/CONSTS(173)))
RATES(3) = ( (1.00000 - STATES(3))*ALGBRC(15) - - ( CONSTS(171)*STATES(3)))/CONSTS(172)
ALGBRC(20) = 1.00000/(1.00000+EXP((STATES(23) - CONSTS(202))/CONSTS(201)))
RATES(4) = ( (1.00000 - STATES(4))*ALGBRC(20) - - ( CONSTS(199)*STATES(4)))/CONSTS(200)
ALGBRC(25) = 1.00000/(1.00000+EXP((STATES(27) - CONSTS(230))/CONSTS(229)))
RATES(5) = ( (1.00000 - STATES(5))*ALGBRC(25) - - ( CONSTS(227)*STATES(5)))/CONSTS(228)
ALGBRC(30) = 1.00000/(1.00000+EXP((STATES(31) - CONSTS(258))/CONSTS(257)))
RATES(6) = ( (1.00000 - STATES(6))*ALGBRC(30) - - ( CONSTS(255)*STATES(6)))/CONSTS(256)
ALGBRC(35) = 1.00000/(1.00000+EXP((STATES(35) - CONSTS(286))/CONSTS(285)))
RATES(7) = ( (1.00000 - STATES(7))*ALGBRC(35) - - ( CONSTS(283)*STATES(7)))/CONSTS(284)
ALGBRC(40) = 1.00000/(1.00000+EXP((STATES(39) - CONSTS(314))/CONSTS(313)))
RATES(8) = ( (1.00000 - STATES(8))*ALGBRC(40) - - ( CONSTS(311)*STATES(8)))/CONSTS(312)
ALGBRC(45) = 1.00000/(1.00000+EXP((STATES(43) - CONSTS(342))/CONSTS(341)))
RATES(9) = ( (1.00000 - STATES(9))*ALGBRC(45) - - ( CONSTS(339)*STATES(9)))/CONSTS(340)
ALGBRC(50) = 1.00000/(1.00000+EXP((STATES(47) - CONSTS(370))/CONSTS(369)))
RATES(10) = ( (1.00000 - STATES(10))*ALGBRC(50) - - ( CONSTS(367)*STATES(10)))/CONSTS(368)
ALGBRC(2) = 1.00000/(1.00000+EXP((STATES(11) - CONSTS(98))/CONSTS(99)))
ALGBRC(52) = CONSTS(97)/cosh((STATES(11) - CONSTS(98))/( 2.00000*CONSTS(99)))
RATES(12) = (ALGBRC(2) - STATES(12))/ALGBRC(52)
ALGBRC(3) = 1.00000/(1.00000+EXP((STATES(11) - CONSTS(103))/CONSTS(104)))
ALGBRC(53) = CONSTS(102)/cosh((STATES(11) - CONSTS(103))/( 2.00000*CONSTS(104)))
RATES(13) = (ALGBRC(3) - STATES(13))/ALGBRC(53)
ALGBRC(4) = 1.00000/(1.00000+EXP((STATES(11) - CONSTS(109))/CONSTS(110)))
ALGBRC(54) = CONSTS(108)/cosh((STATES(11) - CONSTS(109))/( 2.00000*CONSTS(110)))
RATES(14) = (ALGBRC(4) - STATES(14))/ALGBRC(54)
ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(15) - CONSTS(126))/CONSTS(127)))
ALGBRC(56) = CONSTS(125)/cosh((STATES(15) - CONSTS(126))/( 2.00000*CONSTS(127)))
RATES(16) = (ALGBRC(7) - STATES(16))/ALGBRC(56)
ALGBRC(8) = 1.00000/(1.00000+EXP((STATES(15) - CONSTS(131))/CONSTS(132)))
ALGBRC(57) = CONSTS(130)/cosh((STATES(15) - CONSTS(131))/( 2.00000*CONSTS(132)))
RATES(17) = (ALGBRC(8) - STATES(17))/ALGBRC(57)
ALGBRC(9) = 1.00000/(1.00000+EXP((STATES(15) - CONSTS(137))/CONSTS(138)))
ALGBRC(58) = CONSTS(136)/cosh((STATES(15) - CONSTS(137))/( 2.00000*CONSTS(138)))
RATES(18) = (ALGBRC(9) - STATES(18))/ALGBRC(58)
ALGBRC(12) = 1.00000/(1.00000+EXP((STATES(19) - CONSTS(154))/CONSTS(155)))
ALGBRC(60) = CONSTS(153)/cosh((STATES(19) - CONSTS(154))/( 2.00000*CONSTS(155)))
RATES(20) = (ALGBRC(12) - STATES(20))/ALGBRC(60)
ALGBRC(13) = 1.00000/(1.00000+EXP((STATES(19) - CONSTS(159))/CONSTS(160)))
ALGBRC(61) = CONSTS(158)/cosh((STATES(19) - CONSTS(159))/( 2.00000*CONSTS(160)))
RATES(21) = (ALGBRC(13) - STATES(21))/ALGBRC(61)
ALGBRC(14) = 1.00000/(1.00000+EXP((STATES(19) - CONSTS(165))/CONSTS(166)))
ALGBRC(62) = CONSTS(164)/cosh((STATES(19) - CONSTS(165))/( 2.00000*CONSTS(166)))
RATES(22) = (ALGBRC(14) - STATES(22))/ALGBRC(62)
ALGBRC(17) = 1.00000/(1.00000+EXP((STATES(23) - CONSTS(182))/CONSTS(183)))
ALGBRC(64) = CONSTS(181)/cosh((STATES(23) - CONSTS(182))/( 2.00000*CONSTS(183)))
RATES(24) = (ALGBRC(17) - STATES(24))/ALGBRC(64)
ALGBRC(18) = 1.00000/(1.00000+EXP((STATES(23) - CONSTS(187))/CONSTS(188)))
ALGBRC(65) = CONSTS(186)/cosh((STATES(23) - CONSTS(187))/( 2.00000*CONSTS(188)))
RATES(25) = (ALGBRC(18) - STATES(25))/ALGBRC(65)
ALGBRC(19) = 1.00000/(1.00000+EXP((STATES(23) - CONSTS(193))/CONSTS(194)))
ALGBRC(66) = CONSTS(192)/cosh((STATES(23) - CONSTS(193))/( 2.00000*CONSTS(194)))
RATES(26) = (ALGBRC(19) - STATES(26))/ALGBRC(66)
ALGBRC(22) = 1.00000/(1.00000+EXP((STATES(27) - CONSTS(210))/CONSTS(211)))
ALGBRC(68) = CONSTS(209)/cosh((STATES(27) - CONSTS(210))/( 2.00000*CONSTS(211)))
RATES(28) = (ALGBRC(22) - STATES(28))/ALGBRC(68)
ALGBRC(23) = 1.00000/(1.00000+EXP((STATES(27) - CONSTS(215))/CONSTS(216)))
ALGBRC(69) = CONSTS(214)/cosh((STATES(27) - CONSTS(215))/( 2.00000*CONSTS(216)))
RATES(29) = (ALGBRC(23) - STATES(29))/ALGBRC(69)
ALGBRC(24) = 1.00000/(1.00000+EXP((STATES(27) - CONSTS(221))/CONSTS(222)))
ALGBRC(70) = CONSTS(220)/cosh((STATES(27) - CONSTS(221))/( 2.00000*CONSTS(222)))
RATES(30) = (ALGBRC(24) - STATES(30))/ALGBRC(70)
ALGBRC(27) = 1.00000/(1.00000+EXP((STATES(31) - CONSTS(238))/CONSTS(239)))
ALGBRC(72) = CONSTS(237)/cosh((STATES(31) - CONSTS(238))/( 2.00000*CONSTS(239)))
RATES(32) = (ALGBRC(27) - STATES(32))/ALGBRC(72)
ALGBRC(28) = 1.00000/(1.00000+EXP((STATES(31) - CONSTS(243))/CONSTS(244)))
ALGBRC(73) = CONSTS(242)/cosh((STATES(31) - CONSTS(243))/( 2.00000*CONSTS(244)))
RATES(33) = (ALGBRC(28) - STATES(33))/ALGBRC(73)
ALGBRC(29) = 1.00000/(1.00000+EXP((STATES(31) - CONSTS(249))/CONSTS(250)))
ALGBRC(74) = CONSTS(248)/cosh((STATES(31) - CONSTS(249))/( 2.00000*CONSTS(250)))
RATES(34) = (ALGBRC(29) - STATES(34))/ALGBRC(74)
ALGBRC(32) = 1.00000/(1.00000+EXP((STATES(35) - CONSTS(266))/CONSTS(267)))
ALGBRC(76) = CONSTS(265)/cosh((STATES(35) - CONSTS(266))/( 2.00000*CONSTS(267)))
RATES(36) = (ALGBRC(32) - STATES(36))/ALGBRC(76)
ALGBRC(33) = 1.00000/(1.00000+EXP((STATES(35) - CONSTS(271))/CONSTS(272)))
ALGBRC(77) = CONSTS(270)/cosh((STATES(35) - CONSTS(271))/( 2.00000*CONSTS(272)))
RATES(37) = (ALGBRC(33) - STATES(37))/ALGBRC(77)
ALGBRC(34) = 1.00000/(1.00000+EXP((STATES(35) - CONSTS(277))/CONSTS(278)))
ALGBRC(78) = CONSTS(276)/cosh((STATES(35) - CONSTS(277))/( 2.00000*CONSTS(278)))
RATES(38) = (ALGBRC(34) - STATES(38))/ALGBRC(78)
ALGBRC(37) = 1.00000/(1.00000+EXP((STATES(39) - CONSTS(294))/CONSTS(295)))
ALGBRC(80) = CONSTS(293)/cosh((STATES(39) - CONSTS(294))/( 2.00000*CONSTS(295)))
RATES(40) = (ALGBRC(37) - STATES(40))/ALGBRC(80)
ALGBRC(38) = 1.00000/(1.00000+EXP((STATES(39) - CONSTS(299))/CONSTS(300)))
ALGBRC(81) = CONSTS(298)/cosh((STATES(39) - CONSTS(299))/( 2.00000*CONSTS(300)))
RATES(41) = (ALGBRC(38) - STATES(41))/ALGBRC(81)
ALGBRC(39) = 1.00000/(1.00000+EXP((STATES(39) - CONSTS(305))/CONSTS(306)))
ALGBRC(82) = CONSTS(304)/cosh((STATES(39) - CONSTS(305))/( 2.00000*CONSTS(306)))
RATES(42) = (ALGBRC(39) - STATES(42))/ALGBRC(82)
ALGBRC(42) = 1.00000/(1.00000+EXP((STATES(43) - CONSTS(322))/CONSTS(323)))
ALGBRC(84) = CONSTS(321)/cosh((STATES(43) - CONSTS(322))/( 2.00000*CONSTS(323)))
RATES(44) = (ALGBRC(42) - STATES(44))/ALGBRC(84)
ALGBRC(43) = 1.00000/(1.00000+EXP((STATES(43) - CONSTS(327))/CONSTS(328)))
ALGBRC(85) = CONSTS(326)/cosh((STATES(43) - CONSTS(327))/( 2.00000*CONSTS(328)))
RATES(45) = (ALGBRC(43) - STATES(45))/ALGBRC(85)
ALGBRC(44) = 1.00000/(1.00000+EXP((STATES(43) - CONSTS(333))/CONSTS(334)))
ALGBRC(86) = CONSTS(332)/cosh((STATES(43) - CONSTS(333))/( 2.00000*CONSTS(334)))
RATES(46) = (ALGBRC(44) - STATES(46))/ALGBRC(86)
ALGBRC(47) = 1.00000/(1.00000+EXP((STATES(47) - CONSTS(350))/CONSTS(351)))
ALGBRC(88) = CONSTS(349)/cosh((STATES(47) - CONSTS(350))/( 2.00000*CONSTS(351)))
RATES(48) = (ALGBRC(47) - STATES(48))/ALGBRC(88)
ALGBRC(48) = 1.00000/(1.00000+EXP((STATES(47) - CONSTS(355))/CONSTS(356)))
ALGBRC(89) = CONSTS(354)/cosh((STATES(47) - CONSTS(355))/( 2.00000*CONSTS(356)))
RATES(49) = (ALGBRC(48) - STATES(49))/ALGBRC(89)
ALGBRC(49) = 1.00000/(1.00000+EXP((STATES(47) - CONSTS(361))/CONSTS(362)))
ALGBRC(90) = CONSTS(360)/cosh((STATES(47) - CONSTS(361))/( 2.00000*CONSTS(362)))
RATES(50) = (ALGBRC(49) - STATES(50))/ALGBRC(90)
ALGBRC(111) = 1.00000/(1.00000+EXP((STATES(11) - CONSTS(106))/CONSTS(107)))
ALGBRC(121) = CONSTS(105)*ALGBRC(111)*STATES(14)*(STATES(11) - CONSTS(93))
ALGBRC(51) = 1.00000/(1.00000+EXP((STATES(11) - CONSTS(95))/CONSTS(96)))
ALGBRC(91) = CONSTS(94)*ALGBRC(51) ** 3.00000*(1.00000 - STATES(12))*(STATES(11) - CONSTS(93))
ALGBRC(101) = CONSTS(100)*STATES(13) ** 4.00000*(STATES(11) - CONSTS(101))
ALGBRC(131) = CONSTS(111)*(STATES(11) - CONSTS(112))
ALGBRC(141) = CONSTS(114)*(STATES(11) - CONSTS(113))
ALGBRC(1) = CONSTS(1)*STATES(1)+ CONSTS(2)*STATES(1)+ CONSTS(3)*STATES(1)+ CONSTS(4)*STATES(1)+ CONSTS(5)*STATES(1)+ CONSTS(6)*STATES(1)+ CONSTS(7)*STATES(1)+ CONSTS(8)*STATES(1)+ CONSTS(9)*STATES(1)
ALGBRC(151) = ALGBRC(1)*(STATES(11) - CONSTS(113))
RATES(11) = (- (ALGBRC(121)+ALGBRC(91)+ALGBRC(101)+ALGBRC(131)+ALGBRC(141)+ALGBRC(151))+CONSTS(92))/CONSTS(91)
ALGBRC(112) = 1.00000/(1.00000+EXP((STATES(15) - CONSTS(134))/CONSTS(135)))
ALGBRC(122) = CONSTS(133)*ALGBRC(112)*STATES(18)*(STATES(15) - CONSTS(121))
ALGBRC(55) = 1.00000/(1.00000+EXP((STATES(15) - CONSTS(123))/CONSTS(124)))
ALGBRC(92) = CONSTS(122)*ALGBRC(55) ** 3.00000*(1.00000 - STATES(16))*(STATES(15) - CONSTS(121))
ALGBRC(102) = CONSTS(128)*STATES(17) ** 4.00000*(STATES(15) - CONSTS(129))
ALGBRC(132) = CONSTS(139)*(STATES(15) - CONSTS(140))
ALGBRC(142) = CONSTS(142)*(STATES(15) - CONSTS(141))
ALGBRC(6) = CONSTS(10)*STATES(2)+ CONSTS(11)*STATES(2)+ CONSTS(12)*STATES(2)+ CONSTS(13)*STATES(2)+ CONSTS(14)*STATES(2)+ CONSTS(15)*STATES(2)+ CONSTS(16)*STATES(2)+ CONSTS(17)*STATES(2)+ CONSTS(18)*STATES(2)
ALGBRC(152) = ALGBRC(6)*(STATES(15) - CONSTS(141))
RATES(15) = (- (ALGBRC(122)+ALGBRC(92)+ALGBRC(102)+ALGBRC(132)+ALGBRC(142)+ALGBRC(152))+CONSTS(120))/CONSTS(119)
ALGBRC(113) = 1.00000/(1.00000+EXP((STATES(19) - CONSTS(162))/CONSTS(163)))
ALGBRC(123) = CONSTS(161)*ALGBRC(113)*STATES(22)*(STATES(19) - CONSTS(149))
ALGBRC(59) = 1.00000/(1.00000+EXP((STATES(19) - CONSTS(151))/CONSTS(152)))
ALGBRC(93) = CONSTS(150)*ALGBRC(59) ** 3.00000*(1.00000 - STATES(20))*(STATES(19) - CONSTS(149))
ALGBRC(103) = CONSTS(156)*STATES(21) ** 4.00000*(STATES(19) - CONSTS(157))
ALGBRC(133) = CONSTS(167)*(STATES(19) - CONSTS(168))
ALGBRC(143) = CONSTS(170)*(STATES(19) - CONSTS(169))
ALGBRC(11) = CONSTS(19)*STATES(3)+ CONSTS(20)*STATES(3)+ CONSTS(21)*STATES(3)+ CONSTS(22)*STATES(3)+ CONSTS(23)*STATES(3)+ CONSTS(24)*STATES(3)+ CONSTS(25)*STATES(3)+ CONSTS(26)*STATES(3)+ CONSTS(27)*STATES(3)
ALGBRC(153) = ALGBRC(11)*(STATES(19) - CONSTS(169))
RATES(19) = (- (ALGBRC(123)+ALGBRC(93)+ALGBRC(103)+ALGBRC(133)+ALGBRC(143)+ALGBRC(153))+CONSTS(148))/CONSTS(147)
ALGBRC(114) = 1.00000/(1.00000+EXP((STATES(23) - CONSTS(190))/CONSTS(191)))
ALGBRC(124) = CONSTS(189)*ALGBRC(114)*STATES(26)*(STATES(23) - CONSTS(177))
ALGBRC(63) = 1.00000/(1.00000+EXP((STATES(23) - CONSTS(179))/CONSTS(180)))
ALGBRC(94) = CONSTS(178)*ALGBRC(63) ** 3.00000*(1.00000 - STATES(24))*(STATES(23) - CONSTS(177))
ALGBRC(104) = CONSTS(184)*STATES(25) ** 4.00000*(STATES(23) - CONSTS(185))
ALGBRC(134) = CONSTS(195)*(STATES(23) - CONSTS(196))
ALGBRC(144) = CONSTS(198)*(STATES(23) - CONSTS(197))
ALGBRC(16) = CONSTS(28)*STATES(4)+ CONSTS(29)*STATES(4)+ CONSTS(30)*STATES(4)+ CONSTS(31)*STATES(4)+ CONSTS(32)*STATES(4)+ CONSTS(33)*STATES(4)+ CONSTS(34)*STATES(4)+ CONSTS(35)*STATES(4)+ CONSTS(36)*STATES(4)
ALGBRC(154) = ALGBRC(16)*(STATES(23) - CONSTS(197))
RATES(23) = (- (ALGBRC(124)+ALGBRC(94)+ALGBRC(104)+ALGBRC(134)+ALGBRC(144)+ALGBRC(154))+CONSTS(176))/CONSTS(175)
ALGBRC(115) = 1.00000/(1.00000+EXP((STATES(27) - CONSTS(218))/CONSTS(219)))
ALGBRC(125) = CONSTS(217)*ALGBRC(115)*STATES(30)*(STATES(27) - CONSTS(205))
ALGBRC(67) = 1.00000/(1.00000+EXP((STATES(27) - CONSTS(207))/CONSTS(208)))
ALGBRC(95) = CONSTS(206)*ALGBRC(67) ** 3.00000*(1.00000 - STATES(28))*(STATES(27) - CONSTS(205))
ALGBRC(105) = CONSTS(212)*STATES(29) ** 4.00000*(STATES(27) - CONSTS(213))
ALGBRC(135) = CONSTS(223)*(STATES(27) - CONSTS(224))
ALGBRC(145) = CONSTS(226)*(STATES(27) - CONSTS(225))
ALGBRC(21) = CONSTS(37)*STATES(5)+ CONSTS(38)*STATES(5)+ CONSTS(39)*STATES(5)+ CONSTS(40)*STATES(5)+ CONSTS(41)*STATES(5)+ CONSTS(42)*STATES(5)+ CONSTS(43)*STATES(5)+ CONSTS(44)*STATES(5)+ CONSTS(45)*STATES(5)
ALGBRC(155) = ALGBRC(21)*(STATES(27) - CONSTS(225))
RATES(27) = (- (ALGBRC(125)+ALGBRC(95)+ALGBRC(105)+ALGBRC(135)+ALGBRC(145)+ALGBRC(155))+CONSTS(204))/CONSTS(203)
ALGBRC(116) = 1.00000/(1.00000+EXP((STATES(31) - CONSTS(246))/CONSTS(247)))
ALGBRC(126) = CONSTS(245)*ALGBRC(116)*STATES(34)*(STATES(31) - CONSTS(233))
ALGBRC(71) = 1.00000/(1.00000+EXP((STATES(31) - CONSTS(235))/CONSTS(236)))
ALGBRC(96) = CONSTS(234)*ALGBRC(71) ** 3.00000*(1.00000 - STATES(32))*(STATES(31) - CONSTS(233))
ALGBRC(106) = CONSTS(240)*STATES(33) ** 4.00000*(STATES(31) - CONSTS(241))
ALGBRC(136) = CONSTS(251)*(STATES(31) - CONSTS(252))
ALGBRC(146) = CONSTS(254)*(STATES(31) - CONSTS(253))
ALGBRC(26) = CONSTS(46)*STATES(6)+ CONSTS(47)*STATES(6)+ CONSTS(48)*STATES(6)+ CONSTS(49)*STATES(6)+ CONSTS(50)*STATES(6)+ CONSTS(51)*STATES(6)+ CONSTS(52)*STATES(6)+ CONSTS(53)*STATES(6)+ CONSTS(54)*STATES(6)
ALGBRC(156) = ALGBRC(26)*(STATES(31) - CONSTS(253))
RATES(31) = (- (ALGBRC(126)+ALGBRC(96)+ALGBRC(106)+ALGBRC(136)+ALGBRC(146)+ALGBRC(156))+CONSTS(232))/CONSTS(231)
ALGBRC(117) = 1.00000/(1.00000+EXP((STATES(35) - CONSTS(274))/CONSTS(275)))
ALGBRC(127) = CONSTS(273)*ALGBRC(117)*STATES(38)*(STATES(35) - CONSTS(261))
ALGBRC(75) = 1.00000/(1.00000+EXP((STATES(35) - CONSTS(263))/CONSTS(264)))
ALGBRC(97) = CONSTS(262)*ALGBRC(75) ** 3.00000*(1.00000 - STATES(36))*(STATES(35) - CONSTS(261))
ALGBRC(107) = CONSTS(268)*STATES(37) ** 4.00000*(STATES(35) - CONSTS(269))
ALGBRC(137) = CONSTS(279)*(STATES(35) - CONSTS(280))
ALGBRC(147) = CONSTS(282)*(STATES(35) - CONSTS(281))
ALGBRC(31) = CONSTS(55)*STATES(7)+ CONSTS(56)*STATES(7)+ CONSTS(57)*STATES(7)+ CONSTS(58)*STATES(7)+ CONSTS(59)*STATES(7)+ CONSTS(60)*STATES(7)+ CONSTS(61)*STATES(7)+ CONSTS(62)*STATES(7)+ CONSTS(63)*STATES(7)
ALGBRC(157) = ALGBRC(31)*(STATES(35) - CONSTS(281))
RATES(35) = (- (ALGBRC(127)+ALGBRC(97)+ALGBRC(107)+ALGBRC(137)+ALGBRC(147)+ALGBRC(157))+CONSTS(260))/CONSTS(259)
ALGBRC(118) = 1.00000/(1.00000+EXP((STATES(39) - CONSTS(302))/CONSTS(303)))
ALGBRC(128) = CONSTS(301)*ALGBRC(118)*STATES(42)*(STATES(39) - CONSTS(289))
ALGBRC(79) = 1.00000/(1.00000+EXP((STATES(39) - CONSTS(291))/CONSTS(292)))
ALGBRC(98) = CONSTS(290)*ALGBRC(79) ** 3.00000*(1.00000 - STATES(40))*(STATES(39) - CONSTS(289))
ALGBRC(108) = CONSTS(296)*STATES(41) ** 4.00000*(STATES(39) - CONSTS(297))
ALGBRC(138) = CONSTS(307)*(STATES(39) - CONSTS(308))
ALGBRC(148) = CONSTS(310)*(STATES(39) - CONSTS(309))
ALGBRC(36) = CONSTS(64)*STATES(8)+ CONSTS(65)*STATES(8)+ CONSTS(66)*STATES(8)+ CONSTS(67)*STATES(8)+ CONSTS(68)*STATES(8)+ CONSTS(69)*STATES(8)+ CONSTS(70)*STATES(8)+ CONSTS(71)*STATES(8)+ CONSTS(72)*STATES(8)
ALGBRC(158) = ALGBRC(36)*(STATES(39) - CONSTS(309))
RATES(39) = (- (ALGBRC(128)+ALGBRC(98)+ALGBRC(108)+ALGBRC(138)+ALGBRC(148)+ALGBRC(158))+CONSTS(288))/CONSTS(287)
ALGBRC(119) = 1.00000/(1.00000+EXP((STATES(43) - CONSTS(330))/CONSTS(331)))
ALGBRC(129) = CONSTS(329)*ALGBRC(119)*STATES(46)*(STATES(43) - CONSTS(317))
ALGBRC(83) = 1.00000/(1.00000+EXP((STATES(43) - CONSTS(319))/CONSTS(320)))
ALGBRC(99) = CONSTS(318)*ALGBRC(83) ** 3.00000*(1.00000 - STATES(44))*(STATES(43) - CONSTS(317))
ALGBRC(109) = CONSTS(324)*STATES(45) ** 4.00000*(STATES(43) - CONSTS(325))
ALGBRC(139) = CONSTS(335)*(STATES(43) - CONSTS(336))
ALGBRC(149) = CONSTS(338)*(STATES(43) - CONSTS(337))
ALGBRC(41) = CONSTS(73)*STATES(9)+ CONSTS(74)*STATES(9)+ CONSTS(75)*STATES(9)+ CONSTS(76)*STATES(9)+ CONSTS(77)*STATES(9)+ CONSTS(78)*STATES(9)+ CONSTS(79)*STATES(9)+ CONSTS(80)*STATES(9)+ CONSTS(81)*STATES(9)
ALGBRC(159) = ALGBRC(41)*(STATES(43) - CONSTS(337))
RATES(43) = (- (ALGBRC(129)+ALGBRC(99)+ALGBRC(109)+ALGBRC(139)+ALGBRC(149)+ALGBRC(159))+CONSTS(316))/CONSTS(315)
ALGBRC(120) = 1.00000/(1.00000+EXP((STATES(47) - CONSTS(358))/CONSTS(359)))
ALGBRC(130) = CONSTS(357)*ALGBRC(120)*STATES(50)*(STATES(47) - CONSTS(345))
ALGBRC(87) = 1.00000/(1.00000+EXP((STATES(47) - CONSTS(347))/CONSTS(348)))
ALGBRC(100) = CONSTS(346)*ALGBRC(87) ** 3.00000*(1.00000 - STATES(48))*(STATES(47) - CONSTS(345))
ALGBRC(110) = CONSTS(352)*STATES(49) ** 4.00000*(STATES(47) - CONSTS(353))
ALGBRC(140) = CONSTS(363)*(STATES(47) - CONSTS(364))
ALGBRC(150) = CONSTS(366)*(STATES(47) - CONSTS(365))
ALGBRC(46) = CONSTS(82)*STATES(10)+ CONSTS(83)*STATES(10)+ CONSTS(84)*STATES(10)+ CONSTS(85)*STATES(10)+ CONSTS(86)*STATES(10)+ CONSTS(87)*STATES(10)+ CONSTS(88)*STATES(10)+ CONSTS(89)*STATES(10)+ CONSTS(90)*STATES(10)
ALGBRC(160) = ALGBRC(46)*(STATES(47) - CONSTS(365))
RATES(47) = (- (ALGBRC(130)+ALGBRC(100)+ALGBRC(110)+ALGBRC(140)+ALGBRC(150)+ALGBRC(160))+CONSTS(344))/CONSTS(343)
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(5) = 1.00000/(1.00000+EXP((STATES(11) - CONSTS(118))/CONSTS(117)))
ALGBRC(10) = 1.00000/(1.00000+EXP((STATES(15) - CONSTS(146))/CONSTS(145)))
ALGBRC(15) = 1.00000/(1.00000+EXP((STATES(19) - CONSTS(174))/CONSTS(173)))
ALGBRC(20) = 1.00000/(1.00000+EXP((STATES(23) - CONSTS(202))/CONSTS(201)))
ALGBRC(25) = 1.00000/(1.00000+EXP((STATES(27) - CONSTS(230))/CONSTS(229)))
ALGBRC(30) = 1.00000/(1.00000+EXP((STATES(31) - CONSTS(258))/CONSTS(257)))
ALGBRC(35) = 1.00000/(1.00000+EXP((STATES(35) - CONSTS(286))/CONSTS(285)))
ALGBRC(40) = 1.00000/(1.00000+EXP((STATES(39) - CONSTS(314))/CONSTS(313)))
ALGBRC(45) = 1.00000/(1.00000+EXP((STATES(43) - CONSTS(342))/CONSTS(341)))
ALGBRC(50) = 1.00000/(1.00000+EXP((STATES(47) - CONSTS(370))/CONSTS(369)))
ALGBRC(2) = 1.00000/(1.00000+EXP((STATES(11) - CONSTS(98))/CONSTS(99)))
ALGBRC(52) = CONSTS(97)/cosh((STATES(11) - CONSTS(98))/( 2.00000*CONSTS(99)))
ALGBRC(3) = 1.00000/(1.00000+EXP((STATES(11) - CONSTS(103))/CONSTS(104)))
ALGBRC(53) = CONSTS(102)/cosh((STATES(11) - CONSTS(103))/( 2.00000*CONSTS(104)))
ALGBRC(4) = 1.00000/(1.00000+EXP((STATES(11) - CONSTS(109))/CONSTS(110)))
ALGBRC(54) = CONSTS(108)/cosh((STATES(11) - CONSTS(109))/( 2.00000*CONSTS(110)))
ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(15) - CONSTS(126))/CONSTS(127)))
ALGBRC(56) = CONSTS(125)/cosh((STATES(15) - CONSTS(126))/( 2.00000*CONSTS(127)))
ALGBRC(8) = 1.00000/(1.00000+EXP((STATES(15) - CONSTS(131))/CONSTS(132)))
ALGBRC(57) = CONSTS(130)/cosh((STATES(15) - CONSTS(131))/( 2.00000*CONSTS(132)))
ALGBRC(9) = 1.00000/(1.00000+EXP((STATES(15) - CONSTS(137))/CONSTS(138)))
ALGBRC(58) = CONSTS(136)/cosh((STATES(15) - CONSTS(137))/( 2.00000*CONSTS(138)))
ALGBRC(12) = 1.00000/(1.00000+EXP((STATES(19) - CONSTS(154))/CONSTS(155)))
ALGBRC(60) = CONSTS(153)/cosh((STATES(19) - CONSTS(154))/( 2.00000*CONSTS(155)))
ALGBRC(13) = 1.00000/(1.00000+EXP((STATES(19) - CONSTS(159))/CONSTS(160)))
ALGBRC(61) = CONSTS(158)/cosh((STATES(19) - CONSTS(159))/( 2.00000*CONSTS(160)))
ALGBRC(14) = 1.00000/(1.00000+EXP((STATES(19) - CONSTS(165))/CONSTS(166)))
ALGBRC(62) = CONSTS(164)/cosh((STATES(19) - CONSTS(165))/( 2.00000*CONSTS(166)))
ALGBRC(17) = 1.00000/(1.00000+EXP((STATES(23) - CONSTS(182))/CONSTS(183)))
ALGBRC(64) = CONSTS(181)/cosh((STATES(23) - CONSTS(182))/( 2.00000*CONSTS(183)))
ALGBRC(18) = 1.00000/(1.00000+EXP((STATES(23) - CONSTS(187))/CONSTS(188)))
ALGBRC(65) = CONSTS(186)/cosh((STATES(23) - CONSTS(187))/( 2.00000*CONSTS(188)))
ALGBRC(19) = 1.00000/(1.00000+EXP((STATES(23) - CONSTS(193))/CONSTS(194)))
ALGBRC(66) = CONSTS(192)/cosh((STATES(23) - CONSTS(193))/( 2.00000*CONSTS(194)))
ALGBRC(22) = 1.00000/(1.00000+EXP((STATES(27) - CONSTS(210))/CONSTS(211)))
ALGBRC(68) = CONSTS(209)/cosh((STATES(27) - CONSTS(210))/( 2.00000*CONSTS(211)))
ALGBRC(23) = 1.00000/(1.00000+EXP((STATES(27) - CONSTS(215))/CONSTS(216)))
ALGBRC(69) = CONSTS(214)/cosh((STATES(27) - CONSTS(215))/( 2.00000*CONSTS(216)))
ALGBRC(24) = 1.00000/(1.00000+EXP((STATES(27) - CONSTS(221))/CONSTS(222)))
ALGBRC(70) = CONSTS(220)/cosh((STATES(27) - CONSTS(221))/( 2.00000*CONSTS(222)))
ALGBRC(27) = 1.00000/(1.00000+EXP((STATES(31) - CONSTS(238))/CONSTS(239)))
ALGBRC(72) = CONSTS(237)/cosh((STATES(31) - CONSTS(238))/( 2.00000*CONSTS(239)))
ALGBRC(28) = 1.00000/(1.00000+EXP((STATES(31) - CONSTS(243))/CONSTS(244)))
ALGBRC(73) = CONSTS(242)/cosh((STATES(31) - CONSTS(243))/( 2.00000*CONSTS(244)))
ALGBRC(29) = 1.00000/(1.00000+EXP((STATES(31) - CONSTS(249))/CONSTS(250)))
ALGBRC(74) = CONSTS(248)/cosh((STATES(31) - CONSTS(249))/( 2.00000*CONSTS(250)))
ALGBRC(32) = 1.00000/(1.00000+EXP((STATES(35) - CONSTS(266))/CONSTS(267)))
ALGBRC(76) = CONSTS(265)/cosh((STATES(35) - CONSTS(266))/( 2.00000*CONSTS(267)))
ALGBRC(33) = 1.00000/(1.00000+EXP((STATES(35) - CONSTS(271))/CONSTS(272)))
ALGBRC(77) = CONSTS(270)/cosh((STATES(35) - CONSTS(271))/( 2.00000*CONSTS(272)))
ALGBRC(34) = 1.00000/(1.00000+EXP((STATES(35) - CONSTS(277))/CONSTS(278)))
ALGBRC(78) = CONSTS(276)/cosh((STATES(35) - CONSTS(277))/( 2.00000*CONSTS(278)))
ALGBRC(37) = 1.00000/(1.00000+EXP((STATES(39) - CONSTS(294))/CONSTS(295)))
ALGBRC(80) = CONSTS(293)/cosh((STATES(39) - CONSTS(294))/( 2.00000*CONSTS(295)))
ALGBRC(38) = 1.00000/(1.00000+EXP((STATES(39) - CONSTS(299))/CONSTS(300)))
ALGBRC(81) = CONSTS(298)/cosh((STATES(39) - CONSTS(299))/( 2.00000*CONSTS(300)))
ALGBRC(39) = 1.00000/(1.00000+EXP((STATES(39) - CONSTS(305))/CONSTS(306)))
ALGBRC(82) = CONSTS(304)/cosh((STATES(39) - CONSTS(305))/( 2.00000*CONSTS(306)))
ALGBRC(42) = 1.00000/(1.00000+EXP((STATES(43) - CONSTS(322))/CONSTS(323)))
ALGBRC(84) = CONSTS(321)/cosh((STATES(43) - CONSTS(322))/( 2.00000*CONSTS(323)))
ALGBRC(43) = 1.00000/(1.00000+EXP((STATES(43) - CONSTS(327))/CONSTS(328)))
ALGBRC(85) = CONSTS(326)/cosh((STATES(43) - CONSTS(327))/( 2.00000*CONSTS(328)))
ALGBRC(44) = 1.00000/(1.00000+EXP((STATES(43) - CONSTS(333))/CONSTS(334)))
ALGBRC(86) = CONSTS(332)/cosh((STATES(43) - CONSTS(333))/( 2.00000*CONSTS(334)))
ALGBRC(47) = 1.00000/(1.00000+EXP((STATES(47) - CONSTS(350))/CONSTS(351)))
ALGBRC(88) = CONSTS(349)/cosh((STATES(47) - CONSTS(350))/( 2.00000*CONSTS(351)))
ALGBRC(48) = 1.00000/(1.00000+EXP((STATES(47) - CONSTS(355))/CONSTS(356)))
ALGBRC(89) = CONSTS(354)/cosh((STATES(47) - CONSTS(355))/( 2.00000*CONSTS(356)))
ALGBRC(49) = 1.00000/(1.00000+EXP((STATES(47) - CONSTS(361))/CONSTS(362)))
ALGBRC(90) = CONSTS(360)/cosh((STATES(47) - CONSTS(361))/( 2.00000*CONSTS(362)))
ALGBRC(111) = 1.00000/(1.00000+EXP((STATES(11) - CONSTS(106))/CONSTS(107)))
ALGBRC(121) = CONSTS(105)*ALGBRC(111)*STATES(14)*(STATES(11) - CONSTS(93))
ALGBRC(51) = 1.00000/(1.00000+EXP((STATES(11) - CONSTS(95))/CONSTS(96)))
ALGBRC(91) = CONSTS(94)*ALGBRC(51) ** 3.00000*(1.00000 - STATES(12))*(STATES(11) - CONSTS(93))
ALGBRC(101) = CONSTS(100)*STATES(13) ** 4.00000*(STATES(11) - CONSTS(101))
ALGBRC(131) = CONSTS(111)*(STATES(11) - CONSTS(112))
ALGBRC(141) = CONSTS(114)*(STATES(11) - CONSTS(113))
ALGBRC(1) = CONSTS(1)*STATES(1)+ CONSTS(2)*STATES(1)+ CONSTS(3)*STATES(1)+ CONSTS(4)*STATES(1)+ CONSTS(5)*STATES(1)+ CONSTS(6)*STATES(1)+ CONSTS(7)*STATES(1)+ CONSTS(8)*STATES(1)+ CONSTS(9)*STATES(1)
ALGBRC(151) = ALGBRC(1)*(STATES(11) - CONSTS(113))
ALGBRC(112) = 1.00000/(1.00000+EXP((STATES(15) - CONSTS(134))/CONSTS(135)))
ALGBRC(122) = CONSTS(133)*ALGBRC(112)*STATES(18)*(STATES(15) - CONSTS(121))
ALGBRC(55) = 1.00000/(1.00000+EXP((STATES(15) - CONSTS(123))/CONSTS(124)))
ALGBRC(92) = CONSTS(122)*ALGBRC(55) ** 3.00000*(1.00000 - STATES(16))*(STATES(15) - CONSTS(121))
ALGBRC(102) = CONSTS(128)*STATES(17) ** 4.00000*(STATES(15) - CONSTS(129))
ALGBRC(132) = CONSTS(139)*(STATES(15) - CONSTS(140))
ALGBRC(142) = CONSTS(142)*(STATES(15) - CONSTS(141))
ALGBRC(6) = CONSTS(10)*STATES(2)+ CONSTS(11)*STATES(2)+ CONSTS(12)*STATES(2)+ CONSTS(13)*STATES(2)+ CONSTS(14)*STATES(2)+ CONSTS(15)*STATES(2)+ CONSTS(16)*STATES(2)+ CONSTS(17)*STATES(2)+ CONSTS(18)*STATES(2)
ALGBRC(152) = ALGBRC(6)*(STATES(15) - CONSTS(141))
ALGBRC(113) = 1.00000/(1.00000+EXP((STATES(19) - CONSTS(162))/CONSTS(163)))
ALGBRC(123) = CONSTS(161)*ALGBRC(113)*STATES(22)*(STATES(19) - CONSTS(149))
ALGBRC(59) = 1.00000/(1.00000+EXP((STATES(19) - CONSTS(151))/CONSTS(152)))
ALGBRC(93) = CONSTS(150)*ALGBRC(59) ** 3.00000*(1.00000 - STATES(20))*(STATES(19) - CONSTS(149))
ALGBRC(103) = CONSTS(156)*STATES(21) ** 4.00000*(STATES(19) - CONSTS(157))
ALGBRC(133) = CONSTS(167)*(STATES(19) - CONSTS(168))
ALGBRC(143) = CONSTS(170)*(STATES(19) - CONSTS(169))
ALGBRC(11) = CONSTS(19)*STATES(3)+ CONSTS(20)*STATES(3)+ CONSTS(21)*STATES(3)+ CONSTS(22)*STATES(3)+ CONSTS(23)*STATES(3)+ CONSTS(24)*STATES(3)+ CONSTS(25)*STATES(3)+ CONSTS(26)*STATES(3)+ CONSTS(27)*STATES(3)
ALGBRC(153) = ALGBRC(11)*(STATES(19) - CONSTS(169))
ALGBRC(114) = 1.00000/(1.00000+EXP((STATES(23) - CONSTS(190))/CONSTS(191)))
ALGBRC(124) = CONSTS(189)*ALGBRC(114)*STATES(26)*(STATES(23) - CONSTS(177))
ALGBRC(63) = 1.00000/(1.00000+EXP((STATES(23) - CONSTS(179))/CONSTS(180)))
ALGBRC(94) = CONSTS(178)*ALGBRC(63) ** 3.00000*(1.00000 - STATES(24))*(STATES(23) - CONSTS(177))
ALGBRC(104) = CONSTS(184)*STATES(25) ** 4.00000*(STATES(23) - CONSTS(185))
ALGBRC(134) = CONSTS(195)*(STATES(23) - CONSTS(196))
ALGBRC(144) = CONSTS(198)*(STATES(23) - CONSTS(197))
ALGBRC(16) = CONSTS(28)*STATES(4)+ CONSTS(29)*STATES(4)+ CONSTS(30)*STATES(4)+ CONSTS(31)*STATES(4)+ CONSTS(32)*STATES(4)+ CONSTS(33)*STATES(4)+ CONSTS(34)*STATES(4)+ CONSTS(35)*STATES(4)+ CONSTS(36)*STATES(4)
ALGBRC(154) = ALGBRC(16)*(STATES(23) - CONSTS(197))
ALGBRC(115) = 1.00000/(1.00000+EXP((STATES(27) - CONSTS(218))/CONSTS(219)))
ALGBRC(125) = CONSTS(217)*ALGBRC(115)*STATES(30)*(STATES(27) - CONSTS(205))
ALGBRC(67) = 1.00000/(1.00000+EXP((STATES(27) - CONSTS(207))/CONSTS(208)))
ALGBRC(95) = CONSTS(206)*ALGBRC(67) ** 3.00000*(1.00000 - STATES(28))*(STATES(27) - CONSTS(205))
ALGBRC(105) = CONSTS(212)*STATES(29) ** 4.00000*(STATES(27) - CONSTS(213))
ALGBRC(135) = CONSTS(223)*(STATES(27) - CONSTS(224))
ALGBRC(145) = CONSTS(226)*(STATES(27) - CONSTS(225))
ALGBRC(21) = CONSTS(37)*STATES(5)+ CONSTS(38)*STATES(5)+ CONSTS(39)*STATES(5)+ CONSTS(40)*STATES(5)+ CONSTS(41)*STATES(5)+ CONSTS(42)*STATES(5)+ CONSTS(43)*STATES(5)+ CONSTS(44)*STATES(5)+ CONSTS(45)*STATES(5)
ALGBRC(155) = ALGBRC(21)*(STATES(27) - CONSTS(225))
ALGBRC(116) = 1.00000/(1.00000+EXP((STATES(31) - CONSTS(246))/CONSTS(247)))
ALGBRC(126) = CONSTS(245)*ALGBRC(116)*STATES(34)*(STATES(31) - CONSTS(233))
ALGBRC(71) = 1.00000/(1.00000+EXP((STATES(31) - CONSTS(235))/CONSTS(236)))
ALGBRC(96) = CONSTS(234)*ALGBRC(71) ** 3.00000*(1.00000 - STATES(32))*(STATES(31) - CONSTS(233))
ALGBRC(106) = CONSTS(240)*STATES(33) ** 4.00000*(STATES(31) - CONSTS(241))
ALGBRC(136) = CONSTS(251)*(STATES(31) - CONSTS(252))
ALGBRC(146) = CONSTS(254)*(STATES(31) - CONSTS(253))
ALGBRC(26) = CONSTS(46)*STATES(6)+ CONSTS(47)*STATES(6)+ CONSTS(48)*STATES(6)+ CONSTS(49)*STATES(6)+ CONSTS(50)*STATES(6)+ CONSTS(51)*STATES(6)+ CONSTS(52)*STATES(6)+ CONSTS(53)*STATES(6)+ CONSTS(54)*STATES(6)
ALGBRC(156) = ALGBRC(26)*(STATES(31) - CONSTS(253))
ALGBRC(117) = 1.00000/(1.00000+EXP((STATES(35) - CONSTS(274))/CONSTS(275)))
ALGBRC(127) = CONSTS(273)*ALGBRC(117)*STATES(38)*(STATES(35) - CONSTS(261))
ALGBRC(75) = 1.00000/(1.00000+EXP((STATES(35) - CONSTS(263))/CONSTS(264)))
ALGBRC(97) = CONSTS(262)*ALGBRC(75) ** 3.00000*(1.00000 - STATES(36))*(STATES(35) - CONSTS(261))
ALGBRC(107) = CONSTS(268)*STATES(37) ** 4.00000*(STATES(35) - CONSTS(269))
ALGBRC(137) = CONSTS(279)*(STATES(35) - CONSTS(280))
ALGBRC(147) = CONSTS(282)*(STATES(35) - CONSTS(281))
ALGBRC(31) = CONSTS(55)*STATES(7)+ CONSTS(56)*STATES(7)+ CONSTS(57)*STATES(7)+ CONSTS(58)*STATES(7)+ CONSTS(59)*STATES(7)+ CONSTS(60)*STATES(7)+ CONSTS(61)*STATES(7)+ CONSTS(62)*STATES(7)+ CONSTS(63)*STATES(7)
ALGBRC(157) = ALGBRC(31)*(STATES(35) - CONSTS(281))
ALGBRC(118) = 1.00000/(1.00000+EXP((STATES(39) - CONSTS(302))/CONSTS(303)))
ALGBRC(128) = CONSTS(301)*ALGBRC(118)*STATES(42)*(STATES(39) - CONSTS(289))
ALGBRC(79) = 1.00000/(1.00000+EXP((STATES(39) - CONSTS(291))/CONSTS(292)))
ALGBRC(98) = CONSTS(290)*ALGBRC(79) ** 3.00000*(1.00000 - STATES(40))*(STATES(39) - CONSTS(289))
ALGBRC(108) = CONSTS(296)*STATES(41) ** 4.00000*(STATES(39) - CONSTS(297))
ALGBRC(138) = CONSTS(307)*(STATES(39) - CONSTS(308))
ALGBRC(148) = CONSTS(310)*(STATES(39) - CONSTS(309))
ALGBRC(36) = CONSTS(64)*STATES(8)+ CONSTS(65)*STATES(8)+ CONSTS(66)*STATES(8)+ CONSTS(67)*STATES(8)+ CONSTS(68)*STATES(8)+ CONSTS(69)*STATES(8)+ CONSTS(70)*STATES(8)+ CONSTS(71)*STATES(8)+ CONSTS(72)*STATES(8)
ALGBRC(158) = ALGBRC(36)*(STATES(39) - CONSTS(309))
ALGBRC(119) = 1.00000/(1.00000+EXP((STATES(43) - CONSTS(330))/CONSTS(331)))
ALGBRC(129) = CONSTS(329)*ALGBRC(119)*STATES(46)*(STATES(43) - CONSTS(317))
ALGBRC(83) = 1.00000/(1.00000+EXP((STATES(43) - CONSTS(319))/CONSTS(320)))
ALGBRC(99) = CONSTS(318)*ALGBRC(83) ** 3.00000*(1.00000 - STATES(44))*(STATES(43) - CONSTS(317))
ALGBRC(109) = CONSTS(324)*STATES(45) ** 4.00000*(STATES(43) - CONSTS(325))
ALGBRC(139) = CONSTS(335)*(STATES(43) - CONSTS(336))
ALGBRC(149) = CONSTS(338)*(STATES(43) - CONSTS(337))
ALGBRC(41) = CONSTS(73)*STATES(9)+ CONSTS(74)*STATES(9)+ CONSTS(75)*STATES(9)+ CONSTS(76)*STATES(9)+ CONSTS(77)*STATES(9)+ CONSTS(78)*STATES(9)+ CONSTS(79)*STATES(9)+ CONSTS(80)*STATES(9)+ CONSTS(81)*STATES(9)
ALGBRC(159) = ALGBRC(41)*(STATES(43) - CONSTS(337))
ALGBRC(120) = 1.00000/(1.00000+EXP((STATES(47) - CONSTS(358))/CONSTS(359)))
ALGBRC(130) = CONSTS(357)*ALGBRC(120)*STATES(50)*(STATES(47) - CONSTS(345))
ALGBRC(87) = 1.00000/(1.00000+EXP((STATES(47) - CONSTS(347))/CONSTS(348)))
ALGBRC(100) = CONSTS(346)*ALGBRC(87) ** 3.00000*(1.00000 - STATES(48))*(STATES(47) - CONSTS(345))
ALGBRC(110) = CONSTS(352)*STATES(49) ** 4.00000*(STATES(47) - CONSTS(353))
ALGBRC(140) = CONSTS(363)*(STATES(47) - CONSTS(364))
ALGBRC(150) = CONSTS(366)*(STATES(47) - CONSTS(365))
ALGBRC(46) = CONSTS(82)*STATES(10)+ CONSTS(83)*STATES(10)+ CONSTS(84)*STATES(10)+ CONSTS(85)*STATES(10)+ CONSTS(86)*STATES(10)+ CONSTS(87)*STATES(10)+ CONSTS(88)*STATES(10)+ CONSTS(89)*STATES(10)+ CONSTS(90)*STATES(10)
ALGBRC(160) = ALGBRC(46)*(STATES(47) - CONSTS(365))
RETURN
END