Berry64.1.eq 1 Berry et al. Chaos in a bienzymatic cyclic model with two autocatalytic loops, Chaos, Solitons & Fractals 2003, 18:1001-1014 2 Add rate control on p and g, according to paper causes bifurcations (hence alpha and beta parameters). 3 Push me , pull me with p and g, ppp pos and ppg neg gives more interesting results with smaller control 4 Good demonstration of control both limiting (ppg) and accelerating (ppp) to achieve control 5 Changing rm from chaos to multiorbit will cause control on m (see berry6.eq) to lose control but with control on p and g (at fp=fg=1) control will remain throughout the parameter change of rm. 6 Switches ri from 0.0098 to 0.00995 and back to demonstrate the continuous control of the rate control method. 7 Double berry, selective external input to show actvation of different sets, ie dynamic resonance. 8 Same resource 9 different fractions of soluble fragements f but same m 10 proportional increase/decrease sensitivity of control function 11 removed fabs 12 removed external forcing 13 scaling ? 4 cells 14 allow emergent criticality by having 4 separate m matrices with migration of f 15 no switching between modes 16 smaller external input 17 from Berry2.12.eq, 4 by 4 network 18 changed perturbations 19 Step, length 500000, random high, pp's=-1 20 based on berry4.4bnew.eq 21 32 oscillators 22 based on berry32.1.eq, 64 oscillators, simplified parameters Kg=0.1 Kd=1.1 kg=0.05 kd=0.05 %ka=kd/Kd ka=0.0455 Kr=4.5 Ks=1 alpha=0.026 beta=0.00075 n=4 l=4 scm=20 scf=2 ppp=-1 fp=1 adp=0 ppg=-1 fg=1 adg=0 Erim1=7.5 Erim2=1 Erim3=8 Erim4=3.25 Erim5=7.25 Erim6=1.1 Erim7=7.4 Erim8=3.5 Erim9=7.15 Erim10=1.56 Erim11=8.2 Erim12=3.2 Erim13=7.1 Erim14=1.9 Erim15=6.4 Erim16=3 Erim17=4.5 Erim18=4.8 Erim19=8.7 Erim20=3.9 Erim21=2.5 Erim22=5 Erim23=2.6 Erim24=3.66 Erim25=1.34 Erim26=4.1 Erim27=8.2 Erim28=3.125 Erim29=7.85 Erim30=6 Erim31=9.0 Erim32=4.1 Erim33=7.5 Erim34=1 Erim35=8 Erim36=3.25 Erim37=7.25 Erim38=1.1 Erim39=7.4 Erim40=3.5 Erim41=7.15 Erim42=1.56 Erim43=8.2 Erim44=3.2 Erim45=7.1 Erim46=1.9 Erim47=6.4 Erim48=3 Erim49=4.5 Erim50=4.8 Erim51=8.7 Erim52=3.9 Erim53=2.5 Erim54=5 Erim55=2.6 Erim56=3.66 Erim57=1.34 Erim58=4.1 Erim59=8.2 Erim60=3.125 Erim61=7.85 Erim62=6 Erim63=9.0 Erim64=4.1 % No1=0.000316666666667 % No2=0.000326666666667 % No3=0.000316666666667 % No4=0.000326666666667 Nos=0.00005 %rim=0.00098 %rm2=0.00995 %stlen=100000 stlen=50000 stsc=0.00025 rlow=-1 rhigh=1 st=step(t,stlen,0.1) stscalar=abs(stscalar+st*stsc*random(rlow,rhigh)) qf1=f1/(f1+scf) qf2=f2/(f2+scf) qf3=f3/(f3+scf) qf4=f4/(f4+scf) qf5=f5/(f5+scf) qf6=f6/(f6+scf) qf7=f7/(f7+scf) qf8=f8/(f8+scf) qf9=f9/(f9+scf) qf10=f10/(f10+scf) qf11=f11/(f11+scf) qf12=f12/(f12+scf) qf13=f13/(f13+scf) qf14=f14/(f14+scf) qf15=f15/(f15+scf) qf16=f16/(f16+scf) qf17=f17/(f17+scf) qf18=f18/(f18+scf) qf19=f19/(f19+scf) qf20=f20/(f20+scf) qf21=f21/(f21+scf) qf22=f22/(f22+scf) qf23=f23/(f23+scf) qf24=f24/(f24+scf) qf25=f25/(f25+scf) qf26=f26/(f26+scf) qf27=f27/(f27+scf) qf28=f28/(f28+scf) qf29=f29/(f29+scf) qf30=f30/(f30+scf) qf31=f31/(f31+scf) qf32=f32/(f32+scf) qf33=f33/(f33+scf) qf34=f34/(f34+scf) qf35=f35/(f35+scf) qf36=f36/(f36+scf) qf37=f37/(f37+scf) qf38=f38/(f38+scf) qf39=f39/(f39+scf) qf40=f40/(f40+scf) qf41=f41/(f41+scf) qf42=f42/(f42+scf) qf43=f43/(f43+scf) qf44=f44/(f44+scf) qf45=f45/(f45+scf) qf46=f46/(f46+scf) qf47=f47/(f47+scf) qf48=f48/(f47+scf) qf49=f49/(f49+scf) qf50=f50/(f50+scf) qf51=f51/(f51+scf) qf52=f52/(f52+scf) qf53=f53/(f53+scf) qf54=f54/(f54+scf) qf55=f55/(f55+scf) qf56=f56/(f56+scf) qf57=f57/(f57+scf) qf58=f58/(f58+scf) qf59=f59/(f59+scf) qf60=f60/(f60+scf) qf61=f61/(f61+scf) qf62=f62/(f62+scf) qf63=f63/(f63+scf) qf64=f63/(f64+scf) ftotal=f1+f2+f3+f4+f5+f6+f7+f8+f9+f10+f11+f12+f13+f14+f15+f16+f17+f18+f19+f20+f21+f22+f23+f24+f25+f26+f27+f28+f29+f30+f31+f32+f33+f34+f35+f36+f37+f39+f39+f40+f41+f42+f43+f44+f45+f46+f47+f48+f49+f50+f51+f52+f53+f54+f55+f56+f57+f58+f59+f60+f61+f62+f63+f64 mtotal=m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m32+m33+m34+m35+m36+m37+m39+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64 qpp1=qf1+adp sp1=fp*exp(ppp*qpp1) qpg1=qf1+adg sg1=fg*exp(ppg*qpg1) qpp2=qf2+adp sp2=fp*exp(ppp*qpp2) qpg2=qf2+adg sg2=fg*exp(ppg*qpg2) qpp3=qf3+adp sp3=fp*exp(ppp*qpp3) qpg3=qf3+adg sg3=fg*exp(ppg*qpg3) qpp4=qf4+adp sp4=fp*exp(ppp*qpp4) qpg4=qf4+adg sg4=fg*exp(ppg*qpg4) qpp5=qf5+adp sp5=fp*exp(ppp*qpp5) qpg5=qf5+adg sg5=fg*exp(ppg*qpg5) qpp6=qf6+adp sp6=fp*exp(ppp*qpp6) qpg6=qf6+adg sg6=fg*exp(ppg*qpg6) qpp7=qf7+adp sp7=fp*exp(ppp*qpp7) qpg7=qf7+adg sg7=fg*exp(ppg*qpg7) qpp8=qf8+adp sp8=fp*exp(ppp*qpp8) qpg8=qf8+adg sg8=fg*exp(ppg*qpg8) qpp9=qf9+adp sp9=fp*exp(ppp*qpp9) qpg9=qf9+adg sg9=fg*exp(ppg*qpg9) qpp10=qf10+adp sp10=fp*exp(ppp*qpp10) qpg10=qf10+adg sg10=fg*exp(ppg*qpg10) qpp11=qf11+adp sp11=fp*exp(ppp*qpp11) qpg11=qf11+adg sg11=fg*exp(ppg*qpg11) qpp12=qf12+adp sp12=fp*exp(ppp*qpp12) qpg12=qf12+adg sg12=fg*exp(ppg*qpg12) qpp13=qf13+adp sp13=fp*exp(ppp*qpp13) qpg13=qf13+adg sg13=fg*exp(ppg*qpg13) qpp14=qf14+adp sp14=fp*exp(ppp*qpp14) qpg14=qf14+adg sg14=fg*exp(ppg*qpg14) qpp15=qf15+adp sp15=fp*exp(ppp*qpp15) qpg15=qf15+adg sg15=fg*exp(ppg*qpg15) qpp16=qf16+adp sp16=fp*exp(ppp*qpp16) qpg16=qf16+adg sg16=fg*exp(ppg*qpg16) qpp17=qf17+adp sp17=fp*exp(ppp*qpp17) qpg17=qf17+adg sg17=fg*exp(ppg*qpg17) qpp18=qf18+adp sp18=fp*exp(ppp*qpp18) qpg18=qf18+adg sg18=fg*exp(ppg*qpg18) qpp19=qf19+adp sp19=fp*exp(ppp*qpp19) qpg19=qf19+adg sg19=fg*exp(ppg*qpg19) qpp20=qf20+adp sp20=fp*exp(ppp*qpp20) qpg20=qf20+adg sg20=fg*exp(ppg*qpg20) qpp21=qf21+adp sp21=fp*exp(ppp*qpp21) qpg21=qf21+adg sg21=fg*exp(ppg*qpg21) qpp22=qf22+adp sp22=fp*exp(ppp*qpp22) qpg22=qf22+adg sg22=fg*exp(ppg*qpg22) qpp23=qf23+adp sp23=fp*exp(ppp*qpp23) qpg23=qf23+adg sg23=fg*exp(ppg*qpg23) qpp24=qf24+adp sp24=fp*exp(ppp*qpp24) qpg24=qf24+adg sg24=fg*exp(ppg*qpg24) qpp25=qf25+adp sp25=fp*exp(ppp*qpp25) qpg25=qf25+adg sg25=fg*exp(ppg*qpg25) qpp26=qf26+adp sp26=fp*exp(ppp*qpp26) qpg26=qf26+adg sg26=fg*exp(ppg*qpg26) qpp27=qf27+adp sp27=fp*exp(ppp*qpp27) qpg27=qf27+adg sg27=fg*exp(ppg*qpg27) qpp28=qf28+adp sp28=fp*exp(ppp*qpp28) qpg28=qf28+adg sg28=fg*exp(ppg*qpg28) qpp29=qf29+adp sp29=fp*exp(ppp*qpp29) qpg29=qf29+adg sg29=fg*exp(ppg*qpg29) qpp30=qf30+adp sp30=fp*exp(ppp*qpp30) qpg30=qf30+adg sg30=fg*exp(ppg*qpg30) qpp31=qf31+adp sp31=fp*exp(ppp*qpp31) qpg31=qf31+adg sg31=fg*exp(ppg*qpg31) qpp32=qf32+adp sp32=fp*exp(ppp*qpp32) qpg32=qf32+adg sg32=fg*exp(ppg*qpg32) qpp33=qf33+adp sp33=fp*exp(ppp*qpp33) qpg33=qf33+adg sg33=fg*exp(ppg*qpg33) qpp34=qf34+adp sp34=fp*exp(ppp*qpp34) qpg34=qf34+adg sg34=fg*exp(ppg*qpg34) qpp35=qf35+adp sp35=fp*exp(ppp*qpp35) qpg35=qf35+adg sg35=fg*exp(ppg*qpg35) qpp36=qf36+adp sp36=fp*exp(ppp*qpp36) qpg36=qf36+adg sg36=fg*exp(ppg*qpg36) qpp37=qf37+adp sp37=fp*exp(ppp*qpp37) qpg37=qf37+adg sg37=fg*exp(ppg*qpg37) qpp38=qf38+adp sp38=fp*exp(ppp*qpp38) qpg38=qf38+adg sg38=fg*exp(ppg*qpg38) qpp39=qf39+adp sp39=fp*exp(ppp*qpp39) qpg39=qf39+adg sg39=fg*exp(ppg*qpg39) qpp40=qf40+adp sp40=fp*exp(ppp*qpp40) qpg40=qf40+adg sg40=fg*exp(ppg*qpg40) qpp41=qf41+adp sp41=fp*exp(ppp*qpp41) qpg41=qf41+adg sg41=fg*exp(ppg*qpg41) qpp42=qf42+adp sp42=fp*exp(ppp*qpp42) qpg42=qf42+adg sg42=fg*exp(ppg*qpg42) qpp43=qf43+adp sp43=fp*exp(ppp*qpp43) qpg43=qf43+adg sg43=fg*exp(ppg*qpg43) qpp44=qf44+adp sp44=fp*exp(ppp*qpp44) qpg44=qf44+adg sg44=fg*exp(ppg*qpg44) qpp45=qf45+adp sp45=fp*exp(ppp*qpp45) qpg45=qf45+adg sg45=fg*exp(ppg*qpg45) qpp46=qf46+adp sp46=fp*exp(ppp*qpp46) qpg46=qf46+adg sg46=fg*exp(ppg*qpg46) qpp47=qf47+adp sp47=fp*exp(ppp*qpp47) qpg47=qf47+adg sg47=fg*exp(ppg*qpg47) qpp48=qf48+adp sp48=fp*exp(ppp*qpp48) qpg48=qf48+adg sg48=fg*exp(ppg*qpg48) qpp49=qf49+adp sp49=fp*exp(ppp*qpp49) qpg49=qf49+adg sg49=fg*exp(ppg*qpg49) qpp50=qf50+adp sp50=fp*exp(ppp*qpp50) qpg50=qf50+adg sg50=fg*exp(ppg*qpg50) qpp51=qf51+adp sp51=fp*exp(ppp*qpp51) qpg51=qf51+adg sg51=fg*exp(ppg*qpg51) qpp52=qf52+adp sp52=fp*exp(ppp*qpp52) qpg52=qf52+adg sg52=fg*exp(ppg*qpg52) qpp53=qf53+adp sp53=fp*exp(ppp*qpp53) qpg53=qf53+adg sg53=fg*exp(ppg*qpg53) qpp54=qf54+adp sp54=fp*exp(ppp*qpp54) qpg54=qf54+adg sg54=fg*exp(ppg*qpg54) qpp55=qf55+adp sp55=fp*exp(ppp*qpp55) qpg55=qf55+adg sg55=fg*exp(ppg*qpg55) qpp56=qf56+adp sp56=fp*exp(ppp*qpp56) qpg56=qf56+adg sg56=fg*exp(ppg*qpg56) qpp57=qf57+adp sp57=fp*exp(ppp*qpp57) qpg57=qf57+adg sg57=fg*exp(ppg*qpg57) qpp58=qf58+adp sp58=fp*exp(ppp*qpp58) qpg58=qf58+adg sg58=fg*exp(ppg*qpg58) qpp59=qf59+adp sp59=fp*exp(ppp*qpp59) qpg59=qf59+adg sg59=fg*exp(ppg*qpg59) qpp60=qf60+adp sp60=fp*exp(ppp*qpp60) qpg60=qf60+adg sg60=fg*exp(ppg*qpg60) qpp61=qf61+adp sp61=fp*exp(ppp*qpp61) qpg61=qf61+adg sg61=fg*exp(ppg*qpg61) qpp62=qf62+adp sp62=fp*exp(ppp*qpp62) qpg62=qf62+adg sg62=fg*exp(ppg*qpg62) qpp63=qf63+adp sp63=fp*exp(ppp*qpp63) qpg63=qf63+adg sg63=fg*exp(ppg*qpg63) qpp64=qf64+adp sp64=fp*exp(ppp*qpp64) qpg64=qf64+adg sg64=fg*exp(ppg*qpg64) rim1=Nos*(m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m32+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim1*stscalar rim2=Nos*(m1+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m32+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim2*stscalar rim3=Nos*(m1+m2+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m32+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim3*stscalar rim4=Nos*(m1+m2+m3+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m32+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim4*stscalar rim5=Nos*(m1+m2+m3+m4+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m32+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim5*stscalar rim6=Nos*(m1+m2+m3+m4+m5+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m32+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim6*stscalar rim7=Nos*(m1+m2+m3+m4+m5+m6+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m32+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim7*stscalar rim8=Nos*(m1+m2+m3+m4+m5+m6+m7+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m32+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim8*stscalar rim9=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m32+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim9*stscalar rim10=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m32+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim10*stscalar rim11=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m32+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim11*stscalar rim12=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m32+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim12*stscalar rim13=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m32+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim13*stscalar rim14=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m32+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim14*stscalar rim15=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m32+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim15*stscalar rim16=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m32+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim16*stscalar rim17=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m32+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim17*stscalar rim18=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m32+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim18*stscalar rim19=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m32+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim19*stscalar rim20=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m32+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim20*stscalar rim21=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m32+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim21*stscalar rim22=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m23+m24+m25+m26+m27+m28+m29+m30+m31+m32+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim22*stscalar rim23=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m24+m25+m26+m27+m28+m29+m30+m31+m32+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim23*stscalar rim24=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m25+m26+m27+m28+m29+m30+m31+m32+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim24*stscalar rim25=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m26+m27+m28+m29+m30+m31+m32+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim25*stscalar rim26=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m27+m28+m29+m30+m31+m32+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim26*stscalar rim27=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m28+m29+m30+m31+m32+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim27*stscalar rim28=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m29+m30+m31+m32+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim28*stscalar rim29=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m30+m31+m32+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim29*stscalar rim30=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m31+m32+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim30*stscalar rim31=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m32+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim31*stscalar rim32=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim32*stscalar rim33=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m32+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim33*stscalar rim34=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m33+m32+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim34*stscalar rim35=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m33+m34+m32+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim35*stscalar rim36=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m33+m34+m35+m32+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim36*stscalar rim37=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m33+m34+m35+m36+m32+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim37*stscalar rim38=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m33+m34+m35+m36+m37+m32+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim38*stscalar rim39=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m33+m34+m35+m36+m37+m38+m32+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim39*stscalar rim40=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m33+m34+m35+m36+m37+m38+m39+m32+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim40*stscalar rim41=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m33+m34+m35+m36+m37+m38+m39+m40+m32+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim41*stscalar rim42=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m33+m34+m35+m36+m37+m38+m39+m40+m41+m32+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim42*stscalar rim43=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m32+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim43*stscalar rim44=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m32+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim44*stscalar rim45=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m32+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim45*stscalar rim46=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m32+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim46*stscalar rim47=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m32+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim47*stscalar rim48=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m32+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim48*stscalar rim49=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m32+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim49*stscalar rim50=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m32+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim50*stscalar rim51=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m32+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim51*stscalar rim52=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m32+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim52*stscalar rim53=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m32+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim53*stscalar rim54=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m32+m55+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim54*stscalar rim55=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m32+m56+m57+m58+m59+m60+m61+m62+m63+m64)+Erim55*stscalar rim56=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m32+m57+m58+m59+m60+m61+m62+m63+m64)+Erim56*stscalar rim57=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m32+m58+m59+m60+m61+m62+m63+m64)+Erim57*stscalar rim58=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m32+m59+m60+m61+m62+m63+m64)+Erim58*stscalar rim59=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m32+m60+m61+m62+m63+m64)+Erim59*stscalar rim60=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m32+m61+m62+m63+m64)+Erim60*stscalar rim61=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m32+m62+m63+m64)+Erim61*stscalar rim62=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m32+m63+m64)+Erim62*stscalar rim63=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m32+m64)+Erim63*stscalar rim64=Nos*(m1+m2+m3+m4+m5+m6+m7+m8+m9+m10+m11+m12+m13+m14+m15+m16+m17+m18+m19+m20+m21+m22+m23+m24+m25+m26+m27+m28+m29+m30+m31+m33+m34+m35+m36+m37+m38+m39+m40+m41+m42+m43+m44+m45+m46+m47+m48+m49+m50+m51+m52+m53+m54+m55+m56+m57+m58+m59+m60+m61+m62+m63+m32)+Erim64*stscalar m1'=kg*(f1*g1)/(Kg+f1)-(m1*p1)/(1+m1)+rim1 f1'=-kg*(f1*g1)/(Kg+f1)+(m1*p1)/(1+m1)-(f1*p1)/(1+f1) p1'=sp1*alpha*(f1^n)/(Kr^n+f1^n)-ka*p1^2 g1'=sg1*beta*(f1^l)/(Ks^l+f1^l)-kd*(g1*p1)/(Kd+g1) m2'=kg*(f2*g2)/(Kg+f2)-(m2*p2)/(1+m2)+rim2 f2'=-kg*(f2*g2)/(Kg+f2)+(m2*p2)/(1+m2)-(f2*p2)/(1+f2) p2'=sp2*alpha*(f2^n)/(Kr^n+f2^n)-ka*p2^2 g2'=sg2*beta*(f2^l)/(Ks^l+f2^l)-kd*(g2*p2)/(Kd+g2) m3'=kg*(f3*g3)/(Kg+f3)-(m3*p3)/(1+m3)+rim3 f3'=-kg*(f3*g3)/(Kg+f3)+(m3*p3)/(1+m3)-(f3*p3)/(1+f3) p3'=sp3*alpha*(f3^n)/(Kr^n+f3^n)-ka*p3^2 g3'=sg3*beta*(f3^l)/(Ks^l+f3^l)-kd*(g3*p3)/(Kd+g3) m4'=kg*(f4*g4)/(Kg+f4)-(m4*p4)/(1+m4)+rim4 f4'=-kg*(f4*g4)/(Kg+f4)+(m4*p4)/(1+m4)-(f4*p4)/(1+f4) p4'=sp4*alpha*(f4^n)/(Kr^n+f4^n)-ka*p4^2 g4'=sg4*beta*(f4^l)/(Ks^l+f4^l)-kd*(g4*p4)/(Kd+g4) m5'=kg*(f5*g5)/(Kg+f5)-(m5*p5)/(1+m5)+rim5 f5'=-kg*(f5*g5)/(Kg+f5)+(m5*p5)/(1+m5)-(f5*p5)/(1+f5) p5'=sp5*alpha*(f5^n)/(Kr^n+f5^n)-ka*p5^2 g5'=sg5*beta*(f5^l)/(Ks^l+f5^l)-kd*(g5*p5)/(Kd+g5) m6'=kg*(f6*g6)/(Kg+f6)-(m6*p6)/(1+m6)+rim6 f6'=-kg*(f6*g6)/(Kg+f6)+(m6*p6)/(1+m6)-(f6*p6)/(1+f6) p6'=sp6*alpha*(f6^n)/(Kr^n+f6^n)-ka*p6^2 g6'=sg6*beta*(f6^l)/(Ks^l+f6^l)-kd*(g6*p6)/(Kd+g6) m7'=kg*(f7*g7)/(Kg+f7)-(m7*p7)/(1+m7)+rim7 f7'=-kg*(f7*g7)/(Kg+f7)+(m7*p7)/(1+m7)-(f7*p7)/(1+f7) p7'=sp7*alpha*(f7^n)/(Kr^n+f7^n)-ka*p7^2 g7'=sg7*beta*(f7^l)/(Ks^l+f7^l)-kd*(g7*p7)/(Kd+g7) m8'=kg*(f8*g8)/(Kg+f8)-(m8*p8)/(1+m8)+rim8 f8'=-kg*(f8*g8)/(Kg+f8)+(m8*p8)/(1+m8)-(f8*p8)/(1+f8) p8'=sp8*alpha*(f8^n)/(Kr^n+f8^n)-ka*p8^2 g8'=sg8*beta*(f8^l)/(Ks^l+f8^l)-kd*(g8*p8)/(Kd+g8) m9'=kg*(f9*g9)/(Kg+f9)-(m9*p9)/(1+m9)+rim9 f9'=-kg*(f9*g9)/(Kg+f9)+(m9*p9)/(1+m9)-(f9*p9)/(1+f9) p9'=sp9*alpha*(f9^n)/(Kr^n+f9^n)-ka*p9^2 g9'=sg9*beta*(f9^l)/(Ks^l+f9^l)-kd*(g9*p9)/(Kd+g9) m10'=kg*(f10*g10)/(Kg+f10)-(m10*p10)/(1+m10)+rim10 f10'=-kg*(f10*g10)/(Kg+f10)+(m10*p10)/(1+m10)-(f10*p10)/(1+f10) p10'=sp10*alpha*(f10^n)/(Kr^n+f10^n)-ka*p10^2 g10'=sg10*beta*(f10^l)/(Ks^l+f10^l)-kd*(g10*p10)/(Kd+g10) m11'=kg*(f11*g11)/(Kg+f11)-(m11*p11)/(1+m11)+rim11 f11'=-kg*(f11*g11)/(Kg+f11)+(m11*p11)/(1+m11)-(f11*p11)/(1+f11) p11'=sp11*alpha*(f11^n)/(Kr^n+f11^n)-ka*p11^2 g11'=sg11*beta*(f11^l)/(Ks^l+f11^l)-kd*(g11*p11)/(Kd+g11) m12'=kg*(f12*g12)/(Kg+f12)-(m12*p12)/(1+m12)+rim12 f12'=-kg*(f12*g12)/(Kg+f12)+(m12*p12)/(1+m12)-(f12*p12)/(1+f12) p12'=sp12*alpha*(f12^n)/(Kr^n+f12^n)-ka*p12^2 g12'=sg12*beta*(f12^l)/(Ks^l+f12^l)-kd*(g12*p12)/(Kd+g12) m13'=kg*(f13*g13)/(Kg+f13)-(m13*p13)/(1+m13)+rim13 f13'=-kg*(f13*g13)/(Kg+f13)+(m13*p13)/(1+m13)-(f13*p13)/(1+f13) p13'=sp13*alpha*(f13^n)/(Kr^n+f13^n)-ka*p13^2 g13'=sg13*beta*(f13^l)/(Ks^l+f13^l)-kd*(g13*p13)/(Kd+g13) m14'=kg*(f14*g14)/(Kg+f14)-(m14*p14)/(1+m14)+rim14 f14'=-kg*(f14*g14)/(Kg+f14)+(m14*p14)/(1+m14)-(f14*p14)/(1+f14) p14'=sp14*alpha*(f14^n)/(Kr^n+f14^n)-ka*p14^2 g14'=sg14*beta*(f14^l)/(Ks^l+f14^l)-kd*(g14*p14)/(Kd+g14) m15'=kg*(f15*g15)/(Kg+f15)-(m15*p15)/(1+m15)+rim15 f15'=-kg*(f15*g15)/(Kg+f15)+(m15*p15)/(1+m15)-(f15*p15)/(1+f15) p15'=sp15*alpha*(f15^n)/(Kr^n+f15^n)-ka*p15^2 g15'=sg15*beta*(f15^l)/(Ks^l+f15^l)-kd*(g15*p15)/(Kd+g15) m16'=kg*(f16*g16)/(Kg+f16)-(m16*p16)/(1+m16)+rim16 f16'=-kg*(f16*g16)/(Kg+f16)+(m16*p16)/(1+m16)-(f16*p16)/(1+f16) p16'=sp16*alpha*(f16^n)/(Kr^n+f16^n)-ka*p16^2 g16'=sg16*beta*(f16^l)/(Ks^l+f16^l)-kd*(g16*p16)/(Kd+g16) m17'=kg*(f17*g17)/(Kg+f17)-(m17*p17)/(1+m17)+rim17 f17'=-kg*(f17*g17)/(Kg+f17)+(m17*p17)/(1+m17)-(f17*p17)/(1+f17) p17'=sp17*alpha*(f17^n)/(Kr^n+f17^n)-ka*p17^2 g17'=sg17*beta*(f17^l)/(Ks^l+f17^l)-kd*(g17*p17)/(Kd+g17) m18'=kg*(f18*g18)/(Kg+f18)-(m18*p18)/(1+m18)+rim18 f18'=-kg*(f18*g18)/(Kg+f18)+(m18*p18)/(1+m18)-(f18*p18)/(1+f18) p18'=sp18*alpha*(f18^n)/(Kr^n+f18^n)-ka*p18^2 g18'=sg18*beta*(f18^l)/(Ks^l+f18^l)-kd*(g18*p18)/(Kd+g18) m19'=kg*(f19*g19)/(Kg+f19)-(m19*p19)/(1+m19)+rim19 f19'=-kg*(f19*g19)/(Kg+f19)+(m19*p19)/(1+m19)-(f19*p19)/(1+f19) p19'=sp19*alpha*(f19^n)/(Kr^n+f19^n)-ka*p19^2 g19'=sg19*beta*(f19^l)/(Ks^l+f19^l)-kd*(g19*p19)/(Kd+g19) m20'=kg*(f20*g20)/(Kg+f20)-(m20*p20)/(1+m20)+rim20 f20'=-kg*(f20*g20)/(Kg+f20)+(m20*p20)/(1+m20)-(f20*p20)/(1+f20) p20'=sp20*alpha*(f20^n)/(Kr^n+f20^n)-ka*p20^2 g20'=sg20*beta*(f20^l)/(Ks^l+f20^l)-kd*(g20*p20)/(Kd+g20) m21'=kg*(f21*g21)/(Kg+f21)-(m21*p21)/(1+m21)+rim21 f21'=-kg*(f21*g21)/(Kg+f21)+(m21*p21)/(1+m21)-(f21*p21)/(1+f21) p21'=sp21*alpha*(f21^n)/(Kr^n+f21^n)-ka*p21^2 g21'=sg21*beta*(f21^l)/(Ks^l+f21^l)-kd*(g21*p21)/(Kd+g21) m22'=kg*(f22*g22)/(Kg+f22)-(m22*p22)/(1+m22)+rim22 f22'=-kg*(f22*g22)/(Kg+f22)+(m22*p22)/(1+m22)-(f22*p22)/(1+f22) p22'=sp22*alpha*(f22^n)/(Kr^n+f22^n)-ka*p22^2 g22'=sg22*beta*(f22^l)/(Ks^l+f22^l)-kd*(g22*p22)/(Kd+g22) m23'=kg*(f23*g23)/(Kg+f23)-(m23*p23)/(1+m23)+rim23 f23'=-kg*(f23*g23)/(Kg+f23)+(m23*p23)/(1+m23)-(f23*p23)/(1+f23) p23'=sp23*alpha*(f23^n)/(Kr^n+f23^n)-ka*p23^2 g23'=sg23*beta*(f23^l)/(Ks^l+f23^l)-kd*(g23*p23)/(Kd+g23) m24'=kg*(f24*g24)/(Kg+f24)-(m24*p24)/(1+m24)+rim24 f24'=-kg*(f24*g24)/(Kg+f24)+(m24*p24)/(1+m24)-(f24*p24)/(1+f24) p24'=sp24*alpha*(f24^n)/(Kr^n+f24^n)-ka*p24^2 g24'=sg24*beta*(f24^l)/(Ks^l+f24^l)-kd*(g24*p24)/(Kd+g24) m25'=kg*(f25*g25)/(Kg+f25)-(m25*p25)/(1+m25)+rim25 f25'=-kg*(f25*g25)/(Kg+f25)+(m25*p25)/(1+m25)-(f25*p25)/(1+f25) p25'=sp25*alpha*(f25^n)/(Kr^n+f25^n)-ka*p25^2 g25'=sg25*beta*(f25^l)/(Ks^l+f25^l)-kd*(g25*p25)/(Kd+g25) m26'=kg*(f26*g26)/(Kg+f26)-(m26*p26)/(1+m26)+rim26 f26'=-kg*(f26*g26)/(Kg+f26)+(m26*p26)/(1+m26)-(f26*p26)/(1+f26) p26'=sp26*alpha*(f26^n)/(Kr^n+f26^n)-ka*p26^2 g26'=sg26*beta*(f26^l)/(Ks^l+f26^l)-kd*(g26*p26)/(Kd+g26) m27'=kg*(f27*g27)/(Kg+f27)-(m27*p27)/(1+m27)+rim27 f27'=-kg*(f27*g27)/(Kg+f27)+(m27*p27)/(1+m27)-(f27*p27)/(1+f27) p27'=sp27*alpha*(f27^n)/(Kr^n+f27^n)-ka*p27^2 g27'=sg27*beta*(f27^l)/(Ks^l+f27^l)-kd*(g27*p27)/(Kd+g27) m28'=kg*(f28*g28)/(Kg+f28)-(m28*p28)/(1+m28)+rim28 f28'=-kg*(f28*g28)/(Kg+f28)+(m28*p28)/(1+m28)-(f28*p28)/(1+f28) p28'=sp28*alpha*(f28^n)/(Kr^n+f28^n)-ka*p28^2 g28'=sg28*beta*(f28^l)/(Ks^l+f28^l)-kd*(g28*p28)/(Kd+g28) m29'=kg*(f29*g29)/(Kg+f29)-(m29*p29)/(1+m29)+rim29 f29'=-kg*(f29*g29)/(Kg+f29)+(m29*p29)/(1+m29)-(f29*p29)/(1+f29) p29'=sp29*alpha*(f29^n)/(Kr^n+f29^n)-ka*p29^2 g29'=sg29*beta*(f29^l)/(Ks^l+f29^l)-kd*(g29*p29)/(Kd+g29) m30'=kg*(f30*g30)/(Kg+f30)-(m30*p30)/(1+m30)+rim30 f30'=-kg*(f30*g30)/(Kg+f30)+(m30*p30)/(1+m30)-(f30*p30)/(1+f30) p30'=sp30*alpha*(f30^n)/(Kr^n+f30^n)-ka*p30^2 g30'=sg30*beta*(f30^l)/(Ks^l+f30^l)-kd*(g30*p30)/(Kd+g30) m31'=kg*(f31*g31)/(Kg+f31)-(m31*p31)/(1+m31)+rim31 f31'=-kg*(f31*g31)/(Kg+f31)+(m31*p31)/(1+m31)-(f31*p31)/(1+f31) p31'=sp31*alpha*(f31^n)/(Kr^n+f31^n)-ka*p31^2 g31'=sg31*beta*(f31^l)/(Ks^l+f31^l)-kd*(g31*p31)/(Kd+g31) m32'=kg*(f32*g32)/(Kg+f32)-(m32*p32)/(1+m32)+rim32 f32'=-kg*(f32*g32)/(Kg+f32)+(m32*p32)/(1+m32)-(f32*p32)/(1+f32) p32'=sp32*alpha*(f32^n)/(Kr^n+f32^n)-ka*p32^2 g32'=sg32*beta*(f32^l)/(Ks^l+f32^l)-kd*(g32*p32)/(Kd+g32) m33'=kg*(f33*g33)/(Kg+f33)-(m33*p33)/(1+m33)+rim33 f33'=-kg*(f33*g33)/(Kg+f33)+(m33*p33)/(1+m33)-(f33*p33)/(1+f33) p33'=sp33*alpha*(f33^n)/(Kr^n+f33^n)-ka*p33^2 g33'=sg33*beta*(f33^l)/(Ks^l+f33^l)-kd*(g33*p33)/(Kd+g33) m34'=kg*(f34*g34)/(Kg+f34)-(m34*p34)/(1+m34)+rim34 f34'=-kg*(f34*g34)/(Kg+f34)+(m34*p34)/(1+m34)-(f34*p34)/(1+f34) p34'=sp34*alpha*(f34^n)/(Kr^n+f34^n)-ka*p34^2 g34'=sg34*beta*(f34^l)/(Ks^l+f34^l)-kd*(g34*p34)/(Kd+g34) m35'=kg*(f35*g35)/(Kg+f35)-(m35*p35)/(1+m35)+rim35 f35'=-kg*(f35*g35)/(Kg+f35)+(m35*p35)/(1+m35)-(f35*p35)/(1+f35) p35'=sp35*alpha*(f35^n)/(Kr^n+f35^n)-ka*p35^2 g35'=sg35*beta*(f35^l)/(Ks^l+f35^l)-kd*(g35*p35)/(Kd+g35) m36'=kg*(f36*g36)/(Kg+f36)-(m36*p36)/(1+m36)+rim36 f36'=-kg*(f36*g36)/(Kg+f36)+(m36*p36)/(1+m36)-(f36*p36)/(1+f36) p36'=sp36*alpha*(f36^n)/(Kr^n+f36^n)-ka*p36^2 g36'=sg36*beta*(f36^l)/(Ks^l+f36^l)-kd*(g36*p36)/(Kd+g36) m37'=kg*(f37*g37)/(Kg+f37)-(m37*p37)/(1+m37)+rim37 f37'=-kg*(f37*g37)/(Kg+f37)+(m37*p37)/(1+m37)-(f37*p37)/(1+f37) p37'=sp37*alpha*(f37^n)/(Kr^n+f37^n)-ka*p37^2 g37'=sg37*beta*(f37^l)/(Ks^l+f37^l)-kd*(g37*p37)/(Kd+g37) m38'=kg*(f38*g38)/(Kg+f38)-(m38*p38)/(1+m38)+rim38 f38'=-kg*(f38*g38)/(Kg+f38)+(m38*p38)/(1+m38)-(f38*p38)/(1+f38) p38'=sp38*alpha*(f38^n)/(Kr^n+f38^n)-ka*p38^2 g38'=sg38*beta*(f38^l)/(Ks^l+f38^l)-kd*(g38*p38)/(Kd+g38) m39'=kg*(f39*g39)/(Kg+f39)-(m39*p39)/(1+m39)+rim39 f39'=-kg*(f39*g39)/(Kg+f39)+(m39*p39)/(1+m39)-(f39*p39)/(1+f39) p39'=sp39*alpha*(f39^n)/(Kr^n+f39^n)-ka*p39^2 g39'=sg39*beta*(f39^l)/(Ks^l+f39^l)-kd*(g39*p39)/(Kd+g39) m40'=kg*(f40*g40)/(Kg+f40)-(m40*p40)/(1+m40)+rim40 f40'=-kg*(f40*g40)/(Kg+f40)+(m40*p40)/(1+m40)-(f40*p40)/(1+f40) p40'=sp40*alpha*(f40^n)/(Kr^n+f40^n)-ka*p40^2 g40'=sg40*beta*(f40^l)/(Ks^l+f40^l)-kd*(g40*p40)/(Kd+g40) m41'=kg*(f41*g41)/(Kg+f41)-(m41*p41)/(1+m41)+rim41 f41'=-kg*(f41*g41)/(Kg+f41)+(m41*p41)/(1+m41)-(f41*p41)/(1+f41) p41'=sp41*alpha*(f41^n)/(Kr^n+f41^n)-ka*p41^2 g41'=sg41*beta*(f41^l)/(Ks^l+f41^l)-kd*(g41*p41)/(Kd+g41) m42'=kg*(f42*g42)/(Kg+f42)-(m42*p42)/(1+m42)+rim42 f42'=-kg*(f42*g42)/(Kg+f42)+(m42*p42)/(1+m42)-(f42*p42)/(1+f42) p42'=sp42*alpha*(f42^n)/(Kr^n+f42^n)-ka*p42^2 g42'=sg42*beta*(f42^l)/(Ks^l+f42^l)-kd*(g42*p42)/(Kd+g42) m43'=kg*(f43*g43)/(Kg+f43)-(m43*p43)/(1+m43)+rim43 f43'=-kg*(f43*g43)/(Kg+f43)+(m43*p43)/(1+m43)-(f43*p43)/(1+f43) p43'=sp43*alpha*(f43^n)/(Kr^n+f43^n)-ka*p43^2 g43'=sg43*beta*(f43^l)/(Ks^l+f43^l)-kd*(g43*p43)/(Kd+g43) m44'=kg*(f44*g44)/(Kg+f44)-(m44*p44)/(1+m44)+rim44 f44'=-kg*(f44*g44)/(Kg+f44)+(m44*p44)/(1+m44)-(f44*p44)/(1+f44) p44'=sp44*alpha*(f44^n)/(Kr^n+f44^n)-ka*p44^2 g44'=sg44*beta*(f44^l)/(Ks^l+f44^l)-kd*(g44*p44)/(Kd+g44) m45'=kg*(f45*g45)/(Kg+f45)-(m45*p45)/(1+m45)+rim45 f45'=-kg*(f45*g45)/(Kg+f45)+(m45*p45)/(1+m45)-(f45*p45)/(1+f45) p45'=sp45*alpha*(f45^n)/(Kr^n+f45^n)-ka*p45^2 g45'=sg45*beta*(f45^l)/(Ks^l+f45^l)-kd*(g45*p45)/(Kd+g45) m46'=kg*(f46*g46)/(Kg+f46)-(m46*p46)/(1+m46)+rim46 f46'=-kg*(f46*g46)/(Kg+f46)+(m46*p46)/(1+m46)-(f46*p46)/(1+f46) p46'=sp46*alpha*(f46^n)/(Kr^n+f46^n)-ka*p46^2 g46'=sg46*beta*(f46^l)/(Ks^l+f46^l)-kd*(g46*p46)/(Kd+g46) m47'=kg*(f47*g47)/(Kg+f47)-(m47*p47)/(1+m47)+rim47 f47'=-kg*(f47*g47)/(Kg+f47)+(m47*p47)/(1+m47)-(f47*p47)/(1+f47) p47'=sp47*alpha*(f47^n)/(Kr^n+f47^n)-ka*p47^2 g47'=sg47*beta*(f47^l)/(Ks^l+f47^l)-kd*(g47*p47)/(Kd+g47) m48'=kg*(f48*g48)/(Kg+f48)-(m48*p48)/(1+m48)+rim48 f48'=-kg*(f48*g48)/(Kg+f48)+(m48*p48)/(1+m48)-(f48*p48)/(1+f48) p48'=sp48*alpha*(f48^n)/(Kr^n+f48^n)-ka*p48^2 g48'=sg48*beta*(f48^l)/(Ks^l+f48^l)-kd*(g48*p48)/(Kd+g48) m49'=kg*(f49*g49)/(Kg+f49)-(m49*p49)/(1+m49)+rim49 f49'=-kg*(f49*g49)/(Kg+f49)+(m49*p49)/(1+m49)-(f49*p49)/(1+f49) p49'=sp49*alpha*(f49^n)/(Kr^n+f49^n)-ka*p49^2 g49'=sg49*beta*(f49^l)/(Ks^l+f49^l)-kd*(g49*p49)/(Kd+g49) m50'=kg*(f50*g50)/(Kg+f50)-(m50*p50)/(1+m50)+rim50 f50'=-kg*(f50*g50)/(Kg+f50)+(m50*p50)/(1+m50)-(f50*p50)/(1+f50) p50'=sp50*alpha*(f50^n)/(Kr^n+f50^n)-ka*p50^2 g50'=sg50*beta*(f50^l)/(Ks^l+f50^l)-kd*(g50*p50)/(Kd+g50) m51'=kg*(f51*g51)/(Kg+f51)-(m51*p51)/(1+m51)+rim51 f51'=-kg*(f51*g51)/(Kg+f51)+(m51*p51)/(1+m51)-(f51*p51)/(1+f51) p51'=sp51*alpha*(f51^n)/(Kr^n+f51^n)-ka*p51^2 g51'=sg51*beta*(f51^l)/(Ks^l+f51^l)-kd*(g51*p51)/(Kd+g51) m52'=kg*(f52*g52)/(Kg+f52)-(m52*p52)/(1+m52)+rim52 f52'=-kg*(f52*g52)/(Kg+f52)+(m52*p52)/(1+m52)-(f52*p52)/(1+f52) p52'=sp52*alpha*(f52^n)/(Kr^n+f52^n)-ka*p52^2 g52'=sg52*beta*(f52^l)/(Ks^l+f52^l)-kd*(g52*p52)/(Kd+g52) m53'=kg*(f53*g53)/(Kg+f53)-(m53*p53)/(1+m53)+rim53 f53'=-kg*(f53*g53)/(Kg+f53)+(m53*p53)/(1+m53)-(f53*p53)/(1+f53) p53'=sp53*alpha*(f53^n)/(Kr^n+f53^n)-ka*p53^2 g53'=sg53*beta*(f53^l)/(Ks^l+f53^l)-kd*(g53*p53)/(Kd+g53) m54'=kg*(f54*g54)/(Kg+f54)-(m54*p54)/(1+m54)+rim54 f54'=-kg*(f54*g54)/(Kg+f54)+(m54*p54)/(1+m54)-(f54*p54)/(1+f54) p54'=sp54*alpha*(f54^n)/(Kr^n+f54^n)-ka*p54^2 g54'=sg54*beta*(f54^l)/(Ks^l+f54^l)-kd*(g54*p54)/(Kd+g54) m55'=kg*(f55*g55)/(Kg+f55)-(m55*p55)/(1+m55)+rim55 f55'=-kg*(f55*g55)/(Kg+f55)+(m55*p55)/(1+m55)-(f55*p55)/(1+f55) p55'=sp55*alpha*(f55^n)/(Kr^n+f55^n)-ka*p55^2 g55'=sg55*beta*(f55^l)/(Ks^l+f55^l)-kd*(g55*p55)/(Kd+g55) m56'=kg*(f56*g56)/(Kg+f56)-(m56*p56)/(1+m56)+rim56 f56'=-kg*(f56*g56)/(Kg+f56)+(m56*p56)/(1+m56)-(f56*p56)/(1+f56) p56'=sp56*alpha*(f56^n)/(Kr^n+f56^n)-ka*p56^2 g56'=sg56*beta*(f56^l)/(Ks^l+f56^l)-kd*(g56*p56)/(Kd+g56) m57'=kg*(f57*g57)/(Kg+f57)-(m57*p57)/(1+m57)+rim57 f57'=-kg*(f57*g57)/(Kg+f57)+(m57*p57)/(1+m57)-(f57*p57)/(1+f57) p57'=sp57*alpha*(f57^n)/(Kr^n+f57^n)-ka*p57^2 g57'=sg57*beta*(f57^l)/(Ks^l+f57^l)-kd*(g57*p57)/(Kd+g57) m58'=kg*(f58*g58)/(Kg+f58)-(m58*p58)/(1+m58)+rim58 f58'=-kg*(f58*g58)/(Kg+f58)+(m58*p58)/(1+m58)-(f58*p58)/(1+f58) p58'=sp58*alpha*(f58^n)/(Kr^n+f58^n)-ka*p58^2 g58'=sg58*beta*(f58^l)/(Ks^l+f58^l)-kd*(g58*p58)/(Kd+g58) m59'=kg*(f59*g59)/(Kg+f59)-(m59*p59)/(1+m59)+rim59 f59'=-kg*(f59*g59)/(Kg+f59)+(m59*p59)/(1+m59)-(f59*p59)/(1+f59) p59'=sp59*alpha*(f59^n)/(Kr^n+f59^n)-ka*p59^2 g59'=sg59*beta*(f59^l)/(Ks^l+f59^l)-kd*(g59*p59)/(Kd+g59) m60'=kg*(f60*g60)/(Kg+f60)-(m60*p60)/(1+m60)+rim60 f60'=-kg*(f60*g60)/(Kg+f60)+(m60*p60)/(1+m60)-(f60*p60)/(1+f60) p60'=sp60*alpha*(f60^n)/(Kr^n+f60^n)-ka*p60^2 g60'=sg60*beta*(f60^l)/(Ks^l+f60^l)-kd*(g60*p60)/(Kd+g60) m61'=kg*(f61*g61)/(Kg+f61)-(m61*p61)/(1+m61)+rim61 f61'=-kg*(f61*g61)/(Kg+f61)+(m61*p61)/(1+m61)-(f61*p61)/(1+f61) p61'=sp61*alpha*(f61^n)/(Kr^n+f61^n)-ka*p61^2 g61'=sg61*beta*(f61^l)/(Ks^l+f61^l)-kd*(g61*p61)/(Kd+g61) m62'=kg*(f62*g62)/(Kg+f62)-(m62*p62)/(1+m62)+rim62 f62'=-kg*(f62*g62)/(Kg+f62)+(m62*p62)/(1+m62)-(f62*p62)/(1+f62) p62'=sp62*alpha*(f62^n)/(Kr^n+f62^n)-ka*p62^2 g62'=sg62*beta*(f62^l)/(Ks^l+f62^l)-kd*(g62*p62)/(Kd+g62) m63'=kg*(f63*g63)/(Kg+f63)-(m63*p63)/(1+m63)+rim63 f63'=-kg*(f63*g63)/(Kg+f63)+(m63*p63)/(1+m63)-(f63*p63)/(1+f63) p63'=sp63*alpha*(f63^n)/(Kr^n+f63^n)-ka*p63^2 g63'=sg63*beta*(f63^l)/(Ks^l+f63^l)-kd*(g63*p63)/(Kd+g63) m64'=kg*(f64*g64)/(Kg+f64)-(m64*p64)/(1+m64)+rim64 f64'=-kg*(f64*g64)/(Kg+f64)+(m64*p64)/(1+m64)-(f64*p64)/(1+f64) p64'=sp64*alpha*(f64^n)/(Kr^n+f64^n)-ka*p64^2 g64'=sg64*beta*(f64^l)/(Ks^l+f64^l)-kd*(g64*p64)/(Kd+g64) t'=1 m1=7.26 m2=7.26 m3=7.26 m4=7.26 m5=7.26 m6=7.26 m7=7.26 m8=7.26 m9=7.26 m10=7.26 m11=7.26 m12=7.26 m13=7.26 m14=7.26 m15=7.26 m16=7.26 m17=7.26 m18=7.26 m19=7.26 m20=7.26 m21=7.26 m22=7.26 m23=7.26 m24=7.26 m25=7.26 m26=7.26 m27=7.26 m28=7.26 m29=7.26 m30=7.26 m31=7.26 m32=7.26 m33=7.26 m34=7.26 m35=7.26 m36=7.26 m37=7.26 m38=7.26 m39=7.26 m40=7.26 m41=7.26 m42=7.26 m43=7.26 m44=7.26 m45=7.26 m46=7.26 m47=7.26 m48=7.26 m49=7.26 m50=7.26 m51=7.26 m52=7.26 m53=7.26 m54=7.26 m55=7.26 m56=7.26 m57=7.26 m58=7.26 m59=7.26 m60=7.26 m61=7.26 m62=7.26 m63=7.26 m64=7.26 mtotal=32 f1=0.465 f2=0.465 f3=0.465 f4=0.465 f5=0.465 f6=0.465 f7=0.465 f8=0.465 f9=0.465 f10=0.465 f11=0.465 f12=0.465 f13=0.465 f14=0.465 f15=0.465 f16=0.465 f17=0.465 f18=0.465 f19=0.465 f20=0.465 f21=0.465 f22=0.465 f23=0.465 f24=0.465 f25=0.465 f26=0.465 f27=0.465 f28=0.465 f29=0.465 f30=0.465 f31=0.465 f32=0.465 f33=0.465 f34=0.465 f35=0.465 f36=0.465 f37=0.465 f38=0.465 f39=0.465 f40=0.465 f41=0.465 f42=0.465 f43=0.465 f44=0.465 f45=0.465 f46=0.465 f47=0.465 f48=0.465 f49=0.465 f50=0.465 f51=0.465 f52=0.465 f53=0.465 f54=0.465 f55=0.465 f56=0.465 f57=0.465 f58=0.465 f59=0.465 f60=0.465 f61=0.465 f62=0.465 f63=0.465 f64=0.465 ftotal=12 p1=0.016 g1=0.127 p2=0.016 g2=0.127 p3=0.016 g3=0.127 p4=0.016 g4=0.127 p5=0.016 g5=0.127 p6=0.016 g6=0.127 p7=0.016 g7=0.127 p8=0.016 g8=0.127 p9=0.016 g9=0.127 p10=0.016 g10=0.127 p11=0.016 g11=0.127 p12=0.016 g12=0.127 p13=0.016 g13=0.127 p14=0.016 g14=0.127 p15=0.016 g15=0.127 p16=0.016 g16=0.127 p17=0.016 g17=0.127 p18=0.016 g18=0.127 p19=0.016 g19=0.127 p20=0.016 g20=0.127 p21=0.016 g21=0.127 p22=0.016 g22=0.127 p23=0.016 g23=0.127 p24=0.016 g24=0.127 p25=0.016 g25=0.127 p26=0.016 g26=0.127 p27=0.016 g27=0.127 p28=0.016 g28=0.127 p29=0.016 g29=0.127 p30=0.016 g30=0.127 p31=0.016 g31=0.127 p32=0.016 g32=0.127 p33=0.016 g33=0.127 p34=0.016 g34=0.127 p35=0.016 g35=0.127 p36=0.016 g36=0.127 p37=0.016 g37=0.127 p38=0.016 g38=0.127 p39=0.016 g39=0.127 p40=0.016 g40=0.127 p41=0.016 g41=0.127 p42=0.016 g42=0.127 p43=0.016 g43=0.127 p44=0.016 g44=0.127 p45=0.016 g45=0.127 p46=0.016 g46=0.127 p47=0.016 g47=0.127 p48=0.016 g48=0.127 p49=0.016 g49=0.127 p50=0.016 g50=0.127 p51=0.016 g51=0.127 p52=0.016 g52=0.127 p53=0.016 g53=0.127 p54=0.016 g54=0.127 p55=0.016 g55=0.127 p56=0.016 g56=0.127 p57=0.016 g57=0.127 p58=0.016 g58=0.127 p59=0.016 g59=0.127 p60=0.016 g60=0.127 p61=0.016 g61=0.127 p62=0.016 g62=0.127 p63=0.016 g63=0.127 p64=0.016 g64=0.127 qpp1=0 sp1=1 qpg1=0 sg1=1 qpp2=0 sp2=1 qpg2=0 sg2=1 qpp3=0 sp3=1 qpg3=0 sg3=1 qpp4=0 sp4=1 qpg4=0 sg4=1 qpp5=0 sp5=1 qpg5=0 sg5=1 qpp6=0 sp6=1 qpg6=0 sg6=1 qpp7=0 sp7=1 qpg7=0 sg7=1 qpp8=0 sp8=1 qpg8=0 sg8=1 qpp9=0 sp9=1 qpg9=0 sg9=1 qpp10=0 sp10=1 qpg10=0 sg10=1 qpp11=0 sp11=1 qpg11=0 sg11=1 qpp12=0 sp12=1 qpg12=0 sg12=1 qpp13=0 sp13=1 qpg13=0 sg13=1 qpp14=0 sp14=1 qpg14=0 sg14=1 qpp15=0 sp15=1 qpg15=0 sg15=1 qpp16=0 sp16=1 qpg16=0 sg16=1 qpp17=0 sp17=1 qpg17=0 sg17=1 qpp18=0 sp18=1 qpg18=0 sg18=1 qpp19=0 sp19=1 qpg19=0 sg19=1 qpp20=0 sp20=1 qpg20=0 sg20=1 qpp21=0 sp21=1 qpg21=0 sg21=1 qpp22=0 sp22=1 qpg22=0 sg22=1 qpp23=0 sp23=1 qpg23=0 sg23=1 qpp24=0 sp24=1 qpg24=0 sg24=1 qpp25=0 sp25=1 qpg25=0 sg25=1 qpp26=0 sp26=1 qpg26=0 sg26=1 qpp27=0 sp27=1 qpg27=0 sg27=1 qpp28=0 sp28=1 qpg28=0 sg28=1 qpp29=0 sp29=1 qpg29=0 sg29=1 qpp30=0 sp30=1 qpg30=0 sg30=1 qpp31=0 sp31=1 qpg31=0 sg31=1 qpp32=0 sp32=1 qpg32=0 sg32=1 qpp33=0 sp33=1 qpg33=0 sg33=1 qpp34=0 sp34=1 qpg34=0 sg34=1 qpp35=0 sp35=1 qpg35=0 sg35=1 qpp36=0 sp36=1 qpg36=0 sg36=1 qpp37=0 sp37=1 qpg37=0 sg37=1 qpp38=0 sp38=1 qpg38=0 sg38=1 qpp39=0 sp39=1 qpg39=0 sg39=1 qpp40=0 sp40=1 qpg40=0 sg40=1 qpp41=0 sp41=1 qpg41=0 sg41=1 qpp42=0 sp42=1 qpg42=0 sg42=1 qpp43=0 sp43=1 qpg43=0 sg43=1 qpp44=0 sp44=1 qpg44=0 sg44=1 qpp45=0 sp45=1 qpg45=0 sg45=1 qpp46=0 sp46=1 qpg46=0 sg46=1 qpp47=0 sp47=1 qpg47=0 sg47=1 qpp48=0 sp48=1 qpg48=0 sg48=1 qpp49=0 sp49=1 qpg49=0 sg49=1 qpp50=0 sp50=1 qpg50=0 sg50=1 qpp51=0 sp51=1 qpg51=0 sg51=1 qpp52=0 sp52=1 qpg52=0 sg52=1 qpp53=0 sp53=1 qpg53=0 sg53=1 qpp54=0 sp54=1 qpg54=0 sg54=1 qpp55=0 sp55=1 qpg55=0 sg55=1 qpp56=0 sp56=1 qpg56=0 sg56=1 qpp57=0 sp57=1 qpg57=0 sg57=1 qpp58=0 sp58=1 qpg58=0 sg58=1 qpp59=0 sp59=1 qpg59=0 sg59=1 qpp60=0 sp60=1 qpg60=0 sg60=1 qpp61=0 sp61=1 qpg61=0 sg61=1 qpp62=0 sp62=1 qpg62=0 sg62=1 qpp63=0 sp63=1 qpg63=0 sg63=1 qpp64=0 sp64=1 qpg64=0 sg64=1 st=0 t=0 Integrator=Fehlberg-RK %Integrator=Cash-Karp-RK Step=0.5 %Duration=9999999 Duration=4999999 Output=3000 % Graph=Evolution,m1,m2,m3,m4 % Graph=Evolution,f1,f2,f3,f4 % Graph=Evolution,m5,m6,m7,m8 % Graph=Evolution,f5,f6,f7,f8 % Graph=Evolution,m9,m10,m11,m12 % Graph=Evolution,f9,f10,f11,f12 % Graph=Evolution,m13,m14,m15,m16 % Graph=Evolution,f13,f14,f15,f16 Graph=Evolution,stscalar %Graph=ftotal,ptotal %,gtotal % Graph=Evolution,g1,g2,g3,g4 Graph=Evolution,mtotal Graph=Evolution,ftotal % Graph=Evolution,sg1,sp1,sg2,sp2 % Graph=Evolution,sg3,sp3,sg4,sp4 Graph=ftotal,mtotal Lyapunov=false Statistics=false Logfile=berry64.1.log