 %  DOQCS : http://doqcs.ncbs.res.in/ 
 %  Accession Name = AMPAR_CaMKII_weak_coupling 
 %  Accession Number = 65 
 %  Transcriber = Upinder S. Bhalla, NCBS 
 %  Developer = Arnold Hayer and Upinder S. Bhalla, NCBS 
 %  Species = General Mammalian 
 %  Tissue = Neuronal 
 %  Cell Compartment = Synaptic Spine, Postsynaptic Density 
 %  Notes =  This is a model of weak coupling between the AMPAR traffikcing   bistability, and the CaMKII autophosphorylation bistability. In   this model, there are three stable states: Both off, AMPAR on,   or both on. The fourth possible state: CaMKII on but AMPAR off,   is not truly stable, since over the course of hours the AMPAR also   turns on.  
 
   % Note: All the concentrations are in micro Mole (uM)

function moose_convert
%Matlab simulation for chemical kinetics

% Saved on Thu Aug  9 17:15:57 2007
% Automatically generated from kkit model in MOOSE by conv_matlab

global buff;
global t;
global y;

% Define Tables here
tspan=[0, 2000];
% Molecule names
molnames = {
'1 /kinetics/CaM-Ca4'
'2 /kinetics/PP1-active'
'3 /kinetics/PP1-active/Deph-thr286/Deph-thr286_cplx'
'4 /kinetics/PP1-active/Deph-thr286b/Deph-thr286b_cplx'
'5 /kinetics/PP1-active/Deph-thr305/Deph-thr305_cplx'
'6 /kinetics/PP1-active/Deph-thr286c/Deph-thr286c_cplx'
'7 /kinetics/PP1-active/Deph-thr305a/Deph-thr305a_cplx'
'8 /kinetics/cAMP'
'9 /kinetics/Ca'
'10 /kinetics/CaMKII/CaMKII'
'11 /kinetics/CaMKII/CaMKII-CaM'
'12 /kinetics/CaMKII/CaMKII-thr286*-CaM'
'13 /kinetics/CaMKII/CaMKII***'
'14 /kinetics/CaMKII/CaMKII-thr286'
'15 /kinetics/CaMKII/tot_CaM_CaMKII/CaM_act_305/CaM_act_305_cplx'
'16 /kinetics/CaMKII/tot_CaM_CaMKII/CaM_act_286/CaM_act_286_cplx'
'17 /kinetics/CaMKII/tot_autonomous_CaMKII/auton_305/auton_305_cplx'
'18 /kinetics/CaMKII/tot_autonomous_CaMKII/auton_286/auton_286_cplx'
'19 /kinetics/CaMKII/CaMK-thr305'
'20 /kinetics/CaM/CaM'
'21 /kinetics/CaM/neurogranin-CaM'
'22 /kinetics/CaM/neurogranin*'
'23 /kinetics/CaM/neurogranin'
'24 /kinetics/CaM/CaM-PSD'
'25 /kinetics/CaM/neurogranin-CaM_PSD'
'26 /kinetics/CaM/neurogranin_PSD'
'27 /kinetics/CaM/neurogranin*_PSD'
'28 /kinetics/CaM-Ca3'
'29 /kinetics/CaM-TR2-Ca2'
'30 /kinetics/PP1/I1'
'31 /kinetics/PP1/I1*'
'32 /kinetics/PP1/PP1-I1*'
'33 /kinetics/PP1/PP1-I1'
'34 /kinetics/PP2A'
'35 /kinetics/PP2A/PP2A-dephosph-I1/PP2A-dephosph-I1_cplx'
'36 /kinetics/PP2A/PP2A-dephosph-PP1-I*/PP2A-dephosph-PP1-I*_cplx'
'37 /kinetics/PP2A/PP2A-dephosph-I1_PSD/PP2A-dephosph-I1_PSD_cplx'
'38 /kinetics/PP2A/PP2A-dephosph-PP1-I*_PSD/PP2A-dephosph-PP1-I*_PSD_cplx'
'39 /kinetics/PP2B/CaNAB'
'40 /kinetics/PP2B/CaNAB-Ca2'
'41 /kinetics/PP2B/CaM_Ca_n-CaNAB'
'42 /kinetics/PP2B/CaM_Ca_n-CaNAB/dephosph_neurogranin/dephosph_neurogranin_cplx'
'43 /kinetics/PP2B/CaM_Ca_n-CaNAB/dephosph_inhib1/dephosph_inhib1_cplx'
'44 /kinetics/PP2B/CaM_Ca_n-CaNAB/dephosph-PP1-I*/dephosph-PP1-I*_cplx'
'45 /kinetics/PP2B/CaM_Ca_n-CaNAB/dephosph_inhib1_PSD/dephosph_inhib1_PSD_cplx'
'46 /kinetics/PP2B/CaM_Ca_n-CaNAB/dephosph-PP1-I*_PSD/dephosph-PP1-I*_PSD_cplx'
'47 /kinetics/PP2B/CaM_Ca_n-CaNAB/PP2B[17]/PP2B[17]_cplx'
'48 /kinetics/PP2B/CaM_Ca_n-CaNAB/PP2B[18]/PP2B[18]_cplx'
'49 /kinetics/PP2B/CaM_Ca_n-CaNAB/PP2B[19]/PP2B[19]_cplx'
'50 /kinetics/PP2B/CaM_Ca_n-CaNAB/PP2B[20]/PP2B[20]_cplx'
'51 /kinetics/PP2B/CaM_Ca_n-CaNAB/PP2B[21]/PP2B[21]_cplx'
'52 /kinetics/PP2B/CaM_Ca_n-CaNAB/PP2B[22]/PP2B[22]_cplx'
'53 /kinetics/PP2B/CaM_Ca_n-CaNAB/PP2B[23]/PP2B[23]_cplx'
'54 /kinetics/PP2B/CaM_Ca_n-CaNAB/PP2B[24]/PP2B[24]_cplx'
'55 /kinetics/PP2B/CaM_Ca_n-CaNAB/PP2B[25]/PP2B[25]_cplx'
'56 /kinetics/PP2B/CaM_Ca_n-CaNAB/PP2B[26]/PP2B[26]_cplx'
'57 /kinetics/PP2B/CaM_Ca_n-CaNAB/PP2B[27]/PP2B[27]_cplx'
'58 /kinetics/PP2B/CaM_Ca_n-CaNAB/PP2B[28]/PP2B[28]_cplx'
'59 /kinetics/PP2B/CaM_Ca_n-CaNAB/PP2B[29]/PP2B[29]_cplx'
'60 /kinetics/PP2B/CaM_Ca_n-CaNAB/PP2B[30]/PP2B[30]_cplx'
'61 /kinetics/PP2B/CaM_Ca_n-CaNAB/PP2B[31]/PP2B[31]_cplx'
'62 /kinetics/PP2B/CaM_Ca_n-CaNAB/PP2B[32]/PP2B[32]_cplx'
'63 /kinetics/PP2B/CaM_Ca_n-CaNAB/PP2B/PP2B_cplx'
'64 /kinetics/PP2B/CaM_Ca_n-CaNAB/PP2B[1]/PP2B[1]_cplx'
'65 /kinetics/AC/AC1-CaM'
'66 /kinetics/PP2B/CaM_Ca_n-CaNAB/PP2B[2]/PP2B[2]_cplx'
'67 /kinetics/PP2B/CaM_Ca_n-CaNAB/PP2B[3]/PP2B[3]_cplx'
'68 /kinetics/PP2B/CaM_Ca_n-CaNAB/PP2B[4]/PP2B[4]_cplx'
'69 /kinetics/PP2B/CaM_Ca_n-CaNAB/PP2B[5]/PP2B[5]_cplx'
'70 /kinetics/PP2B/CaM_Ca_n-CaNAB/PP2B[6]/PP2B[6]_cplx'
'71 /kinetics/PP2B/CaM_Ca_n-CaNAB/PP2B[7]/PP2B[7]_cplx'
'72 /kinetics/PP2B/CaM_Ca_n-CaNAB/PP2B[8]/PP2B[8]_cplx'
'73 /kinetics/PP2B/CaM_Ca_n-CaNAB/PP2B[9]/PP2B[9]_cplx'
'74 /kinetics/PP2B/CaM_Ca_n-CaNAB/PP2B[10]/PP2B[10]_cplx'
'75 /kinetics/PP2B/CaM_Ca_n-CaNAB/PP2B[11]/PP2B[11]_cplx'
'76 /kinetics/PP2B/CaM_Ca_n-CaNAB/PP2B[12]/PP2B[12]_cplx'
'77 /kinetics/PP2B/CaM_Ca_n-CaNAB/PP2B[13]/PP2B[13]_cplx'
'78 /kinetics/PP2B/CaM_Ca_n-CaNAB/PP2B[14]/PP2B[14]_cplx'
'79 /kinetics/PP2B/CaM_Ca_n-CaNAB/PP2B[15]/PP2B[15]_cplx'
'80 /kinetics/PP2B/CaM_Ca_n-CaNAB/dephosph_neurogranin_PSD/dephosph_neurogranin_PSD_cplx'
'81 /kinetics/PP2B/CaNAB-Ca4'
'82 /kinetics/PP2B/CaNAB-Ca4/dephosph_inhib1_noCaM/dephosph_inhib1_noCaM_cplx'
'83 /kinetics/PP2B/CaNAB-Ca4/dephosph_inhib1_noCaM_PSD/dephosph_inhib1_noCaM_PSD_cplx'
'84 /kinetics/PKA/R2C2'
'85 /kinetics/PKA/R2C2-cAMP'
'86 /kinetics/PKA/R2C2-cAMP2'
'87 /kinetics/PKA/R2C2-cAMP3'
'88 /kinetics/PKA/R2C2-cAMP4'
'89 /kinetics/PKA/R2C-cAMP4'
'90 /kinetics/PKA/R2-cAMP4'
'91 /kinetics/PKA/PKA-inhibitor'
'92 /kinetics/PKA/inhibited-PKA'
'93 /kinetics/PKA/PKA-active'
'94 /kinetics/PKA/PKA-active/phosph-PDE/phosph-PDE_cplx'
'95 /kinetics/PKA/PKA-active/PKA-phosph-I1/PKA-phosph-I1_cplx'
'96 /kinetics/PKA/PKA-active/PKA-phosph-I1_PSD/PKA-phosph-I1_PSD_cplx'
'97 /kinetics/PKA/PKA-active/PKA[17]/PKA[17]_cplx'
'98 /kinetics/PKA/PKA-active/PKA[16]/PKA[16]_cplx'
'99 /kinetics/PKA/PKA-active/PKA[18]/PKA[18]_cplx'
'100 /kinetics/PKA/PKA-active/PKA[19]/PKA[19]_cplx'
'101 /kinetics/PKA/PKA-active/PKA[20]/PKA[20]_cplx'
'102 /kinetics/PKA/PKA-active/PKA[21]/PKA[21]_cplx'
'103 /kinetics/PKA/PKA-active/PKA[22]/PKA[22]_cplx'
'104 /kinetics/PKA/PKA-active/PKA[23]/PKA[23]_cplx'
'105 /kinetics/PKA/PKA-active/PKA[24]/PKA[24]_cplx'
'106 /kinetics/PKA/PKA-active/PKA[25]/PKA[25]_cplx'
'107 /kinetics/PKA/PKA-active/PKA[26]/PKA[26]_cplx'
'108 /kinetics/PKA/PKA-active/PKA[27]/PKA[27]_cplx'
'109 /kinetics/PKA/PKA-active/PKA[28]/PKA[28]_cplx'
'110 /kinetics/PKA/PKA-active/PKA[29]/PKA[29]_cplx'
'111 /kinetics/PKA/PKA-active/PKA[30]/PKA[30]_cplx'
'112 /kinetics/PKA/PKA-active/PKA[31]/PKA[31]_cplx'
'113 /kinetics/PKA/PKA-active/PKA/PKA_cplx'
'114 /kinetics/PKA/PKA-active/PKA[1]/PKA[1]_cplx'
'115 /kinetics/PKA/PKA-active/PKA[2]/PKA[2]_cplx'
'116 /kinetics/PKA/PKA-active/PKA[3]/PKA[3]_cplx'
'117 /kinetics/PKA/PKA-active/PKA[4]/PKA[4]_cplx'
'118 /kinetics/PKA/PKA-active/PKA[5]/PKA[5]_cplx'
'119 /kinetics/PKA/PKA-active/PKA[6]/PKA[6]_cplx'
'120 /kinetics/PKA/PKA-active/PKA[7]/PKA[7]_cplx'
'121 /kinetics/PKA/PKA-active/PKA[8]/PKA[8]_cplx'
'122 /kinetics/PKA/PKA-active/PKA[9]/PKA[9]_cplx'
'123 /kinetics/PKA/PKA-active/PKA[10]/PKA[10]_cplx'
'124 /kinetics/PKA/PKA-active/PKA[11]/PKA[11]_cplx'
'125 /kinetics/PKA/PKA-active/PKA[12]/PKA[12]_cplx'
'126 /kinetics/PKA/PKA-active/PKA[13]/PKA[13]_cplx'
'127 /kinetics/PKA/PKA-active/PKA[14]/PKA[14]_cplx'
'128 /kinetics/PKA/PKA-active/PKA[15]/PKA[15]_cplx'
'129 /kinetics/AC/AC1'
'130 /kinetics/AC/AC2*'
'131 /kinetics/AC/AC2'
'132 /kinetics/AC/AMP'
'133 /kinetics/AC/cAMP-PDE'
'134 /kinetics/AC/cAMP-PDE/PDE/PDE_cplx'
'135 /kinetics/AC/cAMP-PDE*'
'136 /kinetics/AC/cAMP-PDE*/PDE*/PDE*_cplx'
'137 /kinetics/AC/PDE1'
'138 /kinetics/AC/PDE1/PDE1/PDE1_cplx'
'139 /kinetics/AC/CaM.PDE1'
'140 /kinetics/AC/CaM.PDE1/CaM.PDE1/CaM.PDE1_cplx'
'141 /kinetics/AC/cAMP_in_dend'
'142 /kinetics/CaM-TR2-Ca2-PSD'
'143 /kinetics/CaM-Ca3-PSD'
'144 /kinetics/CaM-Ca4-PSD'
'145 /kinetics/Ca-PSD'
'146 /kinetics/PP1_PSD/I1'
'147 /kinetics/PP1_PSD/I1*'
'148 /kinetics/PP1_PSD/PP1-I1*'
'149 /kinetics/PP1_PSD/PP1-I1'
'150 /kinetics/PP1_PSD/PP1-active_PSD'
'151 /kinetics/PP1_PSD/PP1-active_PSD/PP1/PP1_cplx'
'152 /kinetics/PP1_PSD/PP1-active_PSD/PP1[1]/PP1[1]_cplx'
'153 /kinetics/PP1_PSD/PP1-active_PSD/PP1[2]/PP1[2]_cplx'
'154 /kinetics/PP1_PSD/PP1-active_PSD/PP1[3]/PP1[3]_cplx'
'155 /kinetics/PP1_PSD/PP1-active_PSD/PP1[4]/PP1[4]_cplx'
'156 /kinetics/PP1_PSD/PP1-active_PSD/PP1[5]/PP1[5]_cplx'
'157 /kinetics/PP1_PSD/PP1-active_PSD/PP1[6]/PP1[6]_cplx'
'158 /kinetics/PP1_PSD/PP1-active_PSD/PP1[7]/PP1[7]_cplx'
'159 /kinetics/PP1_PSD/PP1-active_PSD/PP1[8]/PP1[8]_cplx'
'160 /kinetics/PP1_PSD/PP1-active_PSD/PP1[9]/PP1[9]_cplx'
'161 /kinetics/PP1_PSD/PP1-active_PSD/PP1[10]/PP1[10]_cplx'
'162 /kinetics/PP1_PSD/PP1-active_PSD/PP1[11]/PP1[11]_cplx'
'163 /kinetics/PP1_PSD/PP1-active_PSD/PP1[12]/PP1[12]_cplx'
'164 /kinetics/PP1_PSD/PP1-active_PSD/PP1[13]/PP1[13]_cplx'
'165 /kinetics/PP1_PSD/PP1-active_PSD/PP1[14]/PP1[14]_cplx'
'166 /kinetics/PP1_PSD/PP1-active_PSD/PP1[15]/PP1[15]_cplx'
'167 /kinetics/PP1_PSD/PP1-active_PSD/PP1[16]/PP1[16]_cplx'
'168 /kinetics/PP1_PSD/PP1-active_PSD/PP1[17]/PP1[17]_cplx'
'169 /kinetics/PP1_PSD/PP1-active_PSD/PP1[18]/PP1[18]_cplx'
'170 /kinetics/PP1_PSD/PP1-active_PSD/PP1[19]/PP1[19]_cplx'
'171 /kinetics/PP1_PSD/PP1-active_PSD/PP1[20]/PP1[20]_cplx'
'172 /kinetics/PP1_PSD/PP1-active_PSD/PP1[21]/PP1[21]_cplx'
'173 /kinetics/PP1_PSD/PP1-active_PSD/PP1[22]/PP1[22]_cplx'
'174 /kinetics/PP1_PSD/PP1-active_PSD/PP1[23]/PP1[23]_cplx'
'175 /kinetics/PP1_PSD/PP1-active_PSD/PP1[24]/PP1[24]_cplx'
'176 /kinetics/PP1_PSD/PP1-active_PSD/PP1[25]/PP1[25]_cplx'
'177 /kinetics/PP1_PSD/PP1-active_PSD/PP1[26]/PP1[26]_cplx'
'178 /kinetics/PP1_PSD/PP1-active_PSD/PP1[27]/PP1[27]_cplx'
'179 /kinetics/PP1_PSD/PP1-active_PSD/PP1[28]/PP1[28]_cplx'
'180 /kinetics/PP1_PSD/PP1-active_PSD/PP1[29]/PP1[29]_cplx'
'181 /kinetics/PP1_PSD/PP1-active_PSD/PP1[30]/PP1[30]_cplx'
'182 /kinetics/PP1_PSD/PP1-active_PSD/PP1[31]/PP1[31]_cplx'
'183 /kinetics/PKC-active'
'184 /kinetics/PKC-active/PKC-phosph-neurogranin/PKC-phosph-neurogranin_cplx'
'185 /kinetics/PKC-active/PKC-phosph-ng-CaM/PKC-phosph-ng-CaM_cplx'
'186 /kinetics/PKC-active/phosph-AC2/phosph-AC2_cplx'
'187 /kinetics/PKC-active/PKC-phosph-neurogranin_PSD/PKC-phosph-neurogranin_PSD_cplx'
'188 /kinetics/PKC-active/PKC-phosph-ng-CaM_PSD/PKC-phosph-ng-CaM_PSD_cplx'
'189 /kinetics/AMPAR/GluR23_M'
'190 /kinetics/AMPAR/GluR23_I'
'191 /kinetics/AMPAR/A_B'
'192 /kinetics/AMPAR/A831*_B'
'193 /kinetics/AMPAR/A845*_B'
'194 /kinetics/AMPAR/A831*845*_B'
'195 /kinetics/AMPAR/A_B845*'
'196 /kinetics/AMPAR/A845*_B845*'
'197 /kinetics/AMPAR/A831*845*_B845*'
'198 /kinetics/AMPAR/A845*_B831*845*'
'199 /kinetics/AMPAR/A835*845*_B835*845*'
'200 /kinetics/AMPAR/A_B831*845*'
'201 /kinetics/AMPAR/A831*_B831*845*'
'202 /kinetics/AMPAR/A_B831*'
'203 /kinetics/AMPAR/A845*_B831*'
'204 /kinetics/AMPAR/A831*_B831*'
'205 /kinetics/AMPAR/A831*845*_B831*'
'206 /kinetics/AMPAR/A831*_B845*'
'207 /kinetics/Anchor'
'208 /kinetics/AMPAR_memb/A_B'
'209 /kinetics/AMPAR_memb/A831*_B'
'210 /kinetics/AMPAR_memb/A_B831*'
'211 /kinetics/AMPAR_memb/A831*_B831*'
'212 /kinetics/AMPAR_memb/A845*_B'
'213 /kinetics/AMPAR_memb/A831*845*_B'
'214 /kinetics/AMPAR_memb/A831*845*_B831*'
'215 /kinetics/AMPAR_memb/A_B845*'
'216 /kinetics/AMPAR_memb/A845*_B831*'
'217 /kinetics/AMPAR_memb/A845*_B845*'
'218 /kinetics/AMPAR_memb/A831*_B845*'
'219 /kinetics/AMPAR_memb/A_B831*845*'
'220 /kinetics/AMPAR_memb/A831*_B831*845*'
'221 /kinetics/AMPAR_memb/A831*845*_B845*'
'222 /kinetics/AMPAR_memb/A845*_B831*845*'
'223 /kinetics/AMPAR_memb/A835*845*_B835*845*'
'224 /kinetics/PP1_CaMKII_PSD/PP1-active_CaMKII_PSD'
'225 /kinetics/PP1_CaMKII_PSD/PP1-active_CaMKII_PSD/Deph-thr286/Deph-thr286_cplx'
'226 /kinetics/PP1_CaMKII_PSD/PP1-active_CaMKII_PSD/Deph-thr286b/Deph-thr286b_cplx'
'227 /kinetics/PP1_CaMKII_PSD/PP1-active_CaMKII_PSD/Deph-thr305/Deph-thr305_cplx'
'228 /kinetics/PP1_CaMKII_PSD/PP1-active_CaMKII_PSD/Deph-thr286c/Deph-thr286c_cplx'
'229 /kinetics/PP1_CaMKII_PSD/PP1-active_CaMKII_PSD/Deph-thr305a/Deph-thr305a_cplx'
'230 /kinetics/CaMKII_PSD/CaMKII-thr305-PSD'
'231 /kinetics/CaMKII_PSD/CaMKII***-PSD'
'232 /kinetics/CaMKII_PSD/CaMKII-PSD'
'233 /kinetics/CaMKII_PSD/CaMKII-thr286-PSD'
'234 /kinetics/CaMKII_PSD/CaMKII-CaM-PSD'
'235 /kinetics/CaMKII_PSD/CaMKII-thr286-CaM-PSD'
'236 /kinetics/CaMKII_PSD/tot-CaM-CaMKII-PSD/CaM-act-305-PSD/CaM-act-305-PSD_cplx'
'237 /kinetics/CaMKII_PSD/tot-CaM-CaMKII-PSD/CaM-act-286-PSD/CaM-act-286-PSD_cplx'
'238 /kinetics/CaMKII_PSD/tot-CaM-CaMKII-PSD/CaM_act_autoph/CaM_act_autoph_cplx'
'239 /kinetics/CaMKII_PSD/tot-auto-PSD/auton-305-PSD/auton-305-PSD_cplx'
'240 /kinetics/CaMKII_PSD/tot-auto-PSD/auton-286-PSD/auton-286-PSD_cplx'
'241 /kinetics/CaMKII_PSD/tot-auto-PSD/auton-autoph/auton-autoph_cplx'
'242 /kinetics/CaMKII_PSD/basal_CaMKII_PSD'
'243 /kinetics/CaMKII_PSD/actCaMKII-PSD/CaMKII/CaMKII_cplx'
'244 /kinetics/CaMKII_PSD/actCaMKII-PSD/CaMKII[1]/CaMKII[1]_cplx'
'245 /kinetics/CaMKII_PSD/actCaMKII-PSD/CaMKII[2]/CaMKII[2]_cplx'
'246 /kinetics/CaMKII_PSD/actCaMKII-PSD/CaMKII[3]/CaMKII[3]_cplx'
'247 /kinetics/CaMKII_PSD/actCaMKII-PSD/CaMKII[4]/CaMKII[4]_cplx'
'248 /kinetics/CaMKII_PSD/actCaMKII-PSD/CaMKII[5]/CaMKII[5]_cplx'
'249 /kinetics/CaMKII_PSD/actCaMKII-PSD/CaMKII[6]/CaMKII[6]_cplx'
'250 /kinetics/CaMKII_PSD/actCaMKII-PSD/CaMKII[7]/CaMKII[7]_cplx'
'251 /kinetics/CaMKII_PSD/actCaMKII-PSD/CaMKII[8]/CaMKII[8]_cplx'
'252 /kinetics/CaMKII_PSD/actCaMKII-PSD/CaMKII[9]/CaMKII[9]_cplx'
'253 /kinetics/CaMKII_PSD/actCaMKII-PSD/CaMKII[10]/CaMKII[10]_cplx'
'254 /kinetics/CaMKII_PSD/actCaMKII-PSD/CaMKII[11]/CaMKII[11]_cplx'
'255 /kinetics/CaMKII_PSD/actCaMKII-PSD/CaMKII[12]/CaMKII[12]_cplx'
'256 /kinetics/CaMKII_PSD/actCaMKII-PSD/CaMKII[13]/CaMKII[13]_cplx'
'257 /kinetics/CaMKII_PSD/actCaMKII-PSD/CaMKII[14]/CaMKII[14]_cplx'
'258 /kinetics/CaMKII_PSD/actCaMKII-PSD/CaMKII[15]/CaMKII[15]_cplx'
'259 /kinetics/CaMKII_PSD/NMDAR'
};

bufnames = {
'1 /kinetics/CaMKII/basal_CaMKII_cyt'
'2 /kinetics/AC/ATP'
'3 /kinetics/temp-PIP2'
'4 /kinetics/AMPAR/AMPAR_deg'
'5 /kinetics/PKC-control'
'6 /kinetics/Ca_control_cyt'
'7 /kinetics/Ca_control_PSD'
'8 /kinetics/AMPAR_bulk'
'9 /kinetics/CaMKII_PSD/basal_CaMKII_PSD_control'
};

sumtotnames = {
'1 /kinetics/CaMKII/tot_CaM_CaMKII'
'2 /kinetics/CaMKII/tot_autonomous_CaMKII'
'3 /kinetics/CaMKII/tot_CaMKII_cyt'
'4 /kinetics/CaMKII/act_CaMKII_cyt'
'5 /kinetics/AMPAR/total_mem'
'6 /kinetics/I_845-P'
'7 /kinetics/tot_I_GluR12'
'8 /kinetics/total_Int'
'9 /kinetics/Ser845'
'10 /kinetics/Ser845-P'
'11 /kinetics/Ser845-PP'
'12 /kinetics/Ser831'
'13 /kinetics/Ser831-P'
'14 /kinetics/Ser831-PP'
'15 /kinetics/tot_mem_GluR12'
'16 /kinetics/I_845'
'17 /kinetics/I_845_PP'
'18 /kinetics/CaMKII_PSD/tot-CaM-CaMKII-PSD'
'19 /kinetics/CaMKII_PSD/tot-auto-PSD'
'20 /kinetics/CaMKII_PSD/tot_CaMKII_PSD'
'21 /kinetics/CaMKII_PSD/actCaMKII-PSD'
'22 /kinetics/CaMKII_PSD/286P-PSD'
};

% Initial concs of 259 molecules
y0(1) = 0;
y0(2) = 1.8;
y0(3) = 0;
y0(4) = 0;
y0(5) = 0;
y0(6) = 0;
y0(7) = 0;
y0(8) = 0;
y0(9) = 0.08;
y0(10) = 20;
y0(11) = 0;
y0(12) = 0;
y0(13) = 0;
y0(14) = 0;
y0(15) = 6.3263e-11;
y0(16) = 2.44537e-08;
y0(17) = 4.90056e-08;
y0(18) = 1.89444e-05;
y0(19) = 0;
y0(20) = 26.3333;
y0(21) = 0;
y0(22) = 0;
y0(23) = 10;
y0(24) = 26.3333;
y0(25) = 0;
y0(26) = 10;
y0(27) = 0;
y0(28) = 0;
y0(29) = 0;
y0(30) = 1.8;
y0(31) = 0;
y0(32) = 0;
y0(33) = 0;
y0(34) = 0.111111;
y0(35) = 1.99333e-06;
y0(36) = 0.00691;
y0(37) = 2.21556e-07;
y0(38) = 0.00127907;
y0(39) = 1;
y0(40) = 0;
y0(41) = 0;
y0(42) = 0;
y0(43) = 0;
y0(44) = 0;
y0(45) = 0;
y0(46) = 0;
y0(47) = 0;
y0(48) = 0;
y0(49) = 0;
y0(50) = 0;
y0(51) = 0;
y0(52) = 0;
y0(53) = 0;
y0(54) = 0;
y0(55) = 0;
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;
y0(69) = 0;
y0(70) = 0;
y0(71) = 0;
y0(72) = 0;
y0(73) = 0;
y0(74) = 0;
y0(75) = 0;
y0(76) = 0;
y0(77) = 0;
y0(78) = 0;
y0(79) = 0;
y0(80) = 0;
y0(81) = 0;
y0(82) = 0;
y0(83) = 0;
y0(84) = 0.5;
y0(85) = 0;
y0(86) = 0;
y0(87) = 0;
y0(88) = 0;
y0(89) = 0;
y0(90) = 0;
y0(91) = 0.259259;
y0(92) = 0;
y0(93) = 0;
y0(94) = 0;
y0(95) = 0;
y0(96) = 0;
y0(97) = 0;
y0(98) = 0;
y0(99) = 0;
y0(100) = 0;
y0(101) = 0;
y0(102) = 0;
y0(103) = 0;
y0(104) = 0;
y0(105) = 0;
y0(106) = 0;
y0(107) = 0;
y0(108) = 0;
y0(109) = 0;
y0(110) = 0;
y0(111) = 0;
y0(112) = 0;
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) = 0;
y0(123) = 0;
y0(124) = 0;
y0(125) = 0;
y0(126) = 0;
y0(127) = 0;
y0(128) = 0;
y0(129) = 0.0740741;
y0(130) = 0;
y0(131) = 0.0740741;
y0(132) = 0;
y0(133) = 0.555556;
y0(134) = 0;
y0(135) = 0;
y0(136) = 0;
y0(137) = 2.59259;
y0(138) = 0;
y0(139) = 0;
y0(140) = 0;
y0(141) = 0;
y0(142) = 0;
y0(143) = 0;
y0(144) = 0;
y0(145) = 0.08;
y0(146) = 4;
y0(147) = 0;
y0(148) = 0;
y0(149) = 0;
y0(150) = 4;
y0(151) = 0;
y0(152) = 0;
y0(153) = 0;
y0(154) = 0;
y0(155) = 0;
y0(156) = 0;
y0(157) = 0;
y0(158) = 0;
y0(159) = 0;
y0(160) = 0;
y0(161) = 0;
y0(162) = 0;
y0(163) = 0;
y0(164) = 0;
y0(165) = 0;
y0(166) = 0;
y0(167) = 0;
y0(168) = 0;
y0(169) = 0;
y0(170) = 0;
y0(171) = 0;
y0(172) = 0;
y0(173) = 0;
y0(174) = 0;
y0(175) = 0;
y0(176) = 0;
y0(177) = 0;
y0(178) = 0;
y0(179) = 0;
y0(180) = 0;
y0(181) = 0;
y0(182) = 0;
y0(183) = 0.1;
y0(184) = 0.00538593;
y0(185) = 0;
y0(186) = 4.00704e-05;
y0(187) = 0;
y0(188) = 0;
y0(189) = 3.5;
y0(190) = 0.0925926;
y0(191) = 0;
y0(192) = 0;
y0(193) = 0;
y0(194) = 0;
y0(195) = 0;
y0(196) = 0;
y0(197) = 0;
y0(198) = 0;
y0(199) = 0;
y0(200) = 0;
y0(201) = 0;
y0(202) = 0;
y0(203) = 0;
y0(204) = 0;
y0(205) = 0;
y0(206) = 0;
y0(207) = 27.3333;
y0(208) = 0;
y0(209) = 0;
y0(210) = 0;
y0(211) = 0;
y0(212) = 0;
y0(213) = 0;
y0(214) = 0;
y0(215) = 0;
y0(216) = 0;
y0(217) = 0;
y0(218) = 0;
y0(219) = 0;
y0(220) = 0;
y0(221) = 0;
y0(222) = 0;
y0(223) = 0;
y0(224) = 2;
y0(225) = 0;
y0(226) = 0;
y0(227) = 0;
y0(228) = 0;
y0(229) = 0;
y0(230) = 0;
y0(231) = 0;
y0(232) = 0;
y0(233) = 0;
y0(234) = 0;
y0(235) = 0;
y0(236) = 0;
y0(237) = 0;
y0(238) = 0;
y0(239) = 0;
y0(240) = 0;
y0(241) = 0;
y0(242) = 2;
y0(243) = 7.1375e-215;
y0(244) = 3.588e-215;
y0(245) = 3.588e-215;
y0(246) = 1.80017e-215;
y0(247) = 7.1375e-215;
y0(248) = 3.588e-215;
y0(249) = 3.588e-215;
y0(250) = 1.80017e-215;
y0(251) = 4.06667e-215;
y0(252) = 2.04433e-215;
y0(253) = 2.04433e-215;
y0(254) = 1.02567e-215;
y0(255) = 4.06667e-215;
y0(256) = 2.04433e-215;
y0(257) = 2.04433e-215;
y0(258) = 1.02567e-215;
y0(259) = 120;
buff(1) = 2;
buff(2) = 2000;
buff(3) = 2.5;
buff(4) = 0;
buff(5) = 0.1;
buff(6) = 0.08;
buff(7) = 0.08;
buff(8) = 0.00925933;
buff(9) = 2;
% Sumtotal eqn for /kinetics/CaMKII/tot_CaM_CaMKII: sumtot(1) =  (y(11) + y(12) )
% Sumtotal eqn for /kinetics/CaMKII/tot_autonomous_CaMKII: sumtot(2) =  (y(14) + y(13) + 1 * buff(1) )
% Sumtotal eqn for /kinetics/CaMKII/tot_CaMKII_cyt: sumtot(3) =  (y(11) + y(12) + y(14) + y(13) + y(19) + y(10) + 1 * buff(1) )
% Sumtotal eqn for /kinetics/CaMKII/act_CaMKII_cyt: sumtot(4) =  (1 *  (y(11) + y(12) ) + 1 *  (y(14) + y(13) + 1 * buff(1) ) )
% Sumtotal eqn for /kinetics/AMPAR/total_mem: sumtot(5) =  (y(189) + 1 *  (1 *  (y(208) + y(209) + y(210) + y(211) ) + 1 *  (y(212) + y(213) + y(216) + y(214) + y(215) + y(218) + y(219) + y(220) ) + 1 *  (y(217) + y(221) + y(222) + y(223) ) ) )
% Sumtotal eqn for /kinetics/I_845-P: sumtot(6) =  (y(205) + y(203) + y(194) + y(193) + y(201) + y(200) + y(206) + y(195) )
% Sumtotal eqn for /kinetics/tot_I_GluR12: sumtot(7) =  (1 *  (y(204) + y(202) + y(192) + y(191) ) + 1 *  (y(205) + y(203) + y(194) + y(193) + y(201) + y(200) + y(206) + y(195) ) + 1 *  (y(199) + y(198) + y(197) + y(196) ) )
% Sumtotal eqn for /kinetics/total_Int: sumtot(8) =  (1 *  (1 *  (y(204) + y(202) + y(192) + y(191) ) + 1 *  (y(205) + y(203) + y(194) + y(193) + y(201) + y(200) + y(206) + y(195) ) + 1 *  (y(199) + y(198) + y(197) + y(196) ) ) + y(190) )
% Sumtotal eqn for /kinetics/Ser845: sumtot(9) =  (y(208) + y(209) + y(210) + y(211) )
% Sumtotal eqn for /kinetics/Ser845-P: sumtot(10) =  (y(212) + y(213) + y(216) + y(214) + y(215) + y(218) + y(219) + y(220) )
% Sumtotal eqn for /kinetics/Ser845-PP: sumtot(11) =  (y(217) + y(221) + y(222) + y(223) )
% Sumtotal eqn for /kinetics/Ser831: sumtot(12) =  (y(217) + y(215) + y(212) + y(208) )
% Sumtotal eqn for /kinetics/Ser831-P: sumtot(13) =  (y(221) + y(222) + y(218) + y(219) + y(213) + y(216) + y(209) + y(210) )
% Sumtotal eqn for /kinetics/Ser831-PP: sumtot(14) =  (y(223) + y(220) + y(214) + y(211) )
% Sumtotal eqn for /kinetics/tot_mem_GluR12: sumtot(15) =  (1 *  (y(208) + y(209) + y(210) + y(211) ) + 1 *  (y(212) + y(213) + y(216) + y(214) + y(215) + y(218) + y(219) + y(220) ) + 1 *  (y(217) + y(221) + y(222) + y(223) ) )
% Sumtotal eqn for /kinetics/I_845: sumtot(16) =  (y(204) + y(202) + y(192) + y(191) )
% Sumtotal eqn for /kinetics/I_845_PP: sumtot(17) =  (y(199) + y(198) + y(197) + y(196) )
% Sumtotal eqn for /kinetics/CaMKII_PSD/tot-CaM-CaMKII-PSD: sumtot(18) =  (y(235) + y(234) )
% Sumtotal eqn for /kinetics/CaMKII_PSD/tot-auto-PSD: sumtot(19) =  (y(233) + y(231) + y(242) )
% Sumtotal eqn for /kinetics/CaMKII_PSD/tot_CaMKII_PSD: sumtot(20) =  (1 *  (y(235) + y(234) + 1 *  (y(233) + y(231) + y(242) ) ) + y(232) + y(230) )
% Sumtotal eqn for /kinetics/CaMKII_PSD/actCaMKII-PSD: sumtot(21) =  (y(235) + y(234) + 1 *  (y(233) + y(231) + y(242) ) )
% Sumtotal eqn for /kinetics/CaMKII_PSD/286P-PSD: sumtot(22) =  (y(235) + y(233) )

% 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, (1 *  (y(:,208) + y(:,209) + y(:,210) + y(:,211) ) + 1 *  (y(:,212) + y(:,213) + y(:,216) + y(:,214) + y(:,215) + y(:,218) + y(:,219) + y(:,220) ) + 1 *  (y(:,217) + y(:,221) + y(:,222) + y(:,223) ) ));
plot(t, (1 *  (y(:,204) + y(:,202) + y(:,192) + y(:,191) ) + 1 *  (y(:,205) + y(:,203) + y(:,194) + y(:,193) + y(:,201) + y(:,200) + y(:,206) + y(:,195) ) + 1 *  (y(:,199) + y(:,198) + y(:,197) + y(:,196) ) ));
plot(t, (y(:,221) + y(:,222) + y(:,218) + y(:,219) + y(:,213) + y(:,216) + y(:,209) + y(:,210) ));
plot(t, (y(:,217) + y(:,215) + y(:,212) + y(:,208) ));
plot(t, (y(:,205) + y(:,203) + y(:,194) + y(:,193) + y(:,201) + y(:,200) + y(:,206) + y(:,195) ));
plot(t, (y(:,204) + y(:,202) + y(:,192) + y(:,191) ));
plot(t, (y(:,199) + y(:,198) + y(:,197) + y(:,196) ));
plot(t, (y(:,223) + y(:,220) + y(:,214) + y(:,211) ));
plot(t,y(:,93));
plot(t,y(:,8));
plot(t,y(:,9));
plot(t,y(:,2));
plot(t,y(:,183));
plot(t,y(:,1));
plot(t,y(:,224));
plot(t,y(:,207));
plot(t,y(:,41));
plot(t,y(:,144));
plot(t,y(:,150));
plot(t, (y(:,208) + y(:,209) + y(:,210) + y(:,211) ));
plot(t, (y(:,235) + y(:,234) + 1 *  (y(:,233) + y(:,231) + y(:,242) ) ));
plot(t, (1 *  (y(:,235) + y(:,234) + 1 *  (y(:,233) + y(:,231) + y(:,242) ) ) + y(:,232) + y(:,230) ));
plot(t, (1 *  (y(:,11) + y(:,12) ) + 1 *  (y(:,14) + y(:,13) + 1 * buff(1) ) ));
plot(t, (y(:,11) + y(:,12) + y(:,14) + y(:,13) + y(:,19) + y(:,10) + 1 * buff(:,1) ));
plot(t, (y(:,217) + y(:,221) + y(:,222) + y(:,223) ));
plot(t, (y(:,212) + y(:,213) + y(:,216) + y(:,214) + y(:,215) + y(:,218) + y(:,219) + y(:,220) ));


% Evaluation function
function dydt = f(t, y)
global buff;
dydt = [
 - 49.9997 * y(1) * y(10) + 5 * y(11) - 1000.19 * y(14) * y(1) + 0.1 * y(12) + 0.464999 * y(28) * y(9) - 10 * y(1) - 599.994 * y(1) * y(81) + 1 * y(41) - 49.9997 * y(1) * y(129) + 1 * y(65) - 719.982 * y(137) * y(1) + 5 * y(139) + 60 * y(144) - 60 * y(1)
 - 499.981 * y(31) * y(2) + 0.1 * y(32) + 1 * y(33) - 2.45144 * y(12) * y(2) + 10 * y(3) + 2.5 * y(3) - 2.45144 * y(14) * y(2) + 10 * y(4) + 2.5 * y(4) - 2.45144 * y(13) * y(2) + 10 * y(5) + 2.5 * y(5) - 2.45144 * y(13) * y(2) + 10 * y(6) + 2.5 * y(6) - 2.45144 * y(19) * y(2) + 10 * y(7) + 2.5 * y(7)
 + 2.45144 * y(12) * y(2) - 10 * y(3) - 2.5 * y(3)
 + 2.45144 * y(14) * y(2) - 10 * y(4) - 2.5 * y(4)
 + 2.45144 * y(13) * y(2) - 10 * y(5) - 2.5 * y(5)
 + 2.45144 * y(13) * y(2) - 10 * y(6) - 2.5 * y(6)
 + 2.45144 * y(19) * y(2) - 10 * y(7) - 2.5 * y(7)
 - 54 * y(84) * y(8) + 33 * y(85) - 54 * y(85) * y(8) + 33 * y(86) - 75.0006 * y(86) * y(8) + 110 * y(87) - 75.0006 * y(8) * y(87) + 32.5 * y(88) + 2 * y(130) * buff(2) / ( 300.002 + buff(2) ) + 4.5 * y(65) * buff(2) / ( 299.998 + buff(2) ) - 5.03998 * y(8) * y(135) + 80 * y(136) - 2.52002 * y(8) * y(133) + 40 * y(134) - 0.210001 * y(8) * y(137) + 6.67 * y(138) - 1.25998 * y(8) * y(139) + 40 * y(140) - 300 * y(8) + 300 * y(141)
 - 143.998 * y(20) * y(9)*y(9) + 144 * y(29) - 3.60002 * y(29) * y(9) + 10 * y(28) - 0.464999 * y(28) * y(9) + 10 * y(1) - 7.20019 * y(9)*y(9) * y(40) + 2 * y(81) - 20016 * y(39) * y(9)*y(9) + 2 * y(40) + 100 * buff(6) - 100 * y(9)
 - 49.9997 * y(1) * y(10) + 5 * y(11) + 0.0333333 * y(232) - 6e-05 * y(10) * y(259) + 2.5 * y(4) + 2.5 * y(7)
 + 49.9997 * y(1) * y(10) - 5 * y(11) - 0.0110133 * y(11) *  (y(11) + y(12) ) + 2 * y(16) - 0.00714312 * y(11) *  (y(14) + y(13) + 1 * buff(1) ) + 2 * y(18) - 0.00012 * y(11) * y(259) + 2.5 * y(3)
 + 1000.19 * y(14) * y(1) - 0.1 * y(12) + 0.5 * y(16) + 0.5 * y(18) - 0.00012 * y(12) * y(259) - 2.45144 * y(12) * y(2) + 10 * y(3)
 + 6 * y(15) + 6 * y(17) - 2.45144 * y(13) * y(2) + 10 * y(5) - 2.45144 * y(13) * y(2) + 10 * y(6)
 - 1000.19 * y(14) * y(1) + 0.1 * y(12) - 0.13216 * y(14) *  (y(11) + y(12) ) + 24 * y(15) - 0.0857142 * y(14) *  (y(14) + y(13) + 1 * buff(1) ) + 24 * y(17) - 2.45144 * y(14) * y(2) + 10 * y(4) + 2.5 * y(5)
 + 0.13216 * y(14) *  (y(11) + y(12) ) - 24 * y(15) - 6 * y(15)
 + 0.0110133 * y(11) *  (y(11) + y(12) ) - 2 * y(16) - 0.5 * y(16)
 + 0.0857142 * y(14) *  (y(14) + y(13) + 1 * buff(1) ) - 24 * y(17) - 6 * y(17)
 + 0.00714312 * y(11) *  (y(14) + y(13) + 1 * buff(1) ) - 2 * y(18) - 0.5 * y(18)
 + 0.0333333 * y(230) - 6e-05 * y(19) * y(259) - 2.45144 * y(19) * y(2) + 10 * y(7) + 2.5 * y(6)
 - 71.999 * y(20) * y(9)*y(9) + 72 * y(29) - 0.300002 * y(23) * y(20) + 1 * y(21) + 0.35 * y(185)
 + 0.300002 * y(23) * y(20) - 1 * y(21) - 0.0611982 * y(21) * y(183) + 1.4 * y(185)
 - 0.333601 * y(22) * y(41) + 2.67 * y(42) - 0.005 * y(22) + 0.58 * y(184) + 0.35 * y(185)
 + 0.67 * y(42) - 0.300002 * y(23) * y(20) + 1 * y(21) + 0.005 * y(22) - 0.102001 * y(23) * y(183) + 2.34 * y(184)
 - 72 * y(24) * y(145)*y(145) + 72 * y(142) + 3.15 * y(188) - 0.3 * y(26) * y(24) + 1 * y(25)
 - 0.0611982 * y(25) * y(183) + 12.6 * y(188) + 0.3 * y(26) * y(24) - 1 * y(25)
 + 6.03 * y(80) - 0.102001 * y(26) * y(183) + 21.06 * y(187) - 0.3 * y(26) * y(24) + 1 * y(25) + 0.005 * y(27)
 + 3.15 * y(188) - 0.333601 * y(27) * y(41) + 24.03 * y(80) + 5.22 * y(187) - 0.005 * y(27)
 + 3.60002 * y(29) * y(9) - 10 * y(28) - 0.464999 * y(28) * y(9) + 10 * y(1)
 + 71.999 * y(20) * y(9)*y(9) - 72 * y(29) - 3.60002 * y(29) * y(9) + 10 * y(28)
 + 0.034 * y(82) + 0.34 * y(43) + 2 * y(35) + 1 * y(33) - 5.99994 * y(30) * y(93) + 36 * y(95)
 - 499.981 * y(31) * y(2) + 0.1 * y(32) - 0.0341998 * y(31) * y(81) + 0.136 * y(82) - 0.341998 * y(31) * y(41) + 1.36 * y(43) - 0.64584 * y(31) * y(34) + 8.3334 * y(35) + 9 * y(95)
 + 499.981 * y(31) * y(2) - 0.1 * y(32) - 0.64584 * y(32) * y(34) + 8.3334 * y(36) - 0.341998 * y(32) * y(41) + 1.36 * y(44)
 - 1 * y(33) + 2 * y(36) + 0.34 * y(44)
 - 0.64584 * y(31) * y(34) + 8.3334 * y(35) + 2 * y(35) - 0.64584 * y(32) * y(34) + 8.3334 * y(36) + 2 * y(36) - 0.07176 * y(147) * y(34) + 8.3334 * y(37) + 2 * y(37) - 0.07176 * y(148) * y(34) + 8.3334 * y(38) + 2 * y(38)
 + 0.64584 * y(31) * y(34) - 8.3334 * y(35) - 2 * y(35)
 + 0.64584 * y(32) * y(34) - 8.3334 * y(36) - 2 * y(36)
 + 0.07176 * y(147) * y(34) - 8.3334 * y(37) - 2 * y(37)
 + 0.07176 * y(148) * y(34) - 8.3334 * y(38) - 2 * y(38)
 - 10008 * y(39) * y(9)*y(9) + 1 * y(40)
 + 10008 * y(39) * y(9)*y(9) - 1 * y(40) - 3.60009 * y(9)*y(9) * y(40) + 1 * y(81)
 - 0.333601 * y(22) * y(41) + 2.67 * y(42) + 0.67 * y(42) - 0.341998 * y(31) * y(41) + 1.36 * y(43) + 0.34 * y(43) - 0.341998 * y(32) * y(41) + 1.36 * y(44) + 0.34 * y(44) - 0.0380004 * y(147) * y(41) + 1.36 * y(45) + 0.34 * y(45) - 0.0380004 * y(148) * y(41) + 1.36 * y(46) + 0.34 * y(46) - 2.01182 * y(193) * y(41) + 8 * y(47) + 2 * y(47) - 2.01182 * y(194) * y(41) + 8 * y(48) + 2 * y(48) - 2.01182 * y(203) * y(41) + 8 * y(49) + 2 * y(49) - 2.01182 * y(205) * y(41) + 8 * y(50) + 2 * y(50) - 2.01182 * y(195) * y(41) + 8 * y(51) + 2 * y(51) - 2.01182 * y(206) * y(41) + 8 * y(52) + 2 * y(52) - 2.01182 * y(200) * y(41) + 8 * y(53) + 2 * y(53) - 2.01182 * y(201) * y(41) + 8 * y(54) + 2 * y(54) - 2.01182 * y(196) * y(41) + 8 * y(55) + 2 * y(55) - 2.01182 * y(197) * y(41) + 8 * y(56) + 2 * y(56) - 2.01182 * y(198) * y(41) + 8 * y(57) + 2 * y(57) - 2.01182 * y(199) * y(41) + 8 * y(58) + 2 * y(58) - 2.01182 * y(196) * y(41) + 8 * y(59) + 2 * y(59) - 2.01182 * y(197) * y(41) + 8 * y(60) + 2 * y(60) - 2.01182 * y(198) * y(41) + 8 * y(61) + 2 * y(61) - 2.01182 * y(199) * y(41) + 8 * y(62) + 2 * y(62) - 0.223536 * y(217) * y(41) + 8 * y(63) + 2 * y(63) - 0.223536 * y(217) * y(41) + 8 * y(64) + 2 * y(64) - 0.223536 * y(215) * y(41) + 8 * y(66) + 2 * y(66) - 0.223536 * y(212) * y(41) + 8 * y(67) + 2 * y(67) - 0.223536 * y(221) * y(41) + 8 * y(68) + 2 * y(68) - 0.223536 * y(221) * y(41) + 8 * y(69) + 2 * y(69) - 0.223536 * y(218) * y(41) + 8 * y(70) + 2 * y(70) - 0.223536 * y(213) * y(41) + 8 * y(71) + 2 * y(71) - 0.223536 * y(222) * y(41) + 8 * y(72) + 2 * y(72) - 0.223536 * y(222) * y(41) + 8 * y(73) + 2 * y(73) - 0.223536 * y(219) * y(41) + 8 * y(74) + 2 * y(74) - 0.223536 * y(216) * y(41) + 8 * y(75) + 2 * y(75) - 0.223536 * y(223) * y(41) + 8 * y(76) + 2 * y(76) - 0.223536 * y(223) * y(41) + 8 * y(77) + 2 * y(77) - 0.223536 * y(214) * y(41) + 8 * y(78) + 2 * y(78) - 0.223536 * y(220) * y(41) + 8 * y(79) + 2 * y(79) - 0.0370668 * y(27) * y(41) + 2.67 * y(80) + 0.67 * y(80) + 599.994 * y(1) * y(81) - 1 * y(41)
 + 0.333601 * y(22) * y(41) - 2.67 * y(42) - 0.67 * y(42)
 + 0.341998 * y(31) * y(41) - 1.36 * y(43) - 0.34 * y(43)
 + 0.341998 * y(32) * y(41) - 1.36 * y(44) - 0.34 * y(44)
 + 0.0380004 * y(147) * y(41) - 1.36 * y(45) - 0.34 * y(45)
 + 0.0380004 * y(148) * y(41) - 1.36 * y(46) - 0.34 * y(46)
 + 2.01182 * y(193) * y(41) - 8 * y(47) - 2 * y(47)
 + 2.01182 * y(194) * y(41) - 8 * y(48) - 2 * y(48)
 + 2.01182 * y(203) * y(41) - 8 * y(49) - 2 * y(49)
 + 2.01182 * y(205) * y(41) - 8 * y(50) - 2 * y(50)
 + 2.01182 * y(195) * y(41) - 8 * y(51) - 2 * y(51)
 + 2.01182 * y(206) * y(41) - 8 * y(52) - 2 * y(52)
 + 2.01182 * y(200) * y(41) - 8 * y(53) - 2 * y(53)
 + 2.01182 * y(201) * y(41) - 8 * y(54) - 2 * y(54)
 + 2.01182 * y(196) * y(41) - 8 * y(55) - 2 * y(55)
 + 2.01182 * y(197) * y(41) - 8 * y(56) - 2 * y(56)
 + 2.01182 * y(198) * y(41) - 8 * y(57) - 2 * y(57)
 + 2.01182 * y(199) * y(41) - 8 * y(58) - 2 * y(58)
 + 2.01182 * y(196) * y(41) - 8 * y(59) - 2 * y(59)
 + 2.01182 * y(197) * y(41) - 8 * y(60) - 2 * y(60)
 + 2.01182 * y(198) * y(41) - 8 * y(61) - 2 * y(61)
 + 2.01182 * y(199) * y(41) - 8 * y(62) - 2 * y(62)
 + 0.223536 * y(217) * y(41) - 8 * y(63) - 2 * y(63)
 + 0.223536 * y(217) * y(41) - 8 * y(64) - 2 * y(64)
 + 49.9997 * y(1) * y(129) - 1 * y(65)
 + 0.223536 * y(215) * y(41) - 8 * y(66) - 2 * y(66)
 + 0.223536 * y(212) * y(41) - 8 * y(67) - 2 * y(67)
 + 0.223536 * y(221) * y(41) - 8 * y(68) - 2 * y(68)
 + 0.223536 * y(221) * y(41) - 8 * y(69) - 2 * y(69)
 + 0.223536 * y(218) * y(41) - 8 * y(70) - 2 * y(70)
 + 0.223536 * y(213) * y(41) - 8 * y(71) - 2 * y(71)
 + 0.223536 * y(222) * y(41) - 8 * y(72) - 2 * y(72)
 + 0.223536 * y(222) * y(41) - 8 * y(73) - 2 * y(73)
 + 0.223536 * y(219) * y(41) - 8 * y(74) - 2 * y(74)
 + 0.223536 * y(216) * y(41) - 8 * y(75) - 2 * y(75)
 + 0.223536 * y(223) * y(41) - 8 * y(76) - 2 * y(76)
 + 0.223536 * y(223) * y(41) - 8 * y(77) - 2 * y(77)
 + 0.223536 * y(214) * y(41) - 8 * y(78) - 2 * y(78)
 + 0.223536 * y(220) * y(41) - 8 * y(79) - 2 * y(79)
 + 0.0370668 * y(27) * y(41) - 2.67 * y(80) - 0.67 * y(80)
 - 0.0341998 * y(31) * y(81) + 0.136 * y(82) + 0.034 * y(82) + 3.60009 * y(9)*y(9) * y(40) - 1 * y(81) - 599.994 * y(1) * y(81) + 1 * y(41) - 0.00380004 * y(147) * y(81) + 0.136 * y(83) + 0.034 * y(83)
 + 0.0341998 * y(31) * y(81) - 0.136 * y(82) - 0.034 * y(82)
 + 0.00380004 * y(147) * y(81) - 0.136 * y(83) - 0.034 * y(83)
 - 54 * y(84) * y(8) + 33 * y(85)
 + 54 * y(84) * y(8) - 33 * y(85) - 54 * y(85) * y(8) + 33 * y(86)
 + 54 * y(85) * y(8) - 33 * y(86) - 75.0006 * y(86) * y(8) + 110 * y(87)
 + 75.0006 * y(86) * y(8) - 110 * y(87) - 75.0006 * y(8) * y(87) + 32.5 * y(88)
 + 75.0006 * y(8) * y(87) - 32.5 * y(88) - 60 * y(88) + 17.9998 * y(93) * y(89)
 + 60 * y(88) - 17.9998 * y(93) * y(89) - 60 * y(89) + 17.9998 * y(93) * y(90)
 + 60 * y(89) - 17.9998 * y(93) * y(90)
 - 59.9994 * y(93) * y(91) + 1 * y(92)
 + 59.9994 * y(93) * y(91) - 1 * y(92)
 + 60 * y(88) - 17.9998 * y(93) * y(89) + 60 * y(89) - 17.9998 * y(93) * y(90) - 59.9994 * y(93) * y(91) + 1 * y(92) - 5.99994 * y(133) * y(93) + 36 * y(94) + 9 * y(94) - 5.99994 * y(30) * y(93) + 36 * y(95) + 9 * y(95) - 0.66666 * y(146) * y(93) + 36 * y(96) + 9 * y(96) - 3.99989 * y(191) * y(93) + 24 * y(97) + 6 * y(97) - 3.99989 * y(192) * y(93) + 24 * y(98) + 6 * y(98) - 3.99989 * y(202) * y(93) + 24 * y(99) + 6 * y(99) - 3.99989 * y(204) * y(93) + 24 * y(100) + 6 * y(100) - 3.99989 * y(191) * y(93) + 24 * y(101) + 6 * y(101) - 3.99989 * y(192) * y(93) + 24 * y(102) + 6 * y(102) - 3.99989 * y(202) * y(93) + 24 * y(103) + 6 * y(103) - 3.99989 * y(204) * y(93) + 24 * y(104) + 6 * y(104) - 3.99989 * y(195) * y(93) + 24 * y(105) + 6 * y(105) - 3.99989 * y(206) * y(93) + 24 * y(106) + 6 * y(106) - 3.99989 * y(200) * y(93) + 24 * y(107) + 6 * y(107) - 3.99989 * y(201) * y(93) + 24 * y(108) + 6 * y(108) - 3.99989 * y(193) * y(93) + 24 * y(109) + 6 * y(109) - 3.99989 * y(194) * y(93) + 24 * y(110) + 6 * y(110) - 3.99989 * y(203) * y(93) + 24 * y(111) + 6 * y(111) - 3.99989 * y(205) * y(93) + 24 * y(112) + 6 * y(112) - 0.444432 * y(215) * y(93) + 24 * y(113) + 6 * y(113) - 0.444432 * y(212) * y(93) + 24 * y(114) + 6 * y(114) - 0.444432 * y(208) * y(93) + 24 * y(115) + 6 * y(115) - 0.444432 * y(208) * y(93) + 24 * y(116) + 6 * y(116) - 0.444432 * y(218) * y(93) + 24 * y(117) + 6 * y(117) - 0.444432 * y(213) * y(93) + 24 * y(118) + 6 * y(118) - 0.444432 * y(209) * y(93) + 24 * y(119) + 6 * y(119) - 0.444432 * y(209) * y(93) + 24 * y(120) + 6 * y(120) - 0.444432 * y(219) * y(93) + 24 * y(121) + 6 * y(121) - 0.444432 * y(216) * y(93) + 24 * y(122) + 6 * y(122) - 0.444432 * y(210) * y(93) + 24 * y(123) + 6 * y(123) - 0.444432 * y(210) * y(93) + 24 * y(124) + 6 * y(124) - 0.444432 * y(220) * y(93) + 24 * y(125) + 6 * y(125) - 0.444432 * y(214) * y(93) + 24 * y(126) + 6 * y(126) - 0.444432 * y(211) * y(93) + 24 * y(127) + 6 * y(127) - 0.444432 * y(211) * y(93) + 24 * y(128) + 6 * y(128)
 + 5.99994 * y(133) * y(93) - 36 * y(94) - 9 * y(94)
 + 5.99994 * y(30) * y(93) - 36 * y(95) - 9 * y(95)
 + 0.66666 * y(146) * y(93) - 36 * y(96) - 9 * y(96)
 + 3.99989 * y(191) * y(93) - 24 * y(97) - 6 * y(97)
 + 3.99989 * y(192) * y(93) - 24 * y(98) - 6 * y(98)
 + 3.99989 * y(202) * y(93) - 24 * y(99) - 6 * y(99)
 + 3.99989 * y(204) * y(93) - 24 * y(100) - 6 * y(100)
 + 3.99989 * y(191) * y(93) - 24 * y(101) - 6 * y(101)
 + 3.99989 * y(192) * y(93) - 24 * y(102) - 6 * y(102)
 + 3.99989 * y(202) * y(93) - 24 * y(103) - 6 * y(103)
 + 3.99989 * y(204) * y(93) - 24 * y(104) - 6 * y(104)
 + 3.99989 * y(195) * y(93) - 24 * y(105) - 6 * y(105)
 + 3.99989 * y(206) * y(93) - 24 * y(106) - 6 * y(106)
 + 3.99989 * y(200) * y(93) - 24 * y(107) - 6 * y(107)
 + 3.99989 * y(201) * y(93) - 24 * y(108) - 6 * y(108)
 + 3.99989 * y(193) * y(93) - 24 * y(109) - 6 * y(109)
 + 3.99989 * y(194) * y(93) - 24 * y(110) - 6 * y(110)
 + 3.99989 * y(203) * y(93) - 24 * y(111) - 6 * y(111)
 + 3.99989 * y(205) * y(93) - 24 * y(112) - 6 * y(112)
 + 0.444432 * y(215) * y(93) - 24 * y(113) - 6 * y(113)
 + 0.444432 * y(212) * y(93) - 24 * y(114) - 6 * y(114)
 + 0.444432 * y(208) * y(93) - 24 * y(115) - 6 * y(115)
 + 0.444432 * y(208) * y(93) - 24 * y(116) - 6 * y(116)
 + 0.444432 * y(218) * y(93) - 24 * y(117) - 6 * y(117)
 + 0.444432 * y(213) * y(93) - 24 * y(118) - 6 * y(118)
 + 0.444432 * y(209) * y(93) - 24 * y(119) - 6 * y(119)
 + 0.444432 * y(209) * y(93) - 24 * y(120) - 6 * y(120)
 + 0.444432 * y(219) * y(93) - 24 * y(121) - 6 * y(121)
 + 0.444432 * y(216) * y(93) - 24 * y(122) - 6 * y(122)
 + 0.444432 * y(210) * y(93) - 24 * y(123) - 6 * y(123)
 + 0.444432 * y(210) * y(93) - 24 * y(124) - 6 * y(124)
 + 0.444432 * y(220) * y(93) - 24 * y(125) - 6 * y(125)
 + 0.444432 * y(214) * y(93) - 24 * y(126) - 6 * y(126)
 + 0.444432 * y(211) * y(93) - 24 * y(127) - 6 * y(127)
 + 0.444432 * y(211) * y(93) - 24 * y(128) - 6 * y(128)
 - 49.9997 * y(1) * y(129) + 1 * y(65)
 - 0.1 * y(130) + 4 * y(186)
 + 0.1 * y(130) - 0.599994 * y(131) * y(183) + 16 * y(186)
 + 20 * y(136) + 10 * y(134) + 10 * y(140) + 1.667 * y(138)
 - 2.52002 * y(8) * y(133) + 40 * y(134) + 10 * y(134) + 0.01 * y(135) - 5.99994 * y(133) * y(93) + 36 * y(94)
 + 2.52002 * y(8) * y(133) - 40 * y(134) - 10 * y(134)
 - 5.03998 * y(8) * y(135) + 80 * y(136) + 20 * y(136) - 0.01 * y(135) + 9 * y(94)
 + 5.03998 * y(8) * y(135) - 80 * y(136) - 20 * y(136)
 - 0.210001 * y(8) * y(137) + 6.67 * y(138) + 1.667 * y(138) - 719.982 * y(137) * y(1) + 5 * y(139)
 + 0.210001 * y(8) * y(137) - 6.67 * y(138) - 1.667 * y(138)
 - 1.25998 * y(8) * y(139) + 40 * y(140) + 10 * y(140) + 719.982 * y(137) * y(1) - 5 * y(139)
 + 1.25998 * y(8) * y(139) - 40 * y(140) - 10 * y(140)
 + 5.4 * y(8) - 5.4 * y(141)
 + 72 * y(24) * y(145)*y(145) - 72 * y(142) - 3.6 * y(142) * y(145) + 10 * y(143)
 + 3.6 * y(142) * y(145) - 10 * y(143) - 0.465012 * y(143) * y(145) + 10 * y(144)
 + 0.465012 * y(143) * y(145) - 10 * y(144) - 1000.02 * y(233) * y(144) + 0.1 * y(235) - 49.9998 * y(232) * y(144) - 540 * y(144) + 540 * y(1) + 5 * y(234)
 - 144 * y(24) * y(145)*y(145) + 144 * y(142) - 3.6 * y(142) * y(145) + 10 * y(143) - 0.465012 * y(143) * y(145) + 10 * y(144) + 100 * buff(7) - 100 * y(145)
 + 1 * y(149) - 5.99994 * y(146) * y(93) + 324 * y(96) + 3.06 * y(45) + 18 * y(37) + 0.306 * y(83)
 - 499.98 * y(147) * y(150) + 0.1 * y(148) + 81 * y(96) - 0.342004 * y(147) * y(41) + 12.24 * y(45) - 0.64584 * y(147) * y(34) + 75.0006 * y(37) - 0.0342004 * y(147) * y(81) + 1.224 * y(83)
 + 499.98 * y(147) * y(150) - 0.1 * y(148) - 0.342004 * y(148) * y(41) + 12.24 * y(46) - 0.64584 * y(148) * y(34) + 75.0006 * y(38)
 + 3.06 * y(46) + 18 * y(38) - 1 * y(149)
 + 1 * y(149) - 499.98 * y(147) * y(150) + 0.1 * y(148) - 0.87498 * y(217) * y(150) + 0.68 * y(151) + 0.17 * y(151) - 0.87498 * y(217) * y(150) + 0.68 * y(152) + 0.17 * y(152) - 0.87498 * y(212) * y(150) + 0.68 * y(153) + 0.17 * y(153) - 0.87498 * y(215) * y(150) + 0.68 * y(154) + 0.17 * y(154) - 0.87498 * y(221) * y(150) + 0.68 * y(155) + 0.17 * y(155) - 0.87498 * y(221) * y(150) + 0.68 * y(156) + 0.17 * y(156) - 0.87498 * y(218) * y(150) + 0.68 * y(157) + 0.17 * y(157) - 0.87498 * y(213) * y(150) + 0.68 * y(158) + 0.17 * y(158) - 0.87498 * y(221) * y(150) + 1.4 * y(159) + 0.35 * y(159) - 0.87498 * y(218) * y(150) + 1.4 * y(160) + 0.35 * y(160) - 0.87498 * y(213) * y(150) + 1.4 * y(161) + 0.35 * y(161) - 0.87498 * y(209) * y(150) + 1.4 * y(162) + 0.35 * y(162) - 0.87498 * y(222) * y(150) + 0.68 * y(163) + 0.17 * y(163) - 0.87498 * y(222) * y(150) + 0.68 * y(164) + 0.17 * y(164) - 0.87498 * y(219) * y(150) + 0.68 * y(165) + 0.17 * y(165) - 0.87498 * y(216) * y(150) + 0.68 * y(166) + 0.17 * y(166) - 0.87498 * y(223) * y(150) + 0.68 * y(167) + 0.17 * y(167) - 0.87498 * y(223) * y(150) + 0.68 * y(168) + 0.17 * y(168) - 0.87498 * y(220) * y(150) + 0.68 * y(169) + 0.17 * y(169) - 0.87498 * y(214) * y(150) + 0.68 * y(170) + 0.17 * y(170) - 0.87498 * y(219) * y(150) + 1.4 * y(171) + 0.35 * y(171) - 0.87498 * y(222) * y(150) + 1.4 * y(172) + 0.35 * y(172) - 0.87498 * y(216) * y(150) + 1.4 * y(173) + 0.35 * y(173) - 0.87498 * y(210) * y(150) + 1.4 * y(174) + 0.35 * y(174) - 0.87498 * y(223) * y(150) + 1.4 * y(175) + 0.35 * y(175) - 0.87498 * y(220) * y(150) + 1.4 * y(176) + 0.35 * y(176) - 0.87498 * y(214) * y(150) + 1.4 * y(177) + 0.35 * y(177) - 0.87498 * y(211) * y(150) + 1.4 * y(178) + 0.35 * y(178) - 0.87498 * y(223) * y(150) + 1.4 * y(179) + 0.35 * y(179) - 0.87498 * y(220) * y(150) + 1.4 * y(180) + 0.35 * y(180) - 0.87498 * y(214) * y(150) + 1.4 * y(181) + 0.35 * y(181) - 0.87498 * y(211) * y(150) + 1.4 * y(182) + 0.35 * y(182)
 + 0.87498 * y(217) * y(150) - 0.68 * y(151) - 0.17 * y(151)
 + 0.87498 * y(217) * y(150) - 0.68 * y(152) - 0.17 * y(152)
 + 0.87498 * y(212) * y(150) - 0.68 * y(153) - 0.17 * y(153)
 + 0.87498 * y(215) * y(150) - 0.68 * y(154) - 0.17 * y(154)
 + 0.87498 * y(221) * y(150) - 0.68 * y(155) - 0.17 * y(155)
 + 0.87498 * y(221) * y(150) - 0.68 * y(156) - 0.17 * y(156)
 + 0.87498 * y(218) * y(150) - 0.68 * y(157) - 0.17 * y(157)
 + 0.87498 * y(213) * y(150) - 0.68 * y(158) - 0.17 * y(158)
 + 0.87498 * y(221) * y(150) - 1.4 * y(159) - 0.35 * y(159)
 + 0.87498 * y(218) * y(150) - 1.4 * y(160) - 0.35 * y(160)
 + 0.87498 * y(213) * y(150) - 1.4 * y(161) - 0.35 * y(161)
 + 0.87498 * y(209) * y(150) - 1.4 * y(162) - 0.35 * y(162)
 + 0.87498 * y(222) * y(150) - 0.68 * y(163) - 0.17 * y(163)
 + 0.87498 * y(222) * y(150) - 0.68 * y(164) - 0.17 * y(164)
 + 0.87498 * y(219) * y(150) - 0.68 * y(165) - 0.17 * y(165)
 + 0.87498 * y(216) * y(150) - 0.68 * y(166) - 0.17 * y(166)
 + 0.87498 * y(223) * y(150) - 0.68 * y(167) - 0.17 * y(167)
 + 0.87498 * y(223) * y(150) - 0.68 * y(168) - 0.17 * y(168)
 + 0.87498 * y(220) * y(150) - 0.68 * y(169) - 0.17 * y(169)
 + 0.87498 * y(214) * y(150) - 0.68 * y(170) - 0.17 * y(170)
 + 0.87498 * y(219) * y(150) - 1.4 * y(171) - 0.35 * y(171)
 + 0.87498 * y(222) * y(150) - 1.4 * y(172) - 0.35 * y(172)
 + 0.87498 * y(216) * y(150) - 1.4 * y(173) - 0.35 * y(173)
 + 0.87498 * y(210) * y(150) - 1.4 * y(174) - 0.35 * y(174)
 + 0.87498 * y(223) * y(150) - 1.4 * y(175) - 0.35 * y(175)
 + 0.87498 * y(220) * y(150) - 1.4 * y(176) - 0.35 * y(176)
 + 0.87498 * y(214) * y(150) - 1.4 * y(177) - 0.35 * y(177)
 + 0.87498 * y(211) * y(150) - 1.4 * y(178) - 0.35 * y(178)
 + 0.87498 * y(223) * y(150) - 1.4 * y(179) - 0.35 * y(179)
 + 0.87498 * y(220) * y(150) - 1.4 * y(180) - 0.35 * y(180)
 + 0.87498 * y(214) * y(150) - 1.4 * y(181) - 0.35 * y(181)
 + 0.87498 * y(211) * y(150) - 1.4 * y(182) - 0.35 * y(182)
 - 0.102001 * y(23) * y(183) + 2.34 * y(184) + 0.58 * y(184) - 0.0611982 * y(21) * y(183) + 1.4 * y(185) + 0.35 * y(185) - 0.599994 * y(131) * y(183) + 16 * y(186) + 4 * y(186) - 0.0067998 * y(25) * y(183) + 1.4 * y(188) + 0.35 * y(188) - 0.0113334 * y(26) * y(183) + 2.34 * y(187) + 0.58 * y(187) + 2.5 * buff(5) - 2.5 * y(183)
 + 0.102001 * y(23) * y(183) - 2.34 * y(184) - 0.58 * y(184)
 + 0.0611982 * y(21) * y(183) - 1.4 * y(185) - 0.35 * y(185)
 + 0.599994 * y(131) * y(183) - 16 * y(186) - 4 * y(186)
 + 0.0113334 * y(26) * y(183) - 2.34 * y(187) - 0.58 * y(187)
 + 0.0067998 * y(25) * y(183) - 1.4 * y(188) - 0.35 * y(188)
 - 0.00035 * y(189) + 0.0126 * y(190)
 + 3.88889e-05 * y(189) - 0.0014 * y(190)
 + 2 * y(47) - 3.99989 * y(191) * y(93) + 24 * y(97) - 3.99989 * y(191) * y(93) + 24 * y(101) + 2 * y(51) + 8.88889e-05 * y(208) + 1 * buff(8) - 1 * y(191)
 + 2 * y(48) - 3.99989 * y(192) * y(93) + 24 * y(98) - 3.99989 * y(192) * y(93) + 24 * y(102) + 2 * y(52) + 8.88889e-05 * y(209)
 - 2.01182 * y(193) * y(41) + 8 * y(47) + 6 * y(97) - 3.99989 * y(193) * y(93) + 24 * y(109) + 2 * y(59)
 - 2.01182 * y(194) * y(41) + 8 * y(48) + 6 * y(98) - 3.99989 * y(194) * y(93) + 24 * y(110) + 2 * y(60)
 + 6 * y(101) - 2.01182 * y(195) * y(41) + 8 * y(51) - 3.99989 * y(195) * y(93) + 24 * y(105) + 2 * y(55)
 + 6 * y(105) - 2.01182 * y(196) * y(41) + 8 * y(55) + 6 * y(109) - 2.01182 * y(196) * y(41) + 8 * y(59) - 0.0012 * y(196) * y(207) + 0.000888889 * y(217) - 3.6e-05 * y(196)
 + 6 * y(106) - 2.01182 * y(197) * y(41) + 8 * y(56) + 6 * y(110) - 2.01182 * y(197) * y(41) + 8 * y(60) - 0.0012 * y(197) * y(207) + 0.000888889 * y(221) - 3.6e-05 * y(197)
 + 6 * y(107) - 2.01182 * y(198) * y(41) + 8 * y(57) + 6 * y(111) - 2.01182 * y(198) * y(41) + 8 * y(61) - 0.0012 * y(198) * y(207) + 0.000888889 * y(222) - 3.6e-05 * y(198)
 + 6 * y(108) - 2.01182 * y(199) * y(41) + 8 * y(58) + 6 * y(112) - 2.01182 * y(199) * y(41) + 8 * y(62) - 0.0012 * y(199) * y(207) + 0.000888889 * y(223) - 3.6e-05 * y(199)
 + 6 * y(103) - 2.01182 * y(200) * y(41) + 8 * y(53) - 3.99989 * y(200) * y(93) + 24 * y(107) + 2 * y(57)
 + 6 * y(104) - 2.01182 * y(201) * y(41) + 8 * y(54) - 3.99989 * y(201) * y(93) + 24 * y(108) + 2 * y(58)
 + 2 * y(49) - 3.99989 * y(202) * y(93) + 24 * y(99) - 3.99989 * y(202) * y(93) + 24 * y(103) + 2 * y(53) + 8.88889e-05 * y(210)
 - 2.01182 * y(203) * y(41) + 8 * y(49) + 6 * y(99) - 3.99989 * y(203) * y(93) + 24 * y(111) + 2 * y(61)
 + 2 * y(50) - 3.99989 * y(204) * y(93) + 24 * y(100) - 3.99989 * y(204) * y(93) + 24 * y(104) + 2 * y(54) + 8.88889e-05 * y(211)
 - 2.01182 * y(205) * y(41) + 8 * y(50) + 6 * y(100) - 3.99989 * y(205) * y(93) + 24 * y(112) + 2 * y(62)
 + 6 * y(102) - 2.01182 * y(206) * y(41) + 8 * y(52) - 3.99989 * y(206) * y(93) + 24 * y(106) + 2 * y(56)
 - 0.0108 * y(199) * y(207) + 0.008 * y(223) - 0.0108 * y(198) * y(207) + 0.008 * y(222) - 0.0108 * y(197) * y(207) + 0.008 * y(221) - 0.0108 * y(196) * y(207) + 0.008 * y(217) + 0.0008 * y(208) + 0.0008 * y(209) + 0.0008 * y(210) + 0.0008 * y(211)
 - 3.99989 * y(208) * y(93) + 216 * y(116) + 0.17 * y(154) - 3.99989 * y(208) * y(93) + 216 * y(115) + 0.17 * y(153) - 0.0008 * y(208) - 0.0277776 * y(208) *  (y(235) + y(234) + 1 *  (y(233) + y(231) + y(242) ) ) + 2 * y(246) + 0.35 * y(162) + 18 * y(66) + 18 * y(67) + 0.35 * y(174) - 0.0277776 * y(208) *  (y(235) + y(234) + 1 *  (y(233) + y(231) + y(242) ) ) + 2 * y(250)
 + 0.17 * y(158) + 0.17 * y(157) - 3.99989 * y(209) * y(93) + 216 * y(120) - 3.99989 * y(209) * y(93) + 216 * y(119) + 0.5 * y(246) - 0.87498 * y(209) * y(150) + 1.4 * y(162) - 0.0008 * y(209) + 18 * y(71) + 18 * y(70) - 0.0277776 * y(209) *  (y(235) + y(234) + 1 *  (y(233) + y(231) + y(242) ) ) + 2 * y(258) + 0.35 * y(182)
 - 3.99989 * y(210) * y(93) + 216 * y(124) - 3.99989 * y(210) * y(93) + 216 * y(123) + 0.17 * y(165) + 0.17 * y(166) - 0.0008 * y(210) + 0.5 * y(250) - 0.0277776 * y(210) *  (y(235) + y(234) + 1 *  (y(233) + y(231) + y(242) ) ) + 2 * y(254) - 0.87498 * y(210) * y(150) + 1.4 * y(174) + 0.35 * y(178) + 18 * y(74) + 18 * y(75)
 - 3.99989 * y(211) * y(93) + 216 * y(128) - 3.99989 * y(211) * y(93) + 216 * y(127) + 0.17 * y(169) + 0.17 * y(170) - 0.0008 * y(211) + 0.5 * y(254) - 0.87498 * y(211) * y(150) + 1.4 * y(178) + 18 * y(78) + 18 * y(79) + 0.5 * y(258) - 0.87498 * y(211) * y(150) + 1.4 * y(182)
 + 0.17 * y(152) - 3.99989 * y(212) * y(93) + 216 * y(114) + 54 * y(115) - 0.87498 * y(212) * y(150) + 0.68 * y(153) - 0.0277776 * y(212) *  (y(235) + y(234) + 1 *  (y(233) + y(231) + y(242) ) ) + 2 * y(245) + 0.35 * y(161) + 18 * y(64) - 2.01182 * y(212) * y(41) + 72 * y(67) + 0.35 * y(173) - 0.0277776 * y(212) *  (y(235) + y(234) + 1 *  (y(233) + y(231) + y(242) ) ) + 2 * y(249)
 + 0.17 * y(156) - 0.87498 * y(213) * y(150) + 0.68 * y(158) + 54 * y(120) - 3.99989 * y(213) * y(93) + 216 * y(118) + 0.5 * y(245) - 0.87498 * y(213) * y(150) + 1.4 * y(161) - 2.01182 * y(213) * y(41) + 72 * y(71) + 18 * y(69) + 0.35 * y(181) - 0.0277776 * y(213) *  (y(235) + y(234) + 1 *  (y(233) + y(231) + y(242) ) ) + 2 * y(257)
 - 3.99989 * y(214) * y(93) + 216 * y(126) + 54 * y(128) + 0.17 * y(168) - 0.87498 * y(214) * y(150) + 0.68 * y(170) + 0.5 * y(253) - 0.87498 * y(214) * y(150) + 1.4 * y(177) + 18 * y(77) - 2.01182 * y(214) * y(41) + 72 * y(78) + 0.5 * y(257) - 0.87498 * y(214) * y(150) + 1.4 * y(181)
 + 0.17 * y(151) - 3.99989 * y(215) * y(93) + 216 * y(113) - 0.87498 * y(215) * y(150) + 0.68 * y(154) + 54 * y(116) - 0.0277776 * y(215) *  (y(235) + y(234) + 1 *  (y(233) + y(231) + y(242) ) ) + 2 * y(244) + 0.35 * y(160) + 18 * y(63) - 2.01182 * y(215) * y(41) + 72 * y(66) + 0.35 * y(171) - 0.0277776 * y(215) *  (y(235) + y(234) + 1 *  (y(233) + y(231) + y(242) ) ) + 2 * y(248)
 - 3.99989 * y(216) * y(93) + 216 * y(122) + 54 * y(124) + 0.17 * y(164) - 0.87498 * y(216) * y(150) + 0.68 * y(166) + 0.5 * y(249) - 0.0277776 * y(216) *  (y(235) + y(234) + 1 *  (y(233) + y(231) + y(242) ) ) + 2 * y(253) - 0.87498 * y(216) * y(150) + 1.4 * y(173) + 0.35 * y(177) + 18 * y(73) - 2.01182 * y(216) * y(41) + 72 * y(75)
 - 0.87498 * y(217) * y(150) + 0.68 * y(151) + 54 * y(113) - 0.87498 * y(217) * y(150) + 0.68 * y(152) + 54 * y(114) + 0.0108 * y(196) * y(207) - 0.008 * y(217) - 0.0277776 * y(217) *  (y(235) + y(234) + 1 *  (y(233) + y(231) + y(242) ) ) + 2 * y(243) + 0.35 * y(159) - 2.01182 * y(217) * y(41) + 72 * y(63) - 2.01182 * y(217) * y(41) + 72 * y(64) - 0.0277776 * y(217) *  (y(235) + y(234) + 1 *  (y(233) + y(231) + y(242) ) ) + 2 * y(247) + 0.35 * y(172)
 + 0.17 * y(155) - 0.87498 * y(218) * y(150) + 0.68 * y(157) + 54 * y(119) - 3.99989 * y(218) * y(93) + 216 * y(117) + 0.5 * y(244) - 0.87498 * y(218) * y(150) + 1.4 * y(160) - 2.01182 * y(218) * y(41) + 72 * y(70) + 18 * y(68) - 0.0277776 * y(218) *  (y(235) + y(234) + 1 *  (y(233) + y(231) + y(242) ) ) + 2 * y(256) + 0.35 * y(180)
 - 3.99989 * y(219) * y(93) + 216 * y(121) + 54 * y(123) + 0.17 * y(163) - 0.87498 * y(219) * y(150) + 0.68 * y(165) + 0.5 * y(248) - 0.0277776 * y(219) *  (y(235) + y(234) + 1 *  (y(233) + y(231) + y(242) ) ) + 2 * y(252) - 0.87498 * y(219) * y(150) + 1.4 * y(171) + 0.35 * y(176) + 18 * y(72) - 2.01182 * y(219) * y(41) + 72 * y(74)
 - 3.99989 * y(220) * y(93) + 216 * y(125) + 54 * y(127) + 0.17 * y(167) - 0.87498 * y(220) * y(150) + 0.68 * y(169) + 0.5 * y(252) - 0.87498 * y(220) * y(150) + 1.4 * y(176) + 18 * y(76) - 2.01182 * y(220) * y(41) + 72 * y(79) - 0.87498 * y(220) * y(150) + 1.4 * y(180) + 0.5 * y(256)
 - 0.87498 * y(221) * y(150) + 0.68 * y(155) - 0.87498 * y(221) * y(150) + 0.68 * y(156) + 54 * y(118) + 54 * y(117) + 0.5 * y(243) - 0.87498 * y(221) * y(150) + 1.4 * y(159) + 0.0108 * y(197) * y(207) - 0.008 * y(221) - 2.01182 * y(221) * y(41) + 72 * y(68) - 2.01182 * y(221) * y(41) + 72 * y(69) - 0.0277776 * y(221) *  (y(235) + y(234) + 1 *  (y(233) + y(231) + y(242) ) ) + 2 * y(255) + 0.35 * y(179)
 + 54 * y(121) + 54 * y(122) - 0.87498 * y(222) * y(150) + 0.68 * y(163) - 0.87498 * y(222) * y(150) + 0.68 * y(164) + 0.0108 * y(198) * y(207) - 0.008 * y(222) + 0.5 * y(247) - 0.0277776 * y(222) *  (y(235) + y(234) + 1 *  (y(233) + y(231) + y(242) ) ) + 2 * y(251) - 0.87498 * y(222) * y(150) + 1.4 * y(172) + 0.35 * y(175) - 2.01182 * y(222) * y(41) + 72 * y(72) - 2.01182 * y(222) * y(41) + 72 * y(73)
 + 54 * y(126) + 54 * y(125) - 0.87498 * y(223) * y(150) + 0.68 * y(167) - 0.87498 * y(223) * y(150) + 0.68 * y(168) + 0.0108 * y(199) * y(207) - 0.008 * y(223) + 0.5 * y(251) - 0.87498 * y(223) * y(150) + 1.4 * y(175) - 2.01182 * y(223) * y(41) + 72 * y(76) - 2.01182 * y(223) * y(41) + 72 * y(77) + 0.5 * y(255) - 0.87498 * y(223) * y(150) + 1.4 * y(179)
 - 0.499998 * y(235) * y(224) + 0.8 * y(225) + 0.2 * y(225) - 0.499998 * y(233) * y(224) + 0.8 * y(226) + 0.2 * y(226) - 0.499998 * y(231) * y(224) + 0.8 * y(227) + 0.2 * y(227) - 0.499998 * y(230) * y(224) + 0.8 * y(229) + 0.2 * y(229) - 0.499998 * y(231) * y(224) + 0.8 * y(228) + 0.2 * y(228)
 + 0.499998 * y(235) * y(224) - 0.8 * y(225) - 0.2 * y(225)
 + 0.499998 * y(233) * y(224) - 0.8 * y(226) - 0.2 * y(226)
 + 0.499998 * y(231) * y(224) - 0.8 * y(227) - 0.2 * y(227)
 + 0.499998 * y(231) * y(224) - 0.8 * y(228) - 0.2 * y(228)
 + 0.499998 * y(230) * y(224) - 0.8 * y(229) - 0.2 * y(229)
 - 0.3 * y(230) + 0.00054 * y(19) * y(259) - 0.499998 * y(230) * y(224) + 0.8 * y(229) + 0.2 * y(228)
 + 6 * y(236) + 6 * y(239) - 0.499998 * y(231) * y(224) + 0.8 * y(227) - 0.499998 * y(231) * y(224) + 0.8 * y(228)
 - 0.3 * y(232) + 0.00054 * y(10) * y(259) - 49.9998 * y(232) * y(144) + 5 * y(234) + 0.2 * y(229) + 0.2 * y(226) - 0.0312498 * y(232) *  (y(235) + y(234) ) + 8 * y(238) - 0.0199998 * y(232) *  (y(233) + y(231) + y(242) ) + 8 * y(241)
 - 0.09375 * y(233) *  (y(235) + y(234) ) + 24 * y(236) - 0.06 * y(233) *  (y(233) + y(231) + y(242) ) + 24 * y(239) - 1000.02 * y(233) * y(144) + 0.1 * y(235) - 0.499998 * y(233) * y(224) + 0.8 * y(226) + 0.2 * y(227) + 2 * y(238) + 2 * y(241)
 + 0.00108 * y(11) * y(259) - 0.0078126 * y(234) *  (y(235) + y(234) ) + 2 * y(237) - 0.00499998 * y(234) *  (y(233) + y(231) + y(242) ) + 2 * y(240) - 5 * y(234) + 49.9998 * y(232) * y(144) + 0.2 * y(225)
 + 0.00108 * y(12) * y(259) + 1000.02 * y(233) * y(144) - 0.1 * y(235) + 0.5 * y(237) + 0.5 * y(240) - 0.499998 * y(235) * y(224) + 0.8 * y(225)
 + 0.09375 * y(233) *  (y(235) + y(234) ) - 24 * y(236) - 6 * y(236)
 + 0.0078126 * y(234) *  (y(235) + y(234) ) - 2 * y(237) - 0.5 * y(237)
 + 0.0312498 * y(232) *  (y(235) + y(234) ) - 8 * y(238) - 2 * y(238)
 + 0.06 * y(233) *  (y(233) + y(231) + y(242) ) - 24 * y(239) - 6 * y(239)
 + 0.00499998 * y(234) *  (y(233) + y(231) + y(242) ) - 2 * y(240) - 0.5 * y(240)
 + 0.0199998 * y(232) *  (y(233) + y(231) + y(242) ) - 8 * y(241) - 2 * y(241)
 + 1 * buff(9) - 1 * y(242)
 + 0.0277776 * y(217) *  (y(235) + y(234) + 1 *  (y(233) + y(231) + y(242) ) ) - 2 * y(243) - 0.5 * y(243)
 + 0.0277776 * y(215) *  (y(235) + y(234) + 1 *  (y(233) + y(231) + y(242) ) ) - 2 * y(244) - 0.5 * y(244)
 + 0.0277776 * y(212) *  (y(235) + y(234) + 1 *  (y(233) + y(231) + y(242) ) ) - 2 * y(245) - 0.5 * y(245)
 + 0.0277776 * y(208) *  (y(235) + y(234) + 1 *  (y(233) + y(231) + y(242) ) ) - 2 * y(246) - 0.5 * y(246)
 + 0.0277776 * y(217) *  (y(235) + y(234) + 1 *  (y(233) + y(231) + y(242) ) ) - 2 * y(247) - 0.5 * y(247)
 + 0.0277776 * y(215) *  (y(235) + y(234) + 1 *  (y(233) + y(231) + y(242) ) ) - 2 * y(248) - 0.5 * y(248)
 + 0.0277776 * y(212) *  (y(235) + y(234) + 1 *  (y(233) + y(231) + y(242) ) ) - 2 * y(249) - 0.5 * y(249)
 + 0.0277776 * y(208) *  (y(235) + y(234) + 1 *  (y(233) + y(231) + y(242) ) ) - 2 * y(250) - 0.5 * y(250)
 + 0.0277776 * y(222) *  (y(235) + y(234) + 1 *  (y(233) + y(231) + y(242) ) ) - 2 * y(251) - 0.5 * y(251)
 + 0.0277776 * y(219) *  (y(235) + y(234) + 1 *  (y(233) + y(231) + y(242) ) ) - 2 * y(252) - 0.5 * y(252)
 + 0.0277776 * y(216) *  (y(235) + y(234) + 1 *  (y(233) + y(231) + y(242) ) ) - 2 * y(253) - 0.5 * y(253)
 + 0.0277776 * y(210) *  (y(235) + y(234) + 1 *  (y(233) + y(231) + y(242) ) ) - 2 * y(254) - 0.5 * y(254)
 + 0.0277776 * y(221) *  (y(235) + y(234) + 1 *  (y(233) + y(231) + y(242) ) ) - 2 * y(255) - 0.5 * y(255)
 + 0.0277776 * y(218) *  (y(235) + y(234) + 1 *  (y(233) + y(231) + y(242) ) ) - 2 * y(256) - 0.5 * y(256)
 + 0.0277776 * y(213) *  (y(235) + y(234) + 1 *  (y(233) + y(231) + y(242) ) ) - 2 * y(257) - 0.5 * y(257)
 + 0.0277776 * y(209) *  (y(235) + y(234) + 1 *  (y(233) + y(231) + y(242) ) ) - 2 * y(258) - 0.5 * y(258)
 - 0.00108 * y(12) * y(259) - 0.00108 * y(11) * y(259) + 0.3 * y(232) - 0.00054 * y(10) * y(259) + 0.3 * y(230) - 0.00054 * y(19) * y(259)
];

% end of model
