From 919bc254074115776a4187e1ca4eac5d8fa49f35 Mon Sep 17 00:00:00 2001 From: Jan Benda Date: Fri, 22 Nov 2024 11:44:18 +0100 Subject: [PATCH] worked on Benjamins comments --- .gitignore | 2 +- model_full.pdf | Bin 271181 -> 271664 bytes model_full.py | 3 +-- susceptibility1.tex | 19 +++++++++---------- 4 files changed, 11 insertions(+), 13 deletions(-) diff --git a/.gitignore b/.gitignore index db23d13..8010db7 100644 --- a/.gitignore +++ b/.gitignore @@ -1,4 +1,4 @@ -__pychache__ +__pycache__/ *.blg *.aux *.bbl diff --git a/model_full.pdf b/model_full.pdf index aa4fd7070edec738b2a6c7d6f891b0c00488bf4d..8830c16cd3e09f75a3f826e601710602aa7ee588 100644 GIT binary patch delta 25924 zcmV)cK&ZdX#t^W|5U`K|e_OL=%TXYOIt($2z6)AuV&4AhO@dyZ+_1 zcUvGexF51RY@!|NCS zOMLz6|9<`5r@wpsf6Hnqy#7gA3wT7@qvTV|Z>`rEU%ANlp1;Spj9!2GC5HXqKmGF8 zpMLRQ<7=RUE%z42?5|({_fIvym6poPy@rHue+xB~m{L#q_3MB9^k@I&&tJd(J${++ z%g|;Uc(7#tvazM!XDVT$BYZf8Admm`uRlp2Kh5%1@+g=uf9sqRd`s^zLMhddYppH+ zvAv7W-o;ek&}Ern?Vgvs@8A2kpet{)_Euwo?1K(I-&ia9m20m1jZX4SZ!GgIygW^CMveC8}t@n=|`r<-%Wh1g3SS2vIJySh2sTPz9NIp^>Pbo1Q|o^Hl8 z-}>E06aC~Ij<=pm4l_mcDdP7yQd4L#xAf6Ub?*vsOlWe=Z6hwSU0j9^SC^|NYZ%URVV9 zQ+!LW-(o8NkvE*XM9$ou>*q|~Km1AF?$5LSJmL05XD_!#)9Y{deRU?netrc0*Y2h) z@9y04W*->P@V>%uamj?`_bYYww~%7YX@r3f&u^Xchn+U6%-Db5o2k8}-U?Rn#30Hp zOa(u~f94DM^%qzorOeV3HkYomukJ#eD>>mW1ONAW_q*TxxBv0mwMM?uY24;`Ausj( zhGQU3bgW^JLVRl>#5{RHyng+IXe~ef5vQs0R6Isa$Zw-$JXq1NYbky8tG?Q) z`4)SNcntf&H|0MRzk6$?^W7;^75i7t+(Y*Fe}D5{_OJi+)iWs9H|)bmg(J=vuRi<> zy8`_~I`B*P_WjEJo}XyE;B>&v5YUgwqkE1P5(W#`+_UYU|I4t-@#4_&9~`d(XXD6Y zR9v09hW_f0`?AVKfIA|+dU>0r=GZua@0$MEKdhyu7!CfF1sdxS~b~dA1boPWsiK^2?w|b?gdA0;jCw4Ctxz2Dfqj>p$Yt&^X2n z!<2B%U^PUn!-*X@A5HT-T#CH46>(%>vzIYi>Z&5tp3yMQBPw4|F8bu+8bg=XU#`Zs z8KeeQcxPbYT%K?EOiF4dITLA_NB@LR>DR~cbMfNqbAMcWNVxFyd30M$ zteJjHZE|mHn5fWA^F6OU;_rMn>^dm^qWm^@oqhCbT3HbXP`+P%#NYX9Xj`}|e?lG; zKX;eTNAIRd1iPR^f4N@nrVciuTKHjD$d1F ziXDsXQia}|LR;&A;SN=^j?#kIe;g6VeHyvK@Z@!W)Q;bq22;MHEu0Hj7f{lAy|kgv zKdo0KEQMOx30MMEme+*F*DwgUhCa}&zT+D?zJ~mGVbR5 zPX3eF0FF;Nyhbh+eX7_|gC@$rruu$(pgfIb*IAgiFsy4VkIo+X_jVQoG^iHM8FK$a zAeXOpmN)sxE#$keNBmPKe^rJ>1SJcnQ^0@WhkQSFAPlYt`1X|s&dwXK9CND}C$xnl zcTnLjP~UQCe&_w;HX@Q3v%i1rH9r@7tGQ41HuFdO;oTDby5anRqL^`3>nN5fZ>w*j zOXJ$HYx|c+@5FWrapEk7QK{PPJK5j=&3hSi-cJdaS;Bo7p^l-Oe|Z0Ba@%*J89B8} z4m18ORcH_f`?s13IsYkAu3cTgF-#$!v- z2G&uZVG_FbPmi}8e_`1_TJGpwfr9kt2jtK;OK4ZGAuX6$SOESDz0{A#EnlPOTRtQ0 z`{C!s_xjAD|H3olAFhwjyd-}4-r|Wb^sW)z$=HK(2I08E*$Ndq*HokpSQ$N+nP$6g z)B66m*yGs6#N@ke+W*NO&M#W8*fFp{O49=i1AyV?IsOipf58BTnsJy-C{VQN-uN8_ zpHc-}irjl)R)ZJfRu-OJ%Gk&S^ET=Bc=x;Cp9AvRXxCr4JAdS#xl!HUpSzoTHWNo{ zGt-Uc_bgAR)SwPsa!PxN#m{x7&IRQr$9rW^UP=zKC-yp`w3QI=HTRJ26?)0{dcIoT zYyawp+bOQke|uM5U$=MQH5P8aQl5HFvLoEn1oykWZa19JwN}Qn4_JSb=kMKPtQ{D-ESj>k?bqJ;Hw$us?VczMa!TtQkpOs?{ne^|(8q8)s#0vHZ3ughb_+|skY zy!@<2-||_xz8`*_`(B?-;lJ?gA3q*`x0Q}gVS4BNbjK=%272#f{Bhz+N%`3if@4^4 z_XzNLGN5k ze<3~8n9{JS>WKZ4Q(EzkeyK5(j!|ZKs4>u3-*8Oze5o;%(tTgAyZ9kJo5jn|*77Z% zZRq>qH?8mW+2Z~S&vD?#!yg`vOSz6d@2BINpo8UZIihzH7x1$m0sR{}rMuTS(`QWg z_ifoV7EZdy&9&*q*!AF=5!YWgf9~7de_VgGGQM_fyx#()$G6FN3urO@LmlSpJA!r% zUL6|U(eHP}|Ng)I+wZPuo(-}IB4(IegYHoJ$UT;vUmcTHAAvO*1M_=q9mdvZTxlFs zn-dOa2knhFbekGTq#lVyBX}_7wO3eGBFFflmnI(}^+vQB>LE0_UPV565HoHee^x?~ zSk=-F@j~cp^afyp!i;pDs(ccM>Jo-EVGc_e7ABh>1kSEF*U?#ld?_DT0GaJ@+c1(( zE}OzpAc8Kcf`#b&^I9M%r@OKXwMz+gwm@_M{u@^JQ0UK;GlE^@TVqi#UX zQ7$1TlNZc)WE2_oI6B@<-a|+ffA_bF{aS~KXX>Ffc^rm+6fa&*5)RlIN&gW^zme~UQ4v0dV1FcV>Alb4jRS@Aq$;?`hxdddf*f7LoIcI*+A zWBrBHFoKk052qPTUaNq6u|snTn9yW$ED_EX#5g=c*Nhx<#gaE1kg-bh z(3f*I@k9e9%Hvp>NmoLK5m^f1K+1986n!4@+{2hzdl?8jn!L^~bH#Q#mavHj)?d+& z#lu<0EMv7-SSv{$$NC<+fAAs4h;Wy2j4{?!UQ~`XElY>Zw>fnwK;E(IbV^`&4IQ&g zC&$1jh#XdLANo}oK&+OK#O;cCZ9SVDYN<>i+bJPh_tRVBX{znR9>&Ty`Difm)O|al zQ3RR8oaI|W7y0PyvdM7IR-EO?@6?5-`(If0oX~W#6L7Q~@=4?87J+c^o*ZbFX^XL0(Fp zd59dk;tcEg)`QI9$oy%aOODHv@&0(>TVPzRto%Zh2cy+Re~uGEw?D+W;vaMvW=Xn# zoOcyZGlp`m$46l7g~)M@7rkhTL{s#nZ=F1zv&jq8P?@YmCqNFvx_;TDER9*-1#lavxmE&4!T?Ypb1U(ed6GGt-Q+ld8n5}momy*La;y;UNqKJw z)*I~3&@m|5OAJSjCdUqoQM-v_iGNf%ZiPv$0Nugntj)w@bQn$=69sKT4;Rr&imrqr z_?6zakn^sUm4H7cdNoa68&2OMa@@Pwz6u&B#*&Zoe`1_dZ$m1q3A?F@9Oea14wEM; za!KSkOS4|)RC@`4=|-9SK{k003XmbPEQbuDWRQk()sN5naRKthYarw;AZr!%D&YOg(x9vxFLfL7su#q1LAI37tdc z7F!RWf2d-Z-O$cklHQM*cLaEaXc6iur z4moZoVtbS`?5LR-=R=7)m<>B(KaSO|bI4|ZoB$bQ^sJLXS$NiHoYo5Z3dZQQA`$6_hF5CL(X#uBqtS zOgK{PyQn)%opr1`*umDH#udqkB(CvH{CBm31E5oDqGk%aK(}JW*k(qBNcdNU!J_I8 z)N$AyDDcx(561+mlD3b_*R~>2 ztf8|zdFq=I2g@ozjffgQgDWd+d1mZcveCy4*i8;2y|-ah zZeqCUVt|n-UrP%;LT+I@A?PImF;-^i$&SsnW}8?bGPj&_jK*pO=6p{&HaK=n4~JP8 z=@Hc>0!@|)ODe~@$h}*R3>6-e!xl*)Tj$aDYI3NtA=$Cnrba?B_6qxUP$TA2YrU&o z068uRy;TM;Ejp%Ej(x&0<_f`xf9YWjYZR8$6s>7etI>cU%3-vPsAJFTOP2!VJcP0y zFwnl5hO|LWUx{cFlFlv zwjuT^=Wr7Ek?cX0SIk@2rBHGD>SfU6?Bf(p=w$dGRWBq7FL1r^Ld-PnObl%cBP;hR zboOR5i$$7o0u^ZnT9|bbJE@MSw|3egkgns*@YrT@*aT6O9iAb@j^0k)9)H!n&$qTE$1Yhj zSc&0kJD49#PZ_oV3Fnv`uL|83RMHbl4DY6Wr_!yeJne35sidokii^%iC*aj(o@Ba* zsunPS15=lD)pQt*db$z?VDzR|Y6GkKe(jch}lO)Kq;L4{#re?yrjOn;FVYSKCD zP(~b|X~}WH$ENM0J59GH@$_ON4Bg3Q97-n?7hO<9eXRLL%d$>=?O^;{H#zLaSL1^D(4xjq20JT@zT+f$rF|? zCy`ogj$zAj)tt{Kc2GChnL+lFOXqOFkD6xgDDJ+Y*8e~bsU#wEw&Y3inp&|y`7 zpaofZ7upz)qURao@A_n?<*QExO`=BU>vH!M3Zpn9ftFd zAXhdq@${}|PQe+XgGrOHgc)?DL48x@5MroY-Lw;&EOkB2E}C)@>S%f_(z?b8Ic(Fm zrWP;6&{R@rf7Y2E_0S!XZXdGUNHdsG#|#(PeO1Q>+Mpp=FwLMzwqXYv4!vw*?8j_6 z3pC=UyM&e*)iQ0SreJpiSr-4>L`_94g7sV1Gr z_B*4QyP@3&+h8!}cT7CKU6>T7x#tTLaG%m{wEl#{4Tz+O*^g^QS5Sg~L$u zgZX#fwSh7BUoOnQ5z#BA#Cy!26IX4UJIo*LfA{9i>!Yd2( zU(Nu!!SE>x*Li0^A^>N%r}qX>rLS2RgL4M3978dk4swQ%^O~VL0clvzZzn;16WeY% z6_%5L@H_UO_GT&)(xv%^8$u6o?>wLnqQ$YEcg_PI*>gD+w)cQ*xf!P8y$6*lX*1;6 ze|ZmZGK~c*dhdZ`{uoBzZTZtp?Bb~0Grc@MbmVq8v$|%A-_8Rj^Ci1yZs$QJKsHRiod-=02F2yP^8l-nH_CECIOhR0i>L>{?Iaj1H?f=t z+etvks4o~FJGyfvsT`A^mp8yW52(s5e{aBc9?&aSi(WOi_h2x}pqhMn4+f9Q9+z)o zdk+MXGTX{~4`@!8p-RW?J?JA8Rn(mKfC;qOARy;G5LB_001N!y|-k5xO4q0)t6zov+%e`|-3QR%wtTr27!8lbe;Ae< z0*AOvI=vQ{AibQtyMd}2W^2i}y8){jQuDcNcY^@)WqX*MTY;Nqtfhas74XuspWki; z0{qz0OYYnX%rWWj+zJHjM{UbdupJ0Izl!N8w*vvzM~TaUkf3=N!@{9t>Fwkl1-#UU zL0h+@py6P(=Pyq|CBi#R?mPv!e@)GUeR~RqxJymGJOx~jZEWPNa>cn8i?X&Xt;Di6*N z8ttZhcyOLDP|kU^v%vWo1-7#3W$5cv`y$aVqP=m@{9qXPEbkzAoIDITe^t+ggWw0D z(-IcezA<>-W;{kszA$(y#mN-88*8T>qfpq5wd2-od13RT_qj~=lGUN%A$Xz{&~dxT zHv})ZUr+1#jlmNq%%Q%+;F-9Rt(!Ln-&ja!Vd@)$@6>)YBm2VO=~{>ujNvhO#+&rn z7uL>F8bMKJoN;N1H#Fa=e@RROfgLnY!?)@mERQ_;O?Cg!t>s3-jLCIqj(xKyYK{)ZOYk?0nLXT2x{47&$#1O*Ib> zkrSH7QMRD-1(DMwFtWLI9YoHJ1$}onj}y1*oM6rj_BzP9I=!ate>aDmbFD#dYER^x zOP1R8Lf{v2PE3Xg!3zpEk3DZfH`WqvMe@SGTlgG$8Y|X$L50cYq7m}W)9wf)}^b}>2JCMA>_R}5eL2{mV zv+1en!0;g>UUbX31Z9E#Vbe~NJCGdRF{6huShTi_oz`*8)sBJZa)7~P)@B@d&QT6) z7#w&$Xd_P>h6_9&IoR=ggXcNK4U!-19CN-PdxxEOMj;w-f8wxnL3XFv(XcX8U;u*6 zIc|(4Um!Ur1!t3DtR4*?5|d_=g)N8AU6TkccP{xt$;AT_hH8o)C8yJFngs?qPFBhP z<(7jsV-A*$p7U{rZ!MU7q31*HosG2{dah%}&CE1X2Iti{K6$gxr5II-p&So-S9a9H$AOy_O4zDdQf%oeZm$@CU>YhzrQTBy+CcO-LT8nRU`9ge>>w&OC3-2oZp1|R%a0MrcOHX z_0W$yP3}f9s#68Z*+Noo^c)Z8x|&(&xxm}#fWTs}u^nhJ zdAtM4scj8&7I>f>!vlr6F8Kz^BhjF$P~{Dj-@O+4_OXp?6>rWx@Ru>H)tRj9X%+yUZT zA=Tj88xY5d6w(KPIO}GGCEv(5L#Ts!lD{IqYb1VDKH=1@`l&JspDa<-2L?x;3M~i_ ze-3$La7g!Zf&8#CT7}n1M>_Np0=B>pMjSqVXu7|T;T)vp= zz(BvJw6Xe+G~vXBJ;koP+QkWnQ?u8*7nq;KNOJ4N?lbGF;Rful>-&gg<6GwI8jUUe zEJ$kOH8<|3smEtQw|&b}twfSu?tV}Qe_?LV-D{c@rBI`AALxXBFfQ>PIoGZmW!H;q zPF$be47%@gv;EQ9SUDDd+DZXFk~I))lSP?n^3qm4Qp%tLX$T&T-N2FduZ1ST)p8f~kGb)8KXAk1 zreGGf?=>SqDvv1tfnMr|a?f2`#7UYkn!{0swUNuQd z{;>2d{$EU9PDnq6A$TV>v?_iUf5HEuPG*%tsUiCQSsVZ}1_Bn8LwTX?M1UR&vw&yN zwis^39d%1#JyL`G{OOY8u1JAM?bR~2Gou;$-4DqRU>3dhWB0ldbxm@VxG!vssm&_f z!?A>VQZj%&=IdEOD5!v5%@SD1v`+3$B+jIi$n#xJL>^#gm+SMq)DOlje_uoATRy|= z`{5_W_xens|G+cdA1;ows{ZiY=L>@rAQ`8vRbj|DbD6Av=z_VjL<1(b=WUERPu zV0(5gTyVOKPaQdOPW;!=f1^N0<~0;csAUUqbV}t9_bd+&zV>RB&QBZ}^Vx;Ql$}1u z-Cvp@m5br}woB?p!pYsZxt84+yMA0Vaysp1&wZYo=#Q4g&RzU^I_&ksDS$6j&BuD$ zaz8oDycYe|O(a7RaPYM8SWgT^VCv<(R!(9$t2v~Y_9ex*(Vcy-f2YLsdFB>e>nSnB z09~3_H|QB5mUJyA1`Z?wdH&myra<(2IUVV}p8V5$Zd|)(9Ak+FJGAv zGeb!AA%FK1Umwx`F@RO%w4mrQmkEdA6~i!js?2Z`IivOk{HA=e_H2^yEXJz~qFSUj zJ{=As$4{#y%2CenBVzZ=g0DkEV<;!kPK=&0@MHF-?=5eMp$=2Ruo!ROn#l>f683-` zR&py#PRbq2f2B_Jb54`b-enO$PISn=i@Xha5Fq694fl04IZLEX5f9dg?Sr-Bg^{Qtr&~04lVc{7@Pi=@!J46tN|;2#p;&~s>8N$acG26z zu&T*r6HQKoSi&ucwT)7y$xBZ$Y4}v8+Sp9a@aL-Ge^WCD^$t>wIq8HBy6D`~+GmaC z#Qmt@Q+3$vGS^@h3B$l_F>IiioKdcb6ih;tq->Unm;Uf33^5nY`NsB(;)kM~IEEQHsphGl%aFxd zSgc<9!-B^8yhMC3Tqp|c9TiWc9GbVxI_(4A2HoV?2ic;kp!`>pwgIb1?UVooY%cl4 zg4JyV(j{2h7=S{-Jgjcf`y?`TS02t zioB4`C=83*V#;(jIsF%`%fhjp9P60b>4ajd0{A#CVXQIe|lJQtd=fGfvf@{L8*{4R*}&u%oJA$TU zvxa~uV}-$rS^Vbmd`ucs_no?8D%|>%O#9ra!g_}IUU}b8_g=wdgd@#ksb=2 zl|=Ry{VHtRCaL<88;z0G+nK?Ne*`@&JdK{pYVwBhikg$17=66j84QKfI4;(8(BAL_ z!gA6uuF71WZgN&(nEKkr$i-puO08yUxEUjcrpcMxTQm#}^Qi$?luzuAuAe)xJ*vL9 zvp&p}KnU~<3WXkyERK-rXB77_bPx@UrJSiX8bbqpSMxe4W|ojexnT7UfBje<9W1J^ z;%1>cl-JgDDKKJSstdLdbj1Wwm74iY_=rr*qa|y?H&iFB6+n3pk&!Ay5OSR2N-T^> zmee`vtSctRWM<3We`H|b(C4uyCl4VF5rLA_lgV-MO$$v)5LXjpYmQV`I8wKd zrE1(-)tIu$pk@6kj2{mFu0datdst153CR|{z=}6@Jr6}US`Y_@R;yAJsT7qpu8Hr( zfM`3+0)QQ2!4&vJy_Q(gfYswIDhrEFrd*bvZBFY#=ix3&1~=G z7R#|)NC{N`Vd6$$e?prd!HUdKjx8`Oz=b-Xo>7H?XV68|#V7H9+B}2zgH(VI5`}Ck zhtzC^8l^t>${Yq0GlV{p+a_B6-}pF32^tl|h? z-RV^$bSZ|`-R5FHsH&$iXtGNTi0s2*#h56^Jm;0~IT#dea$b6U)o*amON1D#i1l2?sk(nh~LKE9H^Os`A0a~~%D?shX*`h_b6efB!lVe{8 z3%6xD$22i3iRhFt14j2I12XY>Q}-Za0QIC`Kqd~2f3%`uSrp0S0-@+8j2Hv^rY@02 z7KF=@QO%hFna~M?<=f#?ZM6=mn_$tluzjYjzaiJLr4eSSYe(^3d5W@A7m?imS)VUS1B;5N!}iOCGo1zxM2 zEeyzhS}?k$m}g!bT0I$uJ;u);m~AU zY1gdbXL1<(!75w8;?BAa8IdU*OfQS_HI#-Cnbb@4y$n%HOHMLYo8CEw#?@0!15%sk ziKf>r12Sg}VTF7b>m~Ll2HE0g8r17hOx{^yf5ZYKJCz%g&mQdz!T!xrJI~;PvF3c> z{86)?>3U2sTPt7B<0jU@T=6pjPehZCj{jW+X<&QnlvE}mNa@ zf1CUf54#1D>=MukaN^XVahwqy{4mpCxRLZonc$q zQv8`ft8>^Tg=g+^8p!|}241oPShD7m-N6jJBqNvJ{bDDs=O@7lD-}?p;av>Kgr1Qs z*s&T`LddO^V0A5+^`iF!BQ>c2>)@hKuGC5wxrBfxT(sT}paSsU!4%Laee`jBhkDy74 zRaQ$^#sBtAz~0eVVajpha7rq_s;%B6T#~HxCUWZiszZ2}1tYB`Yx~RXvWQY>XuC47 zEF@(|7kd_c7k&^a;-@%R&# z|A8m}D9xN;&i}wd%_G=U9sdJOx2g5PF8>3qBOJ{i_#b9oN80Uwf57grGvxL^bgq?T z^6h`1=A^;im-`{JJW4mf|8hSN*WYsIexQA%Xn=P=$f0d>OJC_A%BJpZCmTaXi-Le}LGIhDxYeYr{M|3YoPivdOncf$-lbdVqV6!eAix@{>4^0+X0W zGL-%DDB$$g-2A8SAh5`;$(`?jfLd9ozw;fiwyeSvx9^}2t#agC2P7b484iLFvFe_1Y8)4B_T_`1oLQ-Qt) zT$8p_oby0bzidlWjHKE$9lK1+*Y!>X<|=4j{FlqAFr++IHdKB)6^7=@)b;pwDnuTp zio-Hd9LqO}6xV~Hw$9T)Q?TJ&q3z&Xzhv1 zb3kWA(L}F@e}@3GVBSaPA;6*2*KWH#0kEjeg?xJg*6P0f09dyRi#~n;*+*l$8Jr)0 zi6nZ!dOrY;kkL#_yZrzxtfbb-jrR*4zkqrN_!lO`Hj_KRAGV2R4qkwNlW;B76g=SH z2+G9eH`sW8l{7Lv6FlDEB#qR(Bo6NvG$Z#VUwA*Uf6HX|h6nbUrG`06CeK{1mbb;y zFz*2WDQc^ENgUwMgL+ytfk*s@R`_2If{plh0YVmpzY%}-K4^7)$|a>ovut8n-V+b` z6Si-5(%q0hc9p5N9`e_sfn~`(<{#4nT&Z8Dx>Ejn$sf1^n6+tGy>kW7kWU7L-L8O8 z43zV(e*j4?oP9O3T>(t9$Z;X~+Y>;?M9Jj?PXI0lP1kvP0xG^$+qn>bK~SQ}9pc~V zp1mA)JV8}XNtr{;H@BV?4(+F} zx74s;{SEg|LVm^YbqDuzwH9+Xv*rP{1ogEHf9h^Jk`C5qnr01s$sMeZi#(Xm`FU`S1x<#$;A0O9lWzwAYhaKG!{p8ZK<{%j@y7poLa@NU$v6Jb*^OFJ;==z4 zJofkw{+Ajz7P)Z4{|q)$F!+Z5aX9pF{5cFHiWNMv_ewUD{(>9S&Bu*K^qU z1-Q>x1VhA^U!aiCIh))$24JsSWQ21Jhy%(3DZF2RW>nB>z}_(+^IFCN^tWR`kYlay z|8NX2PO`4&w_||YESC6m3=k%YRz}r32E@UTG~<@?Oeh0cPK*zvVTyKW~Yyhn*{B#(=%5KX6up9=23~>)<6Fa{^B$~1CeCHQn*>#;Q?-$@e z)$!ghz?^c+N51_6)MN5;06@=-e_aQeMy#3{&M_cZ`KV<%2Fw9Kh%pq4(Yzf1f@e&V zZwCO+)7UjUIRPGRc zp;j`dk4Nq`vtIk@ko(9abuEy8A@@;G=)MsBjoc>&+BS0}IOLv@(9uToM(!I`RPCli z?unhuFg}3jdHfs7eM9t9e{sNpII0Lk|DXsZt95iQCsG&~z-7t(Q_fXN38>iQvkqWjJ`UwQ zLxo%?K_@L(KeiXFaj{Yu9j&Z31QOJkNn0kWXVLodQC8b z@yXW*3iwFdi^4y@*|Kh!C$nOa2`LGp_Hu zrDrs}{EVG%`HZvgho2SS=`)A^1J8VaxH|qrve~a_s4C#e6qgYV9^BR8l8af>N!TOO zDJ3FbV!{9hal#}?UA*K(S)|U} z)iC1`NCtoPf&p1jWSBvD^?A}o!<>lg`@{;yrTMx>V@p4)5lMCB@ia>x;i7)D+y%WJ z&u*|Z3~P|MF4I)bAxbLx`8qH4V{yyZ$hmg+7`ra``SATcv**9aka|ok7T{*Six$C|HA4nm9dXia?`t@%p$y+>{QhYbdt#v*LN7 z%;(cmwTm3Au3<*^e7e(MoAuE$0rVJ6e?8qG&mpVKUtOkuT52g>Uw28} zP&m08H`nkPVb_UkI!>4Urp=#R zXNEsq58s@^emEjt>^3VR&;wOig3(P!X9UEg>dk7D_^T#AM%4V7$?SWr}pzbm@GvFLk}%=p^6tMkYJa z|AGN8^@BIY@OQw$z0=KOc{j7)^iHjOP6jCgeed6y z&pW@?1fw8Wm=y-iOs~dTapdJp{QYlUe}DF0|NKQy(rlgS0GPf%``6$7%m4GcfByXc z{rT%JlH@pIPZq|hzW(et|LuSL_Rn9xef|1x-B+@Zh^`ju*Tw1y_kwS-}_%Sule&k$rgccm<)u*y5_^>IUlE|e{|nx zD`Vh#dc@z^Q|5LyCd5vFYhkCHJOz*;h}G_$Sv;5Dt{i!mKnOgVb70nX!fU z?)*vvy!ZZs>U(1FL1sW0{n{_7Q5snDnI76-{~(tp8~n##t`=8Shx^ELU$;&?|GjpW zefo!wj;6g-li`LOsD=0r{n>Y#Xpq-qrrN~422+&=WSDhG#(89Ib}oH zKhtLLfjF}?+7cUqkCg3Tsd^I|e-EF=ikgeQ$7p)!0b*HM?f0&}$P(dQ<8B>CPvcRN ztz!M5=xL0cT{L-zqDN_&W_d95pb-WyY(gI+=xLgQX*0I?dDwJXf2FEV)wEVw&OB^! zA)no(7*~&x*B^2o#^H_`lp*wR^RRG8F?OiK%tL8xi<4#(^U$k=`)n?;@$xj;z8VJL zczJ{%(;T{vmWK>eq)>_`HdY>1Bq5Re5*sNGGfl!UcA3S=L$4HUwYH>2$s^9E613d! zVdP<1?q-1jj*tg0e@xuni+9KPc*qDvKFmeaV{|;if3?WGp~D4iC9Hn(R*}qmiH)R3 zSh40X+Lz7IL$?-ou^tSEqDPMAsZoHAp{HrWXEv$v^K|CzMB_d_^gK<9w;QY3X6NCj zV=4E=*<<89O}J%))+at(JDpaRYVKH!#ZmOLwe}PiA4Gz-|lUcYb7XgnE z?a&29B|Dkecy{=7>c+QnXm-f@#jrj7n9Z`oYdZuxbQs2lCeSVRgwFrPv18mjOiW>9 z)Oy;ecUVkih8Y)=jxq1B?cUamZQ>n{c3w;iG_1&`D^2702unIi3k@+7*;BKoHyKZt z)5U6ZA1Ih&FUxfRmZk~yanuqbiNu$;%Eu+fxk`{oW3MMvQ9;u^~tq}I1c!OMNApGA71hALQY*sSOX!PgjzXu-JKw%|~VoudF5((jF2@ zu`;%Pf8yrsVWn$^1QrwH3t1Ld+`XmhZFTC~d6k`o0YhWKJbAiQO|Eu4K8!|YU7~4X zqw!%(Gp-I}4`Oz=6!T>@7utB)#K!K!NNOA}u~!>7Bp>M)BB zb;LeWPQDut_u=)SxJFYp53LW)Bhr&9<;iC6f8kyVqn1VX17(R9op7fuMdEP!NdIO_ zs}ZC0k?8ek)Zd5Ehr#7>Se^C5=))LtEZeY%eg?nTN@K{Vci{OU@VmvkUOK{8#GNHVNoP2s-AJ(72-e6@m6C3HzBy5$} ze@Ji69~MK2mSgWIe^|%9_a!vGpFwTH3Q;(^AJ)!?LX;N^5E&Mtth~rTipBP0BtSx{ zHjMD(kN|Pzl|`a)IDq&xxA`!!Q2=ocL(C5bASpT&4{?tX0O^DqW%uFv>B3GKsMv-b_2k#H|V635SV6hsngo z1Eif+`|hv+4atOCgSSb5h9vOYBHI)|%%*9{)R4mfG+3S>n6vkiCkrTbyQSc z+rG5afPi#L=L|73ba%HPg3N%Fq;w82z|h?(4HD8NC{ofbr68Rmof1Eu@3+47J+JRs z>zuvLweGd=-2Ytb-nYC|-+eb1oRg}4s^osf`P7a*o8voH+z@^(AUSI_9^Iv`?x~D$ zwsazT2ydo_-L#Y+f|5F3;D&CKLF6hLaB+C|`0OzTorZmb)bUD(6NJ&?0``qF+tM`w zjjEN%qT>1`3s=@6s$3u{;#m&ndnEseEDIy;GZ0 z>62EvwmhZ?)e=0qraUI`Hfq1^?Ep(al;L+bdntLeQmg|JL|xt&-B9ox%Bh;6Nu$qj zAKdhk%6o^uYwWmJ@R18#R#naJMrs(|Rh4Q;qj=B!GdEac7=HiIjXYuFM07>;ox^kw z@AZU+K|S;6)URmqWNES_21v#nWzc$z-W~_7EMtIo1B-&{kywTaTIukpl#tB2pUQ;h zLN|$!CF4gv#2dSAy zMV3TDFlWOfAo=6(*P1R-(OTgk${ zMq*q=vEwH;R%Ij^M;@v0%T(voDg%j-#3wW8dNlIwk3+(=2W=bqRk4=e9Lwh3jr26u zXp{AJb8(B*F?qBp7Q)0{%K7DwS4DW~5E9zOOE{rFv!D(Z=Be)O7A#7xQ5dp3ZzPb6b!lO>AG$SPKw47!-^Nz|zl5zf_r-SUw_6|wBa6Y=d_W{GN^d_inIhbqd2 zQs!D_$ySu~K%kcCV0ci~!(Dozu}m)ICP1zxBVLJqDy_)N&B9U)BvI>@Dy7B-NzY|M zhzid6x5IokUxVllai?(TiLfO^(sSjK_Q-+hk}496-NO>`8xS( zv}dkjF&3wUV!gx8=gE;5Zai(UdD3I%Ot1Wli09$tt8cgclVuMaR*rBwoDu__R*_3< zZcZl=OL|v;iYc59_HO-I^b7!^Bdv~^;b>e5XH)b zFD;hGl<+CNE5(Z%caL#jtMiivXswPBZulY?|7Xy(f_b@JSLlwOh?ap9H)_Y4%>bNr zirlHE1-xm`@Cv$~cGaEm-<7_t6eo#vnd}6jKZy#fY;B_7+d4`l+;K%1Nf|7eFC=38 ziC}`R;Ypf*-oDHSQgn7*8&5FO@EL?a)k|^)Qp>3g{W}9<*e1ItgKmf09cu|mfi-#T zFCgOmLAM&pp@fx7=AZQy1-!ll-G^zXAvAl|j_RgFz{)(RJKoT%z%#Vh z;kxRd$fcf@?m-tLL# zyuBoXRWf~M>`N$e!WDe<7Vl6yd@)y2 z6Yaga4tKFn>7)9^M7EbArH`uox4*cw5>InV^C&APa=L4xiSjDG(sjjqj6Jo_OmscA zf?v%0gFdVn?JLHgF6YyjV|fo95pTXtu8 zjG-Rh}hWcY6Kus-(r-bhg{do3f*y%%()wm;iPr;iZUP9QumaI6sp_hf>1T0J;h zTb^n9se0xsHz7yCbixn+ek*h<5~E{}ap>ZaC;0ZZ7qQ4IphXB}<39H~cxvuI(h-sKJTS`1 zzT6ehJmPFa)<_c#VmfTtk&TXYM^RimX`%rjv$t=VcuE61hf4xYHG@?2xiZo(KeSP72?j~U|$xQG^Gb-v`RaT1sH+*O-596|B)niWE z6b_vntN9}U_4*9*pDFCRYU}*-a14DcOdWXC^IA+Gq{9~1;?ZY!c9jfzm@I6TVk71I z?*b5$YiA@&a<`zTa5K;?`q}(v$r$b29X8S^of?G~oA`TB?L1;yx0XxHvw7vmG<#At zk_PQIlgM2nHYW;Fq7p8Oj*q1YsZKR{svG)dL~*-<@RD@V1MOQL{m;+DmnGEB*m$Ge zkQWVzg|1Tyue{>$S4LWEi|dV)1@G8>SH2)li{D%Q#xHr5F6lG|_P}Z1BK)D=EVk@# z)p}yNe2I2Y{33)Y3ecW*@LR}S9J`HrK6Q97q&JjC&?k{UF_P|9?A8z5kD9$IO`Ip9+N;Y z2>CB|AB9*Uy>|e#z}EuNA9F`dcs}v`)D#t<*+)*+)4!b7^0e_wD4=-%L`29en3p$w zCjW@XT~zIb%}r%{w_D@v5g-Og;s;6d%0FT;Hd=^AuGtnx_d*a~j;6TcMT9>^61nmG z@MrqfgI+-%(Ej9nNx6w>HKryhn8z03FhT?S+C_x#dGT0mMpfWQcCMI5?%Pn7X(+DX z^k+G9#|T2$7JJh|WdgD$iL1;_qVaM{-A6M6Be-9!J(W$b`D?XUZoY1UT?O6&RCIe14Iu~8ww*wy|#_=xmB@YrqR`Xi>e|BLn$!n z*Awc{M+tYOC~af|pQaHuVe^l`xlPkdoC+@^)13@cL#m(@!^y4Exp59lK4z4!QQYJO_v>cUB<+_**bS0PqgW5!6Q4@k3Ev?cOZT#x4(|eh zOBRhLXN7IP&SRawEAy&ujRVYneh75pc(Z&cQ|>)5xgr&=A6%c4tW2k3S3XC)ShQyn z0>bn9j9X!7qjy1P!NA2_Q!?qo>hK7{d``d(P;qc$Ur6Y`-s7dKwBmpDNrCCwvaAu| z`)o{u^IO2ftP0xQ-1Q%bxx6BEj_Bi){Be^@W2PCN;`Ez(nPb)KZ|^N90HOia^wI={ zS<yVxN7k0iz7JKZ))n4=_}& z3NP5r8>!4!zv%_5mL}`#d^8&aRcKnL9+b~HwiMp5i)8${=5@z#&vcXXSqdP$-zm#_ zf1={3D`oWPr*~>GaDEXRi3Yx~bIBS&@8(2(<^U{a=-d_}|j})cPW*(5U%`sZxz27)N~LdT{e> z`WHk0Klpu238wEr^ZIn!_>QJ6r4GHRlyOz3!JAi3Zi!Al76o910LeNvKe<};=Rl8+ zNtI7KDL-(M+((T8>n4#;FKu}a*v_+-2D+j<>QlMBJMZG*SS(a962dj#8pKY(9UP0b zq0z{%7X23q{A;t(r`g-iSpsJ%ujz-+VfmlmIC0arDa=~+S4;NQ<&Cp`YwTKgNR6ZP zk%OV6JdHng({cwV>-P&=15)XccZcHMN}8A*q|HFIMh0;lohRAKN%gdz+WF7hD2Kv~ zN(pStqiJD~wWaxCne>%uL|uZi!|FTL$%v8!wq#8k0hxjoWmwS-0ScnLt46vHuC2wM zSZfOB2Uy2nuXL@R^4l?9w+>Ff6aPxkVQmVh^U-*~-T%IN(wwtaq=Jj_D?gN&Fp;t^#ju<>VYiWDwIc&yY}=wMC?G zh3kfe5hi)OrBUJfSmU_|mDP(>-(R27ZYePImbj52>4#L=?XBQGtDjV^aA5kro?Gw= zXWy3hrkkE+!pglW`-`-JCMD_601imQK3kT+H!5m%KiSM4E*p1aP6p3?OZ8@lN0v!( z3e<4%gKD;i6w-0B`#Mvku1Mr zj?f;+WyS~@3$P|B4U9NPanVK_U`qUb*^NVD%nBbqO_L)q(GwcVEU)~DmXsPv{27G3 zCM$(0J4{COZ8x9cj&U=k9x{uR4A*2dv^-!ld^m=Mr{54EE;wNL#uVGKGY&X4E3N~m z%H;Dln8ka|a3IUJKOLEg$ql9vXM(`0Jj2#&*Xsahn}qO|7fp|CB56cCSX*NF$!#NI zG|HmZTy(e`1ekmZw{K0eW%^jBL+cAE6ZF>J#**No44yW`o-t9vD~;7`D@$Ttg{v%0 z)zQ($5F>h-jIYiWn5TR3=l17E$mLz?#*LN8dz{-oalir3iKEEskIxFo1>V1_HVPDI zq%C2Llxr@SM{$HdF%u5&BcTY>)`+rRaG>-hUS$1Y7ZuXqk9=&wy?fg0ijzR`epErF zo+G%8mfYuaOlKToL0Z8JuVERl3;7m=r2NS6fjp4BJ`8#3W^u#(Q@LZg9KtOpZ|580 zRz&x03OO`R=SI82c-R@WjB>0_jG@^){6owhB4Ks={YEmAMRK9wf}fhf>a5AmvNfw0 zm2V|{)i@^v&XQSx5Hh?m)VqEl9NU}-`)C_D{nGWHE$mWgEem=b^tRggM}JaMjm25> zu~ugBc41q=B-fcZxtPwfxTG{Ro?N@GL0l{I`bDto^E2^jjCa;)x=1*@l-gw%@c1Hk z#_af7aXQ>f@OIo7+f>c?3t)nX?ykijavyheL(FSd3;>J}F<6BvvntQqORgqc2bV;| z45K9bjf5uq60Vf#t!7{F#11Lx^$Y1tyId(t2e(zyB@rPh246bZKVyhEUKYn} zFw!vMAmP9wdFrzotQz%LS;en~=OYWd-_c8nL=I>S=_B|}wXvm2;aFc!-!e0%)nfV% zKqs!KZu8scktyaCO#?LvkxY6R5Mkmcl2wFC0P5#uJ@EV}B0>DR^f_T47BUVR@XwO% zgwC`AnZarSe!Y{;>(r6&a$os=FV=c_MK>$G?aJh&KsBKNi!&ZHiVNNG#Ff(;JrXq(!L^osA((VIU zs5md@GR~;1dUahMWJb_x(3d+`CG<;HFsv#2%@LOTBqF6j#nC=1j#1cu2-PV|!kO~N+HxNjS}Te644=!z^(K>c zNyX{zGm~1}#E6(uCGETWthXT>ZSOPA$5L-C&?8IAdk011J z>8i1gY4}FzHfCV=@21lf@-3!C>-NXHaq0nS=z<0hfipsFkrx6P+C}SRLkL%*K77tI*Tq(IpT*VMQJgNOL@aRR9wY9PV zFY{@I>|436NluwI4@ZY{x|J|9zgd(Ys+QDcqM1IH4tFlbgP1Oo&X~oOKf@`EF2G{Y z%84$Ru6$M}WvZ6nJ|TZzbS{hK_0-WlCa-0LqN6YK4<#e2lQvg1l=Jhb_JB_TbA#if3~W=h+B5~;6r?4A*8 z4quy|2(7Z4%b%VY`r4vTR2++4x}{LoJ(tcSCt6{t>3<|H=GWq}rYkNaHt1+OI1ihr zq7FBAa1x+y`vV;0yooOP$Z}EkFB@{IpPhDne9LHL!nyfG2Vh}c#HJ9$5 zEW*}zCd6DL0{amJqt0s^uSm4c;kT>;tofB5i$<}L)sF{GPScUo4?LmgyDybNhlD=rH2KO!=r83V*H`b@%Jk-_;y|s3WhK=X-=pj8im^ z+fB@fXwG73^YVAf5$T?nptw_oev2M40=e#QRZ2fr785Y{K&Ppa4Pt{ScG0USGZaQ5 zp@Yux>t`?$Cx`~KRJK@=Hld1Zp-L+9Y2Z8~;%o#PnQj;&XIjdwy#XpeI<(Pk=`u%! zDt-Og;!cD=W?mxA;E=FIi5oW!wuuV2aXqY=5^kTKks=eX)1Cd|N$J(H@8z;Rm1Tpn zmYgb2ox61NGFuB9THAx;L;p}6Zg^%aL8pi1#c9*=z!nU(>EeQoc7oxFl`quoQyZKS zGbL(lC9|z61<;e2KV8KZ@WCg&m9GRSif@A}Y){BzUio!?D>@%xGjj%$N-O5!_d;a< zaQ#kc`m!JToHY00=84aB=yNsUHx>99W?tYg0`zBWl7V1o?vr+>SYIx^Z@-LC7En&x z?;369_UgE<*?C@#>(JD8=>wJDJ$&#!1duPlx4UA;Gk)=8qpJcZsc*8co$x~{S+pfV zKZRKg(HJa(!ffOw*Q&=3=Yv(7ZS`%$bJyLz2lB$%-Lf&zay6vzwVkRbO^LlJ2PDk4 zcDnwOPt$5($5DbRd9BkSvXNvY#sTMgiu{Ftjl;{QnTfw@9ft3s>pv(0nqu|_5YhNJ za8hXAn;MGnk%DrO&9)+J_zC?7cl#dSpg96|8bxovpgG~oL`Mu1oKz3M3hdA7AB@J9 z2C`jZI6fcneE#xh(D&i(y%&e4*v+=Iph-jIqwmANAZn3e~mYCi6^r!n4L|c z3<_qo$i6s)aA%#fS}nk5#190>kA1bBGQS18`L}92ky{F$3F2Nsql`dU$o?VOWUXt@CDfLue@XHt*GDq0{nB;L`YB^zOui{^7-gdBa$o zj#?TD`A^bLvL8l3-90bnG~&N;Fi%lLchoqeV>9nU+>M0U-O_6hYMjl4UI4q~R)nxJ zQ&ZdXop1KgM~OB^qT00fPO&>a>N94-U@O=@Q6568UF|AD^sI52J9J?;K3-V^=2sH! zF=tk(Bx{e-3ra*5gn@|7;AW!Tkw1E4w1(7q17{W_hG>4D?EuoZ=_Uv5xQ zp6vPify--Ja!t+MYFk`-{L=mS2D@jyi6ubXDdD?r#ZjpI=wy-)N`&X7%kMUu1AZXS zt$Zi}q`+=HBe&~aimDi41XeaMB)diM3G)-DQD06;(a7Cc01Kj&K6wOfRr+mt$xG-; z2gU$Zee3ceR-%Q&%J{q%Fge>p1@iiL=St#o?P*Q7ht^>>YXSwWBZo+S4|&>b=i17; zoj&;7c|mQ#IMyVZZl=*p$06+VD`&5oO3XNvS|c18Kzm4M+XJda9b@<OEG*C8?pdNC0Rj`~fRmPxoP8mXCFbU3OiUyS=%&{Cxe_kA4`f zWu*kHd~YZY+jJN8%Sxv`+CMw}(#)UF3jBQ)>ZAr{DE#5^2JuOjj1d72o{XC zMEtnfr!eMaxNImLeTm;J(VYf|d+Ckhr_Yf$Ap-52I%tlhA^HL%hzsc^FZHe95(>D+ z4mv$wPiAF&3YasZjKfjSs9M0=!q_OXVG7>zuNnPxht-Olp2ISz^^?jO% z{7PRJi+C@lPb$_g`r71x!ZHSI$-e!i~d8 zcyOArh9Y?Y9O5IW=SsAOXjmhH3r(xf7yqaLsQ~>$V&UaytY_r|i;K;;Zk}QL|)|~`5Q|yk!2#t(C97KCp z*DJbdu|d-|mtR;;_}fB|SpQ{*O|yVUPABW!L^ECrIA`mw+fB)mcxr8(4!J*I8l7m( zh>RDnqB$48%*k69S_Sz~n!E)x>#x|M&0Fh6?>l$lqxTAq_uy1JCM>a$0$W6DJDZLx z#hq=97+xRs&>VbRce2*o1-4T|t5R9pj_fAIX=K{IR@B}$fRevfZ1}ur-=zN#mMEF) z)hJroC;N^9H1c6!Q}S)6(8^@-=r}E}ECdzYAb5X4*Ppqee9VC#Vupx! zXDWfMPy*^U!Arqr?GHu3T)R+1-|rI-aBTFivAx!@fItWPN(K})x0MpUBQ-pYdp2aS zq%d_xw1_iP-^97*j{rhq1;^MQYoyPTCwITgaQGzs2qmiqYtLuw43*s)rNk!%lC{-@ zq?!t%sB%3eNX=8d@~4E3w^g3hB9uz;2Mj5e{Kjox%)7wmadJ`*AJ~t?&uU>VxMAfC z&Y7<9F?zS?r8|QG!BYzF5sr6WYsj>PYW)`ZcF>zDS!%*SsdVrFx3d)Q@gQ`)%heB?Xuvf0jw$D;~SNQ-^h#O+{Do3A7 zq+I=^sTYx!(@|b9=2CBMZ4shySZrUby&oB3jR=}IK1JzT55ev=qO&w z&X`^ESeG`@aj8Zb-?0rVnPE{GTDR$UIbjtR{rwBIAI0|1_7d9^AE>y=>I*gA3YJ#Z z>I!7Jdd<^UcUj>YDph>{_T$*Z)>9;zGx?WG*+~d@?AH*|b$9zSoxPrST=9|ftfWRwJC8%)N(*8RRNZpy{52a%PB&0NnQYgCx& z^}9Ekf1hXbw~2;sR&{g*j!={Wa*DokX|jBXIvRLkN({kfElMN$xX?9bw~#u)Fm|lY z>JiY$bw{=|>`7F|lQ~|!bbOjv2Dvv?@`?$fI5-lmYD^XZ`c&S`*_96drY%%Nzm76 zkZ!$t1fQ{Z7Inp1)GE*huIbTLxlxQ!zY{AaNu2a3g(&?(CUfHt9i0LiQ|Cn7Iv^~}J;gD^1ZckP14 zwa(&)&(`Li-^OOH*`%r%c!KgmI~U70WecTa`xu0Xi)KhV(q-SHVBZu9=}n8)pGNH& zhmfgEOGPCPuMVjeK8iAV{5xX=ZAf_tUt^W$isWrGj-#Q&a=1%HO?;M9%KPSZ=WF`M zLln##Qi$V@DqTct!)V$)ALp;-t8K|Fx5U8buqwFdS-DY5=V(#hK;FUTUKf`?;}sj&wit583DEcZC7U-F`X^QBjV`25=+*0m$Gi-cUsQ20-Xq8J?ojz&Ol|X|kN&4K zS(@!<(I>V~YAcs|=jmC`17%RT^*OqA*0INJjHkRE!rvcAZ`?MgOHj$r4&j^K(yyrS zMZ|Yr3Fegz>mLsj62eW{tFMG5Jd(B~M@LRj6Ne1K>?Bfl>29I2>2f?t(_5*5fdyV#Gkgzh{=Kg3-SLY z`(TeC^ly70A|U8rav%Y)5QP6PJwn1D!N2TzkP{L4%V&@PL_pvldl3I}8YCpl4+Z_T zAyimU}Wk!xcmCWU;gI*{Xf6{=I{Rcmv4XZ_g{W- z|MvCYzT~&xeo1e?eM@iu{+8eV`tA9z3HWzk+M4gZ@6GfBof4t?#9zz1v%z?+fLXJm2%W$hU};lHcb5c7OH>WWSSl{mYx< zb&R(4-}0*Od-#9;^*3+->)SuS{o>CHYlUyd|NdQS_b+&)z1Q-_%G)15DemicZ@>6o z^4nkj-*12We)UD__+_{ER#!uRe_7=M@Kb7s7AmEu+It`M z_pM!itz9nt9aWYd*6Ky%WBoq91yy<9eU6?Ruum%Z`o`W-ui9%rZw$y!dSlz~rPiL- zTx)&u#*bgQ#``QI&3WD#ke~Fz_nLg_jP+i6+j(#%^0hDMkBxxx zK6lF__5!3)uqS&`cj zW?^age;FoUs&9XWA=27zvtV-RIQ!mNh;^kF{AJ<)-#+~A*Z=K*{DY5?@2NC)bG(q7 zdVR+-kT)te&`2r2kCbxV+#ue5^_^HPKl~B5sP^75LGkIGD}Ac(t5rPMQE(`!eD7C% zZ>9Eoo+INi%!gQ1{80SveYDAUmr7IYU)ggnf7svs^+)Vq|LJ?zpxWOt4>Jvpyx+L> z;a^x4*dJ1XUxv4zSDyF$NZ}2u19paldTcJ;dyQ1kSlH%X%l^ZE8Adr?oF@K<<#l0g zTy?FEty9O)U;KVw);S2UN0hhO-goOY4|d>(qW|C@Hp=_jtzo5I=&{Nbcr6&?_*!}S ze_#F{--WXK8gtcz8UB{BMXi+T9wpnJ^p}6kFM}a9F)M&1cG<)lFw5i)ZteP4zsIMc zaI_hksbHJIXvi3c8#8df6wUR}4S6doW68i|Z)=S*)kNsCqF}5?G`?V5%+0~ImJY2y z_rlIGNDo-)!vPEH@_xrB&$Z9#sK6rMe^{LMIUQZ(tnuM@pJDNQ#GT~))7lj?TmBgn zWXd@XANJhM{8yizm#|{Acv_O8Sz>E?Oe%Uz3P&VA9{m$OWj;QRpNkhiKKF;Mhk^}H zpGURjh_%ZPnN8u1109v7S-$tR7yg;=#<>oLzZkzGY-gXnny0Lc1*kq>ec_+^e`;7; z*eg<98$S=5&S&rDkqGC43H#;o@__%$w`1D8bNT5VCZB$K{r=bYj6DXc<4dy#{O7*D zEl%Hs#bS4zaqatGkE2qp8}t6<^%wuSub%@a^oF^F^+d;9|K8VQ1LFGE_RClg_|JTO zN!V7gE~eIOSsV{j=;Kl7SRFXHf5X)5tBoXQj*R8LtlDw#RPW$mi9~Kt&@ZK$^#PaW?pX|0p+r`~Ra?(u_LjY2>LiT0pHC08*TZtC?Akaw6)scH%3l2ER(6rg zSh?>#cKKZ>m+!4CcXBMMe?Xp<4fv0(RTXCb#?_y44=dqppm)G>fj7oW9ImWqC4crgd-~IJR3})|_e}et2pqDZ%GE@`q zUqd$hNkXHR@xXDCe=i-z!vtF2F}C*3f{uilumr`~c9rt>^$T&lZQ+PWe6BvQ@!?pH zZI~<@(1r<~c~*ZxXP;Q^rHp`MzV*7=e0!AKVUFzC%CjB=wM($KtsbY=*o4uG7#WW* z9BAD|d+gC32Sphee@AntP4+rIz20la;r^xMiP|;TNH6sO4ui9$@u)S>hMvVCAbw#) zV|zStZ@pOZZRuE#*e|ixw-@6V-X8ySe0=3T@ym~gp8ST|4H7&Rb5NF{h1D2yZ{QeL znv{Wq#;k3ZJsuls|M*+Xam-@U^3x{m|6~u>Hy&UyV{i;C&%Gvn? z|ICRR{{A|Qn|Sq6bcNw%lT0SDIb0Z)$#y4C1ZTBy-}7Bch|F z{2X&H=oz68fBcN+N6BaGKk5-?N~rUhRgbSbGw{(CVZIVCvlf{V;c4>w(_VK9E~wgQ z>$MITf1B&?v&I-Z7!YMXzb>GUws_#D}&`Y(K4=;_bK#ZRBB_%CdP zG%%-Nz5XYetL68dt3N;r!C2+}7L2C#GSh(1HfMX_f0b8yczLUjz5?1OPF?M#V}b9& zqwr%C;Ou}Cc6%8yPxOwjh~Lo|OTHsF)+3JdSnE3};uqfe6Wb%swl+~IbnkvXJuymQ zgg)mn{)RO=2(eGr-UkpE)1U+y8E-i z@iV3ye=pDqcJhVONBMNRegAxG@4522y3@S=Pp~)s`i(R0&v-aX8x5IP#@T{rkPjH1 z`-it*{mrMJ&6P9vQ96yk{~YT~Txfo0?)yv2>Afw|hR>CVIA1+-#2fY<9C_-~PA7_k zcYd|g8(I~{F}CbDMTpjpP55D_zo2LCJ@7M8e;*~EIsK?dSOFfjKI_Hf3(wl}XpgWG zq2=E7 zf5v|@Ge)jf?C6g_=y_idEFQd|8zq$&Go}rUs=o4k;6yv#F&}0O(Y4x6FJ=rZ)^{vZ zvp&oiqP)THkG=Q>y^|&4cWT9w?=+0{h?6$f`cCoqg?BlK?GYD`!J*vO+|Q@;JG*$b zXUQ40+t`3#>yc2uRZDr+8f*H_<@tRFe|xk=Ku^1cGD91O8a#R=)Hn2B__om3pN)*~ zTQ)vUfj0B|rppDanE9y<^Zhe|=Nh~^4cEWfUhP7RJ?p4+nHntYj46ns_f&KcP z_Zl-Bjsd&qJk#kj*P7J^t%qSgf9WoUpJJy@#c+Ih=Pk#e0~5xumQy%S*wOficTB(8 zis2Z8m@9&(nM$sPFO2G%1t$sGv>dK}m^1h;tOf~nZo}}3p~Maxo#QH5eXt0Ok+Z<` z+OuZBq``>9f|>_r@Wk0!UyCwPOKhdYJFB@3mvwz9o^hlehI0U9HUo|Xe>j!zZumqE z6vh#R_PoNo%ru7JBkw)5^*g#~VkiU0wu`Y~UyF`L#$Lvnj5W#`dj(fLrTHYA@nx~uMWY&%R}C0VWhlslUWY0*%#gFW_jNj6;oRW zvm}*K(K-_wYO38jalMn*f9_Te#&?>Aqn$`ivYeS#+5S9cC_81+e<`WAABuRV{~E_E zipPRtPtz}U%uyZAV;PRsywiJSV$NbYw8q(3FEJv!JYUf?)`u{YP>+@cpXg?JYsI{3 zjoha^jzO2Q;T2WhRg1^EogMFcPP#eq`MioY0&x7Xw;|6rE**+ve+JE!rIAT))sd>< z9mAAk1bv7l(?%?C45QWcwV3Aq^rlPAFu1HDOW0JCmXUzcF&gIf z=Lfe49XUA14jm!D*TBg&>j**J^|h6B9FEZPcqf;h4aZO)m5qY|j@DTIJX(NjgcP56 zHZvTPqbz+b?EhS(e-A6?F#eXM92hiJM^2@C_|T!Su%Qfc#frC%b8ge;*M@m&_~0zl z>mBSH!{MDcQ?wY2`lQ9632S4=5yzZs)>l{~MuzdKt$;PKbi+ISVRZhCHEDw=zA{!1 zfw6=v>q1OhUPV>$7em>7Zo{!vVb%*iOD=O7&RvvO))Rvoe@mK9=Zfi!C3H@~GcVwV z;Vsd7UdqAhhON+W7^)cPz%v)+p$$G=o52dAQ~y&Eo$*$lgPtv~RU=vJcTV3rLToJ|}>(Ck_t zw25xqzzhD?rr>4y+ydTUc1t&Aeqz%dy3xczL#XdA+O_H&;;kO`2wLTf%N?(n49C_q zrnXN%#T4!UJ!$ihPSQ(wnq5t%7GZp?eHFLIn!CPrf6xPK8D2K#Wf2@p!LqNyUoZ-j z`ed>Ht$0jD)R<h zfVn}dhQ;7JwDz>Bn)Z38+gq~QV8)Ron1XNK?K|EF*3%|9y{dF!paqRKn=%A3Rg)Dz z=vlnPf1Mh|V7=?v@IphBje~UxM%nV`8(Vw{z6!>SKfk$2@m?J}(W?5X+{3YW;Zrx> zx?sy?!||y@+;GUhr|j^;ibb|H4=yNO=L4jm&UXA`N!#(!Gv^YgGuCyQe=FaV5>8y`b(+J971svh*jD=1w!!oo zy88pioNQU!rtikBRB;?Mv`cu04r5($Oe>h`Gt2XvzTwsIQWw1sh-3d=I_`nPlp@p0 za;ms=D~t`FOBNr8MPnuzaP}u{9pITBs#*+Y8yY9+YcZsgwkB{)S=cy&!&ots40zvL ze>xFEr_trs%UNY`$p_5ENOzjfb0WTZ=*Ypwp^G4J?1z;$5L=l$nd4{@deKP`&c?qK z!|@?Jd-3@KJFZQ`<-%HEoc@wdpANffp{=yDdgVH7IBff}v~?=)PK)62S3f27y=*00)_1o~B&S1|utw=p&BFan$&6f|l! zc|I6Q24m7r+B9@YoIW%=I*z--5U(@L-(f>!a)l_gUhq8n#@4RRahHH1H3fmAZj46Ye?+8V zVtr1Z-^N{yUrXnn-VDbHr44P(-Y_MW;foH!TX2bYkZR&YRO#EMGa8*Z9LAT=v23EF z3dL~jnRG;wVp=nEVbG6Hj2F%yY!BkyRC)SG{o3&gs%r8L?p;_u?RhK!dc*<_J45OpWQ!V+q$ae@l^SO##OM&pJ~EmSnZna4f{hXwypBJnZ6#oei%W zlUN%cT>dv7^kPf1kUW+zt`P zQjfVvs{?1Unz~s{z%iZDtoFG$zxeZ5G}qR0&Nc2x!{E+Y2MmzzCTP5iOn|sn~qV1I`<0w&bQV z42J35#_$>1lm!9KezIoR^tved_S~-3g`DdGwmT(RfWvZ2G%)C~i7BsUlN>x7o&BUfe`KJt;dCWl`n4N|gw8CU z7&+YK@lv{3zI0a~v@5juvBZbIt#{5dVGJs5{@IqE`Q#kTMO42E_N!V9^yll`ECt8myV>)ZK{C7P zP#fGNis9sNf2Mve4ivV&0Ot{8>k#D1q08Fh^-Gf8mDRcqU6px&U9!kZ@vSw7$?{Jf znuA`x24mOXNgIu8@+l!PeKYf6Qy+V1x{{+~qun6-)aB$Btsv zWXu~zo)7)zeak4nzL&n2KD4?_0Y@LMtnVdHJ{_Mp;|#SD3**}658xgu`-kE6qm$eV z#oy)j1WT6D{5U$%(LiuCh#HTP&M#t)>ai)(`8;}KpLZhu@bfflX`C8<2zqY4DeI0E z7csDVGK} zQ@mg>@R(2p&ou6xL8qq=*9#Og3f9AI_~KGee}a)F0uCzrRalA{I6I{rm}grL-H~zU zGB~%c3x%wIW1(c?8F5=2KTJHQ5ay;c^cxdV;t#pCo0$>_dRh)Z?Ybip@ifD9%Zh|M zZ^tBx`P-9NM<2o{OHR-SVttYN5I7L)1zSY+#{;uY>Iim6HvBN_jT~OD&bS9=9n->y zf3hdD4s)R!iyNr*uFGRL{7~zSR_dCcQ0Q2`>Y>DWZ6-V`%E5}O6u^^OZ$5ScgN{9! zxf`0nheD_Ajj^coXrR#9J@_c6&;wOINO4rd164j50Is!paOLRCsuLx0-;0KX4cW+j z6X#W{;fMR49J2jz;)na@reY)TL47w&f0OO-p}x_bjB(Ri28tW4gAL9Z@xyS__TrQm zeDHElyM%6mbLf&EC~jI5dm5J&DDJssw@9B*+-Uhd3}!SMtrUJ!r>PJ08fiMi%}XcYM#)Cam7K{w%r^rqR#O_}GT*|K%erScK3q=Lf4tP( zcIjBRQY%_~#;GGaeh6`1b|yP?L_(Znl*ua$4}>@wz2(IBA;d8lX-rMS10jwLvVJ7K z#q&4o9dRVSdC0~pH4HxlH|8NP?H!-F=sktFOZ4T0_#wVA7_qZfXR;6R&B0E(*ziMq zH;r{Jh9BmeyjJo0U*?<3e96fff9}J4H(oCmwZtRy-T1lv;D`C_Kby1DCzI;p#pdsch7k*9F_~KxN}+cJuY*vhmtE z{rSUX)4a4B?kAOv`7n}Gf76GuUcAT}*(q+PtQWlnx{|36W}19>`M^ML*}9DLy$m#a z0iEfLGDf9$39N8CL&rs3|Pu_J_%^K z%y*`_l&u}brePKo!=0Fx$yUrd56mbo0YT2bAa7{k`tZKtK#v9sD_3$(QG#&znY{gal) z9<$R2T6)pFCK>Z5gxO(I*S5QbNLSS2<-2or_W7{bTMWP9Y;|C<^Rg?4NMiGrzS@zUM>C#ib#hh)+ITUS+%D`Y zZTd@Y%LtL+hHtOrHp-iNPlFR1wiTn`#T%m<9JuXF8^&O#f3{_kFtFHgtb~2*x6*&W zYICV=){^lP)$qn$(HQ%TbNetnP~KQAz3l_#EdytDr?w;IO^a$9?$ov{U1+pwxKrE2 z$rI`GhuY?A+cu{9GTTC^XQOC=&(wOd9g>W$Kfi0%uEI3oV((P7Gu}x^Vr$F?<2_gb zC&&Y{jgzLie=|LqZM0<40rh0IOBdef+k@E_(ly(qRA9D+OJ4aCVjE}ZY8>q)u*m^o ztdtJ~w)&TsPY7&o>PAln0vpR;Sx!hhfi0|JG0Qm+*fML!(jJUNHEPqGt_wsp={+t4 z%W!9^N$QSf?o3an8jaAFGtz;nE|>zk=+KjEaubpAf3a^BU8sR$!CYod z`9oFnFk8*J5(mGuC+$tWlG=9lbg62lH``SrFTM%I%!O5YWnhuQ5E+Ias+u0G)$TbX zC5>k2e>gU3`48?H^LiYr`f$$~YizOmgUA{4a7{mW!vkrIUmExHB#mLPrRikvA&nVR zxy{Qgki$ZquU1~kVOqn9IX?uyh10@t zRtI$QskaGd4D`#youvHSi;8TzR8rfA<7W(pV<6786M0Li$73_Q)AaexLo3 zGbV?PZDAW8IAi8#(R0X?Ge)1x?p*TWjLF?IWu4)4>udU6Y_`+kLmAVE&Ss;05B^Mi z_*p`Pn~dtm5+i)E(9h^x%LWv??mOk~c5ZV1Rv)e5zHrpoA zf09MatHiAZ@S5U&A+#N~CAoU7>66NPP~PHQ~v4w~RXL3N%o1Qic%jiALYnN!t-mlLaKVgxO3z~$6V2T?RH znT{8f26wMk#KFl*4~7cFg2~0eZko?p@E(D_m@FdJx$y{wP4MmUeU0|ScB}Blb$H zb?C$Jg`vMc8yz+xt6OpTQ$mM!=b#sLr=X$wTgnBolk6 z2mRRTbRcDfTB8TB_cF0<{6d&YotI&6kM=mIE}EMsnf-j4NBk_;ezD}mgbMUic-8}} z$D|cm9WLp+a>d@oEULYvpK8o&hO216Xi7GrOSVp8!?M-**A z1L=X%;O#Kotc-w@GqI;n$>L-TW9Ctf$<8Cu_4P}vGEb$8-dc~3TXfb~?7i3$2i24W zSUzj-&!@S9FbeNcumme{f7^Pg1qK-m&%K^i=DDNuhV1cq4>UAJz(do8l0!3xdOSMi zRyuTG_`ciKpA3r=miYaVu8-fz#Z#X7bvzwB^iClu#zR=_$-?Dg>H-?Oz35eO?~kFw zn|w^e58xaxV~JPzyuEr1r6R$9Fy9{)VE88bc07ktslxfB@fb>lf3aA+|9Z6EMV~_{ zK7D#67JR$#S&Q4@vC4OFgs%(WX@&>p?G?%lO&elaYM_t*J~7amMq@V6?1Fj)KX`8#7Hd!|4n_o_=KJoYbJw#$02Ymm3myjON~R8{T<5?dml(xB3wh zr|ymg9;ebYe>8yjo>NB55T;`bCD$_jV=BXSiH}m1+9S>+XVbHD1+&;`d;PM znl(aaDg-*Vt1}hmuD%uPX(l5lcv8zwajfnt`oa6l?zI^_Uo0(V<2b;@F7 zeo}lg)+38^2iEY_rtopZ$6lLvFm$rw)#osrF_9A2e@Obdr>334e9UDe4R?mMGTJnp zS&Mm%K|+~se2z0Rq-(4!19a`ZyujH!LCJBRj$V!zoFZ#8PoG`qGMpjZYCAKSbykCC zW2dN>752MQ42MO`@Hnx3QGcxtJpk#_p$8)}b8=Tn93L`vwyXg+`ijFnO8CU4EOdPGC@m<6jxtZ(; z*2dP(VuU_}O!VYalLi0~#;Krm_12=_lTIx9e;2|6wUZK#<{IumoGVtc$m5o+Evc|j z8LrDNAG|(phEFCP*TCbe@!qzyxS1rRj zf1J0pMSa=WW>)>la{|Ysc!^cXW1={g4yJ$;oN-yj@G-Y~S&f8&>uCmo&*f8TBPp?*pWy`}35{S`?R-lVJFSHj|gY=t%{P)ZKDHWLl*H&m2~}A z9Dgod%UOHsJqB^i)g*2!5|6C$oy2EOTN{nRSGq1Rq<^$N)LCoLJ=$?tMpc|18lwc& zZ;JTf+Smo>aZ}#zrfKqkG2Iuu>S9b^G(3hRT-b zH*a%>;rOMF1y-S{XE1S?6gn81gtJzyC%1ii0qNMprly+{Js~SE8|BU%b3(5mi4VD@ z)M#1k*~8$)CEvnXF#xQZ^nNG7saO8A=a%7pjBOq*gRXd54lVCAl3vM)f3NH`oXbd6 zM=|l-?uO4q&qz7=GjE2p8a*4i)36v+cvVLSX1U!Jah24*Y7!Gw?1~-c$kZXi!2QJ( zMd}!XiY0@D&g?F3FaO{11 zd-(O#Ksn-=f~%OP5;{jMf5gU1xY{giat(aKFubbcwuWWGy!7Wg%-yoh{|`0JhQp-a z7B*SMsui3BwUF5oH;t)kjmczLbC|~o3~%iaIyDaFGEhe=ETz*Q-&u9Jq%GUf4TvWS zreZXo5E%CMSv1zOaHMpqm{ol_h<9vi7EamrR3|3lXq#$Z)htlCe>E|hQ>)fjST`J} z&&#mq=>d@q-!wYRy{2KW*6L7}Y*cYM5kHOXFoX3y$P%q-Lxhpxdv7hCF(2J=Y)s4J ze>%2&t^JbBEjD&A&^1zqIEQi4>$QY)gl_+h-Zn`RwI}A@VQvb;3vDC4oL%%YDTZT( zXyyhqF{ah@V&YcVe>y{*KRTW9KdCOnS1%`(U z9aD6J#Pja6N*`A(hL>!pf@M6AuolB%R#)>y;R>P0mrkb!f9g2ulvuV;(mk?I^U_J` z{M))DK3f(_iTP_^B;>Vrtfper;mvUVN-0D0m2>lx&Xh2clFmZx8n!FZnQ}7`i@tVK za8B?>7xCR~_Sd92N{D01&801dBa(U;5~pWT()VT#x@9Q*hzf3Bp(UxJ1m(%eGf zQJR8jwld-bf2l@0-7lx){F0VL1)Af`SJ%wt$qoc`zN87R{j|B#Rl^IDvuHnIm&mgj z-j>7}sn_nXa~hp2tA>Gc2^pGBs?p1G@FcB+vi50hh7Ag1%e-nEy!+FIi5)2l`_6R$iec;I-k6`Xtm&I;WYg$ygae;JM%XrGHm{?JKF|4UNS<|qb( z$|lv86M3~}$QIH>oXc>x$m=o3cAVO77O7cl){to?^UIt97Lo}YK>Rp~<8WDqx!3Rj zKMe0Q1oU_d3#sbKi~g6pd*^OuDNsCfb?_qNx^ONo!1PYWhu%md#lvOO_nNu_~oxa~}oj z9A_3sCK@tlQwr)1!zXrwuJ#W`rUrHq$DY%)?K5vi^;;WVTdIXvVnenjZloWi3nhSa zt7Q(D(E;X&yPI5jD2`LiT;>Rg5v8le=nO@Ze>8LAO1NyINY1&^s{4(b;w{3=b~bZ&J;QkXrIVQcUyoZ9nN*N7?MWDdzS08)svQ@3~P-_$JF-ccS+{!hgDDg z$r}UprWnjJK*Uek@By=Xd+cvg-&*1{{E7`Utj0DlnK`VPa~)&Y^jJhkOIUK#yr3B! ze?AT8C1X;Tp~Acr-SElWm1+hSR&5y$gKTtl8KSdD)9P>%_nOo{g&gB@d^EXfqT5=g zwwu8)8Kraa_<#Ya1D_6z{!AdU)RPbe|6jq9+cEFudk7{2kB;6=~QD5jFR?O@jB!@ zE0vkFOAYUOX3)F~jL0y2VG?e0kZeqiD>$Yn*+g@Y#5S#m!^!+cxmq|m`c`9eh%-Nc zIX_KNp)(w3x|G^+>8oBoH#@5i4s2G<#+;C&n`a>YPnRh2P!tyQkvY`^kodq9e@@mq z*q)m|zi4OZYq1RI21IwwrGAsNV3wx-A82lxN~+yLe@${g#|2o!Lx>MGm^J7f7H~0qC#)mpq6d>hHx4pXY(OiY zt5G$aAyuWAXJFR6RKz>Zn{Ao9e<>E_@dLOy%l(1w5?!z?7rKin%MiI-s| zO*mxLCTYaPqMc4*}7^S%MbP5$I!ff5$mFBbfw0EQjrWA=N6a8Ur`*VvDE zQNWQo%!mMvwMfEeh~p$~9CNrpsAA*k0XfTE2w6tr*w@W?i5Cl%Cxun~e_|nAAoShl z{BVIVcw|#w`Iid>G(s= z#@xKMDm$EM4QU$N2lLOse~oVIOJM%99D1iS{DJu|bTa#d`R5LoKQR9!9IG*KyrwQ` znfp^>Mk~k9jR2#trV$ zxDXH_=ha?_e{yl9R)0DW|70q4OL-Fi1&3-i!(WJh7$6^rf8K<3e-_fof;^nPLP* zpxPN4J*)L@Wu;eAm@({njlxU;Fbjr8WbaKcNx>h?JQGxZWab%^-0b=$GV|Q&yRnvm zfv=K6z-s?s;2A-utr?m2BAFbFv42qNjCIs85-D}A|4V)JeT`P zq}0iIJ3}5A^*!hzX82*$`FPXZT#-@d`DX6y#q2t@Yn!-i+dV;G)ESdJ&6Gb=d)91Ue`_7ObK7DKaM>G5|>KME9WB!CuXWoMDWcZ=ZdBZ<6YgDAp zC6c=u+56zmXBY17+ehxaYVMo|cg~%?P5sM z9e;cBf9Y(SwA5UE^6A*WZEcCfIpbBgnH+&Q7oxpS>#9JUOZw=hXL-6Q79K3@tJ1K6 zb_sO}$gVzMF}!BW9tm%r7?)<3r*R6%jyugx90qN5`km=zEaNpVO{c-mvJ0PC%wHt3 z?5G!xcwe>O?OsB0WyuI88D8f1ts)l&F*TNrN-n$#6W-;FY=F9S<}( zHGOqmKRh@od?vOPF!12w?`<9UB*Ghc$7#3|;afL;Lk|tk`-+|poqlL=I^0doBKM%X zX_4%mKhWBniW1A{_rPnjCA(GfL2Kh#JDfyXn|LqTpAWP)cS3~&uT2)*FOi267W zFxM%Em$sIv6Ptw2{J?@suy?k-A+X&-e4SaFR$(IhWe~`Zv(vB=I?a9d{a>pM?+6H6-xhtPf7^#h+j5mn4}wtIgRb z4L-4nIy>$xxI}N$-O6xh!G%ICd3rII1!p4|kNjZ4(Ue98JZW%lCi9gCSJr%0f2Ylz z2d6nh@e1d)od;(mY^4C4;fL>LWA@sdH+Q~UZbXuKB|Q1=#XL_22f7<3Yjc>#FUv=P zr7r7fJ+R&EQ`Qgtls>ZHLI&6^V`RZ4xuvFljchl4Rfmy#f$e7TSV*n% z54IZzy40N=KWsPMZu)K?*ls$yfAne9IjlBgzc6{Sv)8a-M%Gl?WNo_RhrJd#fp(rh z>~&$%9?f+h*=uSrOdGCG_Ii_gce4Ai*F4WPRV8?#uQ6b@liizVX>Bwn>8&4V?JU7^ z2M@G1%A0zWiL^FN(`}6wY3(A%U?mlw^mU~NiPPGVzGe*l2l|@9=522;f5+)-QvSO; z8y(r}-ZkUwgS}>#?^c|pKI}C=oF=2;hrOmRW--6mKwoq9*|yI->1)`R<~sePuZ7xo zTKmw{G|;Cv8pB*A&7>w{NW-ES0)3s;iZ}SM)QQQYG;z=iOFjCstMtfHb5ZE33G-m1 znXIMw{mJL}A0eds;V3Lbe*--szf)R}`6*_^$<8GTJXvOq!fX(N!kn5mCNd-6MYgTT z0ywhWfD}296J_>+s(p!>#>Aa3C@4)yUg4-8Y7{aEs3Hek;`aO=T1f3u}Yu0pLd&l${S z_V~gldfVfPd+WuLZ%fB|#D0mjzP%W~@b>tpt?NM>Mkz>IWXQ0VM<$_|jrmM;W6%{cOi&{p& zc@KdV@=-F=b;0D{f6sa>VcA5MFjiW~25$3@52Vpf!N<2pa~$-mo)@uHygn@vko5F? zluV~)@-)v{Nc1JcnJ>pyP<_VcJ zGfi~*el}hXe?a%7c%&6^a53nZ$na4xbA_^&k6JK8YRwYjAFE7ezflI{@pU{oW$v#p zM4F*rm3X5O^YXO#%u!)96if;bCVAQ~wWy-7$HeZT%2k}dc&ib=E~2q5p18MMEci5X ztVQgRQ03cwfAQ;XPk%Zt{t&U;(KAaN$Y`@kTp=U6Qtl6&2SZrp$plQepHul9;PAejSZaCbI2s4==TYewE`q<>+6Zf2Q+x7G^7!Ob$HD?dh`~7*Z(} zGrdl&hd9;Pf#!T+Kp)NVz)#!Vl4CR9dOSPj);jdy{lM7SpAC(l5YYV?t@OcaQ^P~! zJU@6`66AMlj>q~z>!^$5%X1YgBJ$B?6}}6T+|O>(@(+|h>vVj1E@QyCMPl1~lw3r1 zv8-o3e<-xK=pOc1#~5kdc9AXVejm_+Bb6*dw>2*MT*uei?G~;I)|kx zeBTqjtr79tcCqB!&aob`SHenmyD+X~x3@nX8$UU``w9=dvEIxQz|f6WC7_(hw}@os zA&?wb_jjLR_RAj+?2f-_`u>2RS>L|tlKk%;e?H9XgUvMDbf9MRuy}hEa=N0-Q+Q?G zaD-xJPzDdDqYlnDM+rrLV?ch=8+kVg6Q!ZSTmR&ZIlaE~X2|q+2IMEbvk!4;?i5)a zpS-htcxULx_MvBg(mQ$Co^f~{!XbWqXT9(Iz9Q(606UN*&89+qX9h0joyYHg{q`sS zfAvq_c$Ku7rDO2!;!pndH~;eg{N|s({(pb^_Ge|%>xO=A^y_^4lVAV0|M3rh`t}cR zzxtc-mAUz7p>;ki*1l#D(|p3u@6B(bj>^k0COM`SH`$zyR$f(XXTE(e*quyn$nX6x zN5=d0JBJ9ZQhfwdImyHAwLfl6m*g7Pe-Qsr(+mI1nobFO$X%lKetu2M_tzBr9c45% zwLqw9z<+cUC3+1p*qjkCZ#gmSgA6NGHJ=modrIIo_-C}~3;i1^ehY&OYoq8>QITcw zHuhKF$sWoA|M8bzqNhA~u6p4EebI}bY|Gh~zyEA?dYM%}FVO+hZ3UGVI&|enf1-oG z`~k6dsp`*kO?*pVo;onRSnZwk>_Y~IH>-JB#_lTu&zsV8wBC;fd=pCz=#^#Z?ymrRv&096G8P15T5`l}B|?Mh4B_D#dslR(ClT;& zu0^Aw@bsXB#bwGbM{$OieN@+af020S#bfE3-y!JG@ba#jdlO-lfdemLSShCIdsS#? zUv7qp8f+poUA#TYdD>Lp<%*7uYXTo>8m94PVrkccH)R@LiLL3%?kf1s>=3ZzO&;cm z)rjgnJ$eN`Vl@=KZ^@)>uO9#AE@`1w()8B9$vmnpy>f1w(lGq`fP z9*ybLFuH?LNDi9-Jm1JG>S?g4*H9v8XV+ml;x$wP!pO0nIodVYP%}29=~6T%UBf%G zVxrJjpoTY8j3suZcuQl*3$2{ffmNS|TXox(pY`g*4^aSK*4Rf5o|jxq@%W!(@@S15> z)ABf4HO$L_t~io3CRH;e$*u{JUJ>{)BmrsnLQKGhsi@c3-F0eAzeY1@U!`}GuVDgg zCf{(_gyDH*bk(|D)n-Uug4`T6Q5#-g(e@McNZJhg+tN@m*hJywe-12J=U$;3ap7$m zjT>}s81IpqDIs2&8(xfIj%uFQs|YW(mK6Wrq98oIiRtXn3^q+1dJ_*%D14=Fc(26a zxePYt8z$9G#Zs^zs{rL;QoQiN+tv8@W7dE+2FxZ^Y5LsKZTAyx|=gRu^m5BY49d zzBH4^yGrl8A;Yw?E<1WROw-L2yv<;fyW!JE59MZtlr$#UbXbqr%@pmIJiKTVf2V}h zX#Hrsn*su7Sp+42{B)(VwW=g&>e@h*FLznSKqPL5so)K9LaQ2>XFGYIULpm zTaBsY95oz?7pbo6a>Q^X$xKzJ%%FuM*;RH=t{JaVID8ljyn;svhpAFZvQ*7i2Z#RE zM?&SO;P8G=e|CF=P3%VQ>v}y}H>UDF#ivsP#)1S6(}^%6hfCd<0uEiJVZt_l+7xhD z_M9f)bp>#yCM&YCkN%CyU$qzc8=iQW2znW;>Nk8AwO6H!D}F<#=d?^L`)c2aN@dBW zHGiaUsHF)@NsS^`zF~?IUSO-x%n`n!Dmb>3?n>~cfA}tUN3~kXC#~G$Fu-Jdu(3Eo zB|NVzQ?WTJIZMv#X}yT#(1|!v0$t#vk+a0D$HWe zTn%_h+>$jJMUc&*=ofXTTaz?eST zC2&j{e@9cxl;jFKri?>*urh3RTxA^5A+07LcQxMGB%`kSBJs{9p;Nn-eWcvk?|F}- ze`AX6@^m|>spTxW%PL3SUBz-1cju!!e$?8TO1?4mlEEgSL+^1?@7rKg(Am6NFsIlc zp(83*Wtc$xQP82IJX7(zu#O4nFxAeoJVZf1e}^ryR!c5>WAZtCn(;@?U{lYLLBC8a z@rvj0>FR2xS35__rr3qSDxI?^o12U7>{mKxO2P#bL^{GbO!y<(83vnf4)guwZTbgS zHisRC)?TirW2!kEtz4E~RCA~=fT5rw2Uj(RG6KWZjYvyNxXlu!_?=wI-K24rL~NTa ze{q#@_{PPUP=ievht1c^k!v?u9FCl^O}6HU;%E<;Qu7hT5uG`6;|`iQd|2bEzC}qK zX2PNK`Y_lOaiqJJsUTeC8-8M`rVDU|Z%AdqyqGWQ?fjOKHJB&J-mrJ;@Gn0yHcMSU zw!@}k!>3C(S(hVLLqUbYq!|XARt>*Bf7#+kS2JcPTH3i;{qlo4JF7>#%>`^!*Jg=c zZtgGoHC$Wz*ymxBub~Ja<+WUZ-PCJX4^58}q-Y4@+^>7YRidH1Acs{p*tFN#%I~+lnLy;>N<2a*BJ{#fe?S%glX7ayzU};JU2-*MnEz_pnspUsm{$r@z$8~pV1`R7WaU64PV>F z=~b1H^klV_rz|?#KZ*VFEo#TG?YI}a$hm39h-^qI$(3eII!5$9W^y5zEYfX!cTiJZ z`|T?zs5I#vlmIF%gd`+@bOC`7M0)SNNiPBs>Cyy*5FiwhP^EW}7DPe`z4u;4DHag$ z@_sYFx%cOH_sp4dX3aYLSk2BBS=h;H&*O$yP>T5*mSk3Gzb9CEydkblgOBjg> z!L@_rmjQe=>kzx@MSV)F$RM=goJSND)(JFA6=?CR8{V-`R2yUtW-UqjR`_8>y$Gg}p9s32HO8}=+AD?bPo;u1jtla5E0Q-& z(%9Ro*<5xZWj?IFa&4;+bsl<+;_cgL^%k64iI7T0WW>0W@`WyLq!IGF7l1+ls1(RK z6gAMwn}~_w^SiQm*pG}e?aVA;1D9Up{!2fp!;>sahYwD`CS(NN{Zh$1K-&k9N@R{1s*OA0D$5!-BxDQKnup6+Q_oY)PKWB_qL*EW#TZVi1 zG-tr`-ToRy4dPjPb|SS@&A`~2c<_G;jid)ZP8h`o1T0V%{VX0m@$?xhH3ld)3_rWHG* zZ1Q$ItpML03TxsZ=ke&Q13KzCPG8RqkfQN11DS82+U?6t$!#K%Msk_ehj!YrCw+9M z%VH~&RW=(v;(~r|_oP9)7VVxL_wnH%;~3UG?OXBc*)9@B^Bslsx$r8?4JZ3}hV z%x+cpB%9f>GdoI{mmMCYpjZkSd1_Cx?b0cf z$Z1vFuJU;BQL8Q*(;fu=nP8fw@7s*=gUr<$kuI1t65O_Mz`3{-9d+u=?dHlRd@n{$ zfnJ2jU{04oCDFfD2kR)+H}&TY23P?y9al8LMWj=Y9?w1fnO#aC z@=Vs#(mmLgr|w5y3w^~kK(JG54}Hbt8h0A{DvfbLbx~d3x@`Edf0C(vkFy&YQbDSs z^|35!%F%XAZwp_*Wp|c57;ml3{ylk6N7#t%H~WtDEN`&k2wG;bx;TOP)SK;+0iQ)C zsD%J$>F1MhWgGaX9!#`6@}T|G>sU0>t^UeuUvx0`lE{IqmHj=BLj=aG3g0$mlVf0A zsf7aGfV-40#FQgzUxT}l7f16Z)jllJ0D6>hD$}DHE600i+RYhvwk0|iO9{Q%FKZJB zbQYtRdtqwo113-Y8KIAlK{OFXDWHs~2mD!nJ`+A`F_^6UmaHvCeYA{&cmXhYApcH% z{eD^X*k_xf#zj__ld49N4Lz)`_?dx#!dKm}co}|FR1C&^+j=tmnss@@!3ik7|1Jt2 zvSTgukDx?pnbV#+9cBh}n%pO34YkkEc!kyZwNGp_$6;a=aM(kVi2>H>{qYyYc9GlA z77&#V7+`K^3UvDNOiEksrVf}U`9YwB`4^p8fS!C)LN4I!qI$= zJo7MZ;afYeitYA)r#ROmn%~DuqWRwP&kNB2Lo(yu5Nr494%QioH#5i`nm?@!;c3foAjE_kefS_QJd3gj) zrwqEAlF>k@uyP@*QfP+Va9P1b{5V9qwOP9yoq^llOD9l(s(djY59SzHoHuVG5sYe@R@BtuI>S;2p|t=tpVJ0^ZTLVyTG!yGJnEXUDZD$ z+!xjxd010_JlA7G-?P)%i<<7Yp%`-h{I%bvbZy9cJlzFJGp$0nOVdYFDt=I2)yI0c+#m|%|U;5ra+4ew-O&G3F;n^>%-;~x^#ED z@`oXrcQ9X(bjM0b!7;(c-oZ_;_7xHed@Szw9YUMt2w&(n)PrO{goM&=d1enJQ^hq3 z#ajOwL&m6=)K6qx_^NTt6`O89df-v3Zwh?9x8AMF)+pEFwc5(#u%R1;Wg`1YgEsa< zwpItHqz0>c)Un1Fenh>82Rvs^lO5f+D&TD~q|uScS4NvZ;IQ-+=UyzrQ77Jc*gr^y z>78tl{n9N-Lq9=uKpPY;b%~DdzOJJHQ)*kY>~~8B*u)I~{-IE2U4&Shqp3mKayQtJ z)**-9lrTA4Om1*I&BPPSV&NWh$ULbU>M<_G(Q*bBT?35TldG+6;N^Y`z&z{Y$ILZu zm9z)OJX2B7A6*A?+cFBHwT#?$XoIh0a{zhq168nveeT{WWm3Kz^6ST^e-Kmd?LohJsrG3yDO-cAVRjrI8`;|j1+|?rqV2h0#DAA{gBCI;HnIT>Ajqul6YH&o%50{*LiOwy8*Gt&s2n_slJDnsc}lm6&a+XnF9R>l zF&eKZJ9FoMCarf1t(+EMAbTTUc@sR0dlZpr2l4|)Ki+0jAr)yXhqD5D(;>IH`M`Yi zpLg@wQcq1({giBKndpr{bw$wQjBDeCpR1S_0 zJ+1{qP@hO;shHKM<2iN@x08*xhJ#m*6*N)eBA3>X`5~%*Y};7=m?=4B(uevIlp;n7 ziy{Ldo-e+@6lDJ^ZII6Cqynne^u=aCaCo?#GLY2Nw#w9;jdM zHMo-r^M)gdKVO=sW^a7Z3V76?J39D@_pJbJek?OQk!@V~VY_Ey-HqP3B^J_^77~`s zH{?w731QQ)(pKcefvjyt?7ELxMk?lSae2fgq{BiMRS0Uq`8LR858LK!dNYx zcxd@}6&r#JIBK9d_BzhgmL%r$=OlI16GYjC!|SlNX}y&l=fdz;VOqyUN}{0(sbNhF zGQix#o~s5QJ68=pFpAP#SWm&=-@qWV|1|jc`mWVYt3=*@j$}XosQ&W@<%k2#${dai zZd86}1?}m>YwK|4lDZK#bU8QGAWcd3WzWm=_pme(@^1>kyE$|i?&|V0K~mIF%GEBk zN=!lKL(vd5n~?pUA9>dQR!=!K?Qsts!6er{gJke? z`nUGipcUXNX!{F<|8UXDAWDK#Y5 zQpkkaJuEVD3O*BGPyvP!yI+`TUnE)X>6e+Mx@-%)-4I(dW~k1i7+~3Z37{`XVp+813{Z09U9`lAc1^s`YhgEik5jwcN?t3AXC5UmAQml9{8;!;id& z3su~?vrkt+2VWnpVtLDagQMA{PTGgIkp$j))D-|iH4Vq_E`4*N=+7mDah$Rrx)=x2 z(v7?Sg1C#9bn=G)G+|#V1cUL+@3_>mMQFDgq2V{$Ks+s9{@)D)CzVjq2x(QHpzT#Khr-S~_v@iYkp+ z1Idk}@bJul>FXK+x15m^%0#V>JOGFDs%EC=`P=+lW>0Sq$jxTzsi><$QZ+iJT8gYx zfZBb1w|Zb$pE0U4%XzjBrR2dI?FrFdmB>KBTj#naD=!>V+eYT57R|%gmD&~=^(wx_ z^NJ08XzT4A!&2klIn@6teW7?yNxDH`fa58fpaSs27f+}^f3f}zK0uvUdEBFqBC_XGQvdA~~Ei*KT%ao>OJkiS2e-Xt!%|1E(@Q4%t1f(PbYS zGqlrkOB^3j&%!*lTQwrioe02NPe7_DTxm3*s%+DMi@{^MYK(#EK@nb1bO_+u~KbCNcGIE>$Uz- zB&?Fqy;r^*zn#dafuddCoRIE$f3v(io|n0lt)I@uwK3TuK%6rfz*TYJhZ{io*R2Jp z6A@Gk7AXn#GYIp2wPXHxyj_9E!6;mS14v9|@@Ks6Wk42~+uqxwKH!F|C^=*^k!dMi zO*cwt#$iMkCO##}Ji+z)+IlqF46B|gk(oVW@h7-d5|%1o8I!oQDbfBac=*nYclnc| zMW*Z6M@g^J7EA6<98vQ#_&(p^1mR@lscxz0VFw0Y`wCcs~ zjdIMe$5@&9PPdf6j(VI5%K(7@X>pC+jVFuIvRe?x$x?qs=DSXeEq4-bEmj74U!Eu+ zzAAIHH%lVOn%EEX!;7C$itAooaNxKCJ~{KguGiiBo^M(c8Ar<=qJnCUjQLX4i2y!M zb$%S%b;NsPh3!lV(&}+l!usM}v!-sg?lK)zfIhAz%m_(?-1G}UQ1;05P0QOf`UA@K zINijneJ3-$f-P8XVqj9xriRy|cYn9k8o#Mb!;H3vN%$=g z()vf!U&kYKVrWm9P`WPNP}i@AJXZOgjSD(Grg-H@b2pXb9$pvT`j%PNJ%Q;;2T5gQ zm_~=P^cb#xIX#HV)K}C%mua_T%3=DKyawKs=b%`{rBZe3$y>QYpllyheT>%!{SLIj zcDN;^%x$v4cJ`B*{*z=g)Fi7g7GXZa7T9HZzcxA1Dc&D|qV_SQdIz0PapV#S(`Yi? zPm)fcYjK<$R=xRU!tDAi-=$CCrkXg;KyfJ#q5(mRpq;hM-+s;wmA#ULH62xRri;$C z0<2Q9?Q!u=Gv}MU)+i`-#ak=}GrA*YA{SN053ee{l~fDS>6iUJu80ohwt}}SVtY7q zMfjkgvo-1nI!Zovvq8`ta52#ns8z8tmg%8OmiZ)pzuMUwZ3ByEI~adak^kPUZ;zw1 z#Sxj7YXYcRFjzMCo|rCuhEY}#t*Z$4WJ6(m$2KeM20^6O+sASZ=*E)(!bw_*jT=8Bu<&qtCijcN4rE;%ku$S-KELpN@28!p8%%kth}GqE+GsCA`u{H z`REDOpsPrZ#S1ntw!a|9;zwICrPRS#w1Z+Wy;>U(Mt_X~;V>jQ!YV()0*k?sCA{NDc4%!WFp{XVJF7`aut zUa@WMZAwu$13XK>NY*9;j!ZbTq_~7tcOVvNx_Ft%(DEhH}$6ttb@W7_vew*k*akLf9fIFl6o8yxW zfATfsR9FiKaB^-Ga{Ym4t8PMS)Syo?+7Q5EVo{mW>U}`C>;t@z*}1@qxT5yvfts)9 zyyHlJ_BY-@BiqNnT9i*_m=`atcs^g^8V6P^ql;%`i_CE$;)|fC?T8W>C)CILSvdmy zgC{|T;Ht-bmM{+7u&c?tDz6#II5v;o{~r5Q|7gAQ{299vEQB=VmRh4p)3r=J3rR5#=fM+cW6LqK$FQ+t4j*mTo8Vyee7qn$>Cc}`c1c=#m! z*k2kGKV_jUihgS~*wov1xV zbUcwJQ{kPGgNdQHrY6ruL||fmsyEZ^NTY2;%q8~Yj_-Os^mBnW*_01xQ29|SY$bB z-E&BO=y}VJOIpIwd3Hm=+G~<6=%JZv+cc%`YyI-CbI(x(f8Wh;p5u(U_E`mi*G7Pw$+$pLQH%eh6ijH+;J`JTx z!{|g&pbDkfv;3d;Pf#?uhyv~)(xj_9Bd%AL7H7INKV5`!`<$rLRbY=o?!~Yg9sD-9 z8?#fd``bSWF?BaGZBn;?hg0Nnb=d@1N19$@b5X>RAc3krw*42#nW3W;< z_Q_$$bKpFPkHvI#zWeiADUk#+<)+k?M>T{(6onzy?y?uWHl0fCvbs}^)H9YLy%?yv z!#j!%*&JDX=5XT-!)w*>{wq49o-%mimTh4UtB7#-N`w9xCOL};Ha`ncuMP|a4wI=24Hg4bVj$ z5Ewyln8O#>M9{?suiIv@qUPckZX>*`C#70=C=lQpZ#Rv{91m@)pTnO^j=8+`%&&;> z=tgU~CHGo1m!I0cts)$3R`2lp5d%d>;!?WRAYcTm{b`uw$h`lCJ5$iz3|5SZvPbn9 z-JpZ$u{70Va9GBR3pq#%YO;cup@PHk2L#w=uxpaneC*9p--ZYgTq-&JP9;2T#HLvt z2n;j(UNRAXSVp%x9pf9)VhKv^^pl~)d!^(CS;O6Hk@?)DQL9Vq_l#2e#pH`+rKv{f+N-G}I{+H<-HLuwo8xUOc#M?(_S54d@&9d$*9_u-;3(`Yk__{ClUI z%$xFWoFj1vb;LI82FiuU!r7kQI#w*MfYzx=v^4bxcB1{K|0tH%=?}Yygn6s>9%DNYCyhJ^ZOoQoiepF7khUB|UY99a$c_HqntZeW zy4`_KOr}nX9NBK)EFVmw>}Fq6sTykgkRJE6MtZM6yvdwpPm6T1NUeIh!raO=fQ)rs zbJfph@R#AD75mkjz-7Ao05~LG;&qGsj+(3Ij>M5V)pqyG;%x4^f`HUP^}`_j>x%3q z@Y*S1Xzl)o&U|J!SH;AIueXDL6X3rBZf{5sT-^;dLysuTi`}ghKTf|AfvchNGHVA9 zN_mj2VxZ58F=kQ@#Hb2>Eqf;;Aras$Xewwe-3`bC^J$LXx5qJN zk*EK>AHsE9vEJ9(M^pVLmnQh=^O(qYPKvhXL6(_v?ldw zBNLSX!$45<_`T~06i5twSM=X;*X{1@f6ro~FmW-7|H5FR(Eq~VqLBZ_#Qr}#i1`2J z{ToO8KRk2}&-J_4AreqH8qaf`7X^_3|HZx11Bd^uM;s27_=_hh29bcl|B3^FL7{)? z5fy_%L4W0Zg~1?VSNEc>g#C{Zq7axk9QN0QFqj1B@4yha7+CZ#9s~{%|C RE8Ye}L~rx)DQPO-{y%{k=av8f diff --git a/model_full.py b/model_full.py index 2931e27..c1590a1 100644 --- a/model_full.py +++ b/model_full.py @@ -1,4 +1,3 @@ - import sys import time @@ -572,4 +571,4 @@ def plt_data_full_model(c1, chose_score, detections, devs, dfs, end, grid, mult_ if __name__ == '__main__': - model_full() \ No newline at end of file + model_full() diff --git a/susceptibility1.tex b/susceptibility1.tex index 438f945..8410008 100644 --- a/susceptibility1.tex +++ b/susceptibility1.tex @@ -423,7 +423,7 @@ The transfer function used to describe linear properties of a system is the firs Noise in nonlinear systems, however, linearizes the system's response properties \citep{Yu1989, Chialvo1997}. Vice versa, in the limit to small stimuli, nonlinear systems can be well described by linear response theory \citep{Roddey2000, Doiron2004, Rocha2007, Sharafi2013}. With increasing stimulus amplitude, the contribution of the second-order kernel of the Volterra series becomes more relevant. For these weakly nonlinear responses analytical expressions for the second-order susceptibility have been derived for leaky-integrate-and-fire (LIF) \citep{Voronenko2017} and theta model neurons \citep{Franzen2023}. In the suprathreshold regime, where the LIF generates a baseline firing rate in the absence of an external stimulus, the linear response function has a peak at the baseline firing rate and its harmonics (\subfigrefb{fig:lifresponse}{A}) and the second-order susceptibility shows very distinct ridges of elevated nonlinear responses, exactly where two stimulus frequencies equal or add up to the neuron's baseline firing rate (\subfigrefb{fig:lifresponse}{B}). In experimental data, such structures in the second-order susceptibility have not been reported yet. -Here we search for such weakly nonlinear responses in electroreceptors of the two electrosensory systems of the wave-type electric fish \textit{Apteronotus leptorhynchus}, i.e. the tuberous (active) and the ampullary (passive) electrosensory system. The p-type electroreceptors of the active system (P-units) are driven by the fish's high-frequency, quasi-sinusoidal electric organ discharges (EOD) and encode disturbances of it \citep{Bastian1981a}. The electroreceptors of the passive system are tuned to lower-frequency exogeneous electric fields such as caused by muscle activity of prey \citep{Kalmijn1974}. As different animals have different EOD-frequencies, being exposed to stimuli of multiple distinct frequencies is part of the animal's everyday life \citep{Benda2020,Henninger2020} and weakly non-linear interactions may occur in the electrosensory periphery. In communication contexts \citep{Walz2014, Henninger2018} the EODs of interacting fish superimpose and lead to periodic amplitude modulations (AMs or beats) of the receiver's EOD. Non-linear mechanisms in P-units, enable encoding of AMs in their time-dependent firing rates \citep{Bastian1981a, Walz2014, Middleton2006, Barayeu2023}. When multiple animals interact, the EOD interferences induce second-order amplitude modulations referred to as envelopes \citep{Yu2005, Fotowat2013, Stamper2012Envelope} and saturation non-linearities allow also for the encoding of these in the electrosensory periphery \citep{Savard2011}. Field observations have shown that courting males were able to react to the extremely weak signals of distant intruding males despite the strong foreground EOD of the nearby female \citep{Henninger2018}. Weakly nonlinear interactions at particular combinations of signals can be of immediate relevance in such settings as they could boost detectability of the faint signals \citep{Schlungbaum2023}. +Here we search for such weakly nonlinear responses in electroreceptors of the two electrosensory systems of the wave-type electric fish \textit{Apteronotus leptorhynchus}, i.e. the tuberous (active) and the ampullary (passive) electrosensory system. The p-type electroreceptors of the active system (P-units) are driven by the fish's high-frequency, quasi-sinusoidal electric organ discharges (EOD) and encode disturbances of it \citep{Bastian1981a}. The electroreceptors of the passive system are tuned to lower-frequency exogeneous electric fields such as caused by muscle activity of prey \citep{Kalmijn1974}. As different animals have different EOD-frequencies, being exposed to stimuli of multiple distinct frequencies is part of the animal's everyday life \citep{Benda2020,Henninger2020} and weakly nonlinear interactions may occur in the electrosensory periphery. In communication contexts \citep{Walz2014, Henninger2018} the EODs of interacting fish superimpose and lead to periodic amplitude modulations (AMs or beats) of the receiver's EOD. Nonlinear mechanisms in P-units, enable encoding of AMs in their time-dependent firing rates \citep{Bastian1981a, Walz2014, Middleton2006, Barayeu2023}. When multiple animals interact, the EOD interferences induce second-order amplitude modulations referred to as envelopes \citep{Yu2005, Fotowat2013, Stamper2012Envelope} and saturation nonlinearities allow also for the encoding of these in the electrosensory periphery \citep{Savard2011}. Field observations have shown that courting males were able to react to the extremely weak signals of distant intruding males despite the strong foreground EOD of the nearby female \citep{Henninger2018}. Weakly nonlinear interactions at particular combinations of signals can be of immediate relevance in such settings as they could boost detectability of the faint signals \citep{Schlungbaum2023}. %The population of P-units is heterogeneous with respect to their baseline firing rates as well as in their intrinsic noise levels, as quantified by the coefficient of variation (CV) of the interspike intervals (ISI) \citep{Grewe2017, Hladnik2023}. Low-CV P-units have a less noisy firing pattern that is closer to pacemaker firing, whereas high-CV P-units show a more irregular firing pattern that is closer to a Poisson process. The population of ampullary cells of the passive electrosensory system, on the other hand, is homogeneous in their response properties and CVs are low \citep{Grewe2017}. @@ -449,19 +449,18 @@ When stimulating with both foreign signals simultaneously, additional peaks appe \subsection{Linear and weakly nonlinear regimes} \begin{figure*}[tp] \includegraphics[width=\columnwidth]{nonlin_regime.pdf} - \caption{\label{fig:nonlin_regime} Linear and nonlinear responses of a model P-unit in a three-fish setting in dependence on increasing stimulus amplitudes. The model P-unit was stimulated with two sinewaves of equal amplitude (contrast) at difference frequencies $\bone=30$\,Hz and $\btwo=130$\,Hz relative the receiver's EOD frequency. \btwo{} was set to match the baseline firing rate \fbase{} of the P-unit. The model used has the cell identifier 2013-01-08-aa (table~\ref{modelparams}). \figitem{A--D} Top: the stimulus, an amplitude modulation of the receiver's EOD resulting from the stimulation with the two sine waves. The contrasts of both beats increase from \panel{A} to \panel{D} as indicated. Middle: Spike raster of the model P-unit response. Bottom: power spectrum of the firing rate estimated from the spike raster with a Gaussian kernel ($\sigma=1$\,ms). \figitem{A} At very low stimulus contrasts the response is linear. The only peaks in the response spectrum are at the two stimulating beat frequencies (green and purple marker). \figitem{B} At higher stimulus contrasts, nonlinear responses appear at the sum and the difference of the stimulus frequencies (orange and red marker). \figitem{C} At even stronger stimulation additional peaks appear in the power spectrum. \figitem{D} At a contrast of 10\,\% the response at the sum of the stimulus frequencies almost disappears. \figitem{E} Amplitude of the linear (\bone{}, \btwo{}) and nonlinear (\bdiff{}, \bsum{}) responses of the model P-unit as a function of beat contrast (thick lines). Thin lines indicate the initial linear and quadratic dependence on stimulus amplitude for the linear and non-linear responses, respectively.} + \caption{\label{fig:nonlin_regime} Linear and nonlinear responses of a model P-unit in a three-fish setting in dependence on increasing stimulus amplitudes. The model P-unit was stimulated with two sinewaves of equal amplitude (contrast) at difference frequencies $\bone=30$\,Hz and $\btwo=130$\,Hz relative the receiver's EOD frequency. \btwo{} was set to match the baseline firing rate \fbase{} of the P-unit. The model used has the cell identifier 2013-01-08-aa (table~\ref{modelparams}). \figitem{A--D} Top: the stimulus, an amplitude modulation of the receiver's EOD resulting from the stimulation with the two sine waves. The contrasts of both beats increase from \panel{A} to \panel{D} as indicated. Middle: Spike raster of the model P-unit response. Bottom: power spectrum of the firing rate estimated from the spike raster with a Gaussian kernel ($\sigma=1$\,ms). \figitem{A} At very low stimulus contrasts the response is linear. The only peaks in the response spectrum are at the two stimulating beat frequencies (green and purple marker). \figitem{B} At higher stimulus contrasts, nonlinear responses appear at the sum and the difference of the stimulus frequencies (orange and red marker). \figitem{C} At even stronger stimulation additional peaks appear in the power spectrum. \figitem{D} At a contrast of 10\,\% the response at the sum of the stimulus frequencies almost disappears. \figitem{E} Amplitude of the linear (\bone{}, \btwo{}) and nonlinear (\bdiff{}, \bsum{}) responses of the model P-unit as a function of beat contrast (thick lines). Thin lines indicate the initial linear and quadratic dependence on stimulus amplitude for the linear and nonlinear responses, respectively.} \end{figure*} The stimuli used in \figref{fig:motivation} had the same not-small amplitude. Whether this stimulus conditions falls into the weakly nonlinear regime as in \citet{Voronenko2017} is not clear. In order to illustrate how the responses to two beat frequencies develop over a range of amplitudes we use a stochastic leaky-integrate-and-fire (LIF) based P-unit model fitted to a specific electrophysiologically measured cell \citep{Barayeu2023}. -At very low stimulus contrasts (less than approximately 0.5\,\% relative to the receiver's EOD amplitude) the spectrum has small peaks only at the beat frequencies (note, $\Delta f_2 = f_{base}$, \subfigref{fig:nonlin_regime}{A}, green and purple). The amplitudes of these peaks initially increase linearly with stimulus amplitude (\subfigref{fig:nonlin_regime}{E}, thin lines), an indication of the linear response at lowest stimulus amplitudes. +At very low stimulus contrasts (in the example cell less than approximately 0.5\,\% relative to the receiver's EOD amplitude) the spectrum has small peaks only at the beat frequencies (note, $\Delta f_2 = f_{base}$, \subfigref{fig:nonlin_regime}{A}, green and purple). The amplitudes of these peaks initially increase linearly with stimulus amplitude (\subfigref{fig:nonlin_regime}{E}, thin lines), an indication of the linear response at lowest stimulus amplitudes. -The linear regime is followed by the weakly non-linear regime. In addition to peaks at the stimulus frequencies, peaks at the sum and the difference of the stimulus frequencies appear in the response spectrum (\subfigref{fig:nonlin_regime}{B}, orange and red). The amplitudes of these two peaks initially increase quadratically with stimulus amplitude (\subfigref{fig:nonlin_regime}{E}, thin lines). +This linear regime is followed by the weakly nonlinear regime. In addition to peaks at the stimulus frequencies, peaks at the sum and the difference of the stimulus frequencies appear in the response spectrum (\subfigref{fig:nonlin_regime}{B}, orange and red). The amplitudes of these two peaks initially increase quadratically with stimulus amplitude (\subfigref{fig:nonlin_regime}{E}, thin lines). At higher stimulus amplitudes (\subfigref{fig:nonlin_regime}{C \& D}) additional peaks appear in the response spectrum. The linear response and the weakly-nonlinear response start to deviate from their linear and quadratic dependence on amplitude (\subfigref{fig:nonlin_regime}{E}). The responses may even decrease for intermediate stimulus contrasts (\subfigref{fig:nonlin_regime}{D}). At high stimulus contrasts, additional nonlinearities in the system, in particular clipping of the firing rate, shape the responses. - -For this example, we have chosen two specific stimulus (beat) frequencies. One of these matching the P-unit's baseline firing rate. In the following, however, we are interested in how the non-linear responses depend on different combinations of stimulus frequencies in the weakly nonlinear regime. +For this example, we have chosen two specific stimulus (beat) frequencies. One of these matching the P-unit's baseline firing rate. In the following, however, we are interested in how the nonlinear responses depend on different combinations of stimulus frequencies in the weakly nonlinear regime. \subsection{Nonlinear signal transmission in low-CV P-units} @@ -475,7 +474,7 @@ Weakly nonlinear responses are expected in cells with sufficiently low intrinsic Noise stimuli, here random amplitude modulations (RAM) of the EOD (\subfigref{fig:punit}{C}, top trace, red line), are commonly used to characterize stimulus-driven responses of sensory neurons using transfer functions (first-order susceptibility), spike-triggered averages, or stimulus-response coherences. Here, we additionally estimate the second-order susceptibility to quantify nonlinear encoding. P-unit spikes align more or less clearly with fluctuations in the RAM stimulus. A higher stimulus intensity, here a higher contrast of the RAM relative to the EOD amplitude (see methods), entrains the P-unit response more clearly (light and dark purple for low and high contrast stimuli, \subfigrefb{fig:punit}{C}). Linear encoding, quantified by the transfer function, \eqnref{linearencoding_methods}, is similar for the two RAM contrasts in this low-CV P-unit (\subfigrefb{fig:punit}{D}), as expected for a linear system. The first-order susceptibility is low for low frequencies, peaks in the range below 100\,Hz and then falls off again \citep{Benda2005}. -The second-order susceptibility, \eqnref{eq:susceptibility}, quantifies for each combination of two stimulus frequencies \fone{} and \ftwo{} the amplitude and phase of the stimulus-evoked response at the sum \fsum{} \notebl{but the response spectrum is also multiplied with the stimulus spectra!}\notebl{Also add that this includes responses at f1-f2}. Large values of the second-order susceptibility indicate stimulus-evoked peaks in the response spectrum at the summed frequency that cannot be explained by linear response theory. Similar to the first-order susceptibility, the second-order susceptibility can be estimated directly from the response evoked by a RAM stimulus that stimulates the neuron with a whole range of frequencies simultaneously (\subfigsref{fig:punit}{E, F}). For LIF and theta neuron models driven in the supra-threshold regime, theory predicts nonlinear interactions between the two stimulus frequencies, when the two frequencies \fone{} and \ftwo{} or their sum \fsum{} exactly match the neuron's baseline firing rate \fbase{} \citep{Voronenko2017,Franzen2023}. Only then, additional stimulus-evoked peaks appear in the spectrum of the spiking response that would show up in the second-order susceptibility as a horizontal, a vertical, and an anti-diagonal line (\subfigrefb{fig:lifresponse}{B}, pink triangle in \subfigsref{fig:punit}{E, F}). +The second-order susceptibility, \eqnref{eq:susceptibility}, quantifies for each combination of two stimulus frequencies \fone{} and \ftwo{} the amplitude and phase of the stimulus-evoked response at the sum \fsum{} (and also the difference, see below \subfigrefb{fig:model_full}{A}). Large values of the second-order susceptibility indicate stimulus-evoked peaks in the response spectrum at the summed frequency that cannot be explained by linear response theory. Similar to the first-order susceptibility, the second-order susceptibility can be estimated directly from the response evoked by a RAM stimulus that stimulates the neuron with a whole range of frequencies simultaneously (\subfigsref{fig:punit}{E, F}). For LIF and theta neuron models driven in the supra-threshold regime, theory predicts nonlinear interactions between the two stimulus frequencies, when the two frequencies \fone{} and \ftwo{} or their sum \fsum{} exactly match the neuron's baseline firing rate \fbase{} \citep{Voronenko2017,Franzen2023}. Only then, additional stimulus-evoked peaks appear in the spectrum of the spiking response that would show up in the second-order susceptibility as a horizontal, a vertical, and an anti-diagonal line (\subfigrefb{fig:lifresponse}{B}, pink triangle in \subfigsref{fig:punit}{E, F}). For the low-CV P-unit, we observe a band of slightly elevated second-order susceptibility for the low RAM contrast at \fsumb{} (yellowish anti-diagonal between pink edges, \subfigrefb{fig:punit}{E}). This structure vanishes for the stronger stimulus (\subfigref{fig:punit}{F}). Further, the overall level of the second-order susceptibility is reduced with increasing stimulus strength. To quantify the structural changes in the susceptibility matrices we projected the susceptibility values onto the diagonal by averaging over the anti-diagonals (\subfigrefb{fig:punit}{G}). At low RAM contrast this projected second-order susceptibility indeed has a small peak at \fbase{} (\subfigrefb{fig:punit}{G}, dot on top line). For the higher RAM contrast, however, this peak vanishes and the overall level of the second-order susceptibility is reduced (\subfigrefb{fig:punit}{G}). The reason behind this reduction is that a RAM with a higher contrast is not only a stimulus with an increased amplitude, but also increases the total noise in the system. Increased noise is known to linearize signal transmission \citep{Longtin1993, Chialvo1997, Roddey2000, Voronenko2017} and thus the second-order susceptibility is expected to decrease. @@ -555,12 +554,12 @@ Theoretical work \citep{Voronenko2017,Franzen2023} derived analytical expression % Weakly nonlinear responses versus saturation regime \subsection{Intrinsic noise limits nonlinear responses} -The pattern of elevated second-order susceptibility found in the experimental data matches the theoretical expectations only partially \notebl{too negative. The theory also predicts vanishing patterns for too much noise in the system}. Only P-units with low coefficients of variation (CV $<$ 0.25) of the interspike-interval distribution in their baseline response show the expected nonlinearities (\figref{fig:punit}, \figref{fig:model_full}, \subfigref{fig:dataoverview}{A}). Such low-CV cells are rare among the 221 P-units used in this study. On the other hand, the majority of the ampullary cells have generally lower CVs (median of 0.12) and have an approximately ten-fold higher level of second-order susceptibilities where \fsumb{} (\figref{fig:ampullary}, \subfigrefb{fig:dataoverview}{B}). +The weakly nonlinear regime with its triangular pattern of elevated second-order susceptibility resides between the linear and a stochastic mode-locking regime. Too strong intrinsic noise linearizes the system and wipes out the triangular structure of the second-order susceptibility \citep{Voronenko2017}. Our electrophysiological recordings match this theoretical expectation. Only P-units with low coefficients of variation (CV $<$ 0.25) of the interspike-interval distribution in their baseline response show the expected nonlinearities (\figref{fig:punit}, \figref{fig:model_full}, \subfigref{fig:dataoverview}{A}). Such low-CV cells are rare among the 221 P-units used in this study. On the other hand, the majority of the ampullary cells have generally lower CVs (median of 0.12) and have an approximately ten-fold higher level of second-order susceptibilities where \fsumb{} (\figref{fig:ampullary}, \subfigrefb{fig:dataoverview}{B}). -The CV is a proxy for the intrinsic noise in the cells (\citealp{Vilela2009}, note however the effect of coherence resonance for excitable systems close to a bifurcation, \citealp{Pikovsky1997, Lindner2004}). In both cell types, we observe a negative correlation between the second-order susceptibility at \fsumb{} and the CV, indicating that it is the level of intrinsic noise that shapes nonlinear responses. These findings are in line with previous theoretical and experimental studies showing the linearizing effects of noise in nonlinear systems \citep{Roddey2000, Chialvo1997, Voronenko2017}. Increased intrinsic noise has been demonstrated to increase the CV and to reduce nonlinear phase-locking in vestibular afferents \citep{Schneider2011}. Reduced noise, on the other hand, has been associated with stronger nonlinearity in pyramidal cells of the ELL \citep{Chacron2006}. Only in cells with sufficiently low levels of intrinsic noise, weakly nonlinear responses can be observed. +The CV is a proxy for the intrinsic noise in the cells (\citealp{Vilela2009}, note however the effect of coherence resonance for excitable systems close to a bifurcation, \citealp{Pikovsky1997, Lindner2004}). In both cell types, we observe a negative correlation between the second-order susceptibility at \fsumb{} and the CV (\figrefb{fig:dataoverview}), indicating that it is the level of intrinsic noise that shapes nonlinear responses. These findings are in line with previous theoretical and experimental studies showing the linearizing effects of noise in nonlinear systems \citep{Roddey2000, Chialvo1997, Voronenko2017}. Increased intrinsic noise has been demonstrated to increase the CV and to reduce nonlinear phase-locking in vestibular afferents \citep{Schneider2011}. Reduced noise, on the other hand, has been associated with stronger nonlinearity in pyramidal cells of the ELL \citep{Chacron2006}. Only in cells with sufficiently low levels of intrinsic noise, weakly nonlinear responses can be observed. \subsection{Linearization by white-noise stimulation} -Not only the intrinsic noise but also the stimulation with external white-noise linearizes the cells. This applies to both, the stimulation with AMs in P-units (\subfigrefb{fig:dataoverview}{E}) and direct stimulation in ampullary cells (\subfigrefb{fig:dataoverview}{F}). The stronger the effective stimulus, the less pronounced are the peaks in second-order susceptibility (see \subfigref{fig:punit}{E\&F} for a P-unit example and \subfigref{fig:ampullary}{E\&F} for an ampullary cell). This linearizing effect of noise stimuli limits the weakly nonlinear regime to small stimulus amplitudes. At higher stimulus amplitudes, however, other non-linearities of the system eventually show up in the second-order susceptibility. +Not only the intrinsic noise but also the stimulation with external white-noise linearizes the cells. This applies to both, the stimulation with AMs in P-units (\subfigrefb{fig:dataoverview}{E}) and direct stimulation in ampullary cells (\subfigrefb{fig:dataoverview}{F}). The stronger the effective stimulus, the less pronounced are the peaks in second-order susceptibility (see \subfigref{fig:punit}{E\&F} for a P-unit example and \subfigref{fig:ampullary}{E\&F} for an ampullary cell). This linearizing effect of noise stimuli limits the weakly nonlinear regime to small stimulus amplitudes. At higher stimulus amplitudes, however, other nonlinearities of the system eventually show up in the second-order susceptibility. In order to characterize weakly nonlinear responses of the cells in the limit to vanishing stimulus amplitudes we utilized the Furutsu-Novikov theorem \citep{Novikov1965, Furutsu1963}. Following \citet{Lindner2022}, a substantial part of the intrinsic noise of a P-unit model \citep{Barayeu2023} is treated as signal. Performing this noise-split trick we can estimate the weakly nonlinear response without the linearizing effect of an additional external white noise stimulus. The model of a low-CV P-unit then shows the full nonlinear structure (\figref{model_and_data}) known from analytical derivations and simulations of basic LIF and theta models driven with pairs of sine-wave stimuli \citep{Voronenko2017,Franzen2023}.