
        format none;

        function xdot = f (x, t)
            # Zmena hodnoty x za jednotku casu
            xdot = zeros (210,1);
            xdot(1) = 0; # fixed, cytosol 
xdot(2) = 0; # fixed, external 
xdot(3) = 0; # fixed, K_GOA 
xdot(4) = 0; # fixed, K_Sink_PEP 
xdot(5) = 0; # fixed, KE10 
xdot(6) = 0; # fixed, KE11 
xdot(7) = 0; # fixed, KE12 
xdot(8) = 0; # fixed, KE4 
xdot(9) = 0; # fixed, KE5 
xdot(10) = 0; # fixed, KE7 
xdot(11) = 0; # fixed, KE8 
xdot(12) = 0; # fixed, KE_SS1 
xdot(13) = 0; # fixed, KE_SS2 
xdot(14) = 0; # fixed, Keq_enol 
xdot(15) = 0; # fixed, Keq_gap_dehyd 
xdot(16) = 0; # fixed, Keq_PGM 
xdot(17) = 0; # fixed, kf_CO2 
xdot(18) = 0; # fixed, kf_NADPH 
xdot(19) = 0; # fixed, Kgap 
xdot(20) = 0; # fixed, Kgap_beta 
xdot(21) = 0; # fixed, KM101 
xdot(22) = 0; # fixed, Km112a 
xdot(23) = 0; # fixed, Km112b 
xdot(24) = 0; # fixed, KM132 
xdot(25) = 0; # fixed, KM51 
xdot(26) = 0; # fixed, KM52 
xdot(27) = 0; # fixed, KM53 
xdot(28) = 0; # fixed, KM61 
xdot(29) = 0; # fixed, KM71 
xdot(30) = 0; # fixed, KM71a 
xdot(31) = 0; # fixed, KM72 
xdot(32) = 0; # fixed, KM73 
xdot(33) = 0; # fixed, KM74 
xdot(34) = 0; # fixed, KM81 
xdot(35) = 0; # fixed, KM82 
xdot(36) = 0; # fixed, KM9 
xdot(37) = 0; # fixed, KM_PFK 
xdot(38) = 0; # fixed, Kmp_enol 
xdot(39) = 0; # fixed, Kmp_PGM_alpha 
xdot(40) = 0; # fixed, Kmp_PGM_beta 
xdot(41) = 0; # fixed, Kmp_PGM_gama 
xdot(42) = 0; # fixed, Kms_enol 
xdot(43) = 0; # fixed, Kms_PGM_alpha 
xdot(44) = 0; # fixed, Kms_PGM_beta 
xdot(45) = 0; # fixed, Kms_PGM_gama 
xdot(46) = 0; # fixed, knadph 
xdot(47) = 0; # fixed, knadph_beta 
xdot(48) = 0; # fixed, Knadpp 
xdot(49) = 0; # fixed, Knadpp_beta 
xdot(50) = 0; # fixed, Kpga 
xdot(51) = 0; # fixed, Kpga_beta 
xdot(52) = 0; # fixed, V1 
xdot(53) = 0; # fixed, V10 
xdot(54) = 0; # fixed, V11 
xdot(55) = 0; # fixed, V12 
xdot(56) = 0; # fixed, V13 
xdot(57) = 0; # fixed, V2 
xdot(58) = 0; # fixed, V3 
xdot(59) = 0; # fixed, V4 
xdot(60) = 0; # fixed, V5 
xdot(61) = 0; # fixed, V6 
xdot(62) = 0; # fixed, V7 
xdot(63) = 0; # fixed, V8 
xdot(64) = 0; # fixed, V9 
xdot(65) = 0; # fixed, V_GLY_syn 
xdot(66) = 0; # fixed, V_OXA1 
xdot(67) = 0; # fixed, V_OXA2 
xdot(68) = 0; # fixed, V_PFK 
xdot(69) = 0; # fixed, V_PP1 
xdot(70) = 0; # fixed, V_PP2a 
xdot(71) = 0; # fixed, V_PP2b 
xdot(72) = 0; # fixed, V_PP3 
xdot(73) = 0; # fixed, V_PP4 
xdot(74) = 0; # fixed, V_PP5 
xdot(75) = 0; # fixed, V_PP6 
xdot(76) = 0; # fixed, V_PP7 
xdot(77) = 0; # fixed, V_Sink_DHAP 
xdot(78) = 0; # fixed, V_Sink_E4P 
xdot(79) = 0; # fixed, V_Sink_GLY 
xdot(80) = 0; # fixed, V_Sink_PEP 
xdot(81) = 0; # fixed, V_Sink_Ri5P 
xdot(82) = 0; # fixed, V_Sink_SER 
xdot(83) = 0; # fixed, V_Sink_THR 
xdot(84) = 0; # fixed, V_SS1 
xdot(85) = 0; # fixed, V_synth_SER 
xdot(86) = 0; # fixed, V_synth_THR 
xdot(87) = 0; # fixed, V_TSA1 
xdot(88) = 0; # fixed, V_TSA2 
xdot(89) = 0; # fixed, V_TSA3 
xdot(90) = 0; # fixed, Vf_enol 
xdot(91) = 0; # fixed, Vf_PGM_alpha 
xdot(92) = 0; # fixed, Vf_PGM_beta 
xdot(93) = 0; # fixed, Vf_PGM_gama 
xdot(94) = 0; # fixed, Vgap_dehyd 
xdot(95) = 0; # fixed, Vgap_dehyd_beta 
xdot(96) = 0; # fixed, kf_ATP_LCgene 
xdot(97) = 0; # fixed, kf_NADPH_LCgene 
xdot(98) = 0; # fixed, kf_CO2_LCgene 
xdot(99) = 0; # fixed, V1_LCgene 
xdot(100) = 0; # fixed, V2_LCgene 
xdot(101) = 0; # fixed, V3_LCgene 
xdot(102) = 0; # fixed, V4_LCgene 
xdot(103) = 0; # fixed, V5_LCgene 
xdot(104) = 0; # fixed, V6_LCgene 
xdot(105) = 0; # fixed, V7_LCgene 
xdot(106) = 0; # fixed, V8_LCgene 
xdot(107) = 0; # fixed, V9_LCgene 
xdot(108) = 0; # fixed, V10_LCgene 
xdot(109) = 0; # fixed, V11_LCgene 
xdot(110) = 0; # fixed, V12_LCgene 
xdot(111) = 0; # fixed, V13_LCgene 
xdot(112) = 0; # fixed, V_PFK_LCgene 
xdot(113) = 0; # fixed, Vgap_dehyd_LCgene 
xdot(114) = 0; # fixed, Vgap_dehyd_beta_LCgene 
xdot(115) = 0; # fixed, Vf_PGM_alpha_LCgene 
xdot(116) = 0; # fixed, Vf_PGM_beta_LCgene 
xdot(117) = 0; # fixed, Vf_PGM_gama_LCgene 
xdot(118) = 0; # fixed, Vf_enol_LCgene 
xdot(119) = 0; # fixed, V_PP1_LCgene 
xdot(120) = 0; # fixed, V_PP2a_LCgene 
xdot(121) = 0; # fixed, V_PP2b_LCgene 
xdot(122) = 0; # fixed, V_PP3_LCgene 
xdot(123) = 0; # fixed, V_PP4_LCgene 
xdot(124) = 0; # fixed, V_PP5_LCgene 
xdot(125) = 0; # fixed, V_PP6_LCgene 
xdot(126) = 0; # fixed, V_GLY_syn_LCgene 
xdot(127) = 0; # fixed, V_SS1_LCgene 
xdot(162) = 0; # fixed, Pi 
xdot(163) = 0; # fixed, kf_ATP 
xdot(184) = 0; # fixed, O2 
xdot(210) = 0; # fixed, CO2_ext 
xdot(128) = x(128) = x(163)*x(96); # assignment, light_1  FROM INIT EXPRESSION!!! 
xdot(129) = x(129) = x(18)*x(97); # assignment, light_2  FROM INIT EXPRESSION!!! 
xdot(130) = x(130) = x(17)*x(98); # assignment, co2_diff  FROM INIT EXPRESSION!!! 
xdot(131) = x(131) = x(52)*x(99); # assignment, cc_1  FROM INIT EXPRESSION!!! 
xdot(132) = x(132) = x(57)*x(100); # assignment, cc_2  FROM INIT EXPRESSION!!! 
xdot(133) = x(133) = x(58)*x(101); # assignment, cc_3  FROM INIT EXPRESSION!!! 
xdot(134) = x(134) = x(59)*x(102); # assignment, cc_4  FROM INIT EXPRESSION!!! 
xdot(135) = x(135) = x(60)*x(103); # assignment, cc_5  FROM INIT EXPRESSION!!! 
xdot(136) = x(136) = x(61)*x(104); # assignment, cc_6  FROM INIT EXPRESSION!!! 
xdot(137) = x(137) = x(62)*x(105); # assignment, cc_7  FROM INIT EXPRESSION!!! 
xdot(138) = x(138) = x(63)*x(106); # assignment, cc_8  FROM INIT EXPRESSION!!! 
xdot(139) = x(139) = x(64)*x(107); # assignment, cc_9  FROM INIT EXPRESSION!!! 
xdot(140) = x(140) = x(53)*x(108); # assignment, cc_10  FROM INIT EXPRESSION!!! 
xdot(141) = x(141) = x(54)*x(109); # assignment, cc_11  FROM INIT EXPRESSION!!! 
xdot(142) = x(142) = x(55)*x(110); # assignment, cc_12  FROM INIT EXPRESSION!!! 
xdot(143) = x(143) = x(56)*x(111); # assignment, cc_13  FROM INIT EXPRESSION!!! 
xdot(144) = x(144) = x(68)*x(112); # assignment, gl_1  FROM INIT EXPRESSION!!! 
xdot(145) = x(145) = x(94)*x(113); # assignment, gl_2a  FROM INIT EXPRESSION!!! 
xdot(146) = x(146) = x(95)*x(114); # assignment, gl_2b  FROM INIT EXPRESSION!!! 
xdot(147) = x(147) = x(91)*x(115); # assignment, gl_3a  FROM INIT EXPRESSION!!! 
xdot(148) = x(148) = x(92)*x(116); # assignment, gl_3b  FROM INIT EXPRESSION!!! 
xdot(149) = x(149) = x(93)*x(117); # assignment, gl_3c  FROM INIT EXPRESSION!!! 
xdot(150) = x(150) = x(90)*x(118); # assignment, gl_4  FROM INIT EXPRESSION!!! 
xdot(151) = x(151) = x(69)*x(119); # assignment, pp_1  FROM INIT EXPRESSION!!! 
xdot(152) = x(152) = x(70)*x(120); # assignment, pp_2a  FROM INIT EXPRESSION!!! 
xdot(153) = x(153) = x(71)*x(121); # assignment, pp_2b  FROM INIT EXPRESSION!!! 
xdot(154) = x(154) = x(72)*x(122); # assignment, pp_3  FROM INIT EXPRESSION!!! 
xdot(155) = x(155) = x(73)*x(123); # assignment, pp_4  FROM INIT EXPRESSION!!! 
xdot(156) = x(156) = x(74)*x(124); # assignment, pp_5  FROM INIT EXPRESSION!!! 
xdot(157) = x(157) = x(75)*x(125); # assignment, pp_6  FROM INIT EXPRESSION!!! 
xdot(158) = x(158) = x(65)*x(126); # assignment, gsm_1  FROM INIT EXPRESSION!!! 
xdot(159) = x(159) = x(84)*x(127); # assignment, ss_1  FROM INIT EXPRESSION!!! 
xdot(164) = x(164) = x(161)/(x(161)+x(160)); # assignment, ATPr  FROM INIT EXPRESSION!!! 
xdot(180) = x(180) = x(204)+x(201)+x(202)+x(205)+x(197)+x(198)+x(199)+x(200); # assignment, Mass  FROM INIT EXPRESSION!!! 
xdot(182) = x(182) = x(181)/(x(181)+x(183)); # assignment, NADPHr  FROM INIT EXPRESSION!!! 
xdot(160) =  + ((x(143)*x(161)*x(192))/((x(161)*(1+x(160)/2.5)+x(24)*(1+x(160)/0.4))*(x(192)+0.05*(1+x(188)/2+x(193)/0.7+x(162)/4))))  + ((x(132)*(x(188)*x(161)-(x(165)*x(160))/0.000619))/((x(188)+0.24)*(x(161)+0.39*(1+x(160)/0.23))))  - (x(128)*x(160) * x(162))  + ((x(76)*(x(161)*x(175)-(x(160)*x(188))/300))/((x(175)+0.25)*(x(161)+0.21*(1+x(188)/0.36))))  + ((x(76)*(x(161)*x(175)-(x(160)*x(190))/300))/((x(175)+0.25)*(x(161)+0.21*(1+x(190)/0.36)))) ; # reaction, ADP 
xdot(161) =  - ((x(143)*x(161)*x(192))/((x(161)*(1+x(160)/2.5)+x(24)*(1+x(160)/0.4))*(x(192)+0.05*(1+x(188)/2+x(193)/0.7+x(162)/4))))  - ((x(132)*(x(188)*x(161)-(x(165)*x(160))/0.000619))/((x(188)+0.24)*(x(161)+0.39*(1+x(160)/0.23))))  + (x(128)*x(160) * x(162))  - ((x(76)*(x(161)*x(175)-(x(160)*x(188))/300))/((x(175)+0.25)*(x(161)+0.21*(1+x(188)/0.36))))  - ((x(76)*(x(161)*x(175)-(x(160)*x(190))/300))/((x(175)+0.25)*(x(161)+0.21*(1+x(190)/0.36)))) ; # reaction, ATP 
xdot(165) =  + ((x(132)*(x(188)*x(161)-(x(165)*x(160))/0.000619))/((x(188)+0.24)*(x(161)+0.39*(1+x(160)/0.23))))  - ((x(133)*x(165)*x(181))/((x(165)+0.004)*(x(181)+0.1))) ; # reaction, BPGA 
xdot(166) =  - ((x(193)*(x(131)*x(166)/(x(166)+0.0115*(1+x(184)/0.222))))/(x(193)+0.02*(1+x(188)/0.84+x(170)/0.04+x(195)/0.075+x(162)/0.9+x(181)/0.07)))  + (x(130)*x(210))  - ((x(193)*(x(151)*x(166)/(x(166)+0.0115*(1+x(184)/0.222))))/(x(193)+0.02*(1+x(188)/0.84+x(170)/0.08+x(195)/0.075+x(162)/0.9+x(181)/0.07))) ; # reaction, CO2_cyt 
xdot(167) =  + ((x(134)*(x(173)-x(167)/x(8)))/(x(173)+x(167)))  - ((x(135)*(x(173)*x(167)-x(170)/x(9)))/(x(25)*x(26)*(1+x(173)/x(25)+x(167)/x(26)+x(170)/x(27)+x(173)*x(167)/(x(25)*x(26)))))  - ((x(138)*(x(167)*x(168)-x(195)/x(11)))/((x(168)+x(35))*(x(167)+x(34))))  - ((x(77)*x(167))/(0.6 + x(167))) ; # reaction, DHAP 
xdot(168) =  + 0 * ((x(140)*(x(173)*x(194)-(x(191)*x(209))/x(5)))/((x(30)*x(21)*(1+(1+x(173)/0.27)*(x(169)/0.54+x(194)/0.09)+x(173)/0.09+((x(209)*(1+(x(168)+x(191))/x(30))+x(168)+x(191)))/x(21)))))  + ((x(137)*(x(169)*x(173)-(x(209)*x(168))/x(10)))/((x(169)+x(32)*(1+x(209)/x(29)+x(168)/x(31)))*(x(173)+x(33))))  - ((x(138)*(x(167)*x(168)-x(195)/x(11)))/((x(168)+x(35))*(x(167)+x(34))))  - ((x(78)*x(168))/(0.28 + x(168))) ; # reaction, E4P 
xdot(169) =  + 0 * ((x(140)*(x(173)*x(194)-(x(191)*x(209))/x(5)))/((x(30)*x(21)*(1+(1+x(173)/0.27)*(x(169)/0.54+x(194)/0.09)+x(173)/0.09+((x(209)*(1+(x(168)+x(191))/x(30))+x(168)+x(191)))/x(21)))))  + ((x(136)*x(170))/(x(170)+x(28)*(1+x(169)/0.7+x(162)/12)))  - ((x(137)*(x(169)*x(173)-(x(209)*x(168))/x(10)))/((x(169)+x(32)*(1+x(209)/x(29)+x(168)/x(31)))*(x(173)+x(33))))  - ((x(144)*x(169))/(x(37) + x(169)))  - ((x(159)*(x(12)*x(169)-x(172)))/(x(169)+x(172))) ; # reaction, F6P 
xdot(170) =  + 0 * ((x(193)*(x(131)*x(166)/(x(166)+0.0115*(1+x(184)/0.222))))/(x(193)+0.02*(1+x(188)/0.84+x(170)/0.04+x(195)/0.075+x(162)/0.9+x(181)/0.07)))  + ((x(135)*(x(173)*x(167)-x(170)/x(9)))/(x(25)*x(26)*(1+x(173)/x(25)+x(167)/x(26)+x(170)/x(27)+x(173)*x(167)/(x(25)*x(26)))))  - ((x(136)*x(170))/(x(170)+x(28)*(1+x(169)/0.7+x(162)/12)))  + ((x(144)*x(169))/(x(37) + x(169)))  + 0 * ((x(193)*(x(151)*x(166)/(x(166)+0.0115*(1+x(184)/0.222))))/(x(193)+0.02*(1+x(188)/0.84+x(170)/0.08+x(195)/0.075+x(162)/0.9+x(181)/0.07))) ; # reaction, FBP 
xdot(171) = 0; # not entered !! ; # reaction, G1P 
xdot(172) =  + ((x(159)*(x(12)*x(169)-x(172)))/(x(169)+x(172))) ; # reaction, G6P 
xdot(173) =  - ((x(140)*(x(173)*x(194)-(x(191)*x(209))/x(5)))/((x(30)*x(21)*(1+(1+x(173)/0.27)*(x(169)/0.54+x(194)/0.09)+x(173)/0.09+((x(209)*(1+(x(168)+x(191))/x(30))+x(168)+x(191)))/x(21)))))  + ((x(133)*x(165)*x(181))/((x(165)+0.004)*(x(181)+0.1)))  - ((x(134)*(x(173)-x(167)/x(8)))/(x(173)+x(167)))  - ((x(135)*(x(173)*x(167)-x(170)/x(9)))/(x(25)*x(26)*(1+x(173)/x(25)+x(167)/x(26)+x(170)/x(27)+x(173)*x(167)/(x(25)*x(26)))))  - ((x(137)*(x(169)*x(173)-(x(209)*x(168))/x(10)))/((x(169)+x(32)*(1+x(209)/x(29)+x(168)/x(31)))*(x(173)+x(33))))  - ((x(145)*(x(173)*x(183)-x(15)*x(188)*x(181)))/((1+x(173)/x(19)+x(188)/x(50))*(1+x(183)/x(48)+x(181)/x(46))))  - ((x(146)*(x(173)*x(183)-x(15)*x(188)*x(181)))/((1+x(173)/x(20)+x(188)/x(51))*(1+x(183)/x(49)+x(181)/x(47)))) ; # reaction, GAP 
xdot(174) =  + ((x(152)*x(189))/(x(189)+x(22)*(1+x(174)/94+x(162)/2.55)))  + ((x(153)*x(189))/(x(189)+x(23)*(1+x(174)/94+x(162)/2.55)))  - ((x(154)*x(174))/(0.1 + x(174))) ; # reaction, GCA 
xdot(175) =  + ((x(89)*x(208))/(0.1 + x(208)))  + ((x(157)*(x(179)-x(175)/250000))/(x(179)+0.09*(1+x(179)/12)))  - ((x(76)*(x(161)*x(175)-(x(160)*x(188))/300))/((x(175)+0.25)*(x(161)+0.21*(1+x(188)/0.36))))  - ((x(76)*(x(161)*x(175)-(x(160)*x(190))/300))/((x(175)+0.25)*(x(161)+0.21*(1+x(190)/0.36)))) ; # reaction, GCEA 
xdot(176) =  + ((x(158)*x(178))/(x(3) + x(178)))  + ((x(155)*(x(178)*x(196)-(x(179)*x(176))/607))/((x(178)+0.15)*(x(196)+1.7*(1+x(176)/2))))  -  2 * ((x(156)*x(176))/(6 + x(176)))  - ((x(79)*x(176))/(1 + x(176))) ; # reaction, GLY 
xdot(177) = 0; # not entered !! ; # reaction, Glycogen 
xdot(178) =  -  2 * ((x(87)*x(178))/(0.1 + x(178)))  - ((x(158)*x(178))/(x(3) + x(178)))  - ((x(66)*x(178))/(2 + x(178)))  + ((x(154)*x(174))/(0.1 + x(174)))  - ((x(155)*(x(178)*x(196)-(x(179)*x(176))/607))/((x(178)+0.15)*(x(196)+1.7*(1+x(176)/2)))) ; # reaction, GOA 
xdot(179) =  + ((x(155)*(x(178)*x(196)-(x(179)*x(176))/607))/((x(178)+0.15)*(x(196)+1.7*(1+x(176)/2))))  - ((x(157)*(x(179)-x(175)/250000))/(x(179)+0.09*(1+x(179)/12))) ; # reaction, HPR 
xdot(181) =  + 0 * ((x(193)*(x(131)*x(166)/(x(166)+0.0115*(1+x(184)/0.222))))/(x(193)+0.02*(1+x(188)/0.84+x(170)/0.04+x(195)/0.075+x(162)/0.9+x(181)/0.07)))  - ((x(133)*x(165)*x(181))/((x(165)+0.004)*(x(181)+0.1)))  + ((x(145)*(x(173)*x(183)-x(15)*x(188)*x(181)))/((1+x(173)/x(19)+x(188)/x(50))*(1+x(183)/x(48)+x(181)/x(46))))  + ((x(146)*(x(173)*x(183)-x(15)*x(188)*x(181)))/((1+x(173)/x(20)+x(188)/x(51))*(1+x(183)/x(49)+x(181)/x(47))))  + (x(129)*x(183))  + 0 * ((x(193)*(x(151)*x(166)/(x(166)+0.0115*(1+x(184)/0.222))))/(x(193)+0.02*(1+x(188)/0.84+x(170)/0.08+x(195)/0.075+x(162)/0.9+x(181)/0.07))) ; # reaction, NADPH 
xdot(183) =  + ((x(133)*x(165)*x(181))/((x(165)+0.004)*(x(181)+0.1)))  - ((x(145)*(x(173)*x(183)-x(15)*x(188)*x(181)))/((1+x(173)/x(19)+x(188)/x(50))*(1+x(183)/x(48)+x(181)/x(46))))  - ((x(146)*(x(173)*x(183)-x(15)*x(188)*x(181)))/((1+x(173)/x(20)+x(188)/x(51))*(1+x(183)/x(49)+x(181)/x(47))))  - (x(129)*x(183)) ; # reaction, NADPp 
xdot(185) = 0; # not entered !! ; # reaction, out 
xdot(186) =  + ((x(66)*x(178))/(2 + x(178)))  - ((x(67)*x(186))/(5 + x(186))) ; # reaction, OXA 
xdot(187) =  + ((x(150)*(x(190)-x(187)/x(14)))/(x(42)*(1+x(190)/x(42)+x(187)/x(38))))  - ((x(80)*x(187))/(x(4) + x(187))) ; # reaction, PEP 
xdot(188) =  +  2 * ((x(193)*(x(131)*x(166)/(x(166)+0.0115*(1+x(184)/0.222))))/(x(193)+0.02*(1+x(188)/0.84+x(170)/0.04+x(195)/0.075+x(162)/0.9+x(181)/0.07)))  + 0 * ((x(143)*x(161)*x(192))/((x(161)*(1+x(160)/2.5)+x(24)*(1+x(160)/0.4))*(x(192)+0.05*(1+x(188)/2+x(193)/0.7+x(162)/4))))  - ((x(132)*(x(188)*x(161)-(x(165)*x(160))/0.000619))/((x(188)+0.24)*(x(161)+0.39*(1+x(160)/0.23))))  + ((x(145)*(x(173)*x(183)-x(15)*x(188)*x(181)))/((1+x(173)/x(19)+x(188)/x(50))*(1+x(183)/x(48)+x(181)/x(46))))  + ((x(146)*(x(173)*x(183)-x(15)*x(188)*x(181)))/((1+x(173)/x(20)+x(188)/x(51))*(1+x(183)/x(49)+x(181)/x(47))))  - ((x(147)*(x(188)-x(190)/x(16)))/(x(43)*(1+x(188)/x(43)+x(190)/x(39))))  - ((x(148)*(x(188)-x(190)/x(16)))/(x(44)*(1+x(188)/x(44)+x(190)/x(40))))  - ((x(149)*(x(188)-x(190)/x(16)))/(x(45)*(1+x(188)/x(45)+x(190)/x(41))))  - ((x(85)*x(188))/(2 + x(188)))  + 0 * ((x(193)*(x(151)*x(166)/(x(166)+0.0115*(1+x(184)/0.222))))/(x(193)+0.02*(1+x(188)/0.84+x(170)/0.08+x(195)/0.075+x(162)/0.9+x(181)/0.07)))  + ((x(76)*(x(161)*x(175)-(x(160)*x(188))/300))/((x(175)+0.25)*(x(161)+0.21*(1+x(188)/0.36)))) ; # reaction, PGA 
xdot(189) =  + ((x(193)*(x(151)*x(166)/(x(166)+0.0115*(1+x(184)/0.222))))/(x(193)+0.02*(1+x(188)/0.84+x(170)/0.08+x(195)/0.075+x(162)/0.9+x(181)/0.07)))  - ((x(152)*x(189))/(x(189)+x(22)*(1+x(174)/94+x(162)/2.55)))  - ((x(153)*x(189))/(x(189)+x(23)*(1+x(174)/94+x(162)/2.55))) ; # reaction, PGCA 
xdot(190) =  + ((x(147)*(x(188)-x(190)/x(16)))/(x(43)*(1+x(188)/x(43)+x(190)/x(39))))  + ((x(148)*(x(188)-x(190)/x(16)))/(x(44)*(1+x(188)/x(44)+x(190)/x(40))))  + ((x(149)*(x(188)-x(190)/x(16)))/(x(45)*(1+x(188)/x(45)+x(190)/x(41))))  - ((x(150)*(x(190)-x(187)/x(14)))/(x(42)*(1+x(190)/x(42)+x(187)/x(38))))  + ((x(76)*(x(161)*x(175)-(x(160)*x(190))/300))/((x(175)+0.25)*(x(161)+0.21*(1+x(190)/0.36)))) ; # reaction, ppGA 
xdot(191) =  + ((x(140)*(x(173)*x(194)-(x(191)*x(209))/x(5)))/((x(30)*x(21)*(1+(1+x(173)/0.27)*(x(169)/0.54+x(194)/0.09)+x(173)/0.09+((x(209)*(1+(x(168)+x(191))/x(30))+x(168)+x(191)))/x(21)))))  - ((x(141)*(x(191)-x(192)/x(6)))/(x(191)+x(192)))  - ((x(81)*x(191))/(0.61 + x(191))) ; # reaction, Ri5P 
xdot(192) =  + ((x(141)*(x(191)-x(192)/x(6)))/(x(191)+x(192)))  + ((x(142)*(x(209)-x(192)/x(7)))/(x(209)+x(192)))  - ((x(143)*x(161)*x(192))/((x(161)*(1+x(160)/2.5)+x(24)*(1+x(160)/0.4))*(x(192)+0.05*(1+x(188)/2+x(193)/0.7+x(162)/4)))) ; # reaction, Ru5P 
xdot(193) =  - ((x(193)*(x(131)*x(166)/(x(166)+0.0115*(1+x(184)/0.222))))/(x(193)+0.02*(1+x(188)/0.84+x(170)/0.04+x(195)/0.075+x(162)/0.9+x(181)/0.07)))  + ((x(143)*x(161)*x(192))/((x(161)*(1+x(160)/2.5)+x(24)*(1+x(160)/0.4))*(x(192)+0.05*(1+x(188)/2+x(193)/0.7+x(162)/4))))  - ((x(193)*(x(151)*x(166)/(x(166)+0.0115*(1+x(184)/0.222))))/(x(193)+0.02*(1+x(188)/0.84+x(170)/0.08+x(195)/0.075+x(162)/0.9+x(181)/0.07))) ; # reaction, RuBP 
xdot(194) =  - ((x(140)*(x(173)*x(194)-(x(191)*x(209))/x(5)))/((x(30)*x(21)*(1+(1+x(173)/0.27)*(x(169)/0.54+x(194)/0.09)+x(173)/0.09+((x(209)*(1+(x(168)+x(191))/x(30))+x(168)+x(191)))/x(21)))))  + ((x(139)*x(195))/(x(195)+x(36)*(1+x(162)/12))) ; # reaction, S7P 
xdot(195) =  + 0 * ((x(193)*(x(131)*x(166)/(x(166)+0.0115*(1+x(184)/0.222))))/(x(193)+0.02*(1+x(188)/0.84+x(170)/0.04+x(195)/0.075+x(162)/0.9+x(181)/0.07)))  + ((x(138)*(x(167)*x(168)-x(195)/x(11)))/((x(168)+x(35))*(x(167)+x(34))))  - ((x(139)*x(195))/(x(195)+x(36)*(1+x(162)/12)))  + 0 * ((x(193)*(x(151)*x(166)/(x(166)+0.0115*(1+x(184)/0.222))))/(x(193)+0.02*(1+x(188)/0.84+x(170)/0.08+x(195)/0.075+x(162)/0.9+x(181)/0.07))) ; # reaction, SBP 
xdot(196) =  + ((x(85)*x(188))/(2 + x(188)))  - ((x(155)*(x(178)*x(196)-(x(179)*x(176))/607))/((x(178)+0.15)*(x(196)+1.7*(1+x(176)/2))))  + ((x(156)*x(176))/(6 + x(176)))  - ((x(82)*x(196))/(1 + x(196))) ; # reaction, SER 
xdot(197) = 0; # not entered !! ; # reaction, Sink1 
xdot(198) =  + ((x(79)*x(176))/(1 + x(176)))  + ((x(82)*x(196))/(1 + x(196)))  + ((x(83)*x(207))/(1 + x(207))) ; # reaction, Sink2 
xdot(199) = 0; # not entered !! ; # reaction, Sink3 
xdot(200) =  + ((x(67)*x(186))/(5 + x(186))) ; # reaction, Sink4 
xdot(201) =  +  3 * ((x(77)*x(167))/(0.6 + x(167))) ; # reaction, Sink_DHAP 
xdot(202) =  +  4 * ((x(78)*x(168))/(0.28 + x(168))) ; # reaction, Sink_E4P 
xdot(203) = 0; # not entered !! ; # reaction, Sink_GAP 
xdot(204) =  +  3 * ((x(80)*x(187))/(x(4) + x(187))) ; # reaction, Sink_PEP 
xdot(205) =  +  5 * ((x(81)*x(191))/(0.61 + x(191))) ; # reaction, Sink_Ri5P 
xdot(206) = 0; # not entered !! ; # reaction, Sucrose 
xdot(207) =  - ((x(83)*x(207))/(1 + x(207))) ; # reaction, THR 
xdot(208) =  + ((x(87)*x(178))/(0.1 + x(178)))  - ((x(89)*x(208))/(0.1 + x(208))) ; # reaction, TSA 
xdot(209) =  + ((x(140)*(x(173)*x(194)-(x(191)*x(209))/x(5)))/((x(30)*x(21)*(1+(1+x(173)/0.27)*(x(169)/0.54+x(194)/0.09)+x(173)/0.09+((x(209)*(1+(x(168)+x(191))/x(30))+x(168)+x(191)))/x(21)))))  - ((x(142)*(x(209)-x(192)/x(7)))/(x(209)+x(192)))  + ((x(137)*(x(169)*x(173)-(x(209)*x(168))/x(10)))/((x(169)+x(32)*(1+x(209)/x(29)+x(168)/x(31)))*(x(173)+x(33)))) ; # reaction, Xu5P 
xdot(128) = 0; # end
xdot(129) = 0; # end
xdot(130) = 0; # end
xdot(131) = 0; # end
xdot(132) = 0; # end
xdot(133) = 0; # end
xdot(134) = 0; # end
xdot(135) = 0; # end
xdot(136) = 0; # end
xdot(137) = 0; # end
xdot(138) = 0; # end
xdot(139) = 0; # end
xdot(140) = 0; # end
xdot(141) = 0; # end
xdot(142) = 0; # end
xdot(143) = 0; # end
xdot(144) = 0; # end
xdot(145) = 0; # end
xdot(146) = 0; # end
xdot(147) = 0; # end
xdot(148) = 0; # end
xdot(149) = 0; # end
xdot(150) = 0; # end
xdot(151) = 0; # end
xdot(152) = 0; # end
xdot(153) = 0; # end
xdot(154) = 0; # end
xdot(155) = 0; # end
xdot(156) = 0; # end
xdot(157) = 0; # end
xdot(158) = 0; # end
xdot(159) = 0; # end
xdot(164) = 0; # end
xdot(180) = 0; # end
xdot(182) = 0; # end

        endfunction

        # time scale
        t = linspace(0, 846000, 200);

        # default values
        x0 = zeros (210,1);
        x0(1) = 1; # fixed compartment, cytosol
x0(2) = 1; # fixed compartment, external
x0(3) = 0.1; # fixed parameter, K_GOA
x0(4) = 7.81; # fixed parameter, K_Sink_PEP
x0(5) = 1.54; # fixed parameter, KE10
x0(6) = 2.35; # fixed parameter, KE11
x0(7) = 0.32; # fixed parameter, KE12
x0(8) = 0.1; # fixed parameter, KE4
x0(9) = 1.57; # fixed parameter, KE5
x0(10) = 0.3; # fixed parameter, KE7
x0(11) = 0.83; # fixed parameter, KE8
x0(12) = 1.24; # fixed parameter, KE_SS1
x0(13) = 0.0001; # fixed parameter, KE_SS2
x0(14) = 1.56; # fixed parameter, Keq_enol
x0(15) = 0.548; # fixed parameter, Keq_gap_dehyd
x0(16) = 0.035; # fixed parameter, Keq_PGM
x0(17) = 1.91; # fixed parameter, kf_CO2
x0(18) = 0.68; # fixed parameter, kf_NADPH
x0(19) = 1; # fixed parameter, Kgap
x0(20) = 0.54; # fixed parameter, Kgap_beta
x0(21) = 0.118; # fixed parameter, KM101
x0(22) = 0.8; # fixed parameter, Km112a
x0(23) = 0.15; # fixed parameter, Km112b
x0(24) = 0.05; # fixed parameter, KM132
x0(25) = 0.3; # fixed parameter, KM51
x0(26) = 0.4; # fixed parameter, KM52
x0(27) = 0.02; # fixed parameter, KM53
x0(28) = 0.033; # fixed parameter, KM61
x0(29) = 0.1; # fixed parameter, KM71
x0(30) = 0.616; # fixed parameter, KM71a
x0(31) = 0.1; # fixed parameter, KM72
x0(32) = 0.1; # fixed parameter, KM73
x0(33) = 0.1; # fixed parameter, KM74
x0(34) = 0.4; # fixed parameter, KM81
x0(35) = 0.2; # fixed parameter, KM82
x0(36) = 0.05; # fixed parameter, KM9
x0(37) = 0.8; # fixed parameter, KM_PFK
x0(38) = 0.95; # fixed parameter, Kmp_enol
x0(39) = 0.32; # fixed parameter, Kmp_PGM_alpha
x0(40) = 3.36; # fixed parameter, Kmp_PGM_beta
x0(41) = 0.6; # fixed parameter, Kmp_PGM_gama
x0(42) = 0.39; # fixed parameter, Kms_enol
x0(43) = 1.2; # fixed parameter, Kms_PGM_alpha
x0(44) = 0.57; # fixed parameter, Kms_PGM_beta
x0(45) = 2.58; # fixed parameter, Kms_PGM_gama
x0(46) = 0.927; # fixed parameter, knadph
x0(47) = 1.184; # fixed parameter, knadph_beta
x0(48) = 2; # fixed parameter, Knadpp
x0(49) = 1.3; # fixed parameter, Knadpp_beta
x0(50) = 0.357; # fixed parameter, Kpga
x0(51) = 1.01; # fixed parameter, Kpga_beta
x0(52) = 0.991; # fixed parameter, V1
x0(53) = 0.23; # fixed parameter, V10
x0(54) = 1.3; # fixed parameter, V11
x0(55) = 1.685; # fixed parameter, V12
x0(56) = 0.43; # fixed parameter, V13
x0(57) = 0.602; # fixed parameter, V2
x0(58) = 0.536; # fixed parameter, V3
x0(59) = 0.1; # fixed parameter, V4
x0(60) = 1.613; # fixed parameter, V5
x0(61) = 0.11; # fixed parameter, V6
x0(62) = 0.3; # fixed parameter, V7
x0(63) = 1.11; # fixed parameter, V8
x0(64) = 0.159; # fixed parameter, V9
x0(65) = 0.1; # fixed parameter, V_GLY_syn
x0(66) = 0.25; # fixed parameter, V_OXA1
x0(67) = 0.0002; # fixed parameter, V_OXA2
x0(68) = 0.0145; # fixed parameter, V_PFK
x0(69) = 0.009; # fixed parameter, V_PP1
x0(70) = 0.001; # fixed parameter, V_PP2a
x0(71) = 0.002; # fixed parameter, V_PP2b
x0(72) = 0.0042; # fixed parameter, V_PP3
x0(73) = 0.43; # fixed parameter, V_PP4
x0(74) = 0.002; # fixed parameter, V_PP5
x0(75) = 0.004; # fixed parameter, V_PP6
x0(76) = 0.05; # fixed parameter, V_PP7
x0(77) = 0; # fixed parameter, V_Sink_DHAP
x0(78) = 0.0009; # fixed parameter, V_Sink_E4P
x0(79) = 0.0004; # fixed parameter, V_Sink_GLY
x0(80) = 0.5; # fixed parameter, V_Sink_PEP
x0(81) = 0.006; # fixed parameter, V_Sink_Ri5P
x0(82) = 0.0021; # fixed parameter, V_Sink_SER
x0(83) = 0.001; # fixed parameter, V_Sink_THR
x0(84) = 0.052; # fixed parameter, V_SS1
x0(85) = 0.0018; # fixed parameter, V_synth_SER
x0(86) = 0.0005; # fixed parameter, V_synth_THR
x0(87) = 0.01; # fixed parameter, V_TSA1
x0(88) = 0.023; # fixed parameter, V_TSA2
x0(89) = 0.1; # fixed parameter, V_TSA3
x0(90) = 1.69; # fixed parameter, Vf_enol
x0(91) = 0.64; # fixed parameter, Vf_PGM_alpha
x0(92) = 1.9; # fixed parameter, Vf_PGM_beta
x0(93) = 1.59; # fixed parameter, Vf_PGM_gama
x0(94) = 0.824; # fixed parameter, Vgap_dehyd
x0(95) = 2.45; # fixed parameter, Vgap_dehyd_beta
x0(96) = 1; # fixed parameter, kf_ATP_LCgene
x0(97) = 1; # fixed parameter, kf_NADPH_LCgene
x0(98) = 1; # fixed parameter, kf_CO2_LCgene
x0(99) = 1; # fixed parameter, V1_LCgene
x0(100) = 1; # fixed parameter, V2_LCgene
x0(101) = 1; # fixed parameter, V3_LCgene
x0(102) = 1; # fixed parameter, V4_LCgene
x0(103) = 1; # fixed parameter, V5_LCgene
x0(104) = 1; # fixed parameter, V6_LCgene
x0(105) = 1; # fixed parameter, V7_LCgene
x0(106) = 1; # fixed parameter, V8_LCgene
x0(107) = 1; # fixed parameter, V9_LCgene
x0(108) = 1; # fixed parameter, V10_LCgene
x0(109) = 1; # fixed parameter, V11_LCgene
x0(110) = 1; # fixed parameter, V12_LCgene
x0(111) = 1; # fixed parameter, V13_LCgene
x0(112) = 1; # fixed parameter, V_PFK_LCgene
x0(113) = 1; # fixed parameter, Vgap_dehyd_LCgene
x0(114) = 1; # fixed parameter, Vgap_dehyd_beta_LCgene
x0(115) = 1; # fixed parameter, Vf_PGM_alpha_LCgene
x0(116) = 1; # fixed parameter, Vf_PGM_beta_LCgene
x0(117) = 1; # fixed parameter, Vf_PGM_gama_LCgene
x0(118) = 1; # fixed parameter, Vf_enol_LCgene
x0(119) = 1; # fixed parameter, V_PP1_LCgene
x0(120) = 1; # fixed parameter, V_PP2a_LCgene
x0(121) = 1; # fixed parameter, V_PP2b_LCgene
x0(122) = 1; # fixed parameter, V_PP3_LCgene
x0(123) = 1; # fixed parameter, V_PP4_LCgene
x0(124) = 1; # fixed parameter, V_PP5_LCgene
x0(125) = 1; # fixed parameter, V_PP6_LCgene
x0(126) = 1; # fixed parameter, V_GLY_syn_LCgene
x0(127) = 1; # fixed parameter, V_SS1_LCgene
x0(160) = 0.3447986077; # fixed specie, ADP
x0(161) = 1.15520132; # fixed specie, ATP
x0(162) = 4.999999759; # fixed specie, Pi
x0(163) = 0.19; # fixed parameter, kf_ATP
x0(165) = 0.007482485571; # fixed specie, BPGA
x0(166) = 0.02; # fixed specie, CO2_cyt
x0(167) = 0.1; # fixed specie, DHAP
x0(168) = 0.2288202267; # fixed specie, E4P
x0(169) = 0.4; # fixed specie, F6P
x0(170) = 0.07617467533; # fixed specie, FBP
x0(171) = 0.01359336237; # fixed specie, G1P
x0(172) = 1.94057631; # fixed specie, G6P
x0(173) = 0.3683626365; # fixed specie, GAP
x0(174) = 0.01235938354; # fixed specie, GCA
x0(175) = 0.004000971532; # fixed specie, GCEA
x0(176) = 0.5; # fixed specie, GLY
x0(177) = 20; # fixed specie, Glycogen
x0(178) = 0.032736383; # fixed specie, GOA
x0(179) = 0.1; # fixed specie, HPR
x0(181) = 0.1489955134; # fixed specie, NADPH
x0(183) = 0.2710044663; # fixed specie, NADPp
x0(184) = 0.02599999875; # fixed specie, O2
x0(185) = 0; # fixed specie, out
x0(186) = 0.1; # fixed specie, OXA
x0(187) = 4.19091919; # fixed specie, PEP
x0(188) = 2.07; # fixed specie, PGA
x0(189) = 0.01; # fixed specie, PGCA
x0(190) = 4.19091918; # fixed specie, ppGA
x0(191) = 0.0744045363; # fixed specie, Ri5P
x0(192) = 0.02689822284; # fixed specie, Ru5P
x0(193) = 0.3105938683; # fixed specie, RuBP
x0(194) = 0.5113194231; # fixed specie, S7P
x0(195) = 1.362618878; # fixed specie, SBP
x0(196) = 0.5; # fixed specie, SER
x0(197) = 0; # fixed specie, Sink1
x0(198) = 0; # fixed specie, Sink2
x0(199) = 0; # fixed specie, Sink3
x0(200) = 0; # fixed specie, Sink4
x0(201) = 0; # fixed specie, Sink_DHAP
x0(202) = 0; # fixed specie, Sink_E4P
x0(203) = 0; # fixed specie, Sink_GAP
x0(204) = 0; # fixed specie, Sink_PEP
x0(205) = 0; # fixed specie, Sink_Ri5P
x0(206) = 1.5; # fixed specie, Sucrose
x0(207) = 0.3; # fixed specie, THR
x0(208) = 0.03243314432; # fixed specie, TSA
x0(209) = 0.05580778704; # fixed specie, Xu5P
x0(210) = 0.1; # fixed specie, CO2_ext
x0(128) = x0(163)*x0(96); # assignment parameter, light_1
x0(129) = x0(18)*x0(97); # assignment parameter, light_2
x0(130) = x0(17)*x0(98); # assignment parameter, co2_diff
x0(131) = x0(52)*x0(99); # assignment parameter, cc_1
x0(132) = x0(57)*x0(100); # assignment parameter, cc_2
x0(133) = x0(58)*x0(101); # assignment parameter, cc_3
x0(134) = x0(59)*x0(102); # assignment parameter, cc_4
x0(135) = x0(60)*x0(103); # assignment parameter, cc_5
x0(136) = x0(61)*x0(104); # assignment parameter, cc_6
x0(137) = x0(62)*x0(105); # assignment parameter, cc_7
x0(138) = x0(63)*x0(106); # assignment parameter, cc_8
x0(139) = x0(64)*x0(107); # assignment parameter, cc_9
x0(140) = x0(53)*x0(108); # assignment parameter, cc_10
x0(141) = x0(54)*x0(109); # assignment parameter, cc_11
x0(142) = x0(55)*x0(110); # assignment parameter, cc_12
x0(143) = x0(56)*x0(111); # assignment parameter, cc_13
x0(144) = x0(68)*x0(112); # assignment parameter, gl_1
x0(145) = x0(94)*x0(113); # assignment parameter, gl_2a
x0(146) = x0(95)*x0(114); # assignment parameter, gl_2b
x0(147) = x0(91)*x0(115); # assignment parameter, gl_3a
x0(148) = x0(92)*x0(116); # assignment parameter, gl_3b
x0(149) = x0(93)*x0(117); # assignment parameter, gl_3c
x0(150) = x0(90)*x0(118); # assignment parameter, gl_4
x0(151) = x0(69)*x0(119); # assignment parameter, pp_1
x0(152) = x0(70)*x0(120); # assignment parameter, pp_2a
x0(153) = x0(71)*x0(121); # assignment parameter, pp_2b
x0(154) = x0(72)*x0(122); # assignment parameter, pp_3
x0(155) = x0(73)*x0(123); # assignment parameter, pp_4
x0(156) = x0(74)*x0(124); # assignment parameter, pp_5
x0(157) = x0(75)*x0(125); # assignment parameter, pp_6
x0(158) = x0(65)*x0(126); # assignment parameter, gsm_1
x0(159) = x0(84)*x0(127); # assignment parameter, ss_1
x0(164) = x0(161)/(x0(161)+x0(160)); # assignment specie, ATPr
x0(180) = x0(204)+x0(201)+x0(202)+x0(205)+x0(197)+x0(198)+x0(199)+x0(200); # assignment specie, Mass
x0(182) = x0(181)/(x0(181)+x0(183)); # assignment specie, NADPHr


        #lsode_options ( 'absolute tolerance', max (50*eps, 0.5e-28) );
        lsode_options ( 'relative tolerance', max (50*eps, 0.5e-28) ); 

        y = lsode ("f", x0, t);
        
        for i = 1:length(t)
            y(i, 128) = y(i, 163)*y(i, 96); # assignment, light_1  FROM INIT EXPRESSION!!! 
y(i, 129) = y(i, 18)*y(i, 97); # assignment, light_2  FROM INIT EXPRESSION!!! 
y(i, 130) = y(i, 17)*y(i, 98); # assignment, co2_diff  FROM INIT EXPRESSION!!! 
y(i, 131) = y(i, 52)*y(i, 99); # assignment, cc_1  FROM INIT EXPRESSION!!! 
y(i, 132) = y(i, 57)*y(i, 100); # assignment, cc_2  FROM INIT EXPRESSION!!! 
y(i, 133) = y(i, 58)*y(i, 101); # assignment, cc_3  FROM INIT EXPRESSION!!! 
y(i, 134) = y(i, 59)*y(i, 102); # assignment, cc_4  FROM INIT EXPRESSION!!! 
y(i, 135) = y(i, 60)*y(i, 103); # assignment, cc_5  FROM INIT EXPRESSION!!! 
y(i, 136) = y(i, 61)*y(i, 104); # assignment, cc_6  FROM INIT EXPRESSION!!! 
y(i, 137) = y(i, 62)*y(i, 105); # assignment, cc_7  FROM INIT EXPRESSION!!! 
y(i, 138) = y(i, 63)*y(i, 106); # assignment, cc_8  FROM INIT EXPRESSION!!! 
y(i, 139) = y(i, 64)*y(i, 107); # assignment, cc_9  FROM INIT EXPRESSION!!! 
y(i, 140) = y(i, 53)*y(i, 108); # assignment, cc_10  FROM INIT EXPRESSION!!! 
y(i, 141) = y(i, 54)*y(i, 109); # assignment, cc_11  FROM INIT EXPRESSION!!! 
y(i, 142) = y(i, 55)*y(i, 110); # assignment, cc_12  FROM INIT EXPRESSION!!! 
y(i, 143) = y(i, 56)*y(i, 111); # assignment, cc_13  FROM INIT EXPRESSION!!! 
y(i, 144) = y(i, 68)*y(i, 112); # assignment, gl_1  FROM INIT EXPRESSION!!! 
y(i, 145) = y(i, 94)*y(i, 113); # assignment, gl_2a  FROM INIT EXPRESSION!!! 
y(i, 146) = y(i, 95)*y(i, 114); # assignment, gl_2b  FROM INIT EXPRESSION!!! 
y(i, 147) = y(i, 91)*y(i, 115); # assignment, gl_3a  FROM INIT EXPRESSION!!! 
y(i, 148) = y(i, 92)*y(i, 116); # assignment, gl_3b  FROM INIT EXPRESSION!!! 
y(i, 149) = y(i, 93)*y(i, 117); # assignment, gl_3c  FROM INIT EXPRESSION!!! 
y(i, 150) = y(i, 90)*y(i, 118); # assignment, gl_4  FROM INIT EXPRESSION!!! 
y(i, 151) = y(i, 69)*y(i, 119); # assignment, pp_1  FROM INIT EXPRESSION!!! 
y(i, 152) = y(i, 70)*y(i, 120); # assignment, pp_2a  FROM INIT EXPRESSION!!! 
y(i, 153) = y(i, 71)*y(i, 121); # assignment, pp_2b  FROM INIT EXPRESSION!!! 
y(i, 154) = y(i, 72)*y(i, 122); # assignment, pp_3  FROM INIT EXPRESSION!!! 
y(i, 155) = y(i, 73)*y(i, 123); # assignment, pp_4  FROM INIT EXPRESSION!!! 
y(i, 156) = y(i, 74)*y(i, 124); # assignment, pp_5  FROM INIT EXPRESSION!!! 
y(i, 157) = y(i, 75)*y(i, 125); # assignment, pp_6  FROM INIT EXPRESSION!!! 
y(i, 158) = y(i, 65)*y(i, 126); # assignment, gsm_1  FROM INIT EXPRESSION!!! 
y(i, 159) = y(i, 84)*y(i, 127); # assignment, ss_1  FROM INIT EXPRESSION!!! 
y(i, 164) = y(i, 161)/(y(i, 161)+y(i, 160)); # assignment, ATPr  FROM INIT EXPRESSION!!! 
y(i, 180) = y(i, 204)+y(i, 201)+y(i, 202)+y(i, 205)+y(i, 197)+y(i, 198)+y(i, 199)+y(i, 200); # assignment, Mass  FROM INIT EXPRESSION!!! 
y(i, 182) = y(i, 181)/(y(i, 181)+y(i, 183)); # assignment, NADPHr  FROM INIT EXPRESSION!!! 

        endfor
        
        [t',y]
        