% M-file for model CCP created by BioNetGen 2.0.46 function [t_out,obs_out,x_out]= CCP(tend) Nspecies=355; Nreactions=1043; % Parameters v_R00299=1; v_R00771=0.7; v_R00756=0.85; v_R01068=0.85; v_R01015_f=0.15; v_R01015_b=1.0; v_R01061_f=2.0; v_R01061_b=0.25; v_R01512_f=2.0; v_R01512_b=0.35; v_R01518_f=0.35; v_R01518_b=2.0; v_R00658_f=2.0; v_R00658_b=0.35; v_R00200=1.27; v_R00345=0.97; v_R00341=0.66; v_R00210=1.09; v_R00351=0.89; v_R01324_f=1.0; v_R01324_b=0.11; v_R00267=0.89; v_R01197=0.8; v_R00405=0.8; v_R00412=0.8; v_R01082=0.8; v_R00342=0.72; v_R00214=0.07; v_R00835=0.29; v_R02035=0.29; v_R01528=0.29; v_R01529_f=0.30; v_R01529_b=0.20; v_R01056_f=0.20; v_R01056_b=0.30; v_R01641_f=0.18; v_R01641_b=0.09; v_R01827_f=0.18; v_R01827_b=0.09; v_R01067_f=0.06; v_R01067_b=0.12; v_vs_G6P=0.01; v_vs_F6P=0.01; v_vs_T3P=0.01; v_vs_G3P=0.1; v_vs_PEP=0.06; v_vs_PYR=0.26; v_vs_ACCOA=0.2; v_vs_Ru5P=0.04; v_vs_E4P=0.03; v_vs_OAA=0.14; v_vs_AKG=0.1; v_vi_GLU_unlabeled=1; v_vi_GLU_labeled=0; % Intial concentrations x0= [ 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 1; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0; 0;]; % Stoichiometry matrix % Reaction flux function function f= flux(t,x) f(1)=v_R00299*x(27); f(2)=v_R00771*x(26); f(3)=v_R00756*x(10); f(4)=v_R01068*x(16); f(5)=v_R01015_f*x(28); f(6)=v_R01015_b*x(13); f(7)=v_R01061_f*x(28); f(8)=v_R01061_b*x(19); f(9)=v_R01512_f*x(5); f(10)=v_R01512_b*x(19); f(11)=v_R01518_f*x(12); f(12)=v_R01518_b*x(5); f(13)=v_R00658_f*x(12); f(14)=v_R00658_b*x(21); f(15)=v_R00200*x(21); f(16)=v_R00345*x(21)*x(7); f(17)=v_R00341*x(9); f(18)=v_R00210*x(20)*x(17); f(19)=v_R00351*x(6)*x(9); f(20)=v_R01324_f*x(22); f(21)=v_R01324_b*x(29); f(22)=v_R00267*x(29); f(23)=v_R01197*x(25)*x(17); f(24)=v_R00405*x(15); f(25)=v_R00412*x(23); f(26)=4*v_R01082*x(2); f(27)=v_R00342*x(8); f(28)=v_R00214*x(8); f(29)=v_R00835*x(26); f(30)=v_R02035*x(1); f(31)=v_R01528*x(3); f(32)=v_R01529_f*x(24); f(33)=v_R01529_b*x(11); f(34)=v_R01056_f*x(14); f(35)=v_R01056_b*x(24); f(36)=v_R01641_f*x(4)*x(28); f(37)=v_R01641_b*x(14)*x(11); f(38)=v_R01827_f*x(4)*x(28); f(39)=v_R01827_b*x(18)*x(10); f(40)=v_R01067_f*x(10)*x(28); f(41)=v_R01067_b*x(18)*x(11); f(42)=v_vs_G6P*x(26); f(43)=v_vs_F6P*x(10); f(44)=v_vs_T3P*x(28); f(45)=v_vs_G3P*x(5); f(46)=v_vs_PEP*x(21); f(47)=v_vs_PYR*x(20); f(48)=v_vs_ACCOA*x(6); f(49)=v_vs_Ru5P*x(24); f(50)=v_vs_E4P*x(18); f(51)=v_vs_OAA*x(9); f(52)=v_vs_AKG*x(25); f(53)=v_vi_GLU_unlabeled*x(30); f(54)=v_vi_GLU_labeled*x(30); f(55)=v_R00299*x(32); f(56)=v_R00771*x(33); f(57)=v_R00835*x(33); f(58)=v_vs_G6P*x(33); f(59)=v_R00756*x(34); f(60)=v_R02035*x(35); f(61)=v_R01827_b*x(18)*x(34); f(62)=v_R01067_f*x(34)*x(28); f(63)=v_vs_F6P*x(34); f(64)=v_R01068*x(36); f(65)=v_R01015_f*x(38); f(66)=v_R01061_f*x(38); f(67)=v_R01528*x(37); f(68)=v_R01641_f*x(4)*x(38); f(69)=v_R01827_f*x(4)*x(38); f(70)=v_R01827_b*x(39)*x(10); f(71)=v_R01827_b*x(39)*x(34); f(72)=v_R01067_f*x(10)*x(38); f(73)=v_R01067_f*x(34)*x(38); f(74)=v_R01067_b*x(39)*x(11); f(75)=v_vs_T3P*x(38); f(76)=v_vs_E4P*x(39); f(77)=v_R01015_b*x(40); f(78)=v_R01061_b*x(41); f(79)=v_R01512_b*x(41); f(80)=v_R01529_f*x(42); f(81)=v_R01529_b*x(43); f(82)=v_R01056_b*x(42); f(83)=v_R01641_f*x(44)*x(28); f(84)=v_R01641_f*x(44)*x(38); f(85)=v_R01641_b*x(14)*x(43); f(86)=v_R01827_f*x(44)*x(28); f(87)=v_R01827_f*x(44)*x(38); f(88)=v_R01067_b*x(18)*x(43); f(89)=v_R01067_b*x(39)*x(43); f(90)=v_vs_Ru5P*x(42); f(91)=v_R01512_f*x(45); f(92)=v_R01518_b*x(45); f(93)=v_R01056_f*x(46); f(94)=v_R01641_b*x(46)*x(11); f(95)=v_R01641_b*x(46)*x(43); f(96)=v_vs_G3P*x(45); f(97)=v_R01518_f*x(47); f(98)=v_R00658_f*x(47); f(99)=v_R00658_b*x(48); f(100)=v_R00200*x(48); f(101)=v_R00345*x(48)*x(7); f(102)=v_vs_PEP*x(48); f(103)=v_R00341*x(50); f(104)=v_R00210*x(49)*x(17); f(105)=v_R00351*x(6)*x(50); f(106)=v_vs_PYR*x(49); f(107)=v_vs_OAA*x(50); f(108)=v_R00351*x(51)*x(9); f(109)=v_R00351*x(51)*x(50); f(110)=v_R01324_f*x(52); f(111)=v_vs_ACCOA*x(51); f(112)=v_R01324_f*x(53); f(113)=v_R01324_f*x(54); f(114)=v_R01324_b*x(55); f(115)=v_R00267*x(55); f(116)=v_R01324_b*x(56); f(117)=v_R01324_b*x(57); f(118)=v_R00267*x(56); f(119)=v_R00267*x(57); f(120)=v_R01197*x(58)*x(17); f(121)=v_vs_AKG*x(58); f(122)=v_R01197*x(59)*x(17); f(123)=v_R01197*x(60)*x(17); f(124)=v_R00405*x(61); f(125)=v_vs_AKG*x(59); f(126)=v_vs_AKG*x(60); f(127)=v_R00405*x(62); f(128)=v_R00405*x(63); f(129)=v_R00412*x(64); f(130)=v_R00412*x(65); f(131)=v_R00412*x(66); f(132)=2*v_R01082*x(67); f(133)=2*v_R01082*x(67); f(134)=2*v_R01082*x(68); f(135)=2*v_R01082*x(68); f(136)=v_R01082*x(69); f(137)=v_R01082*x(69); f(138)=v_R01082*x(69); f(139)=v_R01082*x(69); f(140)=v_R00342*x(70); f(141)=v_R00342*x(71); f(142)=v_R00214*x(70); f(143)=v_R00214*x(71); f(144)=v_R00345*x(21)*x(80); f(145)=v_R00345*x(48)*x(80); f(146)=v_R00341*x(78); f(147)=v_R00341*x(79); f(148)=v_R00210*x(81)*x(17); f(149)=v_R00351*x(6)*x(78); f(150)=v_R00351*x(6)*x(79); f(151)=v_R00351*x(51)*x(78); f(152)=v_R00351*x(51)*x(79); f(153)=v_R00342*x(72); f(154)=v_R00342*x(73); f(155)=v_R00342*x(74); f(156)=v_R00342*x(75); f(157)=v_R00342*x(76); f(158)=v_R00342*x(77); f(159)=v_R00214*x(72); f(160)=v_R00214*x(73); f(161)=v_R00214*x(74); f(162)=v_R00214*x(75); f(163)=v_R00214*x(76); f(164)=v_R00214*x(77); f(165)=v_vs_PYR*x(81); f(166)=v_vs_OAA*x(78); f(167)=v_vs_OAA*x(79); f(168)=v_R00658_b*x(83); f(169)=v_R00200*x(83); f(170)=v_R00345*x(83)*x(7); f(171)=v_R00345*x(83)*x(80); f(172)=v_R00341*x(82); f(173)=v_R00341*x(88); f(174)=v_R00341*x(89); f(175)=v_R00341*x(90); f(176)=v_R00341*x(91); f(177)=v_R00210*x(92)*x(17); f(178)=v_R00210*x(93)*x(17); f(179)=v_R00210*x(94)*x(17); f(180)=v_R00351*x(6)*x(82); f(181)=v_R00351*x(6)*x(88); f(182)=v_R00351*x(6)*x(89); f(183)=v_R00351*x(6)*x(90); f(184)=v_R00351*x(6)*x(91); f(185)=v_R00351*x(51)*x(82); f(186)=v_R00351*x(51)*x(88); f(187)=v_R00351*x(51)*x(89); f(188)=v_R00351*x(51)*x(90); f(189)=v_R00351*x(51)*x(91); f(190)=v_R01324_f*x(84); f(191)=v_R01324_f*x(85); f(192)=v_R01324_f*x(86); f(193)=v_R01324_f*x(87); f(194)=v_vs_PEP*x(83); f(195)=v_vs_PYR*x(92); f(196)=v_vs_PYR*x(93); f(197)=v_vs_PYR*x(94); f(198)=v_vs_OAA*x(82); f(199)=v_vs_OAA*x(88); f(200)=v_vs_OAA*x(89); f(201)=v_vs_OAA*x(90); f(202)=v_vs_OAA*x(91); f(203)=v_R01518_f*x(95); f(204)=v_R00658_f*x(95); f(205)=v_R00658_b*x(97); f(206)=v_R00658_b*x(98); f(207)=v_R00658_b*x(99); f(208)=v_R00200*x(97); f(209)=v_R00200*x(98); f(210)=v_R00200*x(99); f(211)=v_R00345*x(97)*x(7); f(212)=v_R00345*x(97)*x(80); f(213)=v_R00345*x(98)*x(7); f(214)=v_R00345*x(98)*x(80); f(215)=v_R00345*x(99)*x(7); f(216)=v_R00345*x(99)*x(80); f(217)=v_R00341*x(96); f(218)=v_R00351*x(6)*x(96); f(219)=v_R00351*x(51)*x(96); f(220)=v_R00351*x(100)*x(9); f(221)=v_R00351*x(100)*x(50); f(222)=v_R00351*x(100)*x(78); f(223)=v_R00351*x(100)*x(79); f(224)=v_R00351*x(100)*x(82); f(225)=v_R00351*x(100)*x(88); f(226)=v_R00351*x(100)*x(89); f(227)=v_R00351*x(100)*x(90); f(228)=v_R00351*x(100)*x(91); f(229)=v_R00351*x(100)*x(96); f(230)=v_R01324_f*x(101); f(231)=v_R01324_f*x(102); f(232)=v_R01324_f*x(103); f(233)=v_R01324_f*x(104); f(234)=v_R01324_f*x(105); f(235)=v_R01324_f*x(106); f(236)=v_R01324_f*x(107); f(237)=v_R01324_f*x(108); f(238)=v_R01324_f*x(109); f(239)=v_R01324_f*x(110); f(240)=v_R01324_b*x(111); f(241)=v_R01324_b*x(112); f(242)=v_R01324_b*x(113); f(243)=v_R01324_b*x(114); f(244)=v_R00267*x(111); f(245)=v_R00267*x(112); f(246)=v_R00267*x(113); f(247)=v_R00267*x(114); f(248)=v_vs_PEP*x(97); f(249)=v_vs_PEP*x(98); f(250)=v_vs_PEP*x(99); f(251)=v_vs_ACCOA*x(100); f(252)=v_vs_OAA*x(96); f(253)=v_R01512_f*x(115); f(254)=v_R01518_f*x(116); f(255)=v_R01518_f*x(117); f(256)=v_R01518_f*x(118); f(257)=v_R01518_b*x(115); f(258)=v_R00658_f*x(116); f(259)=v_R00658_f*x(117); f(260)=v_R00658_f*x(118); f(261)=v_R00341*x(119); f(262)=v_R00341*x(120); f(263)=v_R00351*x(6)*x(119); f(264)=v_R00351*x(6)*x(120); f(265)=v_R00351*x(51)*x(119); f(266)=v_R00351*x(51)*x(120); f(267)=v_R00351*x(100)*x(119); f(268)=v_R00351*x(100)*x(120); f(269)=v_R01324_f*x(121); f(270)=v_R01324_f*x(122); f(271)=v_R01324_f*x(123); f(272)=v_R01324_f*x(124); f(273)=v_R01324_f*x(125); f(274)=v_R01324_f*x(126); f(275)=v_R01324_f*x(127); f(276)=v_R01324_f*x(128); f(277)=v_R01324_f*x(129); f(278)=v_R01324_f*x(130); f(279)=v_R01324_f*x(131); f(280)=v_R01324_f*x(132); f(281)=v_R01324_b*x(133); f(282)=v_R01324_b*x(134); f(283)=v_R01324_b*x(135); f(284)=v_R01324_b*x(136); f(285)=v_R01324_b*x(137); f(286)=v_R01324_b*x(138); f(287)=v_R01324_b*x(139); f(288)=v_R01324_b*x(140); f(289)=v_R01324_b*x(141); f(290)=v_R01324_b*x(142); f(291)=v_R00267*x(133); f(292)=v_R00267*x(134); f(293)=v_R00267*x(135); f(294)=v_R00267*x(136); f(295)=v_R00267*x(137); f(296)=v_R00267*x(138); f(297)=v_R00267*x(139); f(298)=v_R00267*x(140); f(299)=v_R00267*x(141); f(300)=v_R00267*x(142); f(301)=v_R01197*x(143)*x(17); f(302)=v_R01197*x(144)*x(17); f(303)=v_vs_G3P*x(115); f(304)=v_vs_OAA*x(119); f(305)=v_vs_OAA*x(120); f(306)=v_vs_AKG*x(143); f(307)=v_vs_AKG*x(144); f(308)=v_R01061_b*x(145); f(309)=v_R01512_f*x(146); f(310)=v_R01512_f*x(147); f(311)=v_R01512_f*x(148); f(312)=v_R01512_b*x(145); f(313)=v_R01518_b*x(146); f(314)=v_R01518_b*x(147); f(315)=v_R01518_b*x(148); f(316)=v_R01324_f*x(149); f(317)=v_R01324_f*x(150); f(318)=v_R01324_f*x(151); f(319)=v_R01324_f*x(152); f(320)=v_R01324_f*x(153); f(321)=v_R01324_f*x(154); f(322)=v_R01324_b*x(155); f(323)=v_R01324_b*x(156); f(324)=v_R01324_b*x(157); f(325)=v_R01324_b*x(158); f(326)=v_R01324_b*x(159); f(327)=v_R01324_b*x(160); f(328)=v_R01324_b*x(161); f(329)=v_R01324_b*x(162); f(330)=v_R01324_b*x(163); f(331)=v_R01324_b*x(164); f(332)=v_R01324_b*x(165); f(333)=v_R01324_b*x(166); f(334)=v_R00267*x(155); f(335)=v_R00267*x(156); f(336)=v_R00267*x(157); f(337)=v_R00267*x(158); f(338)=v_R00267*x(159); f(339)=v_R00267*x(160); f(340)=v_R00267*x(161); f(341)=v_R00267*x(162); f(342)=v_R00267*x(163); f(343)=v_R00267*x(164); f(344)=v_R00267*x(165); f(345)=v_R00267*x(166); f(346)=v_R01197*x(167)*x(17); f(347)=v_R01197*x(168)*x(17); f(348)=v_R01197*x(169)*x(17); f(349)=v_R01197*x(170)*x(17); f(350)=v_R01197*x(171)*x(17); f(351)=v_R01197*x(172)*x(17); f(352)=v_vs_G3P*x(146); f(353)=v_vs_G3P*x(147); f(354)=v_vs_G3P*x(148); f(355)=v_vs_AKG*x(167); f(356)=v_vs_AKG*x(168); f(357)=v_vs_AKG*x(169); f(358)=v_vs_AKG*x(170); f(359)=v_vs_AKG*x(171); f(360)=v_vs_AKG*x(172); f(361)=v_R01015_f*x(173); f(362)=v_R01061_f*x(173); f(363)=v_R01061_b*x(174); f(364)=v_R01061_b*x(175); f(365)=v_R01061_b*x(176); f(366)=v_R01512_b*x(174); f(367)=v_R01512_b*x(175); f(368)=v_R01512_b*x(176); f(369)=v_R01324_b*x(177); f(370)=v_R01324_b*x(178); f(371)=v_R01324_b*x(179); f(372)=v_R01324_b*x(180); f(373)=v_R01324_b*x(181); f(374)=v_R01324_b*x(182); f(375)=v_R00267*x(177); f(376)=v_R00267*x(178); f(377)=v_R00267*x(179); f(378)=v_R00267*x(180); f(379)=v_R00267*x(181); f(380)=v_R00267*x(182); f(381)=v_R01197*x(183)*x(17); f(382)=v_R01197*x(184)*x(17); f(383)=v_R01197*x(185)*x(17); f(384)=v_R01197*x(186)*x(17); f(385)=v_R01197*x(187)*x(17); f(386)=v_R01197*x(188)*x(17); f(387)=v_R00405*x(189); f(388)=v_R00405*x(190); f(389)=v_R01641_f*x(4)*x(173); f(390)=v_R01641_f*x(44)*x(173); f(391)=v_R01827_f*x(4)*x(173); f(392)=v_R01827_f*x(44)*x(173); f(393)=v_R01067_f*x(10)*x(173); f(394)=v_R01067_f*x(34)*x(173); f(395)=v_vs_T3P*x(173); f(396)=v_vs_AKG*x(183); f(397)=v_vs_AKG*x(184); f(398)=v_vs_AKG*x(185); f(399)=v_vs_AKG*x(186); f(400)=v_vs_AKG*x(187); f(401)=v_vs_AKG*x(188); f(402)=v_R00756*x(200); f(403)=v_R01015_f*x(192); f(404)=v_R01015_f*x(193); f(405)=v_R01015_f*x(194); f(406)=v_R01015_b*x(191); f(407)=v_R01061_f*x(192); f(408)=v_R01061_f*x(193); f(409)=v_R01061_f*x(194); f(410)=v_R00405*x(195); f(411)=v_R00405*x(196); f(412)=v_R00405*x(197); f(413)=v_R00412*x(198); f(414)=v_R01529_b*x(199); f(415)=v_R01641_f*x(4)*x(192); f(416)=v_R01641_f*x(4)*x(193); f(417)=v_R01641_f*x(4)*x(194); f(418)=v_R01641_f*x(44)*x(192); f(419)=v_R01641_f*x(44)*x(193); f(420)=v_R01641_f*x(44)*x(194); f(421)=v_R01641_b*x(14)*x(199); f(422)=v_R01641_b*x(46)*x(199); f(423)=v_R01827_f*x(4)*x(192); f(424)=v_R01827_f*x(4)*x(193); f(425)=v_R01827_f*x(4)*x(194); f(426)=v_R01827_f*x(44)*x(192); f(427)=v_R01827_f*x(44)*x(193); f(428)=v_R01827_f*x(44)*x(194); f(429)=v_R01827_b*x(18)*x(200); f(430)=v_R01827_b*x(39)*x(200); f(431)=v_R01067_f*x(10)*x(192); f(432)=v_R01067_f*x(10)*x(193); f(433)=v_R01067_f*x(10)*x(194); f(434)=v_R01067_f*x(34)*x(192); f(435)=v_R01067_f*x(34)*x(193); f(436)=v_R01067_f*x(34)*x(194); f(437)=v_R01067_f*x(200)*x(28); f(438)=v_R01067_f*x(200)*x(38); f(439)=v_R01067_f*x(200)*x(173); f(440)=v_R01067_f*x(200)*x(192); f(441)=v_R01067_f*x(200)*x(193); f(442)=v_R01067_f*x(200)*x(194); f(443)=v_R01067_b*x(18)*x(199); f(444)=v_R01067_b*x(39)*x(199); f(445)=v_vs_F6P*x(200); f(446)=v_vs_T3P*x(192); f(447)=v_vs_T3P*x(193); f(448)=v_vs_T3P*x(194); f(449)=v_R00756*x(211); f(450)=v_R00756*x(212); f(451)=v_R00756*x(213); f(452)=v_R01068*x(201); f(453)=v_R01015_b*x(202); f(454)=v_R01015_b*x(203); f(455)=v_R01015_b*x(204); f(456)=v_R00412*x(205); f(457)=4*v_R01082*x(206); f(458)=v_R01529_f*x(207); f(459)=v_R01529_b*x(208); f(460)=v_R01529_b*x(209); f(461)=v_R01529_b*x(210); f(462)=v_R01056_b*x(207); f(463)=v_R01641_b*x(14)*x(208); f(464)=v_R01641_b*x(14)*x(209); f(465)=v_R01641_b*x(14)*x(210); f(466)=v_R01641_b*x(46)*x(208); f(467)=v_R01641_b*x(46)*x(209); f(468)=v_R01641_b*x(46)*x(210); f(469)=v_R01827_b*x(18)*x(211); f(470)=v_R01827_b*x(18)*x(212); f(471)=v_R01827_b*x(18)*x(213); f(472)=v_R01827_b*x(39)*x(211); f(473)=v_R01827_b*x(39)*x(212); f(474)=v_R01827_b*x(39)*x(213); f(475)=v_R01827_b*x(214)*x(10); f(476)=v_R01827_b*x(214)*x(34); f(477)=v_R01827_b*x(214)*x(200); f(478)=v_R01827_b*x(214)*x(211); f(479)=v_R01827_b*x(214)*x(212); f(480)=v_R01827_b*x(214)*x(213); f(481)=v_R01067_f*x(211)*x(28); f(482)=v_R01067_f*x(211)*x(38); f(483)=v_R01067_f*x(211)*x(173); f(484)=v_R01067_f*x(211)*x(192); f(485)=v_R01067_f*x(211)*x(193); f(486)=v_R01067_f*x(211)*x(194); f(487)=v_R01067_f*x(212)*x(28); f(488)=v_R01067_f*x(212)*x(38); f(489)=v_R01067_f*x(212)*x(173); f(490)=v_R01067_f*x(212)*x(192); f(491)=v_R01067_f*x(212)*x(193); f(492)=v_R01067_f*x(212)*x(194); f(493)=v_R01067_f*x(213)*x(28); f(494)=v_R01067_f*x(213)*x(38); f(495)=v_R01067_f*x(213)*x(173); f(496)=v_R01067_f*x(213)*x(192); f(497)=v_R01067_f*x(213)*x(193); f(498)=v_R01067_f*x(213)*x(194); f(499)=v_R01067_b*x(18)*x(208); f(500)=v_R01067_b*x(18)*x(209); f(501)=v_R01067_b*x(18)*x(210); f(502)=v_R01067_b*x(39)*x(208); f(503)=v_R01067_b*x(39)*x(209); f(504)=v_R01067_b*x(39)*x(210); f(505)=v_R01067_b*x(214)*x(11); f(506)=v_R01067_b*x(214)*x(43); f(507)=v_R01067_b*x(214)*x(199); f(508)=v_R01067_b*x(214)*x(208); f(509)=v_R01067_b*x(214)*x(209); f(510)=v_R01067_b*x(214)*x(210); f(511)=v_vs_F6P*x(211); f(512)=v_vs_F6P*x(212); f(513)=v_vs_F6P*x(213); f(514)=v_vs_Ru5P*x(207); f(515)=v_vs_E4P*x(214); f(516)=v_R01068*x(215); f(517)=v_R01068*x(216); f(518)=v_R01068*x(217); f(519)=4*v_R01082*x(218); f(520)=v_R00342*x(219); f(521)=v_R00214*x(219); f(522)=v_R01529_f*x(220); f(523)=v_R01529_f*x(221); f(524)=v_R01529_f*x(222); f(525)=v_R01056_f*x(223); f(526)=v_R01056_b*x(220); f(527)=v_R01056_b*x(221); f(528)=v_R01056_b*x(222); f(529)=v_R01641_f*x(224)*x(28); f(530)=v_R01641_f*x(224)*x(38); f(531)=v_R01641_f*x(224)*x(173); f(532)=v_R01641_f*x(224)*x(192); f(533)=v_R01641_f*x(224)*x(193); f(534)=v_R01641_f*x(224)*x(194); f(535)=v_R01641_b*x(223)*x(11); f(536)=v_R01641_b*x(223)*x(43); f(537)=v_R01641_b*x(223)*x(199); f(538)=v_R01641_b*x(223)*x(208); f(539)=v_R01641_b*x(223)*x(209); f(540)=v_R01641_b*x(223)*x(210); f(541)=v_R01827_f*x(224)*x(28); f(542)=v_R01827_f*x(224)*x(38); f(543)=v_R01827_f*x(224)*x(173); f(544)=v_R01827_f*x(224)*x(192); f(545)=v_R01827_f*x(224)*x(193); f(546)=v_R01827_f*x(224)*x(194); f(547)=v_R01827_b*x(225)*x(10); f(548)=v_R01827_b*x(225)*x(34); f(549)=v_R01827_b*x(225)*x(200); f(550)=v_R01827_b*x(225)*x(211); f(551)=v_R01827_b*x(225)*x(212); f(552)=v_R01827_b*x(225)*x(213); f(553)=v_R01827_b*x(226)*x(10); f(554)=v_R01827_b*x(226)*x(34); f(555)=v_R01827_b*x(226)*x(200); f(556)=v_R01827_b*x(226)*x(211); f(557)=v_R01827_b*x(226)*x(212); f(558)=v_R01827_b*x(226)*x(213); f(559)=v_R01827_b*x(227)*x(10); f(560)=v_R01827_b*x(227)*x(34); f(561)=v_R01827_b*x(227)*x(200); f(562)=v_R01827_b*x(227)*x(211); f(563)=v_R01827_b*x(227)*x(212); f(564)=v_R01827_b*x(227)*x(213); f(565)=v_R01067_b*x(225)*x(11); f(566)=v_R01067_b*x(225)*x(43); f(567)=v_R01067_b*x(225)*x(199); f(568)=v_R01067_b*x(225)*x(208); f(569)=v_R01067_b*x(225)*x(209); f(570)=v_R01067_b*x(225)*x(210); f(571)=v_R01067_b*x(226)*x(11); f(572)=v_R01067_b*x(226)*x(43); f(573)=v_R01067_b*x(226)*x(199); f(574)=v_R01067_b*x(226)*x(208); f(575)=v_R01067_b*x(226)*x(209); f(576)=v_R01067_b*x(226)*x(210); f(577)=v_R01067_b*x(227)*x(11); f(578)=v_R01067_b*x(227)*x(43); f(579)=v_R01067_b*x(227)*x(199); f(580)=v_R01067_b*x(227)*x(208); f(581)=v_R01067_b*x(227)*x(209); f(582)=v_R01067_b*x(227)*x(210); f(583)=v_vs_Ru5P*x(220); f(584)=v_vs_Ru5P*x(221); f(585)=v_vs_Ru5P*x(222); f(586)=v_vs_E4P*x(225); f(587)=v_vs_E4P*x(226); f(588)=v_vs_E4P*x(227); f(589)=v_R00341*x(229); f(590)=v_R00210*x(230)*x(17); f(591)=v_R00351*x(6)*x(229); f(592)=v_R00351*x(51)*x(229); f(593)=v_R00351*x(100)*x(229); f(594)=v_R00342*x(228); f(595)=v_R00214*x(228); f(596)=v_R01056_f*x(231); f(597)=v_R01056_f*x(232); f(598)=v_R01056_f*x(233); f(599)=v_R01641_f*x(234)*x(28); f(600)=v_R01641_f*x(234)*x(38); f(601)=v_R01641_f*x(234)*x(173); f(602)=v_R01641_f*x(234)*x(192); f(603)=v_R01641_f*x(234)*x(193); f(604)=v_R01641_f*x(234)*x(194); f(605)=v_R01641_f*x(235)*x(28); f(606)=v_R01641_f*x(235)*x(38); f(607)=v_R01641_f*x(235)*x(173); f(608)=v_R01641_f*x(235)*x(192); f(609)=v_R01641_f*x(235)*x(193); f(610)=v_R01641_f*x(235)*x(194); f(611)=v_R01641_f*x(236)*x(28); f(612)=v_R01641_f*x(236)*x(38); f(613)=v_R01641_f*x(236)*x(173); f(614)=v_R01641_f*x(236)*x(192); f(615)=v_R01641_f*x(236)*x(193); f(616)=v_R01641_f*x(236)*x(194); f(617)=v_R01641_b*x(231)*x(11); f(618)=v_R01641_b*x(231)*x(43); f(619)=v_R01641_b*x(231)*x(199); f(620)=v_R01641_b*x(231)*x(208); f(621)=v_R01641_b*x(231)*x(209); f(622)=v_R01641_b*x(231)*x(210); f(623)=v_R01641_b*x(232)*x(11); f(624)=v_R01641_b*x(232)*x(43); f(625)=v_R01641_b*x(232)*x(199); f(626)=v_R01641_b*x(232)*x(208); f(627)=v_R01641_b*x(232)*x(209); f(628)=v_R01641_b*x(232)*x(210); f(629)=v_R01641_b*x(233)*x(11); f(630)=v_R01641_b*x(233)*x(43); f(631)=v_R01641_b*x(233)*x(199); f(632)=v_R01641_b*x(233)*x(208); f(633)=v_R01641_b*x(233)*x(209); f(634)=v_R01641_b*x(233)*x(210); f(635)=v_R01827_f*x(234)*x(28); f(636)=v_R01827_f*x(234)*x(38); f(637)=v_R01827_f*x(234)*x(173); f(638)=v_R01827_f*x(234)*x(192); f(639)=v_R01827_f*x(234)*x(193); f(640)=v_R01827_f*x(234)*x(194); f(641)=v_R01827_f*x(235)*x(28); f(642)=v_R01827_f*x(235)*x(38); f(643)=v_R01827_f*x(235)*x(173); f(644)=v_R01827_f*x(235)*x(192); f(645)=v_R01827_f*x(235)*x(193); f(646)=v_R01827_f*x(235)*x(194); f(647)=v_R01827_f*x(236)*x(28); f(648)=v_R01827_f*x(236)*x(38); f(649)=v_R01827_f*x(236)*x(173); f(650)=v_R01827_f*x(236)*x(192); f(651)=v_R01827_f*x(236)*x(193); f(652)=v_R01827_f*x(236)*x(194); f(653)=v_vs_PYR*x(230); f(654)=v_vs_OAA*x(229); f(655)=v_R00658_b*x(237); f(656)=v_R00200*x(237); f(657)=v_R00345*x(237)*x(7); f(658)=v_R00345*x(237)*x(80); f(659)=v_R00351*x(238)*x(9); f(660)=v_R00351*x(238)*x(50); f(661)=v_R00351*x(238)*x(78); f(662)=v_R00351*x(238)*x(79); f(663)=v_R00351*x(238)*x(82); f(664)=v_R00351*x(238)*x(88); f(665)=v_R00351*x(238)*x(89); f(666)=v_R00351*x(238)*x(90); f(667)=v_R00351*x(238)*x(91); f(668)=v_R00351*x(238)*x(96); f(669)=v_R00351*x(238)*x(119); f(670)=v_R00351*x(238)*x(120); f(671)=v_R00351*x(238)*x(229); f(672)=v_R01324_f*x(239); f(673)=v_R01324_f*x(240); f(674)=v_R01324_f*x(241); f(675)=v_vs_PEP*x(237); f(676)=v_vs_ACCOA*x(238); f(677)=v_R01518_f*x(242); f(678)=v_R00658_f*x(242); f(679)=v_R00341*x(243); f(680)=v_R00351*x(6)*x(243); f(681)=v_R00351*x(51)*x(243); f(682)=v_R00351*x(100)*x(243); f(683)=v_R00351*x(238)*x(243); f(684)=v_R01324_f*x(244); f(685)=v_R01324_f*x(245); f(686)=v_R01324_f*x(246); f(687)=v_R01324_f*x(247); f(688)=v_R01324_f*x(248); f(689)=v_R01324_f*x(249); f(690)=v_R01324_f*x(250); f(691)=v_R01324_f*x(251); f(692)=v_R01324_f*x(252); f(693)=v_R01324_f*x(253); f(694)=v_R01324_f*x(254); f(695)=v_R01324_f*x(255); f(696)=v_R01324_f*x(256); f(697)=v_R01324_b*x(257); f(698)=v_R01324_b*x(258); f(699)=v_R01324_b*x(259); f(700)=v_R00267*x(257); f(701)=v_R00267*x(258); f(702)=v_R00267*x(259); f(703)=v_vs_OAA*x(243); f(704)=v_R01512_f*x(260); f(705)=v_R01518_b*x(260); f(706)=v_R01324_f*x(261); f(707)=v_R01324_f*x(262); f(708)=v_R01324_f*x(263); f(709)=v_R01324_f*x(264); f(710)=v_R01324_b*x(265); f(711)=v_R01324_b*x(266); f(712)=v_R01324_b*x(267); f(713)=v_R01324_b*x(268); f(714)=v_R01324_b*x(269); f(715)=v_R01324_b*x(270); f(716)=v_R01324_b*x(271); f(717)=v_R01324_b*x(272); f(718)=v_R01324_b*x(273); f(719)=v_R01324_b*x(274); f(720)=v_R01324_b*x(275); f(721)=v_R01324_b*x(276); f(722)=v_R01324_b*x(277); f(723)=v_R00267*x(265); f(724)=v_R00267*x(266); f(725)=v_R00267*x(267); f(726)=v_R00267*x(268); f(727)=v_R00267*x(269); f(728)=v_R00267*x(270); f(729)=v_R00267*x(271); f(730)=v_R00267*x(272); f(731)=v_R00267*x(273); f(732)=v_R00267*x(274); f(733)=v_R00267*x(275); f(734)=v_R00267*x(276); f(735)=v_R00267*x(277); f(736)=v_R01197*x(278)*x(17); f(737)=v_R01197*x(279)*x(17); f(738)=v_R01197*x(280)*x(17); f(739)=v_vs_G3P*x(260); f(740)=v_vs_AKG*x(278); f(741)=v_vs_AKG*x(279); f(742)=v_vs_AKG*x(280); f(743)=v_R01061_b*x(281); f(744)=v_R01512_b*x(281); f(745)=v_R01324_b*x(282); f(746)=v_R01324_b*x(283); f(747)=v_R01324_b*x(284); f(748)=v_R01324_b*x(285); f(749)=v_R00267*x(282); f(750)=v_R00267*x(283); f(751)=v_R00267*x(284); f(752)=v_R00267*x(285); f(753)=v_R01197*x(286)*x(17); f(754)=v_R01197*x(287)*x(17); f(755)=v_R01197*x(288)*x(17); f(756)=v_R01197*x(289)*x(17); f(757)=v_R01197*x(290)*x(17); f(758)=v_R01197*x(291)*x(17); f(759)=v_R01197*x(292)*x(17); f(760)=v_R00405*x(293); f(761)=v_R00405*x(294); f(762)=v_R00405*x(295); f(763)=v_vs_AKG*x(286); f(764)=v_vs_AKG*x(287); f(765)=v_vs_AKG*x(288); f(766)=v_vs_AKG*x(289); f(767)=v_vs_AKG*x(290); f(768)=v_vs_AKG*x(291); f(769)=v_vs_AKG*x(292); f(770)=v_R01015_f*x(296); f(771)=v_R01061_f*x(296); f(772)=v_R01197*x(297)*x(17); f(773)=v_R01197*x(298)*x(17); f(774)=v_R01197*x(299)*x(17); f(775)=v_R01197*x(300)*x(17); f(776)=v_R00405*x(301); f(777)=v_R00405*x(302); f(778)=v_R00405*x(303); f(779)=v_R00405*x(304); f(780)=v_R00412*x(305); f(781)=v_R00412*x(306); f(782)=v_R01641_f*x(4)*x(296); f(783)=v_R01641_f*x(44)*x(296); f(784)=v_R01641_f*x(224)*x(296); f(785)=v_R01641_f*x(234)*x(296); f(786)=v_R01641_f*x(235)*x(296); f(787)=v_R01641_f*x(236)*x(296); f(788)=v_R01827_f*x(4)*x(296); f(789)=v_R01827_f*x(44)*x(296); f(790)=v_R01827_f*x(224)*x(296); f(791)=v_R01827_f*x(234)*x(296); f(792)=v_R01827_f*x(235)*x(296); f(793)=v_R01827_f*x(236)*x(296); f(794)=v_R01067_f*x(10)*x(296); f(795)=v_R01067_f*x(34)*x(296); f(796)=v_R01067_f*x(200)*x(296); f(797)=v_R01067_f*x(211)*x(296); f(798)=v_R01067_f*x(212)*x(296); f(799)=v_R01067_f*x(213)*x(296); f(800)=v_vs_T3P*x(296); f(801)=v_vs_AKG*x(297); f(802)=v_vs_AKG*x(298); f(803)=v_vs_AKG*x(299); f(804)=v_vs_AKG*x(300); f(805)=v_R00756*x(312); f(806)=v_R01015_b*x(307); f(807)=v_R00412*x(308); f(808)=2*v_R01082*x(309); f(809)=2*v_R01082*x(309); f(810)=2*v_R01082*x(310); f(811)=2*v_R01082*x(310); f(812)=v_R01529_b*x(311); f(813)=v_R01641_b*x(14)*x(311); f(814)=v_R01641_b*x(46)*x(311); f(815)=v_R01641_b*x(223)*x(311); f(816)=v_R01641_b*x(231)*x(311); f(817)=v_R01641_b*x(232)*x(311); f(818)=v_R01641_b*x(233)*x(311); f(819)=v_R01827_b*x(18)*x(312); f(820)=v_R01827_b*x(39)*x(312); f(821)=v_R01827_b*x(214)*x(312); f(822)=v_R01827_b*x(225)*x(312); f(823)=v_R01827_b*x(226)*x(312); f(824)=v_R01827_b*x(227)*x(312); f(825)=v_R01067_f*x(312)*x(28); f(826)=v_R01067_f*x(312)*x(38); f(827)=v_R01067_f*x(312)*x(173); f(828)=v_R01067_f*x(312)*x(192); f(829)=v_R01067_f*x(312)*x(193); f(830)=v_R01067_f*x(312)*x(194); f(831)=v_R01067_f*x(312)*x(296); f(832)=v_R01067_b*x(18)*x(311); f(833)=v_R01067_b*x(39)*x(311); f(834)=v_R01067_b*x(214)*x(311); f(835)=v_R01067_b*x(225)*x(311); f(836)=v_R01067_b*x(226)*x(311); f(837)=v_R01067_b*x(227)*x(311); f(838)=v_vs_F6P*x(312); f(839)=v_R01068*x(313); f(840)=4*v_R01082*x(314); f(841)=v_R00342*x(315); f(842)=v_R00342*x(316); f(843)=v_R00342*x(317); f(844)=v_R00342*x(318); f(845)=v_R00214*x(315); f(846)=v_R00214*x(316); f(847)=v_R00214*x(317); f(848)=v_R00214*x(318); f(849)=v_R01529_f*x(319); f(850)=v_R01056_b*x(319); f(851)=v_R01827_b*x(320)*x(10); f(852)=v_R01827_b*x(320)*x(34); f(853)=v_R01827_b*x(320)*x(200); f(854)=v_R01827_b*x(320)*x(211); f(855)=v_R01827_b*x(320)*x(212); f(856)=v_R01827_b*x(320)*x(213); f(857)=v_R01827_b*x(320)*x(312); f(858)=v_R01067_b*x(320)*x(11); f(859)=v_R01067_b*x(320)*x(43); f(860)=v_R01067_b*x(320)*x(199); f(861)=v_R01067_b*x(320)*x(208); f(862)=v_R01067_b*x(320)*x(209); f(863)=v_R01067_b*x(320)*x(210); f(864)=v_R01067_b*x(320)*x(311); f(865)=v_vs_Ru5P*x(319); f(866)=v_vs_E4P*x(320); f(867)=v_R00341*x(322); f(868)=v_R00210*x(323)*x(17); f(869)=v_R00351*x(6)*x(322); f(870)=v_R00351*x(51)*x(322); f(871)=v_R00351*x(100)*x(322); f(872)=v_R00351*x(238)*x(322); f(873)=v_R00342*x(321); f(874)=v_R00214*x(321); f(875)=v_R01056_f*x(324); f(876)=v_R01641_f*x(325)*x(28); f(877)=v_R01641_f*x(325)*x(38); f(878)=v_R01641_f*x(325)*x(173); f(879)=v_R01641_f*x(325)*x(192); f(880)=v_R01641_f*x(325)*x(193); f(881)=v_R01641_f*x(325)*x(194); f(882)=v_R01641_f*x(325)*x(296); f(883)=v_R01641_b*x(324)*x(11); f(884)=v_R01641_b*x(324)*x(43); f(885)=v_R01641_b*x(324)*x(199); f(886)=v_R01641_b*x(324)*x(208); f(887)=v_R01641_b*x(324)*x(209); f(888)=v_R01641_b*x(324)*x(210); f(889)=v_R01641_b*x(324)*x(311); f(890)=v_R01827_f*x(325)*x(28); f(891)=v_R01827_f*x(325)*x(38); f(892)=v_R01827_f*x(325)*x(173); f(893)=v_R01827_f*x(325)*x(192); f(894)=v_R01827_f*x(325)*x(193); f(895)=v_R01827_f*x(325)*x(194); f(896)=v_R01827_f*x(325)*x(296); f(897)=v_vs_PYR*x(323); f(898)=v_vs_OAA*x(322); f(899)=v_R00658_b*x(326); f(900)=v_R00200*x(326); f(901)=v_R00345*x(326)*x(7); f(902)=v_R00345*x(326)*x(80); f(903)=v_R00341*x(331); f(904)=v_R00351*x(6)*x(331); f(905)=v_R00351*x(51)*x(331); f(906)=v_R00351*x(100)*x(331); f(907)=v_R00351*x(238)*x(331); f(908)=v_R01324_f*x(327); f(909)=v_R01324_f*x(328); f(910)=v_R01324_f*x(329); f(911)=v_R01324_f*x(330); f(912)=v_vs_PEP*x(326); f(913)=v_vs_OAA*x(331); f(914)=v_R01518_f*x(332); f(915)=v_R00658_f*x(332); f(916)=v_R01324_f*x(333); f(917)=v_R01324_f*x(334); f(918)=v_R01324_f*x(335); f(919)=v_R01324_f*x(336); f(920)=v_R01324_b*x(337); f(921)=v_R01324_b*x(338); f(922)=v_R01324_b*x(339); f(923)=v_R01324_b*x(340); f(924)=v_R00267*x(337); f(925)=v_R00267*x(338); f(926)=v_R00267*x(339); f(927)=v_R00267*x(340); f(928)=v_R01512_f*x(341); f(929)=v_R01518_b*x(341); f(930)=v_R01324_b*x(342); f(931)=v_R01324_b*x(343); f(932)=v_R01324_b*x(344); f(933)=v_R01324_b*x(345); f(934)=v_R00267*x(342); f(935)=v_R00267*x(343); f(936)=v_R00267*x(344); f(937)=v_R00267*x(345); f(938)=v_vs_G3P*x(341); f(939)=v_R01061_b*x(346); f(940)=v_R01512_b*x(346); f(941)=v_R01015_f*x(347); f(942)=v_R01061_f*x(347); f(943)=v_R01641_f*x(4)*x(347); f(944)=v_R01641_f*x(44)*x(347); f(945)=v_R01641_f*x(224)*x(347); f(946)=v_R01641_f*x(234)*x(347); f(947)=v_R01641_f*x(235)*x(347); f(948)=v_R01641_f*x(236)*x(347); f(949)=v_R01641_f*x(325)*x(347); f(950)=v_R01827_f*x(4)*x(347); f(951)=v_R01827_f*x(44)*x(347); f(952)=v_R01827_f*x(224)*x(347); f(953)=v_R01827_f*x(234)*x(347); f(954)=v_R01827_f*x(235)*x(347); f(955)=v_R01827_f*x(236)*x(347); f(956)=v_R01827_f*x(325)*x(347); f(957)=v_R01067_f*x(10)*x(347); f(958)=v_R01067_f*x(34)*x(347); f(959)=v_R01067_f*x(200)*x(347); f(960)=v_R01067_f*x(211)*x(347); f(961)=v_R01067_f*x(212)*x(347); f(962)=v_R01067_f*x(213)*x(347); f(963)=v_R01067_f*x(312)*x(347); f(964)=v_vs_T3P*x(347); f(965)=v_R00756*x(350); f(966)=v_R01015_b*x(348); f(967)=v_R01529_b*x(349); f(968)=v_R01641_b*x(14)*x(349); f(969)=v_R01641_b*x(46)*x(349); f(970)=v_R01641_b*x(223)*x(349); f(971)=v_R01641_b*x(231)*x(349); f(972)=v_R01641_b*x(232)*x(349); f(973)=v_R01641_b*x(233)*x(349); f(974)=v_R01641_b*x(324)*x(349); f(975)=v_R01827_b*x(18)*x(350); f(976)=v_R01827_b*x(39)*x(350); f(977)=v_R01827_b*x(214)*x(350); f(978)=v_R01827_b*x(225)*x(350); f(979)=v_R01827_b*x(226)*x(350); f(980)=v_R01827_b*x(227)*x(350); f(981)=v_R01827_b*x(320)*x(350); f(982)=v_R01067_f*x(350)*x(28); f(983)=v_R01067_f*x(350)*x(38); f(984)=v_R01067_f*x(350)*x(173); f(985)=v_R01067_f*x(350)*x(192); f(986)=v_R01067_f*x(350)*x(193); f(987)=v_R01067_f*x(350)*x(194); f(988)=v_R01067_f*x(350)*x(296); f(989)=v_R01067_f*x(350)*x(347); f(990)=v_R01067_b*x(18)*x(349); f(991)=v_R01067_b*x(39)*x(349); f(992)=v_R01067_b*x(214)*x(349); f(993)=v_R01067_b*x(225)*x(349); f(994)=v_R01067_b*x(226)*x(349); f(995)=v_R01067_b*x(227)*x(349); f(996)=v_R01067_b*x(320)*x(349); f(997)=v_vs_F6P*x(350); f(998)=v_R01068*x(351); f(999)=v_R01529_f*x(352); f(1000)=v_R01056_b*x(352); f(1001)=v_R01827_b*x(353)*x(10); f(1002)=v_R01827_b*x(353)*x(34); f(1003)=v_R01827_b*x(353)*x(200); f(1004)=v_R01827_b*x(353)*x(211); f(1005)=v_R01827_b*x(353)*x(212); f(1006)=v_R01827_b*x(353)*x(213); f(1007)=v_R01827_b*x(353)*x(312); f(1008)=v_R01827_b*x(353)*x(350); f(1009)=v_R01067_b*x(353)*x(11); f(1010)=v_R01067_b*x(353)*x(43); f(1011)=v_R01067_b*x(353)*x(199); f(1012)=v_R01067_b*x(353)*x(208); f(1013)=v_R01067_b*x(353)*x(209); f(1014)=v_R01067_b*x(353)*x(210); f(1015)=v_R01067_b*x(353)*x(311); f(1016)=v_R01067_b*x(353)*x(349); f(1017)=v_vs_Ru5P*x(352); f(1018)=v_vs_E4P*x(353); f(1019)=v_R01056_f*x(354); f(1020)=v_R01641_f*x(355)*x(28); f(1021)=v_R01641_f*x(355)*x(38); f(1022)=v_R01641_f*x(355)*x(173); f(1023)=v_R01641_f*x(355)*x(192); f(1024)=v_R01641_f*x(355)*x(193); f(1025)=v_R01641_f*x(355)*x(194); f(1026)=v_R01641_f*x(355)*x(296); f(1027)=v_R01641_f*x(355)*x(347); f(1028)=v_R01641_b*x(354)*x(11); f(1029)=v_R01641_b*x(354)*x(43); f(1030)=v_R01641_b*x(354)*x(199); f(1031)=v_R01641_b*x(354)*x(208); f(1032)=v_R01641_b*x(354)*x(209); f(1033)=v_R01641_b*x(354)*x(210); f(1034)=v_R01641_b*x(354)*x(311); f(1035)=v_R01641_b*x(354)*x(349); f(1036)=v_R01827_f*x(355)*x(28); f(1037)=v_R01827_f*x(355)*x(38); f(1038)=v_R01827_f*x(355)*x(173); f(1039)=v_R01827_f*x(355)*x(192); f(1040)=v_R01827_f*x(355)*x(193); f(1041)=v_R01827_f*x(355)*x(194); f(1042)=v_R01827_f*x(355)*x(296); f(1043)=v_R01827_f*x(355)*x(347); end % Derivative function function d= xdot(t,x) f=flux(t,x); d(1,1)= +f(29) -f(30); d(2,1)= +f(25) -f(26); d(3,1)= +f(30) -f(31); d(4,1)= -f(36) +f(37) -f(38) +f(39) +f(61) -f(68) -f(69) +f(85) -f(389) -f(391) -f(415) -f(416) -f(417) +f(421) -f(423) -f(424) -f(425) +f(429) +f(463) +f(464) +f(465) +f(469) +f(470) +f(471) -f(782) -f(788) +f(813) +f(819) -f(943) -f(950) +f(968) +f(975); d(5,1)= -f(9) +f(10) +f(11) -f(12) -f(45); d(6,1)= +f(18) -f(19) -f(48) -f(105) +f(148) -f(149) -f(150) -f(180) -f(181) -f(182) -f(183) -f(184) -f(218) -f(263) -f(264) -f(591) -f(680) -f(869) -f(904); d(7,1)= -f(16) +f(17) +f(18) +f(22) +f(23) +f(28) +f(31) +f(67) -f(101) +f(103) +f(104) +f(115) +f(118) +f(119) +f(120) +f(122) +f(123) +f(143) +f(147) +f(159) +f(160) +f(162) +f(164) -f(170) +f(173) +f(174) +f(176) +f(177) -f(211) -f(213) -f(215) +f(244) +f(246) +f(291) +f(292) +f(294) +f(296) +f(297) +f(299) +f(336) +f(337) +f(338) +f(340) +f(341) +f(343) +f(347) +f(350) +f(381) +f(382) +f(385) +f(521) +f(589) +f(590) -f(657) +f(700) +f(701) +f(702) +f(723) +f(724) +f(725) +f(727) +f(728) +f(730) +f(735) +f(736) +f(737) +f(738) +f(749) +f(750) +f(751) +f(752) +f(753) +f(754) +f(757) +f(759) +f(846) +f(867) -f(901); d(8,1)= +f(26) -f(27) -f(28); d(9,1)= +f(16) -f(17) -f(19) +f(27) -f(51) -f(108) -f(220) -f(659); d(10,1)= +f(2) -f(3) +f(38) -f(39) -f(40) +f(41) -f(43) -f(70) -f(72) +f(86) +f(88) -f(393) -f(431) -f(432) -f(433) +f(443) -f(475) +f(499) +f(500) +f(501) +f(541) -f(547) -f(553) -f(559) +f(635) +f(641) +f(647) -f(794) +f(832) -f(851) +f(890) -f(957) +f(990) -f(1001) +f(1036); d(11,1)= +f(32) -f(33) +f(36) -f(37) +f(40) -f(41) +f(62) -f(74) +f(83) -f(94) +f(437) +f(481) +f(487) +f(493) -f(505) +f(529) -f(535) -f(565) -f(571) -f(577) +f(599) +f(605) +f(611) -f(617) -f(623) -f(629) +f(825) -f(858) +f(876) -f(883) +f(982) -f(1009) +f(1020) -f(1028); d(12,1)= -f(11) +f(12) -f(13) +f(14); d(13,1)= +f(4) +f(5) -f(6) +f(64) +f(452) +f(516) +f(517) +f(518) +f(839) +f(998); d(14,1)= -f(34) +f(35) +f(36) -f(37) +f(68) -f(85) +f(389) +f(415) +f(416) +f(417) -f(421) -f(463) -f(464) -f(465) +f(782) -f(813) +f(943) -f(968); d(15,1)= +f(23) -f(24) +f(301); d(16,1)= +f(3) -f(4); d(17,1)= -f(18) +f(19) -f(23) +f(24) -f(104) +f(105) +f(108) +f(109) -f(120) -f(122) -f(123) +f(124) +f(127) +f(128) -f(148) +f(149) +f(150) +f(151) +f(152) -f(177) -f(178) -f(179) +f(180) +f(181) +f(182) +f(183) +f(184) +f(185) +f(186) +f(187) +f(188) +f(189) +f(218) +f(219) +f(220) +f(221) +f(222) +f(223) +f(224) +f(225) +f(226) +f(227) +f(228) +f(229) +f(263) +f(264) +f(265) +f(266) +f(267) +f(268) -f(301) -f(302) -f(346) -f(347) -f(348) -f(349) -f(350) -f(351) -f(381) -f(382) -f(383) -f(384) -f(385) -f(386) +f(387) +f(388) +f(410) +f(411) +f(412) -f(590) +f(591) +f(592) +f(593) +f(659) +f(660) +f(661) +f(662) +f(663) +f(664) +f(665) +f(666) +f(667) +f(668) +f(669) +f(670) +f(671) +f(680) +f(681) +f(682) +f(683) -f(736) -f(737) -f(738) -f(753) -f(754) -f(755) -f(756) -f(757) -f(758) -f(759) +f(760) +f(761) +f(762) -f(772) -f(773) -f(774) -f(775) +f(776) +f(777) +f(778) +f(779) -f(868) +f(869) +f(870) +f(871) +f(872) +f(904) +f(905) +f(906) +f(907); d(18,1)= +f(38) -f(39) +f(40) -f(41) -f(50) -f(61) +f(69) +f(72) -f(88) +f(391) +f(393) +f(423) +f(424) +f(425) -f(429) +f(431) +f(432) +f(433) -f(443) -f(469) -f(470) -f(471) -f(499) -f(500) -f(501) +f(788) +f(794) -f(819) -f(832) +f(950) +f(957) -f(975) -f(990); d(19,1)= +f(7) -f(8) +f(9) -f(10); d(20,1)= +f(15) -f(18) +f(28) -f(47) +f(142); d(21,1)= +f(13) -f(14) -f(15) -f(16) +f(17) -f(46) -f(144) +f(146); d(22,1)= +f(19) -f(20) +f(21); d(23,1)= +f(24) -f(25); d(24,1)= +f(31) -f(32) +f(33) +f(34) -f(35) -f(49); d(25,1)= +f(22) -f(23) -f(52) +f(245); d(26,1)= +f(1) -f(2) -f(29) -f(42); d(27,1)= -f(1) +f(53); d(28,1)= +f(4) -f(5) +f(6) -f(7) +f(8) -f(36) +f(37) -f(38) +f(39) -f(40) +f(41) -f(44) -f(62) +f(70) +f(74) -f(83) -f(86) +f(94) -f(437) +f(475) -f(481) -f(487) -f(493) +f(505) -f(529) +f(535) -f(541) +f(547) +f(553) +f(559) +f(565) +f(571) +f(577) -f(599) -f(605) -f(611) +f(617) +f(623) +f(629) -f(635) -f(641) -f(647) -f(825) +f(851) +f(858) -f(876) +f(883) -f(890) -f(982) +f(1001) +f(1009) -f(1020) +f(1028) -f(1036); d(29,1)= +f(20) -f(21) -f(22); d(30,1)= 0*f(53) 0*f(54); d(31,1)= +f(42) +f(43) +f(44) +f(45) +f(46) +f(47) +f(48) +f(49) +f(50) +f(51) +f(52) +f(58) +f(63) +f(75) +f(76) +f(90) +f(96) +f(102) +f(106) +f(107) +f(111) +f(121) +f(125) +f(126) +f(165) +f(166) +f(167) +f(194) +f(195) +f(196) +f(197) +f(198) +f(199) +f(200) +f(201) +f(202) +f(248) +f(249) +f(250) +f(251) +f(252) +f(303) +f(304) +f(305) +f(306) +f(307) +f(352) +f(353) +f(354) +f(355) +f(356) +f(357) +f(358) +f(359) +f(360) +f(395) +f(396) +f(397) +f(398) +f(399) +f(400) +f(401) +f(445) +f(446) +f(447) +f(448) +f(511) +f(512) +f(513) +f(514) +f(515) +f(583) +f(584) +f(585) +f(586) +f(587) +f(588) +f(653) +f(654) +f(675) +f(676) +f(703) +f(739) +f(740) +f(741) +f(742) +f(763) +f(764) +f(765) +f(766) +f(767) +f(768) +f(769) +f(800) +f(801) +f(802) +f(803) +f(804) +f(838) +f(865) +f(866) +f(897) +f(898) +f(912) +f(913) +f(938) +f(964) +f(997) +f(1017) +f(1018); d(32,1)= +f(54) -f(55); d(33,1)= +f(55) -f(56) -f(57) -f(58); d(34,1)= +f(56) -f(59) -f(61) -f(62) -f(63) +f(69) -f(71) -f(73) +f(74) +f(87) +f(89) -f(394) -f(434) -f(435) -f(436) +f(444) -f(476) +f(502) +f(503) +f(504) +f(542) -f(548) -f(554) -f(560) +f(636) +f(642) +f(648) -f(795) +f(833) -f(852) +f(891) -f(958) +f(991) -f(1002) +f(1037); d(35,1)= +f(57) -f(60); d(36,1)= +f(59) -f(64); d(37,1)= +f(60) -f(67); d(38,1)= +f(61) +f(64) -f(65) -f(66) -f(68) -f(69) +f(71) -f(72) -f(73) -f(75) +f(77) +f(78) -f(84) +f(85) -f(87) +f(88) +f(89) +f(95) -f(438) +f(476) -f(482) -f(488) -f(494) +f(506) -f(530) +f(536) -f(542) +f(548) +f(554) +f(560) +f(566) +f(572) +f(578) -f(600) -f(606) -f(612) +f(618) +f(624) +f(630) -f(636) -f(642) -f(648) -f(826) +f(852) +f(859) -f(877) +f(884) -f(891) -f(983) +f(1002) +f(1010) -f(1021) +f(1029) -f(1037); d(39,1)= +f(62) -f(70) -f(71) +f(73) -f(74) -f(76) +f(86) +f(87) -f(89) +f(392) +f(394) +f(426) +f(427) +f(428) -f(430) +f(434) +f(435) +f(436) -f(444) -f(472) -f(473) -f(474) -f(502) -f(503) -f(504) +f(789) +f(795) -f(820) -f(833) +f(951) +f(958) -f(976) -f(991); d(40,1)= +f(65) -f(77); d(41,1)= +f(66) -f(78) -f(79) +f(91); d(42,1)= +f(67) -f(80) +f(81) -f(82) -f(90) +f(93); d(43,1)= +f(68) +f(72) +f(73) +f(80) -f(81) +f(84) -f(85) -f(88) -f(89) -f(95) +f(438) +f(482) +f(488) +f(494) -f(506) +f(530) -f(536) -f(566) -f(572) -f(578) +f(600) +f(606) +f(612) -f(618) -f(624) -f(630) +f(826) -f(859) +f(877) -f(884) +f(983) -f(1010) +f(1021) -f(1029); d(44,1)= +f(70) +f(71) -f(83) -f(84) -f(86) -f(87) +f(94) +f(95) -f(390) -f(392) -f(418) -f(419) -f(420) +f(422) -f(426) -f(427) -f(428) +f(430) +f(466) +f(467) +f(468) +f(472) +f(473) +f(474) -f(783) -f(789) +f(814) +f(820) -f(944) -f(951) +f(969) +f(976); d(45,1)= +f(79) -f(91) -f(92) -f(96) +f(97); d(46,1)= +f(82) +f(83) +f(84) -f(93) -f(94) -f(95) +f(390) +f(418) +f(419) +f(420) -f(422) -f(466) -f(467) -f(468) +f(783) -f(814) +f(944) -f(969); d(47,1)= +f(92) -f(97) -f(98) +f(99); d(48,1)= +f(98) -f(99) -f(100) -f(101) -f(102) +f(103) -f(145) +f(172); d(49,1)= +f(100) -f(104) -f(106) +f(159) +f(161); d(50,1)= +f(101) -f(103) -f(105) -f(107) -f(109) +f(153) -f(221) -f(660); d(51,1)= +f(104) -f(108) -f(109) -f(111) -f(151) -f(152) +f(178) -f(185) -f(186) -f(187) -f(188) -f(189) -f(219) -f(265) -f(266) -f(592) -f(681) -f(870) -f(905); d(52,1)= +f(105) -f(110) +f(114); d(53,1)= +f(108) -f(112) +f(116); d(54,1)= +f(109) -f(113) +f(117); d(55,1)= +f(110) -f(114) -f(115); d(56,1)= +f(112) -f(116) -f(118); d(57,1)= +f(113) -f(117) -f(119); d(58,1)= +f(115) -f(120) -f(121) +f(293); d(59,1)= +f(118) -f(122) -f(125) +f(247); d(60,1)= +f(119) -f(123) -f(126) +f(298); d(61,1)= +f(120) -f(124) +f(346); d(62,1)= +f(122) -f(127) +f(302); d(63,1)= +f(123) -f(128) +f(349); d(64,1)= +f(124) -f(129) +f(410); d(65,1)= +f(127) -f(130) +f(387); d(66,1)= +f(128) -f(131) +f(412) +f(760) +f(776); d(67,1)= +f(129) -f(132) -f(133); d(68,1)= +f(130) -f(134) -f(135); d(69,1)= +f(131) -f(136) -f(137) -f(138) -f(139); d(70,1)= +f(132) -f(140) -f(142); d(71,1)= +f(133) -f(141) -f(143); d(72,1)= +f(134) -f(153) -f(159); d(73,1)= +f(135) -f(154) -f(160); d(74,1)= +f(136) -f(155) -f(161); d(75,1)= +f(137) -f(156) -f(162); d(76,1)= +f(138) -f(157) -f(163); d(77,1)= +f(139) -f(158) -f(164); d(78,1)= +f(140) +f(144) -f(146) -f(149) -f(151) -f(166) -f(222) -f(661); d(79,1)= +f(141) -f(147) -f(150) -f(152) -f(167) +f(170) -f(223) -f(662); d(80,1)= +f(142) -f(144) -f(145) +f(146) +f(148) +f(161) +f(163) -f(171) +f(172) +f(175) +f(178) +f(179) -f(212) -f(214) -f(216) +f(217) +f(245) +f(247) +f(261) +f(262) +f(293) +f(295) +f(298) +f(300) +f(301) +f(302) +f(334) +f(335) +f(339) +f(342) +f(344) +f(345) +f(346) +f(348) +f(349) +f(351) +f(375) +f(376) +f(377) +f(378) +f(379) +f(380) +f(383) +f(384) +f(386) +f(595) -f(658) +f(679) +f(726) +f(729) +f(731) +f(732) +f(733) +f(734) +f(755) +f(756) +f(758) +f(772) +f(773) +f(774) +f(775) +f(845) +f(847) +f(848) +f(868) +f(874) -f(902) +f(903) +f(924) +f(925) +f(926) +f(927) +f(934) +f(935) +f(936) +f(937); d(81,1)= +f(143) -f(148) -f(165) +f(169) +f(595); d(82,1)= +f(145) +f(155) -f(172) -f(180) -f(185) -f(198) -f(224) -f(663); d(83,1)= +f(147) -f(168) -f(169) -f(170) -f(171) -f(194) +f(204) +f(217); d(84,1)= +f(149) -f(190) +f(240); d(85,1)= +f(150) -f(191) +f(241); d(86,1)= +f(151) -f(192) +f(242); d(87,1)= +f(152) -f(193) +f(243); d(88,1)= +f(154) -f(173) -f(181) -f(186) -f(199) +f(211) -f(225) -f(664); d(89,1)= +f(156) -f(174) -f(182) -f(187) -f(200) +f(213) -f(226) -f(665); d(90,1)= +f(157) -f(175) -f(183) -f(188) -f(201) +f(212) -f(227) -f(666); d(91,1)= +f(158) -f(176) -f(184) -f(189) -f(202) +f(215) -f(228) -f(667); d(92,1)= +f(160) +f(163) -f(177) -f(195) +f(208); d(93,1)= +f(162) -f(178) -f(196) +f(209) +f(847); d(94,1)= +f(164) -f(179) -f(197) +f(210) +f(848); d(95,1)= +f(168) -f(203) -f(204) +f(257); d(96,1)= +f(171) -f(217) -f(218) -f(219) -f(229) -f(252) +f(594) -f(668); d(97,1)= +f(173) +f(175) -f(205) -f(208) -f(211) -f(212) -f(248) +f(258); d(98,1)= +f(174) -f(206) -f(209) -f(213) -f(214) -f(249) +f(259) +f(261); d(99,1)= +f(176) -f(207) -f(210) -f(215) -f(216) -f(250) +f(260) +f(262); d(100,1)= +f(177) +f(179) -f(220) -f(221) -f(222) -f(223) -f(224) -f(225) -f(226) -f(227) -f(228) -f(229) -f(251) -f(267) -f(268) -f(593) -f(682) -f(871) -f(906); d(101,1)= +f(180) -f(230) +f(281); d(102,1)= +f(181) -f(231) +f(282); d(103,1)= +f(182) -f(232) +f(283); d(104,1)= +f(183) -f(233) +f(284); d(105,1)= +f(184) -f(234) +f(285); d(106,1)= +f(185) -f(235) +f(286); d(107,1)= +f(186) -f(236) +f(287); d(108,1)= +f(187) -f(237) +f(288); d(109,1)= +f(188) -f(238) +f(289); d(110,1)= +f(189) -f(239) +f(290); d(111,1)= +f(190) -f(240) -f(244); d(112,1)= +f(191) -f(241) -f(245); d(113,1)= +f(192) -f(242) -f(246); d(114,1)= +f(193) -f(243) -f(247); d(115,1)= +f(203) -f(253) -f(257) -f(303) +f(312); d(116,1)= +f(205) -f(254) -f(258) +f(313); d(117,1)= +f(206) -f(255) -f(259) +f(314); d(118,1)= +f(207) -f(256) -f(260) +f(315); d(119,1)= +f(214) -f(261) -f(263) -f(265) -f(267) -f(304) -f(669) +f(843); d(120,1)= +f(216) -f(262) -f(264) -f(266) -f(268) -f(305) -f(670) +f(844); d(121,1)= +f(218) -f(269) +f(322); d(122,1)= +f(219) -f(270) +f(323); d(123,1)= +f(220) -f(271) +f(324); d(124,1)= +f(221) -f(272) +f(325); d(125,1)= +f(222) -f(273) +f(326); d(126,1)= +f(223) -f(274) +f(327); d(127,1)= +f(224) -f(275) +f(328); d(128,1)= +f(225) -f(276) +f(329); d(129,1)= +f(226) -f(277) +f(330); d(130,1)= +f(227) -f(278) +f(331); d(131,1)= +f(228) -f(279) +f(332); d(132,1)= +f(229) -f(280) +f(333); d(133,1)= +f(230) -f(281) -f(291); d(134,1)= +f(231) -f(282) -f(292); d(135,1)= +f(232) -f(283) -f(293); d(136,1)= +f(233) -f(284) -f(294); d(137,1)= +f(234) -f(285) -f(295); d(138,1)= +f(235) -f(286) -f(296); d(139,1)= +f(236) -f(287) -f(297); d(140,1)= +f(237) -f(288) -f(298); d(141,1)= +f(238) -f(289) -f(299); d(142,1)= +f(239) -f(290) -f(300); d(143,1)= +f(244) -f(301) -f(306) +f(334); d(144,1)= +f(246) -f(302) -f(307) +f(335); d(145,1)= +f(253) -f(308) -f(312) +f(362); d(146,1)= +f(254) -f(309) -f(313) -f(352) +f(366); d(147,1)= +f(255) -f(310) -f(314) -f(353) +f(367); d(148,1)= +f(256) -f(311) -f(315) -f(354) +f(368); d(149,1)= +f(263) -f(316) +f(369); d(150,1)= +f(264) -f(317) +f(370); d(151,1)= +f(265) -f(318) +f(371); d(152,1)= +f(266) -f(319) +f(372); d(153,1)= +f(267) -f(320) +f(373); d(154,1)= +f(268) -f(321) +f(374); d(155,1)= +f(269) -f(322) -f(334); d(156,1)= +f(270) -f(323) -f(335); d(157,1)= +f(271) -f(324) -f(336); d(158,1)= +f(272) -f(325) -f(337); d(159,1)= +f(273) -f(326) -f(338); d(160,1)= +f(274) -f(327) -f(339); d(161,1)= +f(275) -f(328) -f(340); d(162,1)= +f(276) -f(329) -f(341); d(163,1)= +f(277) -f(330) -f(342); d(164,1)= +f(278) -f(331) -f(343); d(165,1)= +f(279) -f(332) -f(344); d(166,1)= +f(280) -f(333) -f(345); d(167,1)= +f(291) -f(346) -f(355) +f(375); d(168,1)= +f(292) +f(295) -f(347) -f(356); d(169,1)= +f(294) -f(348) -f(357) +f(376); d(170,1)= +f(296) -f(349) -f(358) +f(377); d(171,1)= +f(297) +f(300) -f(350) -f(359); d(172,1)= +f(299) -f(351) -f(360) +f(378); d(173,1)= +f(308) -f(361) -f(362) -f(389) -f(390) -f(391) -f(392) -f(393) -f(394) -f(395) +f(406) +f(421) +f(422) +f(429) +f(430) -f(439) +f(443) +f(444) +f(452) +f(477) -f(483) -f(489) -f(495) +f(507) -f(531) +f(537) -f(543) +f(549) +f(555) +f(561) +f(567) +f(573) +f(579) -f(601) -f(607) -f(613) +f(619) +f(625) +f(631) -f(637) -f(643) -f(649) -f(827) +f(853) +f(860) -f(878) +f(885) -f(892) -f(984) +f(1003) +f(1011) -f(1022) +f(1030) -f(1038); d(174,1)= +f(309) -f(363) -f(366) +f(407); d(175,1)= +f(310) -f(364) -f(367) +f(408); d(176,1)= +f(311) -f(365) -f(368) +f(409); d(177,1)= +f(316) -f(369) -f(375); d(178,1)= +f(317) -f(370) -f(376); d(179,1)= +f(318) -f(371) -f(377); d(180,1)= +f(319) -f(372) -f(378); d(181,1)= +f(320) -f(373) -f(379); d(182,1)= +f(321) -f(374) -f(380); d(183,1)= +f(336) +f(339) -f(381) -f(396); d(184,1)= +f(337) +f(342) -f(382) -f(397); d(185,1)= +f(338) +f(345) -f(383) -f(398); d(186,1)= +f(340) +f(379) -f(384) -f(399); d(187,1)= +f(341) +f(344) -f(385) -f(400); d(188,1)= +f(343) +f(380) -f(386) -f(401); d(189,1)= +f(347) +f(348) -f(387); d(190,1)= +f(350) +f(351) -f(388); d(191,1)= +f(361) -f(406); d(192,1)= +f(363) -f(403) -f(407) -f(415) -f(418) -f(423) -f(426) -f(431) -f(434) -f(440) -f(446) +f(453) +f(463) +f(466) +f(469) +f(472) +f(478) -f(484) -f(490) -f(496) +f(499) +f(502) +f(508) +f(516) -f(532) +f(538) -f(544) +f(550) +f(556) +f(562) +f(568) +f(574) +f(580) -f(602) -f(608) -f(614) +f(620) +f(626) +f(632) -f(638) -f(644) -f(650) -f(828) +f(854) +f(861) -f(879) +f(886) -f(893) -f(985) +f(1004) +f(1012) -f(1023) +f(1031) -f(1039); d(193,1)= +f(364) -f(404) -f(408) -f(416) -f(419) -f(424) -f(427) -f(432) -f(435) -f(441) -f(447) +f(454) +f(464) +f(467) +f(470) +f(473) +f(479) -f(485) -f(491) -f(497) +f(500) +f(503) +f(509) +f(517) -f(533) +f(539) -f(545) +f(551) +f(557) +f(563) +f(569) +f(575) +f(581) -f(603) -f(609) -f(615) +f(621) +f(627) +f(633) -f(639) -f(645) -f(651) -f(829) +f(855) +f(862) -f(880) +f(887) -f(894) -f(986) +f(1005) +f(1013) -f(1024) +f(1032) -f(1040); d(194,1)= +f(365) -f(405) -f(409) -f(417) -f(420) -f(425) -f(428) -f(433) -f(436) -f(442) -f(448) +f(455) +f(465) +f(468) +f(471) +f(474) +f(480) -f(486) -f(492) -f(498) +f(501) +f(504) +f(510) +f(518) -f(534) +f(540) -f(546) +f(552) +f(558) +f(564) +f(570) +f(576) +f(582) -f(604) -f(610) -f(616) +f(622) +f(628) +f(634) -f(640) -f(646) -f(652) -f(830) +f(856) +f(863) -f(881) +f(888) -f(895) -f(987) +f(1006) +f(1014) -f(1025) +f(1033) -f(1041); d(195,1)= +f(381) +f(383) -f(410); d(196,1)= +f(382) +f(384) -f(411); d(197,1)= +f(385) +f(386) -f(412); d(198,1)= +f(388) -f(413); d(199,1)= +f(389) +f(390) +f(393) +f(394) -f(414) -f(421) -f(422) +f(439) -f(443) -f(444) +f(458) +f(483) +f(489) +f(495) -f(507) +f(531) -f(537) -f(567) -f(573) -f(579) +f(601) +f(607) +f(613) -f(619) -f(625) -f(631) +f(827) -f(860) +f(878) -f(885) +f(984) -f(1011) +f(1022) -f(1030); d(200,1)= +f(391) +f(392) -f(402) -f(429) -f(430) -f(437) -f(438) -f(439) -f(440) -f(441) -f(442) -f(445) -f(477) +f(505) +f(506) +f(507) +f(508) +f(509) +f(510) +f(543) -f(549) -f(555) -f(561) +f(637) +f(643) +f(649) -f(796) +f(834) -f(853) +f(892) -f(959) +f(992) -f(1003) +f(1038); d(201,1)= +f(402) -f(452); d(202,1)= +f(403) -f(453); d(203,1)= +f(404) -f(454); d(204,1)= +f(405) -f(455); d(205,1)= +f(411) -f(456); d(206,1)= +f(413) -f(457); d(207,1)= +f(414) -f(458) -f(462) -f(514) +f(525); d(208,1)= +f(415) +f(418) +f(431) +f(434) +f(440) -f(459) -f(463) -f(466) +f(484) +f(490) +f(496) -f(499) -f(502) -f(508) +f(522) +f(532) -f(538) -f(568) -f(574) -f(580) +f(602) +f(608) +f(614) -f(620) -f(626) -f(632) +f(828) -f(861) +f(879) -f(886) +f(985) -f(1012) +f(1023) -f(1031); d(209,1)= +f(416) +f(419) +f(432) +f(435) +f(441) -f(460) -f(464) -f(467) +f(485) +f(491) +f(497) -f(500) -f(503) -f(509) +f(523) +f(533) -f(539) -f(569) -f(575) -f(581) +f(603) +f(609) +f(615) -f(621) -f(627) -f(633) +f(829) -f(862) +f(880) -f(887) +f(986) -f(1013) +f(1024) -f(1032); d(210,1)= +f(417) +f(420) +f(433) +f(436) +f(442) -f(461) -f(465) -f(468) +f(486) +f(492) +f(498) -f(501) -f(504) -f(510) +f(524) +f(534) -f(540) -f(570) -f(576) -f(582) +f(604) +f(610) +f(616) -f(622) -f(628) -f(634) +f(830) -f(863) +f(881) -f(888) +f(987) -f(1014) +f(1025) -f(1033); d(211,1)= +f(423) +f(426) -f(449) -f(469) -f(472) -f(478) -f(481) -f(482) -f(483) -f(484) -f(485) -f(486) -f(511) +f(544) -f(550) -f(556) -f(562) +f(565) +f(566) +f(567) +f(568) +f(569) +f(570) +f(638) +f(644) +f(650) -f(797) +f(835) -f(854) +f(893) -f(960) +f(993) -f(1004) +f(1039); d(212,1)= +f(424) +f(427) -f(450) -f(470) -f(473) -f(479) -f(487) -f(488) -f(489) -f(490) -f(491) -f(492) -f(512) +f(545) -f(551) -f(557) -f(563) +f(571) +f(572) +f(573) +f(574) +f(575) +f(576) +f(639) +f(645) +f(651) -f(798) +f(836) -f(855) +f(894) -f(961) +f(994) -f(1005) +f(1040); d(213,1)= +f(425) +f(428) -f(451) -f(471) -f(474) -f(480) -f(493) -f(494) -f(495) -f(496) -f(497) -f(498) -f(513) +f(546) -f(552) -f(558) -f(564) +f(577) +f(578) +f(579) +f(580) +f(581) +f(582) +f(640) +f(646) +f(652) -f(799) +f(837) -f(856) +f(895) -f(962) +f(995) -f(1006) +f(1041); d(214,1)= +f(437) +f(438) +f(439) +f(440) +f(441) +f(442) -f(475) -f(476) -f(477) -f(478) -f(479) -f(480) -f(505) -f(506) -f(507) -f(508) -f(509) -f(510) -f(515) +f(541) +f(542) +f(543) +f(544) +f(545) +f(546) +f(790) +f(796) -f(821) -f(834) +f(952) +f(959) -f(977) -f(992); d(215,1)= +f(449) -f(516); d(216,1)= +f(450) -f(517); d(217,1)= +f(451) -f(518); d(218,1)= +f(456) -f(519); d(219,1)= +f(457) -f(520) -f(521); d(220,1)= +f(459) -f(522) -f(526) -f(583) +f(596); d(221,1)= +f(460) -f(523) -f(527) -f(584) +f(597); d(222,1)= +f(461) -f(524) -f(528) -f(585) +f(598); d(223,1)= +f(462) -f(525) +f(529) +f(530) +f(531) +f(532) +f(533) +f(534) -f(535) -f(536) -f(537) -f(538) -f(539) -f(540) +f(784) -f(815) +f(945) -f(970); d(224,1)= +f(475) +f(476) +f(477) +f(478) +f(479) +f(480) -f(529) -f(530) -f(531) -f(532) -f(533) -f(534) +f(535) +f(536) +f(537) +f(538) +f(539) +f(540) -f(541) -f(542) -f(543) -f(544) -f(545) -f(546) -f(784) -f(790) +f(815) +f(821) -f(945) -f(952) +f(970) +f(977); d(225,1)= +f(481) +f(482) +f(483) +f(484) +f(485) +f(486) -f(547) -f(548) -f(549) -f(550) -f(551) -f(552) -f(565) -f(566) -f(567) -f(568) -f(569) -f(570) -f(586) +f(635) +f(636) +f(637) +f(638) +f(639) +f(640) +f(791) +f(797) -f(822) -f(835) +f(953) +f(960) -f(978) -f(993); d(226,1)= +f(487) +f(488) +f(489) +f(490) +f(491) +f(492) -f(553) -f(554) -f(555) -f(556) -f(557) -f(558) -f(571) -f(572) -f(573) -f(574) -f(575) -f(576) -f(587) +f(641) +f(642) +f(643) +f(644) +f(645) +f(646) +f(792) +f(798) -f(823) -f(836) +f(954) +f(961) -f(979) -f(994); d(227,1)= +f(493) +f(494) +f(495) +f(496) +f(497) +f(498) -f(559) -f(560) -f(561) -f(562) -f(563) -f(564) -f(577) -f(578) -f(579) -f(580) -f(581) -f(582) -f(588) +f(647) +f(648) +f(649) +f(650) +f(651) +f(652) +f(793) +f(799) -f(824) -f(837) +f(955) +f(962) -f(980) -f(995); d(228,1)= +f(519) -f(594) -f(595); d(229,1)= +f(520) -f(589) -f(591) -f(592) -f(593) -f(654) +f(657) -f(671); d(230,1)= +f(521) -f(590) -f(653) +f(656) +f(845); d(231,1)= +f(526) -f(596) +f(599) +f(600) +f(601) +f(602) +f(603) +f(604) -f(617) -f(618) -f(619) -f(620) -f(621) -f(622) +f(785) -f(816) +f(946) -f(971); d(232,1)= +f(527) -f(597) +f(605) +f(606) +f(607) +f(608) +f(609) +f(610) -f(623) -f(624) -f(625) -f(626) -f(627) -f(628) +f(786) -f(817) +f(947) -f(972); d(233,1)= +f(528) -f(598) +f(611) +f(612) +f(613) +f(614) +f(615) +f(616) -f(629) -f(630) -f(631) -f(632) -f(633) -f(634) +f(787) -f(818) +f(948) -f(973); d(234,1)= +f(547) +f(548) +f(549) +f(550) +f(551) +f(552) -f(599) -f(600) -f(601) -f(602) -f(603) -f(604) +f(617) +f(618) +f(619) +f(620) +f(621) +f(622) -f(635) -f(636) -f(637) -f(638) -f(639) -f(640) -f(785) -f(791) +f(816) +f(822) -f(946) -f(953) +f(971) +f(978); d(235,1)= +f(553) +f(554) +f(555) +f(556) +f(557) +f(558) -f(605) -f(606) -f(607) -f(608) -f(609) -f(610) +f(623) +f(624) +f(625) +f(626) +f(627) +f(628) -f(641) -f(642) -f(643) -f(644) -f(645) -f(646) -f(786) -f(792) +f(817) +f(823) -f(947) -f(954) +f(972) +f(979); d(236,1)= +f(559) +f(560) +f(561) +f(562) +f(563) +f(564) -f(611) -f(612) -f(613) -f(614) -f(615) -f(616) +f(629) +f(630) +f(631) +f(632) +f(633) +f(634) -f(647) -f(648) -f(649) -f(650) -f(651) -f(652) -f(787) -f(793) +f(818) +f(824) -f(948) -f(955) +f(973) +f(980); d(237,1)= +f(589) -f(655) -f(656) -f(657) -f(658) -f(675) +f(678) +f(679); d(238,1)= +f(590) -f(659) -f(660) -f(661) -f(662) -f(663) -f(664) -f(665) -f(666) -f(667) -f(668) -f(669) -f(670) -f(671) -f(676) -f(683) +f(868) -f(872) -f(907); d(239,1)= +f(591) -f(672) +f(697); d(240,1)= +f(592) -f(673) +f(698); d(241,1)= +f(593) -f(674) +f(699); d(242,1)= +f(655) -f(677) -f(678) +f(705); d(243,1)= +f(658) -f(679) -f(680) -f(681) -f(682) -f(683) -f(703) +f(841); d(244,1)= +f(659) -f(684) +f(710); d(245,1)= +f(660) -f(685) +f(711); d(246,1)= +f(661) -f(686) +f(712); d(247,1)= +f(662) -f(687) +f(713); d(248,1)= +f(663) -f(688) +f(714); d(249,1)= +f(664) -f(689) +f(715); d(250,1)= +f(665) -f(690) +f(716); d(251,1)= +f(666) -f(691) +f(717); d(252,1)= +f(667) -f(692) +f(718); d(253,1)= +f(668) -f(693) +f(719); d(254,1)= +f(669) -f(694) +f(720); d(255,1)= +f(670) -f(695) +f(721); d(256,1)= +f(671) -f(696) +f(722); d(257,1)= +f(672) -f(697) -f(700); d(258,1)= +f(673) -f(698) -f(701); d(259,1)= +f(674) -f(699) -f(702); d(260,1)= +f(677) -f(704) -f(705) -f(739) +f(744); d(261,1)= +f(680) -f(706) +f(745); d(262,1)= +f(681) -f(707) +f(746); d(263,1)= +f(682) -f(708) +f(747); d(264,1)= +f(683) -f(709) +f(748); d(265,1)= +f(684) -f(710) -f(723); d(266,1)= +f(685) -f(711) -f(724); d(267,1)= +f(686) -f(712) -f(725); d(268,1)= +f(687) -f(713) -f(726); d(269,1)= +f(688) -f(714) -f(727); d(270,1)= +f(689) -f(715) -f(728); d(271,1)= +f(690) -f(716) -f(729); d(272,1)= +f(691) -f(717) -f(730); d(273,1)= +f(692) -f(718) -f(731); d(274,1)= +f(693) -f(719) -f(732); d(275,1)= +f(694) -f(720) -f(733); d(276,1)= +f(695) -f(721) -f(734); d(277,1)= +f(696) -f(722) -f(735); d(278,1)= +f(700) -f(736) -f(740) +f(924); d(279,1)= +f(701) -f(737) -f(741) +f(925); d(280,1)= +f(702) -f(738) -f(742) +f(926); d(281,1)= +f(704) -f(743) -f(744) +f(771); d(282,1)= +f(706) -f(745) -f(749); d(283,1)= +f(707) -f(746) -f(750); d(284,1)= +f(708) -f(747) -f(751); d(285,1)= +f(709) -f(748) -f(752); d(286,1)= +f(723) +f(726) -f(753) -f(763); d(287,1)= +f(724) +f(729) -f(754) -f(764); d(288,1)= +f(725) +f(732) -f(755) -f(765); d(289,1)= +f(727) +f(733) -f(756) -f(766); d(290,1)= +f(728) +f(731) -f(757) -f(767); d(291,1)= +f(730) +f(734) -f(758) -f(768); d(292,1)= +f(735) -f(759) -f(769) +f(927); d(293,1)= +f(736) -f(760) +f(772); d(294,1)= +f(737) -f(761) +f(773); d(295,1)= +f(738) -f(762) +f(774); d(296,1)= +f(743) -f(770) -f(771) -f(782) -f(783) -f(784) -f(785) -f(786) -f(787) -f(788) -f(789) -f(790) -f(791) -f(792) -f(793) -f(794) -f(795) -f(796) -f(797) -f(798) -f(799) -f(800) +f(806) +f(813) +f(814) +f(815) +f(816) +f(817) +f(818) +f(819) +f(820) +f(821) +f(822) +f(823) +f(824) -f(831) +f(832) +f(833) +f(834) +f(835) +f(836) +f(837) +f(839) +f(857) +f(864) -f(882) +f(889) -f(896) -f(988) +f(1007) +f(1015) -f(1026) +f(1034) -f(1042); d(297,1)= +f(749) -f(772) -f(801) +f(934); d(298,1)= +f(750) -f(773) -f(802) +f(935); d(299,1)= +f(751) -f(774) -f(803) +f(936); d(300,1)= +f(752) -f(775) -f(804) +f(937); d(301,1)= +f(753) +f(755) -f(776); d(302,1)= +f(754) +f(756) -f(777); d(303,1)= +f(757) +f(758) -f(778); d(304,1)= +f(759) +f(775) -f(779); d(305,1)= +f(761) +f(778) -f(780); d(306,1)= +f(762) +f(777) -f(781); d(307,1)= +f(770) -f(806); d(308,1)= +f(779) -f(807); d(309,1)= +f(780) -f(808) -f(809); d(310,1)= +f(781) -f(810) -f(811); d(311,1)= +f(782) +f(783) +f(784) +f(785) +f(786) +f(787) +f(794) +f(795) +f(796) +f(797) +f(798) +f(799) -f(812) -f(813) -f(814) -f(815) -f(816) -f(817) -f(818) +f(831) -f(832) -f(833) -f(834) -f(835) -f(836) -f(837) +f(849) -f(864) +f(882) -f(889) +f(988) -f(1015) +f(1026) -f(1034); d(312,1)= +f(788) +f(789) +f(790) +f(791) +f(792) +f(793) -f(805) -f(819) -f(820) -f(821) -f(822) -f(823) -f(824) -f(825) -f(826) -f(827) -f(828) -f(829) -f(830) -f(831) -f(838) -f(857) +f(858) +f(859) +f(860) +f(861) +f(862) +f(863) +f(864) +f(896) -f(963) +f(996) -f(1007) +f(1042); d(313,1)= +f(805) -f(839); d(314,1)= +f(807) -f(840); d(315,1)= +f(808) -f(841) -f(845); d(316,1)= +f(809) -f(842) -f(846); d(317,1)= +f(810) -f(843) -f(847); d(318,1)= +f(811) -f(844) -f(848); d(319,1)= +f(812) -f(849) -f(850) -f(865) +f(875); d(320,1)= +f(825) +f(826) +f(827) +f(828) +f(829) +f(830) +f(831) -f(851) -f(852) -f(853) -f(854) -f(855) -f(856) -f(857) -f(858) -f(859) -f(860) -f(861) -f(862) -f(863) -f(864) -f(866) +f(890) +f(891) +f(892) +f(893) +f(894) +f(895) +f(896) +f(956) +f(963) -f(981) -f(996); d(321,1)= +f(840) -f(873) -f(874); d(322,1)= +f(842) -f(867) -f(869) -f(870) -f(871) -f(872) -f(898) +f(901); d(323,1)= +f(846) -f(868) +f(874) -f(897) +f(900); d(324,1)= +f(850) -f(875) +f(876) +f(877) +f(878) +f(879) +f(880) +f(881) +f(882) -f(883) -f(884) -f(885) -f(886) -f(887) -f(888) -f(889) +f(949) -f(974); d(325,1)= +f(851) +f(852) +f(853) +f(854) +f(855) +f(856) +f(857) -f(876) -f(877) -f(878) -f(879) -f(880) -f(881) -f(882) +f(883) +f(884) +f(885) +f(886) +f(887) +f(888) +f(889) -f(890) -f(891) -f(892) -f(893) -f(894) -f(895) -f(896) -f(949) -f(956) +f(974) +f(981); d(326,1)= +f(867) -f(899) -f(900) -f(901) -f(902) +f(903) -f(912) +f(915); d(327,1)= +f(869) -f(908) +f(920); d(328,1)= +f(870) -f(909) +f(921); d(329,1)= +f(871) -f(910) +f(922); d(330,1)= +f(872) -f(911) +f(923); d(331,1)= +f(873) +f(902) -f(903) -f(904) -f(905) -f(906) -f(907) -f(913); d(332,1)= +f(899) -f(914) -f(915) +f(929); d(333,1)= +f(904) -f(916) +f(930); d(334,1)= +f(905) -f(917) +f(931); d(335,1)= +f(906) -f(918) +f(932); d(336,1)= +f(907) -f(919) +f(933); d(337,1)= +f(908) -f(920) -f(924); d(338,1)= +f(909) -f(921) -f(925); d(339,1)= +f(910) -f(922) -f(926); d(340,1)= +f(911) -f(923) -f(927); d(341,1)= +f(914) -f(928) -f(929) -f(938) +f(940); d(342,1)= +f(916) -f(930) -f(934); d(343,1)= +f(917) -f(931) -f(935); d(344,1)= +f(918) -f(932) -f(936); d(345,1)= +f(919) -f(933) -f(937); d(346,1)= +f(928) -f(939) -f(940) +f(942); d(347,1)= +f(939) -f(941) -f(942) -f(943) -f(944) -f(945) -f(946) -f(947) -f(948) -f(949) -f(950) -f(951) -f(952) -f(953) -f(954) -f(955) -f(956) -f(957) -f(958) -f(959) -f(960) -f(961) -f(962) -f(963) -f(964) +f(966) +f(968) +f(969) +f(970) +f(971) +f(972) +f(973) +f(974) +f(975) +f(976) +f(977) +f(978) +f(979) +f(980) +f(981) -f(989) +f(990) +f(991) +f(992) +f(993) +f(994) +f(995) +f(996) +f(998) +f(1008) +f(1016) -f(1027) +f(1035) -f(1043); d(348,1)= +f(941) -f(966); d(349,1)= +f(943) +f(944) +f(945) +f(946) +f(947) +f(948) +f(949) +f(957) +f(958) +f(959) +f(960) +f(961) +f(962) +f(963) -f(967) -f(968) -f(969) -f(970) -f(971) -f(972) -f(973) -f(974) +f(989) -f(990) -f(991) -f(992) -f(993) -f(994) -f(995) -f(996) +f(999) -f(1016) +f(1027) -f(1035); d(350,1)= +f(950) +f(951) +f(952) +f(953) +f(954) +f(955) +f(956) -f(965) -f(975) -f(976) -f(977) -f(978) -f(979) -f(980) -f(981) -f(982) -f(983) -f(984) -f(985) -f(986) -f(987) -f(988) -f(989) -f(997) -f(1008) +f(1009) +f(1010) +f(1011) +f(1012) +f(1013) +f(1014) +f(1015) +f(1016) +f(1043); d(351,1)= +f(965) -f(998); d(352,1)= +f(967) -f(999) -f(1000) -f(1017) +f(1019); d(353,1)= +f(982) +f(983) +f(984) +f(985) +f(986) +f(987) +f(988) +f(989) -f(1001) -f(1002) -f(1003) -f(1004) -f(1005) -f(1006) -f(1007) -f(1008) -f(1009) -f(1010) -f(1011) -f(1012) -f(1013) -f(1014) -f(1015) -f(1016) -f(1018) +f(1036) +f(1037) +f(1038) +f(1039) +f(1040) +f(1041) +f(1042) +f(1043); d(354,1)= +f(1000) -f(1019) +f(1020) +f(1021) +f(1022) +f(1023) +f(1024) +f(1025) +f(1026) +f(1027) -f(1028) -f(1029) -f(1030) -f(1031) -f(1032) -f(1033) -f(1034) -f(1035); d(355,1)= +f(1001) +f(1002) +f(1003) +f(1004) +f(1005) +f(1006) +f(1007) +f(1008) -f(1020) -f(1021) -f(1022) -f(1023) -f(1024) -f(1025) -f(1026) -f(1027) +f(1028) +f(1029) +f(1030) +f(1031) +f(1032) +f(1033) +f(1034) +f(1035) -f(1036) -f(1037) -f(1038) -f(1039) -f(1040) -f(1041) -f(1042) -f(1043); end % Integrate ODEs [t_out,x_out]= ode15s(@xdot, [0 tend], x0); % Plot results plot(t_out,x_out); end