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