 %  DOQCS : http://doqcs.ncbs.res.in/ 
 %  Accession Name = Ajay_Bhalla_2004_Feedback_Tuning 
 %  Accession Number = 78 
 %  Transcriber = Upinder S. Bhalla, NCBS 
 %  Developer = Sriram M. Ajay and Upinder S. Bhalla, NCBS 
 %  Species = Rat 
 %  Tissue = Hippocampal CA1 
 %  Cell Compartment = Dendrites 
 %  Notes =  This model is taken from <a href=http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=pubmed&cmd=Retrieve&dopt=AbstractPlus&list_uids=15548210&query_hl=1&itool=pubmed_docsum>Ajay SM, Bhalla US. Eur J Neurosci. 2004 Nov;20(10):2671-80</a>. This is the feedback model from Figure 8a. 
 
   % Note: All the concentrations are in micro Mole (uM)

function moose_convert
%Matlab simulation for chemical kinetics

% Saved on Wed Jul 25 16:14:42 2007
% Automatically generated from kkit model in MOOSE by conv_matlab

global buff;
global t;
global y;

% Define Tables here
tspan=[0, 2500];
% Molecule names
molnames = {
'1 /kinetics/PKC/PKC-Ca'
'2 /kinetics/PKC/PKC-DAG-AA*'
'3 /kinetics/PKC/PKC-Ca-AA*'
'4 /kinetics/PKC/PKC-Ca-memb*'
'5 /kinetics/PKC/PKC-DAG-memb*'
'6 /kinetics/PKC/PKC-basal*'
'7 /kinetics/PKC/PKC-AA*'
'8 /kinetics/PKC/PKC-Ca-DAG'
'9 /kinetics/PKC/PKC-DAG'
'10 /kinetics/PKC/PKC-DAG-AA'
'11 /kinetics/PKC/PKC-cytosolic'
'12 /kinetics/DAG'
'13 /kinetics/AA'
'14 /kinetics/PKC-active/PKC-act-raf/PKC-act-raf_cplx'
'15 /kinetics/PKC-active/PKC-inact-GAP/PKC-inact-GAP_cplx'
'16 /kinetics/PKC-active/PKC-act-GEF/PKC-act-GEF_cplx'
'17 /kinetics/PKC-active/PKC-phosph-neurogranin/PKC-phosph-neurogranin_cplx'
'18 /kinetics/PKC-active/PKC-phosph-ng-CaM/PKC-phosph-ng-CaM_cplx'
'19 /kinetics/PKC-active/phosph-AC2/phosph-AC2_cplx'
'20 /kinetics/PLA2/PLA2-cytosolic'
'21 /kinetics/PLA2/PLA2-Ca*'
'22 /kinetics/PLA2/PLA2-Ca*/kenz/kenz_cplx'
'23 /kinetics/PLA2/PIP2-PLA2*'
'24 /kinetics/PLA2/PIP2-PLA2*/kenz/kenz_cplx'
'25 /kinetics/PLA2/PIP2-Ca-PLA2*'
'26 /kinetics/PLA2/PIP2-Ca-PLA2*/kenz/kenz_cplx'
'27 /kinetics/PLA2/DAG-Ca-PLA2*'
'28 /kinetics/PLA2/DAG-Ca-PLA2*/kenz/kenz_cplx'
'29 /kinetics/PLA2/PLA2*-Ca'
'30 /kinetics/PLA2/PLA2*-Ca/kenz/kenz_cplx'
'31 /kinetics/PLA2/PLA2*'
'32 /kinetics/MAPK*'
'33 /kinetics/MAPK*/MAPK*/MAPK*_cplx'
'34 /kinetics/MAPK*/MAPK*-feedback/MAPK*-feedback_cplx'
'35 /kinetics/MAPK*/phosph_Sos/phosph_Sos_cplx'
'36 /kinetics/IP3'
'37 /kinetics/PLCbeta/PLC'
'38 /kinetics/PLCbeta/PLC-Ca'
'39 /kinetics/PLCbeta/PLC-Ca/PLC-Ca/PLC-Ca_cplx'
'40 /kinetics/PLCbeta/PLC-Ca-Gq'
'41 /kinetics/PLCbeta/PLC-Ca-Gq/PLCb-Ca-Gq/PLCb-Ca-Gq_cplx'
'42 /kinetics/PLCbeta/PLC-Gq'
'43 /kinetics/BetaGamma'
'44 /kinetics/G*GTP'
'45 /kinetics/G*GDP'
'46 /kinetics/Gq/G-GDP'
'47 /kinetics/Gq/mGluR'
'48 /kinetics/Gq/Rec-Glu'
'49 /kinetics/Gq/Rec-Gq'
'50 /kinetics/Gq/Rec-Glu-Gq'
'51 /kinetics/Gq/Blocked-rec-Gq'
'52 /kinetics/MAPK/craf-1'
'53 /kinetics/MAPK/craf-1*'
'54 /kinetics/MAPK/MAPKK'
'55 /kinetics/MAPK/MAPK'
'56 /kinetics/MAPK/craf-1**'
'57 /kinetics/MAPK/MAPK-tyr'
'58 /kinetics/MAPK/MAPKK*'
'59 /kinetics/MAPK/MAPKK*/MAPKKtyr/MAPKKtyr_cplx'
'60 /kinetics/MAPK/MAPKK*/MAPKKthr/MAPKKthr_cplx'
'61 /kinetics/MAPK/MAPKK-ser'
'62 /kinetics/MAPK/Raf-GTP-Ras'
'63 /kinetics/MAPK/Raf-GTP-Ras/Raf-GTP-Ras.1/Raf-GTP-Ras.1_cplx'
'64 /kinetics/MAPK/Raf-GTP-Ras/Raf-GTP-Ras.2/Raf-GTP-Ras.2_cplx'
'65 /kinetics/MAPK/Raf*-GTP-Ras'
'66 /kinetics/MAPK/Raf*-GTP-Ras/Raf*-GTP-Ras.1/Raf*-GTP-Ras.1_cplx'
'67 /kinetics/MAPK/Raf*-GTP-Ras/Raf*-GTP-Ras.2/Raf*-GTP-Ras.2_cplx'
'68 /kinetics/MKP-1'
'69 /kinetics/MKP-1/MKP1-tyr-deph/MKP1-tyr-deph_cplx'
'70 /kinetics/MKP-1/MKP1-thr-deph/MKP1-thr-deph_cplx'
'71 /kinetics/PPhosphatase2A'
'72 /kinetics/PPhosphatase2A/craf-deph/craf-deph_cplx'
'73 /kinetics/PPhosphatase2A/MAPKK-deph/MAPKK-deph_cplx'
'74 /kinetics/PPhosphatase2A/MAPKK-deph-ser/MAPKK-deph-ser_cplx'
'75 /kinetics/PPhosphatase2A/craf**-deph/craf**-deph_cplx'
'76 /kinetics/Ras/GEF-Gprot-bg'
'77 /kinetics/Ras/GEF-Gprot-bg/GEF-bg_act-ras/GEF-bg_act-ras_cplx'
'78 /kinetics/Ras/inact-GEF'
'79 /kinetics/Ras/GEF*'
'80 /kinetics/Ras/GEF*/GEF*-act-ras/GEF*-act-ras_cplx'
'81 /kinetics/Ras/GTP-Ras'
'82 /kinetics/Ras/GDP-Ras'
'83 /kinetics/Ras/GAP*'
'84 /kinetics/Ras/GAP'
'85 /kinetics/Ras/GAP/GAP-inact-ras/GAP-inact-ras_cplx'
'86 /kinetics/Ras/inact-GEF*'
'87 /kinetics/Ras/CaM-GEF'
'88 /kinetics/Ras/CaM-GEF/CaM-GEF-act-ras/CaM-GEF-act-ras_cplx'
'89 /kinetics/PKA-active'
'90 /kinetics/PKA-active/PKA-phosph-GEF/PKA-phosph-GEF_cplx'
'91 /kinetics/PKA-active/PKA-phosph-I1/PKA-phosph-I1_cplx'
'92 /kinetics/PKA-active/phosph-PDE/phosph-PDE_cplx'
'93 /kinetics/CaM-Ca4'
'94 /kinetics/Shc*.Sos.Grb2'
'95 /kinetics/Shc*.Sos.Grb2/Sos.Ras_GEF/Sos.Ras_GEF_cplx'
'96 /kinetics/EGFR/EGFR'
'97 /kinetics/EGFR/L.EGFR'
'98 /kinetics/EGFR/L.EGFR/phosph_PLC_g/phosph_PLC_g_cplx'
'99 /kinetics/EGFR/L.EGFR/phosph_Shc/phosph_Shc_cplx'
'100 /kinetics/EGFR/SHC'
'101 /kinetics/EGFR/SHC*'
'102 /kinetics/EGFR/Internal_L.EGFR'
'103 /kinetics/Sos/Sos*.Grb2'
'104 /kinetics/Sos/Grb2'
'105 /kinetics/Sos/Sos.Grb2'
'106 /kinetics/Sos/Sos*'
'107 /kinetics/Sos/Sos'
'108 /kinetics/PLC_g/PLC_g'
'109 /kinetics/PLC_g/PLC_G*'
'110 /kinetics/PLC_g/Ca.PLC_g'
'111 /kinetics/PLC_g/Ca.PLC_g*'
'112 /kinetics/CaMKII/CaMKII'
'113 /kinetics/CaMKII/CaMKII-CaM'
'114 /kinetics/CaMKII/CaMKII-thr286*-CaM'
'115 /kinetics/CaMKII/CaMKII***'
'116 /kinetics/CaMKII/CaMKII-thr286'
'117 /kinetics/CaMKII/CaMK-thr306'
'118 /kinetics/CaMKII/tot_CaM_CaMKII/CaM_act_305/CaM_act_305_cplx'
'119 /kinetics/CaMKII/tot_CaM_CaMKII/CaM_act_286/CaM_act_286_cplx'
'120 /kinetics/CaMKII/tot_autonomous_CaMKII/auton_305/auton_305_cplx'
'121 /kinetics/CaMKII/tot_autonomous_CaMKII/auton_286/auton_286_cplx'
'122 /kinetics/PP1-active'
'123 /kinetics/PP1-active/Deph-thr286/Deph-thr286_cplx'
'124 /kinetics/PP1-active/Deph-thr305/Deph-thr305_cplx'
'125 /kinetics/PP1-active/Deph-thr306/Deph-thr306_cplx'
'126 /kinetics/PP1-active/Deph-thr286c/Deph-thr286c_cplx'
'127 /kinetics/PP1-active/Deph_thr286b/Deph_thr286b_cplx'
'128 /kinetics/CaM/CaM'
'129 /kinetics/CaM/neurogranin-CaM'
'130 /kinetics/CaM/neurogranin*'
'131 /kinetics/CaM/neurogranin'
'132 /kinetics/CaM-Ca3'
'133 /kinetics/CaM-TR2-Ca2'
'134 /kinetics/CaM(Ca)n-CaNAB/dephosph_neurogranin/dephosph_neurogranin_cplx'
'135 /kinetics/CaM(Ca)n-CaNAB/dephosph_inhib1/dephosph_inhib1_cplx'
'136 /kinetics/CaM(Ca)n-CaNAB/dephosph-PP1-I*/dephosph-PP1-I*_cplx'
'137 /kinetics/PP1/I1'
'138 /kinetics/PP1/I1*'
'139 /kinetics/PP1/PP1-I1*'
'140 /kinetics/PP1/PP1-I1'
'141 /kinetics/PP2A'
'142 /kinetics/PP2A/PP2A-dephosph-I1/PP2A-dephosph-I1_cplx'
'143 /kinetics/PP2A/PP2A-dephosph-PP1-I*/PP2A-dephosph-PP1-I*_cplx'
'144 /kinetics/CaNAB-Ca4'
'145 /kinetics/CaNAB-Ca4/dephosph_inhib1_noCaM/dephosph_inhib1_noCaM_cplx'
'146 /kinetics/PP2B/CaNAB'
'147 /kinetics/PP2B/CaNAB-Ca2'
'148 /kinetics/PP2B/CaMCa3-CaNAB'
'149 /kinetics/PP2B/CaMCa2-CANAB'
'150 /kinetics/PP2B/CaMCa4-CaNAB'
'151 /kinetics/PKA/R2C2'
'152 /kinetics/PKA/R2C2-cAMP'
'153 /kinetics/PKA/R2C2-cAMP2'
'154 /kinetics/PKA/R2C2-cAMP3'
'155 /kinetics/PKA/R2C2-cAMP4'
'156 /kinetics/PKA/R2C-cAMP4'
'157 /kinetics/PKA/R2-cAMP4'
'158 /kinetics/PKA/PKA-inhibitor'
'159 /kinetics/PKA/inhibited-PKA'
'160 /kinetics/cAMP'
'161 /kinetics/AC/AC1-CaM'
'162 /kinetics/AC/AC1'
'163 /kinetics/AC/AC2*'
'164 /kinetics/AC/AC2-Gs'
'165 /kinetics/AC/AC2'
'166 /kinetics/AC/AC1-Gs'
'167 /kinetics/AC/AMP'
'168 /kinetics/AC/AC2*-Gs'
'169 /kinetics/AC/cAMP-PDE'
'170 /kinetics/AC/cAMP-PDE/PDE/PDE_cplx'
'171 /kinetics/AC/cAMP-PDE*'
'172 /kinetics/AC/cAMP-PDE*/PDE*/PDE*_cplx'
'173 /kinetics/AC/PDE1'
'174 /kinetics/AC/PDE1/PDE1/PDE1_cplx'
'175 /kinetics/AC/CaM.PDE1'
'176 /kinetics/AC/CaM.PDE1/CaM.PDE1/CaM.PDE1_cplx'
'177 /kinetics/Gs-alpha'
'178 /kinetics/Ca'
'179 /kinetics/Ca_input'
};

bufnames = {
'1 /kinetics/PLA2/APC'
'2 /kinetics/temp-PIP2'
'3 /kinetics/Glu'
'4 /kinetics/PLCbeta/Inositol'
'5 /kinetics/PLCbeta/PC'
'6 /kinetics/PIP2'
'7 /kinetics/Gq/mGluRAntag'
'8 /kinetics/EGFR/EGF'
'9 /kinetics/AC/ATP'
};

sumtotnames = {
'1 /kinetics/PKC-active'
'2 /kinetics/CaMKII/tot_CaM_CaMKII'
'3 /kinetics/CaMKII/tot_autonomous_CaMKII'
'4 /kinetics/CaM(Ca)n-CaNAB'
};

% Initial concs of 179 molecules
y0(1) = 3.72078e-17;
y0(2) = 4.91367e-18;
y0(3) = 1.75e-16;
y0(4) = 1.38956e-17;
y0(5) = 9.43511e-21;
y0(6) = 0;
y0(7) = 1.81333e-17;
y0(8) = 8.46311e-23;
y0(9) = 1.161e-16;
y0(10) = 2.51889e-19;
y0(11) = 1;
y0(12) = 0;
y0(13) = 0;
y0(14) = 0;
y0(15) = 0;
y0(16) = 0;
y0(17) = 0;
y0(18) = 0;
y0(19) = 0;
y0(20) = 0.4;
y0(21) = 0;
y0(22) = 0;
y0(23) = 0;
y0(24) = 0;
y0(25) = 0;
y0(26) = 0;
y0(27) = 0;
y0(28) = 0;
y0(29) = 0;
y0(30) = 0;
y0(31) = 0;
y0(32) = 0;
y0(33) = 0;
y0(34) = 0;
y0(35) = 0;
y0(36) = 0.73;
y0(37) = 0.8;
y0(38) = 0;
y0(39) = 0;
y0(40) = 0;
y0(41) = 0;
y0(42) = 0;
y0(43) = 0;
y0(44) = 0;
y0(45) = 0;
y0(46) = 1;
y0(47) = 0.3;
y0(48) = 0;
y0(49) = 0;
y0(50) = 0;
y0(51) = 0;
y0(52) = 0.5;
y0(53) = 0;
y0(54) = 0.5;
y0(55) = 3.6;
y0(56) = 0;
y0(57) = 0;
y0(58) = 0;
y0(59) = 0;
y0(60) = 0;
y0(61) = 0;
y0(62) = 0;
y0(63) = 0;
y0(64) = 0;
y0(65) = 0;
y0(66) = 0;
y0(67) = 0;
y0(68) = 0.024;
y0(69) = 0;
y0(70) = 0;
y0(71) = 1;
y0(72) = 0;
y0(73) = 0;
y0(74) = 0;
y0(75) = 0;
y0(76) = 0;
y0(77) = 0;
y0(78) = 0.1;
y0(79) = 0;
y0(80) = 0;
y0(81) = 0;
y0(82) = 0.5;
y0(83) = 0;
y0(84) = 0.01;
y0(85) = 0;
y0(86) = 0;
y0(87) = 0;
y0(88) = 0;
y0(89) = 0;
y0(90) = 0;
y0(91) = 0;
y0(92) = 0;
y0(93) = 0;
y0(94) = 0;
y0(95) = 0;
y0(96) = 0.166667;
y0(97) = 0;
y0(98) = 0;
y0(99) = 0;
y0(100) = 0.5;
y0(101) = 0;
y0(102) = 0;
y0(103) = 0;
y0(104) = 1;
y0(105) = 0;
y0(106) = 0;
y0(107) = 0.1;
y0(108) = 0.82;
y0(109) = 0;
y0(110) = 0;
y0(111) = 0;
y0(112) = 70;
y0(113) = 0;
y0(114) = 0;
y0(115) = 0;
y0(116) = 0;
y0(117) = 0;
y0(118) = 0;
y0(119) = 0;
y0(120) = 0;
y0(121) = 0;
y0(122) = 1.8;
y0(123) = 0;
y0(124) = 0;
y0(125) = 0;
y0(126) = 0;
y0(127) = 0;
y0(128) = 20;
y0(129) = 0;
y0(130) = 0;
y0(131) = 10;
y0(132) = 0;
y0(133) = 0;
y0(134) = 0;
y0(135) = 0;
y0(136) = 0;
y0(137) = 1.8;
y0(138) = 0.001;
y0(139) = 0;
y0(140) = 0;
y0(141) = 0.12;
y0(142) = 0;
y0(143) = 0;
y0(144) = 0;
y0(145) = 0;
y0(146) = 1;
y0(147) = 0;
y0(148) = 0;
y0(149) = 0;
y0(150) = 0;
y0(151) = 0.5;
y0(152) = 0;
y0(153) = 0;
y0(154) = 0;
y0(155) = 0;
y0(156) = 0;
y0(157) = 0;
y0(158) = 0.25;
y0(159) = 0;
y0(160) = 0;
y0(161) = 0;
y0(162) = 0.02;
y0(163) = 0;
y0(164) = 0;
y0(165) = 0.015;
y0(166) = 0;
y0(167) = 325493;
y0(168) = 0;
y0(169) = 0.45;
y0(170) = 0;
y0(171) = 0;
y0(172) = 0;
y0(173) = 2;
y0(174) = 0;
y0(175) = 0;
y0(176) = 0;
y0(177) = 0;
y0(178) = 0.08;
y0(179) = 0;
buff(1) = 30;
buff(2) = 2.5;
buff(3) = 0;
buff(4) = 0;
buff(5) = 0;
buff(6) = 10;
buff(7) = 0;
buff(8) = 0;
buff(9) = 5000;
% Sumtotal eqn for /kinetics/PKC-active: sumtot(1) =  (y(2) + y(4) + y(3) + y(5) + y(6) + y(7) )
% Sumtotal eqn for /kinetics/CaMKII/tot_CaM_CaMKII: sumtot(2) =  (y(113) + y(114) )
% Sumtotal eqn for /kinetics/CaMKII/tot_autonomous_CaMKII: sumtot(3) =  (y(116) + y(115) )
% Sumtotal eqn for /kinetics/CaM(Ca)n-CaNAB: sumtot(4) =  (y(150) + y(148) + y(149) )

% Call the ODE solver and display results
[t, y] = ode23s(@f, tspan, y0);

hold on;
do_plots;

function do_plots()
global t;
global y;
global buff;
plot(t,y(:,113));
plot(t,y(:,178));
plot(t,y(:,13));
plot(t,y(:,36));
plot(t,y(:,114));
plot(t,y(:,116));
plot(t,y(:,115));
plot(t,y(:,32));
plot(t,y(:,89));
plot(t, (y(:,2) + y(:,4) + y(:,3) + y(:,5) + y(:,6) + y(:,7) ));
plot(t,y(:,122));
plot(t,y(:,68));
plot(t,y(:,93));

% Evaluation function
function dydt = f(t, y)
global buff;
dydt = [
 + 0.600003 * y(11) * y(178) - 0.5 * y(1) - 0.00799983 * y(12) * y(1) + 8.6348 * y(8) - 1.2705 * y(1) + 3.5026 * y(4) - 0.00119997 * y(1) * y(13) + 0.1 * y(3)
 + 2 * y(10) - 0.2 * y(2)
 + 0.00119997 * y(1) * y(13) - 0.1 * y(3)
 + 1.2705 * y(1) - 3.5026 * y(4)
 + 1 * y(8) - 0.1 * y(5)
 + 1 * y(11) - 50 * y(6)
 + 0.000119997 * y(13) * y(11) - 0.1 * y(7)
 + 0.00799983 * y(12) * y(1) - 8.6348 * y(8) - 1 * y(8) + 0.1 * y(5)
 + 0.000600003 * y(11) * y(12) - 0.1 * y(9) - 0.018 * y(9) * y(13) + 2 * y(10)
 + 0.018 * y(9) * y(13) - 2 * y(10) - 2 * y(10) + 0.2 * y(2)
 - 0.600003 * y(11) * y(178) + 0.5 * y(1) - 1 * y(11) + 50 * y(6) - 0.000119997 * y(13) * y(11) + 0.1 * y(7) - 0.000600003 * y(11) * y(12) + 0.1 * y(9)
 - 0.00799983 * y(12) * y(1) + 8.6348 * y(8) - 0.000600003 * y(11) * y(12) + 0.1 * y(9) - 0.00299997 * y(12) * y(21) + 4 * y(27) + 48 * y(41) - 0.15 * y(12) + 10 * y(39) + 14 * y(110) * buff(6) / ( 97.2222 + buff(6) ) + 57 * y(111) * buff(6) / ( 19.7917 + buff(6) )
 - 0.00119997 * y(1) * y(13) + 0.1 * y(3) - 0.000119997 * y(13) * y(11) + 0.1 * y(7) - 0.018 * y(9) * y(13) + 2 * y(10) + 5.4 * y(22) + 11.04 * y(24) + 36 * y(26) + 60 * y(28) + 120 * y(30) - 0.4 * y(13)
 + 0.299997 * y(52) *  (y(2) + y(4) + y(3) + y(5) + y(6) + y(7) ) - 16 * y(14) - 4 * y(14)
 + 6.00003 * y(84) *  (y(2) + y(4) + y(3) + y(5) + y(6) + y(7) ) - 16 * y(15) - 4 * y(15)
 + 6.00003 * y(78) *  (y(2) + y(4) + y(3) + y(5) + y(6) + y(7) ) - 16 * y(16) - 4 * y(16)
 + 0.101997 * y(131) *  (y(2) + y(4) + y(3) + y(5) + y(6) + y(7) ) - 2.34 * y(17) - 0.58 * y(17)
 + 0.0612 * y(129) *  (y(2) + y(4) + y(3) + y(5) + y(6) + y(7) ) - 1.4 * y(18) - 0.35 * y(18)
 + 0.600003 * y(165) *  (y(2) + y(4) + y(3) + y(5) + y(6) + y(7) ) - 16 * y(19) - 4 * y(19)
 - 0.99999 * y(20) * y(178) + 0.1 * y(21) - 0.00119997 * buff(2) * y(20) + 0.5 * y(23) - 1.95003 * y(20) * y(32) + 40 * y(33) + 0.17 * y(31)
 + 0.99999 * y(20) * y(178) - 0.1 * y(21) - 1.35 * buff(1) * y(21) + 21.6 * y(22) + 5.4 * y(22) - 0.0119997 * buff(2) * y(21) + 0.1 * y(25) - 0.00299997 * y(12) * y(21) + 4 * y(27)
 + 1.35 * buff(1) * y(21) - 21.6 * y(22) - 5.4 * y(22)
 + 0.00119997 * buff(2) * y(20) - 0.5 * y(23) - 2.76003 * buff(1) * y(23) + 44.16 * y(24) + 11.04 * y(24)
 + 2.76003 * buff(1) * y(23) - 44.16 * y(24) - 11.04 * y(24)
 + 0.0119997 * buff(2) * y(21) - 0.1 * y(25) - 9 * buff(1) * y(25) + 144 * y(26) + 36 * y(26)
 + 9 * buff(1) * y(25) - 144 * y(26) - 36 * y(26)
 + 0.00299997 * y(12) * y(21) - 4 * y(27) - 15.0003 * buff(1) * y(27) + 240 * y(28) + 60 * y(28)
 + 15.0003 * buff(1) * y(27) - 240 * y(28) - 60 * y(28)
 - 29.9997 * buff(1) * y(29) + 480 * y(30) + 120 * y(30) + 6.00003 * y(31) * y(178) - 0.1 * y(29)
 + 29.9997 * buff(1) * y(29) - 480 * y(30) - 120 * y(30)
 + 10 * y(33) - 6.00003 * y(31) * y(178) + 0.1 * y(29) - 0.17 * y(31)
 - 1.95003 * y(20) * y(32) + 40 * y(33) + 10 * y(33) - 1.95003 * y(53) * y(32) + 40 * y(34) + 10 * y(34) + 0.3 * y(60) - 150.003 * y(32) * y(68) + 16 * y(70) - 19.5003 * y(107) * y(32) + 40 * y(35) + 10 * y(35)
 + 1.95003 * y(20) * y(32) - 40 * y(33) - 10 * y(33)
 + 1.95003 * y(53) * y(32) - 40 * y(34) - 10 * y(34)
 + 19.5003 * y(107) * y(32) - 40 * y(35) - 10 * y(35)
 + 48 * y(41) - 2.5 * y(36) + 10 * y(39) + 14 * y(110) * buff(6) / ( 97.2222 + buff(6) ) + 57 * y(111) * buff(6) / ( 19.7917 + buff(6) )
 - 2.99997 * y(178) * y(37) + 1 * y(38) - 2.52 * y(37) * y(44) + 1 * y(42)
 + 2.99997 * y(178) * y(37) - 1 * y(38) - 25.2 * y(44) * y(38) + 1 * y(40) + 0.0133 * y(40) - 2.52 * buff(6) * y(38) + 40 * y(39) + 10 * y(39)
 + 2.52 * buff(6) * y(38) - 40 * y(39) - 10 * y(39)
 - 47.9997 * buff(6) * y(40) + 192 * y(41) + 48 * y(41) + 25.2 * y(44) * y(38) - 1 * y(40) - 0.0133 * y(40) + 29.9997 * y(178) * y(42) - 1 * y(40)
 + 47.9997 * buff(6) * y(40) - 192 * y(41) - 48 * y(41)
 + 2.52 * y(37) * y(44) - 1 * y(42) - 29.9997 * y(178) * y(42) + 1 * y(40)
 - 6.00003 * y(45) * y(43) + 0.0001 * y(46) + 0.01 * y(50) - 6.00003 * y(43) * y(78) + 1 * y(76)
 - 25.2 * y(44) * y(38) + 1 * y(40) - 2.52 * y(37) * y(44) + 1 * y(42) - 0.0133 * y(44) + 0.0001 * y(46) + 0.01 * y(50)
 + 0.0133 * y(40) + 0.0133 * y(44) - 6.00003 * y(45) * y(43)
 + 6.00003 * y(45) * y(43) - 0.0001 * y(46) - 0.00600003 * y(46) * y(48) + 0.0001 * y(50) - 0.600003 * y(46) * y(47) + 1 * y(49)
 - 16.8003 * y(47) * buff(3) + 10 * y(48) - 0.600003 * y(46) * y(47) + 1 * y(49)
 + 16.8003 * y(47) * buff(3) - 10 * y(48) - 0.00600003 * y(46) * y(48) + 0.0001 * y(50) + 0.01 * y(50)
 - 16.8003 * buff(3) * y(49) + 0.1 * y(50) + 0.600003 * y(46) * y(47) - 1 * y(49) - 60.0003 * y(49) * buff(7) + 0.01 * y(51)
 + 0.00600003 * y(46) * y(48) - 0.0001 * y(50) + 16.8003 * buff(3) * y(49) - 0.1 * y(50) - 0.01 * y(50)
 + 60.0003 * y(49) * buff(7) - 0.01 * y(51)
 - 0.299997 * y(52) *  (y(2) + y(4) + y(3) + y(5) + y(6) + y(7) ) + 16 * y(14) + 6 * y(72) - 6.00003 * y(52) * y(81) + 1 * y(62)
 + 4 * y(14) - 1.95003 * y(53) * y(32) + 40 * y(34) - 1.9161 * y(53) * y(71) + 24 * y(72) + 6 * y(75) - 9.9999 * y(53) * y(81) + 0.5 * y(65)
 + 6 * y(74) - 9.4284 * y(54) * y(65) + 1.2 * y(66) - 9.4284 * y(54) * y(62) + 1.2 * y(63)
 - 32.4 * y(55) * y(58) + 1.2 * y(59) + 4 * y(69)
 + 10 * y(34) - 1.9161 * y(56) * y(71) + 24 * y(75)
 + 0.3 * y(59) - 32.4 * y(57) * y(58) + 1.2 * y(60) - 150.003 * y(57) * y(68) + 16 * y(69) + 4 * y(70)
 - 32.4 * y(55) * y(58) + 1.2 * y(59) + 0.3 * y(59) - 32.4 * y(57) * y(58) + 1.2 * y(60) + 0.3 * y(60) - 1.9161 * y(58) * y(71) + 24 * y(73) + 0.3 * y(67) + 0.3 * y(64)
 + 32.4 * y(55) * y(58) - 1.2 * y(59) - 0.3 * y(59)
 + 32.4 * y(57) * y(58) - 1.2 * y(60) - 0.3 * y(60)
 + 6 * y(73) - 1.9161 * y(61) * y(71) + 24 * y(74) + 0.3 * y(66) - 9.4284 * y(61) * y(65) + 1.2 * y(67) + 0.3 * y(63) - 9.4284 * y(61) * y(62) + 1.2 * y(64)
 + 6.00003 * y(52) * y(81) - 1 * y(62) - 9.4284 * y(54) * y(62) + 1.2 * y(63) + 0.3 * y(63) - 9.4284 * y(61) * y(62) + 1.2 * y(64) + 0.3 * y(64)
 + 9.4284 * y(54) * y(62) - 1.2 * y(63) - 0.3 * y(63)
 + 9.4284 * y(61) * y(62) - 1.2 * y(64) - 0.3 * y(64)
 + 9.9999 * y(53) * y(81) - 0.5 * y(65) - 9.4284 * y(54) * y(65) + 1.2 * y(66) + 0.3 * y(66) - 9.4284 * y(61) * y(65) + 1.2 * y(67) + 0.3 * y(67)
 + 9.4284 * y(54) * y(65) - 1.2 * y(66) - 0.3 * y(66)
 + 9.4284 * y(61) * y(65) - 1.2 * y(67) - 0.3 * y(67)
 - 150.003 * y(57) * y(68) + 16 * y(69) + 4 * y(69) - 150.003 * y(32) * y(68) + 16 * y(70) + 4 * y(70)
 + 150.003 * y(57) * y(68) - 16 * y(69) - 4 * y(69)
 + 150.003 * y(32) * y(68) - 16 * y(70) - 4 * y(70)
 - 1.9161 * y(53) * y(71) + 24 * y(72) + 6 * y(72) - 1.9161 * y(58) * y(71) + 24 * y(73) + 6 * y(73) - 1.9161 * y(61) * y(71) + 24 * y(74) + 6 * y(74) - 1.9161 * y(56) * y(71) + 24 * y(75) + 6 * y(75)
 + 1.9161 * y(53) * y(71) - 24 * y(72) - 6 * y(72)
 + 1.9161 * y(58) * y(71) - 24 * y(73) - 6 * y(73)
 + 1.9161 * y(61) * y(71) - 24 * y(74) - 6 * y(74)
 + 1.9161 * y(56) * y(71) - 24 * y(75) - 6 * y(75)
 - 0.198 * y(82) * y(76) + 0.08 * y(77) + 0.02 * y(77) + 6.00003 * y(43) * y(78) - 1 * y(76)
 + 0.198 * y(82) * y(76) - 0.08 * y(77) - 0.02 * y(77)
 - 6.00003 * y(43) * y(78) + 1 * y(76) - 6.00003 * y(78) *  (y(2) + y(4) + y(3) + y(5) + y(6) + y(7) ) + 16 * y(16) + 1 * y(79) - 6.00003 * y(78) * y(89) + 36 * y(90) - 60.0003 * y(78) * y(93) + 1 * y(87) + 1 * y(86)
 + 4 * y(16) - 1 * y(79) - 0.198 * y(82) * y(79) + 0.08 * y(80) + 0.02 * y(80)
 + 0.198 * y(82) * y(79) - 0.08 * y(80) - 0.02 * y(80)
 + 0.02 * y(77) - 49.4856 * y(81) * y(84) + 40 * y(85) - 0.0001 * y(81) + 0.02 * y(80) + 0.02 * y(88) - 9.9999 * y(53) * y(81) + 0.5 * y(65) + 0.02 * y(95) - 6.00003 * y(52) * y(81) + 1 * y(62)
 - 0.198 * y(82) * y(76) + 0.08 * y(77) + 10 * y(85) + 0.0001 * y(81) - 0.198 * y(82) * y(79) + 0.08 * y(80) - 0.198 * y(82) * y(87) + 0.08 * y(88) - 0.198 * y(82) * y(94) + 0.08 * y(95)
 + 4 * y(15) - 0.1 * y(83)
 - 49.4856 * y(81) * y(84) + 40 * y(85) + 10 * y(85) - 6.00003 * y(84) *  (y(2) + y(4) + y(3) + y(5) + y(6) + y(7) ) + 16 * y(15) + 0.1 * y(83)
 + 49.4856 * y(81) * y(84) - 40 * y(85) - 10 * y(85)
 + 9 * y(90) - 1 * y(86)
 + 60.0003 * y(78) * y(93) - 1 * y(87) - 0.198 * y(82) * y(87) + 0.08 * y(88) + 0.02 * y(88)
 + 0.198 * y(82) * y(87) - 0.08 * y(88) - 0.02 * y(88)
 - 6.00003 * y(78) * y(89) + 36 * y(90) + 9 * y(90) - 6.00003 * y(137) * y(89) + 36 * y(91) + 9 * y(91) + 60 * y(155) - 18 * y(89) * y(156) + 60 * y(156) - 18 * y(89) * y(157) - 60.0003 * y(89) * y(158) + 1 * y(159) - 6.00003 * y(169) * y(89) + 36 * y(92) + 9 * y(92)
 + 6.00003 * y(78) * y(89) - 36 * y(90) - 9 * y(90)
 + 6.00003 * y(137) * y(89) - 36 * y(91) - 9 * y(91)
 + 6.00003 * y(169) * y(89) - 36 * y(92) - 9 * y(92)
 - 60.0003 * y(78) * y(93) + 1 * y(87) - 49.9995 * y(93) * y(112) + 5 * y(113) - 1000.17 * y(116) * y(93) + 0.1 * y(114) + 0.465003 * y(132) * y(178) - 10 * y(93) - 600.003 * y(93) * y(144) + 1 * y(150) - 49.9995 * y(93) * y(162) + 1 * y(161) - 720 * y(173) * y(93) + 5 * y(175)
 + 0.499977 * y(105) * y(101) - 0.1 * y(94) - 0.198 * y(82) * y(94) + 0.08 * y(95) + 0.02 * y(95)
 + 0.198 * y(82) * y(94) - 0.08 * y(95) - 0.02 * y(95)
 - 4.20003 * y(96) * buff(8) + 0.25 * y(97)
 + 4.20003 * y(96) * buff(8) - 0.25 * y(97) - 2.99997 * y(110) * y(97) + 0.8 * y(98) + 0.2 * y(98) - 1.19997 * y(100) * y(97) + 0.8 * y(99) + 0.2 * y(99) - 0.002 * y(97) + 0.00033 * y(102)
 + 2.99997 * y(110) * y(97) - 0.8 * y(98) - 0.2 * y(98)
 + 1.19997 * y(100) * y(97) - 0.8 * y(99) - 0.2 * y(99)
 + 0.0016667 * y(101) - 1.19997 * y(100) * y(97) + 0.8 * y(99)
 - 0.0016667 * y(101) - 0.499977 * y(105) * y(101) + 0.1 * y(94) + 0.2 * y(99)
 + 0.002 * y(97) - 0.00033 * y(102)
 + 0.0250002 * y(106) * y(104) - 0.0168 * y(103)
 - 0.0250002 * y(104) * y(107) + 0.0168 * y(105) - 0.0250002 * y(106) * y(104) + 0.0168 * y(103)
 + 0.0250002 * y(104) * y(107) - 0.0168 * y(105) - 0.499977 * y(105) * y(101) + 0.1 * y(94)
 + 10 * y(35) - 0.0250002 * y(106) * y(104) + 0.0168 * y(103) - 0.001 * y(106)
 - 0.0250002 * y(104) * y(107) + 0.0168 * y(105) - 19.5003 * y(107) * y(32) + 40 * y(35) + 0.001 * y(106)
 - 180 * y(108) * y(178) + 10 * y(110)
 - 11.9997 * y(178) * y(109) + 10 * y(111)
 + 180 * y(108) * y(178) - 10 * y(110) - 2.99997 * y(110) * y(97) + 0.8 * y(98) + 0.05 * y(111)
 + 11.9997 * y(178) * y(109) - 10 * y(111) + 0.2 * y(98) - 0.05 * y(111)
 - 49.9995 * y(93) * y(112) + 5 * y(113) + 0.35 * y(125) + 0.35 * y(127) - 0.003 * y(112)
 + 49.9995 * y(93) * y(112) - 5 * y(113) + 0.35 * y(123) - 0.0220068 * y(113) *  (y(113) + y(114) ) + 2 * y(119) - 0.0142857 * y(113) *  (y(116) + y(115) ) + 2 * y(121)
 + 1000.17 * y(116) * y(93) - 0.1 * y(114) - 0.343197 * y(114) * y(122) + 1.4 * y(123) + 0.5 * y(119) + 0.5 * y(121)
 - 0.343197 * y(115) * y(122) + 1.4 * y(124) - 0.343197 * y(115) * y(122) + 1.4 * y(126) + 6 * y(118) + 6 * y(120)
 - 1000.17 * y(116) * y(93) + 0.1 * y(114) + 0.35 * y(124) - 0.343197 * y(116) * y(122) + 1.4 * y(127) + 0.003 * y(112) - 0.264087 * y(116) *  (y(113) + y(114) ) + 24 * y(118) - 0.171423 * y(116) *  (y(116) + y(115) ) + 24 * y(120)
 - 0.343197 * y(117) * y(122) + 1.4 * y(125) + 0.35 * y(126)
 + 0.264087 * y(116) *  (y(113) + y(114) ) - 24 * y(118) - 6 * y(118)
 + 0.0220068 * y(113) *  (y(113) + y(114) ) - 2 * y(119) - 0.5 * y(119)
 + 0.171423 * y(116) *  (y(116) + y(115) ) - 24 * y(120) - 6 * y(120)
 + 0.0142857 * y(113) *  (y(116) + y(115) ) - 2 * y(121) - 0.5 * y(121)
 - 0.343197 * y(114) * y(122) + 1.4 * y(123) + 0.35 * y(123) - 0.343197 * y(115) * y(122) + 1.4 * y(124) + 0.35 * y(124) - 0.343197 * y(117) * y(122) + 1.4 * y(125) + 0.35 * y(125) - 0.343197 * y(116) * y(122) + 1.4 * y(127) + 0.35 * y(127) - 0.343197 * y(115) * y(122) + 1.4 * y(126) + 0.35 * y(126) - 499.977 * y(138) * y(122) + 0.1 * y(139) + 1 * y(140)
 + 0.343197 * y(114) * y(122) - 1.4 * y(123) - 0.35 * y(123)
 + 0.343197 * y(115) * y(122) - 1.4 * y(124) - 0.35 * y(124)
 + 0.343197 * y(117) * y(122) - 1.4 * y(125) - 0.35 * y(125)
 + 0.343197 * y(115) * y(122) - 1.4 * y(126) - 0.35 * y(126)
 + 0.343197 * y(116) * y(122) - 1.4 * y(127) - 0.35 * y(127)
 - 72.0001 * y(128) * y(178)*y(178) + 72 * y(133) - 0.299997 * y(131) * y(128) + 1 * y(129) + 0.35 * y(18)
 + 0.299997 * y(131) * y(128) - 1 * y(129) - 0.0612 * y(129) *  (y(2) + y(4) + y(3) + y(5) + y(6) + y(7) ) + 1.4 * y(18)
 - 0.333603 * y(130) *  (y(150) + y(148) + y(149) ) + 2.67 * y(134) + 0.58 * y(17) + 0.35 * y(18) - 0.005 * y(130)
 + 0.67 * y(134) - 0.299997 * y(131) * y(128) + 1 * y(129) - 0.101997 * y(131) *  (y(2) + y(4) + y(3) + y(5) + y(6) + y(7) ) + 2.34 * y(17) + 0.005 * y(130)
 + 3.6 * y(133) * y(178) - 10 * y(132) - 0.465003 * y(132) * y(178) + 10 * y(93) - 2.23803 * y(132) * y(144) + 1 * y(148)
 + 72.0001 * y(128) * y(178)*y(178) - 72 * y(133) - 3.6 * y(133) * y(178) + 10 * y(132) - 0.240003 * y(144) * y(133) + 1 * y(149)
 + 0.333603 * y(130) *  (y(150) + y(148) + y(149) ) - 2.67 * y(134) - 0.67 * y(134)
 + 0.342 * y(138) *  (y(150) + y(148) + y(149) ) - 1.36 * y(135) - 0.34 * y(135)
 + 0.342 * y(139) *  (y(150) + y(148) + y(149) ) - 1.36 * y(136) - 0.34 * y(136)
 - 6.00003 * y(137) * y(89) + 36 * y(91) + 0.034 * y(145) + 0.34 * y(135) + 6 * y(142) + 1 * y(140)
 - 499.977 * y(138) * y(122) + 0.1 * y(139) + 9 * y(91) - 0.0342 * y(138) * y(144) + 0.136 * y(145) - 0.342 * y(138) *  (y(150) + y(148) + y(149) ) + 1.36 * y(135) - 3.96 * y(138) * y(141) + 25 * y(142)
 + 499.977 * y(138) * y(122) - 0.1 * y(139) - 3.96 * y(139) * y(141) + 25 * y(143) - 0.342 * y(139) *  (y(150) + y(148) + y(149) ) + 1.36 * y(136)
 - 1 * y(140) + 6 * y(143) + 0.34 * y(136)
 - 3.96 * y(138) * y(141) + 25 * y(142) + 6 * y(142) - 3.96 * y(139) * y(141) + 25 * y(143) + 6 * y(143)
 + 3.96 * y(138) * y(141) - 25 * y(142) - 6 * y(142)
 + 3.96 * y(139) * y(141) - 25 * y(143) - 6 * y(143)
 - 0.0342 * y(138) * y(144) + 0.136 * y(145) + 0.034 * y(145) + 3.59996 * y(178)*y(178) * y(147) - 1 * y(144) - 0.240003 * y(144) * y(133) + 1 * y(149) - 600.003 * y(93) * y(144) + 1 * y(150) - 2.23803 * y(132) * y(144) + 1 * y(148)
 + 0.0342 * y(138) * y(144) - 0.136 * y(145) - 0.034 * y(145)
 - 10008.4 * y(146) * y(178)*y(178) + 1 * y(147)
 + 10008.4 * y(146) * y(178)*y(178) - 1 * y(147) - 3.59996 * y(178)*y(178) * y(147) + 1 * y(144)
 + 2.23803 * y(132) * y(144) - 1 * y(148)
 + 0.240003 * y(144) * y(133) - 1 * y(149)
 + 600.003 * y(93) * y(144) - 1 * y(150)
 - 54 * y(151) * y(160) + 33 * y(152)
 + 54 * y(151) * y(160) - 33 * y(152) - 54 * y(152) * y(160) + 33 * y(153)
 + 54 * y(152) * y(160) - 33 * y(153) - 74.9997 * y(153) * y(160) + 110 * y(154)
 + 74.9997 * y(153) * y(160) - 110 * y(154) - 74.9997 * y(160) * y(154) + 32.5 * y(155)
 + 74.9997 * y(160) * y(154) - 32.5 * y(155) - 60 * y(155) + 18 * y(89) * y(156)
 + 60 * y(155) - 18 * y(89) * y(156) - 60 * y(156) + 18 * y(89) * y(157)
 + 60 * y(156) - 18 * y(89) * y(157)
 - 60.0003 * y(89) * y(158) + 1 * y(159)
 + 60.0003 * y(89) * y(158) - 1 * y(159)
 - 54 * y(151) * y(160) + 33 * y(152) - 54 * y(152) * y(160) + 33 * y(153) - 74.9997 * y(153) * y(160) + 110 * y(154) - 74.9997 * y(160) * y(154) + 32.5 * y(155) + 7 * y(163) * buff(9) / ( 20.1153 + buff(9) ) + 18 * y(164) * buff(9) / ( 20 + buff(9) ) + 18 * y(161) * buff(9) / ( 20 + buff(9) ) + 18 * y(166) * buff(9) / ( 20 + buff(9) ) + 54 * y(168) * buff(9) / ( 60 + buff(9) ) - 5.04 * y(160) * y(171) + 80 * y(172) - 2.52 * y(160) * y(169) + 40 * y(170) - 0.209997 * y(160) * y(173) + 6.67 * y(174) - 1.26 * y(160) * y(175) + 40 * y(176)
 + 49.9995 * y(93) * y(162) - 1 * y(161)
 - 49.9995 * y(93) * y(162) + 1 * y(161) - 126 * y(177) * y(162) + 1 * y(166)
 + 4 * y(19) - 0.1 * y(163) - 833.283 * y(177) * y(163) + 1 * y(168)
 + 499.995 * y(165) * y(177) - 1 * y(164)
 - 499.995 * y(165) * y(177) + 1 * y(164) - 0.600003 * y(165) *  (y(2) + y(4) + y(3) + y(5) + y(6) + y(7) ) + 16 * y(19) + 0.1 * y(163)
 + 126 * y(177) * y(162) - 1 * y(166)
 + 1.2e+07 * y(172) + 6e+06 * y(170) + 6e+06 * y(176) + 1.0002e+06 * y(174)
 + 833.283 * y(177) * y(163) - 1 * y(168)
 - 2.52 * y(160) * y(169) + 40 * y(170) + 10 * y(170) + 0.1 * y(171) - 6.00003 * y(169) * y(89) + 36 * y(92)
 + 2.52 * y(160) * y(169) - 40 * y(170) - 10 * y(170)
 - 5.04 * y(160) * y(171) + 80 * y(172) + 20 * y(172) - 0.1 * y(171) + 9 * y(92)
 + 5.04 * y(160) * y(171) - 80 * y(172) - 20 * y(172)
 - 0.209997 * y(160) * y(173) + 6.67 * y(174) + 1.667 * y(174) - 720 * y(173) * y(93) + 5 * y(175)
 + 0.209997 * y(160) * y(173) - 6.67 * y(174) - 1.667 * y(174)
 - 1.26 * y(160) * y(175) + 40 * y(176) + 10 * y(176) + 720 * y(173) * y(93) - 5 * y(175)
 + 1.26 * y(160) * y(175) - 40 * y(176) - 10 * y(176)
 - 499.995 * y(165) * y(177) + 1 * y(164) - 126 * y(177) * y(162) + 1 * y(166) - 833.283 * y(177) * y(163) + 1 * y(168)
 - 0.600003 * y(11) * y(178) + 0.5 * y(1) - 0.99999 * y(20) * y(178) + 0.1 * y(21) - 6.00003 * y(31) * y(178) + 0.1 * y(29) - 2.99997 * y(178) * y(37) + 1 * y(38) - 29.9997 * y(178) * y(42) + 1 * y(40) - 180 * y(108) * y(178) + 10 * y(110) - 11.9997 * y(178) * y(109) + 10 * y(111) - 144 * y(128) * y(178)*y(178) + 144 * y(133) - 3.6 * y(133) * y(178) + 10 * y(132) - 0.465003 * y(132) * y(178) + 10 * y(93) - 7.19993 * y(178)*y(178) * y(147) + 2 * y(144) - 20016.7 * y(146) * y(178)*y(178) + 2 * y(147) + 100 * y(179) - 100 * y(178)
 - 100 * y(179) + 100 * y(178)
];

% end of model
