From bbc1945609db69acfef4b7797190e2cc2e18a364 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Gi=C3=B9=20Marcer?= Date: Wed, 10 Jun 2020 16:18:56 +0200 Subject: [PATCH] ex-1: change name from slides.py to misc.py Also add the possibility to create the median.pdf plot. --- ex-1/plots/{slides.py => misc.py} | 17 +++++++++++++++-- slides/images/median.pdf | Bin 0 -> 17811 bytes 2 files changed, 15 insertions(+), 2 deletions(-) rename ex-1/plots/{slides.py => misc.py} (82%) create mode 100644 slides/images/median.pdf diff --git a/ex-1/plots/slides.py b/ex-1/plots/misc.py similarity index 82% rename from ex-1/plots/slides.py rename to ex-1/plots/misc.py index 95342ba..85710dd 100755 --- a/ex-1/plots/slides.py +++ b/ex-1/plots/misc.py @@ -11,11 +11,12 @@ import sys # # -n → agrs.note → do the Landau plot with the notes # -b → args.both → do the Landau-Moyal plot +# -m → args.area → do the Landau median plot def moyal(x, μ, σ): N = (1)/(np.sqrt(2 * np.pi) * σ) - return N * np.exp(- 0.5 * ((x - μ)/σ + np.exp( - (x - μ)/σ))) + return N * np.exp(- 0.5 * ((x - μ)/σ + np.exp(- (x - μ)/σ))) def main(args): @@ -25,7 +26,7 @@ def main(args): elif args.save: plt.rcParams['font.size'] = 8 - if args.both: + if (args.both or args.area): plt.figure(figsize=(3, 2)) elif args.note: plt.figure(figsize=(5, 3)) @@ -69,6 +70,13 @@ def main(args): x, y = np.loadtxt(sys.stdin, unpack=True) plt.plot(x, y, color='#92182b', label='Landau') + # Color the area under the function + if args.area: + median = 1.3557804 + x0 = np.argmin(abs(median - x)) + plt.fill_between(x[:x0], y[:x0], color='#ff99a8') + plt.fill_between(x[x0:], y[x0:], color='#b28e94') + # do Moyal plot if args.both: μ = -0.22278298 @@ -86,6 +94,8 @@ def main(args): plt.savefig('notes/images/1-notes.pdf') if args.both: plt.savefig('slides/images/both-pdf.pdf', transparent=True) + if args.area: + plt.savefig('slides/images/median.pdf', transparent=True) if __name__ == '__main__': @@ -103,5 +113,8 @@ if __name__ == '__main__': parser.add_argument('-n', '--note', action='store_true', default=False, help='do the note plot') + parser.add_argument('-a', '--area', + action='store_true', default=False, + help='do the area plot') args = parser.parse_args() main(args) diff --git a/slides/images/median.pdf b/slides/images/median.pdf new file mode 100644 index 0000000000000000000000000000000000000000..30bd305cf7ecad5c2b2accf922d72f98f22947cf GIT binary patch literal 17811 zcmb`v2UJtd+6GFIB1Ms25{MumO%g&2RRK}yf?()`5PI)T1Q8H~pdeMiLKBEUXoB<_ z=>k$jIs(!`lp>&a2S3mEopb*6-+$e8!=Za$u9z=T8_{mX@W=dSpK|1_{rKtSnbr86Q+NIA zvUn@LeQf2szi&J@AKw^q)IGl!@AgyaYfS6>KJLhESuv0CTFBA1%=78_#=8fyn`o~O z{f(2qCgbBlN?r(wyA5}+3EMmq?bx5#3hw%jFBQY%KiuW4xpSu+juZlEVDfQ@wsRSI zFwbQjO$RB~e6c6e2}4tZrJC2DU5hy{osR+qUsN~Ps8P>fcqMfvzjAT+`)FWnwm^P% z)T6b?9lPs%*bPZ+vXe*58v_Gx6V)GI_e^QCrSti~=RY2^)@7a9kX%e(bx)5sGw?BS zT)XBuJ_j35+1ZtDFmr1^bXdNqL;myq%b{EPy=OG3kL`7w!2XYPg4WJc0G`_9#JMUn7pC!4Y&JKo+4h7lkcvNQJMO0NWYIu(uXrN@lXpFqSHx$DmriTMpj?=fcaTnyzNgB~je~~1#Qi@wSuqRmBYqHSh z8PsJt8kID4B##ZbbE-Dcd8V&~UYn_COa7X&Gh)#-kZDbWww)AmUON|v)lpSELbS9| zYRxV(qeU0inCUu39k}r617zUpF=rV%t{^&`rq1@Bs$VMl#sRIEfBg$t7laj!+&bF( zS!*di_in&uDV}shUS3;lekmjQ)WeC)5szCkCZ+hAJ!|54NuGG}5-K|v9aMTkfaCL; z*3=LLxowb?E`0^+re_fqm+1;B%OoPM-O(Y7)82VYsCFn8GZD*rG5Nk@7UxB?jGqhd zJ)X2GRMK#n;(3|#p~jsGQtv;g!+;5LjCu{oUuS#6O`g|(KIy_`4b`QWtl;~|OP$VZ z;qDppy$Deha1G^qcIU(8YxDz=*CV3o`O@+!IZZ_t3&N~FO9SoQ*C(A>?M1eMq2BlT z(({erQ@Xmq-HtC#a0qL>Er6>N|3;ip#Rszeb~%V+006^8n@z-PpPQOb8AFkq)=Sn) zZtz+@;<0j%A^AELDE14p=5u8ja_qrfHyw0*FMv5?HjXMEL$@b?+0oW`kTPrmSx!< z+prO`?39e%N6z&WZN53F`9>&f{6`60dgfLb= zmQDm@ls@1XyZg-@ zze*V}o!%NRaF<@1)z4X-RuaqIYe*guo&_A+)hkkqD925{*%`360;r|-tL}h^>cvKb zGZwRLfWXLYd0L=@dYTT*UVk7^+}tLBJ3Z;ZE|tDOq0PUvzd*-c+YZQ#86E^w+4X4! z&IdltqEAUnr`QAX!F0eJHP(XZ86(i~6jjsk#OaI{HUwbFEf=Z~5xF=DxV>yIsX!GX zqS79aTt(d|A>01LCpCWrg@#uR4GZlCeFWyb)sezSTcrl{Xn8;fH+~Kdk|H|T@nwJ` zSb(+mI&d*GQ1nssy_pjni0lr9o=uYWJSo6@a#mzLwm$n+=+t_W`}{at>npH3riY%Q zad~N?dZHJb0rT_>+qV|9wk>7>u8Bk1^%#(AO97_6JWE9f9G%WKV&?SaR@!b6r`q09!g24 z4<|eqnvsq)pyTQ%4;@JY!N%3Q77KiS2E9`9lv8B~#$u@304leeM&(sK1qobl!>@I6mL3eNc! zVY?5uEx#;MM>GHz-)^f=hq;Wd%c)3cvD)Pqh*ajtEsawO4jmtq6e|VWyE1g$KHr#q zWgZ=>k4g{#`?JlS;ZtRw|MtN#qf_lYJ6cw3GqEM9nvt~Bl9qKN!>7?@=C)w_>g)XM z4vgCX_IMoej#A;$U>kKu8dXzaEvVH}U;^(Ff= z(PVW8u;TKm>InGRPz+?5k<#om!m;OOuV!d;yRc|t`^3k4k7>u=TE?y=u(s<5(`KffYkalV&QgR?Yl(i1D_Ov_abbq-xC+zZdsgkS(71qF_vb1l4_lT?R^RJ9tg;d#^pLt6n5%LO(-|)`H1Z}}@Eh*dmoN_)dq&4t?WGPE&h^RX zX)ZVK@|2@&Nhnp#-N#rX(28&pcRsGgioW+EUWPgBZ_fA0bAd#)NG~AEu~8L}Y4%{r zrk>CPq$I!6HGQ)iCi!4G3vUI1_s#4XeR`T0N(d}DTIMkuCXjIl8o^Z^?|m zW2|+buj3K@tU_R0lfvbDOs7Pqy4py^9iZn(Cc<_YI-~|D%YMGF;S4jbp+^a+ zeTh#J(rB-PdL0z?t#RY-cNx7X*f-3XYco&Ll4xxjEQRLA)z91%$W!$iN{sP0O*#l| zo4n35I#A6ba4Hlfk|iiaHPZ6q9hLc~4{xFMH-QDYN77rPWsD6j_-I>XtJ9G+G7haD zNV)~&z4@XSZMfvXtE1ABZ!zx`SdhU{7dQvn|_7#o~;~BDI+IF13sq!L^&- zUD#PQuAYT8Cpmd1Y&G;{!@kX4Q8sz5O&ElG%n!DhU1@E_x)&t4_hfdGsfU*`5-w(( z9gE(w50))YM$eK(QM~DZsyA2VP`_bM6Su0ive7y_+jk{}dm^+|3A}1woB+;_=Q7ty z4;h5Ii$W*Un?%fIXLIM>Jfn=p55`eXP~Ee`sQw|3GT*FIRZPm?Q-f%gC}+>?vAfmw zFLB0Sq#F$lAop))57zcQlZ@~>QvMR0$Bkyxp_0CKd`0kVhi)rrRn6da#TgkXR+{L! zWta~%lhw@kTt{uBPre_C>ibi(pcz(fv-{l&&?XLUTz(*3q`RX@3T4j%vsa7=eU+-y z8bvWAj7{xO!_LNgTYr`?LMY4eK58?yIo*jifpUflPsy)}b7y>yE&dh|K{|qHW#6T; zJHMgehf?T#0BahU-ZDF6&`r9rLq&y-avRB(?cpf;C@TTIgg(~+w(e41e7MpBe~7N% z#JS7c>gC1l?88z8gL`N@5=B-$trTe!o#PK~qwbIRlj~ABqcizE{gcsYSB}ol_cm>| zme2D;e^$wY}22e2br#9ldsWzd94{Ppx*tsol zH=fGgFSWYdWMN17s0k55CD=a~tAg&Z2Dgt1u?x(V%JxV-nqR&zRd8^fckg(n`rKPn zZICIwbGjg*B6*u^lS79_oGHJ}(LZC9{Fhi>^ft4nsW?UC(|C~X`9Tvx6|11TbCp%F z0Aca-dxO@e_7e}n;MNc8hXR7Nzv&ASD%%3Tjx2u;lTmui5GA`yzBMR7nDXpJ*Xb#A zE1*OQG8mZRv;3a7PHmNz%VucWLf#oqJ-)wO#O*OxXig3TS!#_$rwxp5wLamfxHwa) z#!lyAj=1#3omJ<|P?Z)9bdA&~E*k&Vr+zfj<;RuLbkM^m84T{uPN{K_O@A_>hiCy9 z&lN0-!JnUB0V3a; z7!TEj1*G$NdngDk3J5Gk&M7EMmrN`d&+nLQhEeUYp-`^RNBh5?w!Z$&2QjAcMbiDy zqqZWp?)AVu=N9jyw479sKw!qw>~8hyY)%MwTqZ^E*mqI#Gdgqg)LOTC?BcJw_B>rs z5DOo1ZxE5ic3Z{<^N#7oe|KCa7<#zMv;^LMEQKmMpT^*>l_-%M%{j(W^5XOk zWz?(TXxd&6w#e6Z8sQnPy`C#g1_YYq6qHR>O?2zc=pksCKRqSn!EK!7Kr^qDNM~7Z2i^IQF6-c%r4E#>AsWbX9n{yqwv`Gj@`(TvqhdL}Ja# zcEfjOt5K#Qcw+0<#maRO0-bdWIvt|axN{_mL=2i41~e}%<9BjcPQ%kiQogtNhs|OJ z(dn}pfoJDnz9tP0pO@>^J@(HcA1)3~8{ZjgMo+y~Se;nSvlVjQWHT;KBrO$Q7+fMcJD}~7j#Jjq3dKc3O}HIq zb|W-3Lh!8+(if>TuKw1FeWm6>D;gtwKVyNudnr$>sf;0cI$|$#lt!6PJo{j547Y=U@`sA8- zeKCY_v__i4E1-==GA}XqO#uOcR$@2w@l>xKi6x1mw~KXj#n&X#&2b!7RA+3}b!oAc z&{6FzSHy?7wlxR6L%$vNU#(sW{>l-ik;qf9>rW_wZSXY9)u7EEqRro;CkxPu?1Vt9 z%u(INjT2Cq4=m2I*kF`tpNsp!{!e!*h# z*d>z{J=zlV>mh61U~(`h0M8M^V-7FZTuFWiH4G{p=^#KMno+rJ-wFAU2f?~T*To@Yt*E`%s`$7aqb+ZE;S-*08p|TkI4@GTfxO~Ri1{dQ z^k5p#NQO;Ps-CsWYu@hCkTs-f#eW53lfWK0vo5e@8XKW&NY==RfcLzNO@evijuPgR zV&4vlcr73UE)lw75JqQu@e;1sq-VABB2Sl)HOho8W;r8^fo%;HIm3>x*pj0>G)5IH z;Fp$_AMV>*H2*xfq^I5h za{w+9n8B9NkP!-H+;JnvRL}aQ5p`SZ!IXR8w%Up!dz81yk~Q3SVQIq&iU<7 zIn9bre@ubAVo}kC_b_VmqcQ%w#wWu`pL{$1vCJ$(v|pf6>g2L-QHM~#J?_>L+|2<&9G~x^GodSV7FZx@C*!K0af2dfl7+ZePTwev6 z6YIa!=k4CGTy=Z)tMm=3ARl5SrvT8a-~81;dILT&93#m+=I=o{XPC)pEg48 z@HH1NhyMCmd1cld-rBH@n!bQ{do?CwxyC8As||Mhp}hN=_pmDA4sBPlIc9qNYcaj; zCs?^#@kg1twRGQIdZo?uu-!P-LxJo&s)NM=4%1)0y2C8Cz<7St{=X=cMsK@{j?(vDckOp&lSJ+m?o)l zNS=@}YH(?pka6C4D!&)^yf$y(PWW=I?FX0Xm4Lv@QmH4GnHr}yl)-{SVQ_<(wBA9l0X-yR)Gre0a7H2ar zv{@)wi%a2Ca>q20Yo9Oe4zh3BNLIyg1$kiW9}(2Qp~V@#8fn>xE0nPbeWXoaYt!9T z5pt6>-2EyNRFGVc>$V`X6=W+n>f02nx&evsMV2&8mObZ^2+M z$Xzka5h!@0B1L>$6TFUjyvwK|Uz1xeIV8`KU;lZ?g4?m^LM2n?IPI-^rt)!G%k`^e zyZ21ww_n%mjmqCFsV5o3u3v56C9{#om(`bkls{KdZ#8aV<|5u!#pLpdwyu8Q-ly)4 zb^gU&Pj9Z}xAi@frZ?*9gQqMU-Nd`!q-=e9b-W%yIvq{DK|!@QcR%QSYvz|(3&?h6 z>YRnYhxlOiu-C*5#K!P0AF)mic^TldRROH#&E;X0R<#$4E{x5Gd<|sI z3$&dj-z6yu6OQ=!DDG1i=`4dZcc8%kjv!%6#qo4Qv8GnEjqWHfA=X|HPaiN-_ci3LNSJIic6%&kYGU%)Jl zH)n#>=^+*}j(nyF`{|Lh!$t)%s(z%kXq{;pND2%6*G$0c ztvUNOy{X948f(_yOKa%B&Lt~ea9CILl{7=U$a6cEd+#WUVHDHGh-+YO)V}FynI5|H zNddG*@(uCmHO(>l%+`be5%$-^HDgVx^O9AZy;nK$2N!%pSP(svD0 zbJdRjcHt6&W)}yQ$b3BVT<2I#{Pc#Y;2HLERpzDKqETaKmhisV+i{%_dqszI{Oi;G z!mno3eF?uBs>-JQ=FAzp?7r(R1*1Elh;`}+_Op=H+?~{S>8!Tl>|~1Z9co=j($sy# zxMSw2_BtF!_-QFF#7D>O?dzlnDXk?dJyBGs7LQChO%#Gtt76zWYnfr5OwJ?r8dk95 zxNJWp7Fe0UJ@!aTxSUs1%bV1x%`srRsnOuv^MA?7$ z<#2l%)+85|U&rbSnx(M%)x!X%3k2dSL)f1uR-hHKynXY&o1WTB3Ct+}g`|`z?pU+T zT)Ugmxm!$zD znQ^PR!;>#RX-6f4z7;#4)ZI6MQDpTE%Wb30>&1kb?XqX&Z$ubAZu3XoIrvN`1ztOZEVe!^Olb=HB&y4Vk@p}vV<74>a%=PkPQ^kZt8z{@kBi&iT`qjqQ5{Tou(fxdT ziRmH2|0;cza7;@Y3eG65oP)W|;o28hOEu-xJm;D9k|*Qf$wcck=l0PR=LbmUo+av0 z?f&q9_9j+gB_DRY*?m3%S#JL6n=^eH>{Z;2S3yPJl4;7HPVgi|W{XkiSg)=6qrFNm zxbPq>X1JTd8LG6KUW!VK$YrD#DB1W`&BX~j+Zc|#2>xNT_17vqDo!A$8z3;AB zr?Lz!e!{4f1=zV2^3N`Ae-%%ic?GXAZ_!$CRJh6OTiTDso7gt$f@X{jW*)U$)Qh|2 zXu%ABsKk!=QhG$}q@E7o3|A0oWXpPKW*`SE6IxbCU*FuZ61~J2{c`&`nQr#c-1DDb z7O-O5XC`1ddzRPHX676Rm%;m@#oMR{^NCw6(V8>Ddd9}~eNv(sl z*E}pvuL-zvZ;tAh)_h66N&-#hBwqvMjL*=fDIXp=g zY;NN>vjcv^k6!zsvozqoTLW@0aL&KLaa&H--hY&i&Qe0ijE&r``@_Kf#>xIfD>vHy z)aue5M{w_4ig8kFC8*q8aji>*Lx87dt6V{6SqU!xZWGNPc8do7oRL z@IdSZFIw@~YrP?MfsoQkC+gcHq^J<1xq^!G8U>g3RlHrtu9)2)vUje_6FD`Sm%%Q7 zzRCgqF81I~hG3YO6O?=CX<@?o*Ela4TxRa3Vfrm}YAL}Q(I2KbB0)0rVr5P2*#=j3 ze@m)zX;~Z6t|I{2>g(U&O*+dF$#k2ty+dEMJ$$*vAh3cF{qu%Pgasr^Udrs_Zt!)H z7$vkdf6q6)+A468k6lFlK~10BmOFO_eq%R4Qaw$JxBk{zpU`X9Leuo3q}4eG_yO z_wX4zZ4n)1H~izsH3eGXr_e|aTJxI2mS2`Su`g7cV>6cOOEN`8QppHxe(uI1+8D*< zTY(F{XAKmy2-a_JR+F>C3zjR(^`-O{!-G1`%|~$Lb)4!SGDb{`x0iF0lftva6`ki) zYGSn}zT%V5xC)^2q)*F~M9`dVNtHyNeQ;$*f=yJYn@UD%f;7HDawxjy{$Q7fal8q{ z7)R;!>AJ^lLg2^*`A|)otTlO=bZ4ATM>J>9MD%>eXwYYMPN!Q>-}HH!HxE$xSW%;z zM&rE1v$!g|stJU*wcd@$D!-lt7xhysFNCD`cA{ZswrqiY>upl%E29iN<^zV56f=fa z11KD+(q$*A4FTWjlM{yjl<79fyOo;Lcj|gv#8D zU&fJEx*k)}@r=cZPh=y4cbrpG-g3=LL8+ija4G#9hVEd>D#@LhE~6rW7|q;qM%Wp1 zcx?M?%%$;I*iDb>)98y(n~n#}7Gs=AO=D~Nd(rIZ#>{K|vNq-~KbjYQh-2cxe)=$o zUHhzfm$4t$SCnTlXY$y1Ip|c!FF2oE)aU0Ml%_F{L>>u9bx608G4o%@v%8aZ|KK8x z1rP2ym(!$;hi6yM{m(ULZ(YAfiqM$1nJWO}q4Jt%rnc#Jq#|bQ^xSp{3{eJ&<#@YX z7=1$ij7jg(4J%|C<%>GmFs|j)AUWe9T@sp9@M`e8-Q8GCkMl(LBWgVQjPSL|)I64oa%yeR#aA z%VI5RdLT!!^O?Zxzr=+s68%K|DoT=K3T1_H{q*7Wd|a;a)gX{yBvX^`OnFJD9Df+qA+8{g&S zFI2+pv4)4KdrJHxtUYl&^$$ha+d&Vn zpDUkq+dVIk{A(^cgWXAbFRRU?gNv=>ZTAK0)?XETSJ{K;P1vj{e>{ac^CUvF&Ns0X zzRYmhVczmup?O~_LZ!oOc${yqW+4z|OdDg>+uce7$#t9gVEabUi&dkM7rw>W(PDm>q6sVCKYQUGb)uuXI|3omWjcWwCtva zi`EKwv~hBu#i!-73|>|~v$0MQ%I_-9RlVWVTbrY2?mmxQ`ti=UqDKX%#uAsKO|oM$ zcH_wb2c<77j>ksC12YyLPrF>B8)Aw8# zV{4%FDm}*{Ir>_>ctcv1kVd+sO`&E^C?dqSl#jl)P#M zJ6cyCQmzFQUepxrsQcc=tj@6E|IKGAsO5=;Z&Cl2>^wdEFv_lcBa(g@$+GmlJt4V$ z+1cWm_IC#wffSB>U7m6YDgl_5|e8_t}Evd4B}^4=t)#MB#;JNBt(@~)&R zwjf#Y^0^dLtfhb6LLqH7MA8Tg zg6I3ZyqfKI$V964y9ZEjIW=!*S+c*t#r3Su1bOU;RX;d?{_0sqd zU7^J-Ta^d+H zEi1RJK*cZKQn%N{eb492N+|Q&sDffWZ-$P5G#b=PZUJ-+n{$92e!0%C?n{Lk+b*HJ z!FNi&+^?GxpO!_@?hwLEDIUmi+32zrc3_Yc9`5<`LU%!+<+De#cEtjD2^u+Pm_<$p z8l3WnYIFFzwb3a_q13|(IUjnyyi9t{6USaW#rYNY$*Np4i`jL_$CkA}HXy*6&%)x< z^9n61feQmMQicMhy!XY-eCEcnl|huYR!%*pQBHh3PlRbBuJ{gS=TxRW=OsU6T*@>6 zo@He9hdh&*rIq~l-nkH zRGVSBu{*0u*=A2H&$91)=#j(jfQ}?IGyM|B92+ZYl*7~s60U5Qfz^o0KEo>F3ofK#TAO} zS6|J;Hx){=4;472BZMN1&oi0I#$STJM0liqJx#ll#=SH0(LWDA&hCd2xSz^fzUpf` zx0}j?ya)WEqqrDlyx;s(GU#1ObFKhw?n+k<>Y|hyDO=39_5&#a+w>UsT$R~W_Qcy> znT6(SBj(et`uBbG=88b8s=Apb*{c?*WjY>rzMN*8_8!oBj3>WiEgF-aTiO@gvh(9r zMcv4XV94vH!*e(<`H`fFZJ$Z3Pa^jFbi!~|9PKF@zKC4x$Hy*!E%S!) zbohBpWoB2K{K&l-WH$k9`7jPGrj4z{A}k(e8N7akt<0L*uKXc4VzQ1DsYM@JcUrWfTm*|6J>_W6$@^oNpO5*A* zZan&+AG_cppX|1Q%&J2l9Qy^YGlHCSXOZ#@gz*wti>F-q;wPLLDT z67o=yFb)C`Bmu8;^3g5^;1b2x1}>psr`k2-(^dI~syS?>hab+R0Zl;x61hKvB#n}9 zc*%sleX+G*jvLFu?}1B=L0`4n*pC^NHS#8G(@-D`C*I*-YGGL(d*b$o!IlG;2r+R8 zqYJ(GIae%8t`{zK3jB=pay|{KoRR$ip6oG}rOF!zSpr+0ze%8ZE@xCdjK5coWqIdQ z_p(8Ug|}hW?(#=`OFxzc=Zkw^g{Ei$i#Uf6ZaH2GyTUW`y?BF83Ot)BN? z1c5|UXOYWv1dR|`qrx^kqljN{E(#>#Fz=7?&kzQGMt1TMG^7wl)+^W&$z(KtJxJv0 zBJ#8pL8B32RJMWb@5wUAB?XH(d_fAS6Ev>N8T|lW>)7PaCvo@`N*1E|&A}p9SCOa9 z2^znE*E*uUxDRD${+D187vQyyD?x-WB>(o034FpjHYwc?r`CYZJq8zD2VUzWkP(EH z-Er0w1YsBVx)n%9p?wQpL)3NO_|YtKfeUN6I*Uv?k2SnDhcp!Q3m$$YYjpMDyy5ZO zty~GLA#xsxk@cJCdkH28Z+qbGK?uS;o;Yhhg0QwHP9;xXsVFO>P^}A3QwScYoV zf{H3uu_Y={f-sjCju%c4zTt%d?bO@|TK!OyM2BgNw?hw;jo= zhl=c)@N#sZBA#Y^Ivc2HxEasMx3TY-?T*W%B{%{vFGJ1{iVj)}3A&O--Uz)7wsd)i9j;hIg8Q(y z&P4BSE9dTrniUN*O7|T*x~aV{*{UHhMP5LpF*H@bOUL=g&l@(BEh@ZVeBTd z_Ye^HLp~WAIxoaKqj!%S%db0ww5G)J-*G3T=-oqi3d<%_hof&u%O*F4qdTumYzu*B zZ%Hs=`88&djI3CG|5@Z$c9|pQA$QRQdcxqcyQtw=!XSx9-6_=!4fD@faWO^$<0TK# zG$z8JnujO_3t`Z-X5P23TsArRAzG+W)*)m~jXpM4dI33p3!7`RfP8x!o15T6$j?bOSaw4so1~x}fLGy+sp#cMpu5B(GTjEt zk6J{=*kSn}Esk9svB;U%!Bi6lgXNVDi&j3(oBmYhDW4YVJKLufSHoj*J(R^J4cy1Z zefMCW{O6wk(PUWp_6xg^wxDF}O$Kf^96AS3)n7s1|Dcv>O6diDUNdT&*HhVFCt!bbhqUDcagETV{A6J4+GEd1ofQ@@>+Qj zxAgsM;UZJ;OyM^!@vxj7nXE5_kZdjd+rAPb>#)XSb^Nq{LAb)*`RI@}nT{hBkD;cz zWA#$5$6d2hF?aW`J5RrjTkX|#7s_0U#mR^vXBV&b3P#+nJf!b1Sy}rme|P>y$dOF- zLVVIV{TL*p@uqKMbIh3#(XR(wx?krT2l59IETLr~s9!l3Nk0a0ZFTlP^S`3}lj|tv zQf@t*l6|T^Wby6eUk|_71Ro70Bp!+C4yiPEvKHKa5D{u%YeTo(#gFx&%(wa6-*`Oc zX!&#NN805L*#d`O{>LBA2S4j`i%@i(wld?I{&+Ru4u$8Isi!9G=L4{Y!lmT(QNn6I z)F72{u}E|7;A&qR_X5N1lIfjK&&%H_Yvj!wJl2=Hmwk42Rix1QwXK%@?CN|8bLWg^ zs$qY>_=LpExZ@?^iMt^Pe^&i1gUPo3P09FWn~q`Ttt0mDN4p;Rj{BcS36-OEA2cl{ z`}*VWp>KyC?_Jk)S^MD_@Br82nsCNHL>Vj{e%C;Qy#H(sRYOEh*~-$Jc@L12NNrz8 zNU@Upj3*RVUP_DQWKM)1<72wWbKz%fb>6jVY>S^+O6ugy;^zK>JO-8G$MSa!E`O#v zoLKzS3ZApqi1frQaEzbFg80x=qD+`u&r%h(19l$j9^J-yIqt_sSjO1AeR4Z=TBzu@ zNZGREuG9|GQV(0Z!5Odopqks(6Hcca3|c>R*VGU2eWG=&KX*Sc+pZ#I>f^)ax*@>; zzIc3!mH%knue^nz>Au6_U!K$Kd&T)*w8h(hjojRdQ4ZN&NgnXMltNcLMks1Oo|XJs zj>a6fQ)@oT%DVq^Kd|kX?#;aA&%4KQ>l`7w+sdD|eV2%gx4*boe{tVVI9WvU77Q*e z3;V-}1NdP7H%ga;q|CqZxwLQAx~ReER;QWGW=K|@0tsBXd$(2K#{rQ>Z zI!M$`ePnyoYn2^~PkPku$^D8Zfx^UUbeZes!+Ro6MoxRb9BF&oGgGv*>no!w{&nWn zTmL(BtD2=PUvh&)gNs#odcJAN``zHH-^yb;XM2`HKP~fdlANOToi~?%J;R9XA9nWd z9MiH~G!yf3GSXz_o=u{BUVISR3lFWX*bI@7TGG(` z4@mz>e*E8{$U^_k(N|+@>6@U`YC(Bcr$Pj5pIuBUzsm*Dm~hv98&oaFC}A*#ii|py zMJF0*uuYkkZ6#H$t{<1VmR;&<=@YgXXi)Thk4jkA*VjX_x0~bJ|A6cNN-C5@{2MN~ zbcG&N4=M?zN>GoyocuUiCFVTab&|H8g?pzLLyajfbH0?NDi(+mp(pA1LND;&f+haR zzWm<|A_4F;|7m_Hy>5C{71|QRjZCAm2Tq4k+mC6I#|( zlG8d{b*T~0OZ(b~f2}JWk|nbFBCQNRFg9o9&M9CK$|_^^#UYMU%6s-+-lZg!uy+2f zv|yj0FAcgJ{*D};zNgNqc3Pcwe9NkhtWbFNNQZRpS4iuI*yk4sJ6%~^2V5EJTk_?u zvx1JkyaS#r5+BMX3GCENSGg4euI9PVm5=xhOEP}2rPWC$E!I(vJwsx}YhfI^R;-55 z9K_~6bE_lloxMzAQ1RcUfe|UO|B7mq7v_W;3lL8oQAR!n0N)fKqS~W8iL_WR2<$JC z>{VBQ14?AQ5_zg>4z{)!A}1AKl$t^$r6J<(7*7W`8;GPVL_84V=>`#ZbHzZU;1F?J z2VV^F@%T;J{nso+_Nuyrrh~fY32?M&QVZ6|u4({G=o)Fjxn-`b{Xx2t~ zdpexZU@t+TM0zUl_pdWRx)$m$NfJ_05OE!p3&3XmTi<^?oJ5!a>fdO8+vT={jki5u zh~MmE0QjGjOO*e|j90wSM7}LT0s$Bg<$fLGU}x`5)Qbq1ct=bE4wwm`)7p7K;1a)4 z{~z@IPpAM78K6Bops%>vIb$GDfC7y2KLH4sLs9~u40~f-jEHjo*5Dr|6G#3F*gy0j zDn``96e10U{{K;E{=d!xhsi=D0sl!#Lm{$oS%@@Tng$^Sg~$L|OF~8p0!KhkiWGp8 zc%RrGE)Dx#!)1v5|F=T(=RO?hODx1W{uC)V3<84^4InEC6ksBREJQ{IxQ0tWqyXVE z#3CgDkw8e%00kl~B@L010ww`g0vRCiWn=*WP~eA2`Wh$7*XvL_{2E)H>gCNPJ#&tDB`Gp%|ZkPh0>hd`)dM(Gy?Lc5Ch?N z3=res-^+=ufw&;X#6QYOJP>2zPdSN|zhZ(I5C1GgTqoDS9K_i8Q;3>NL!|&80$7PL z@IMrUEZ`GjTmWhLpXEmw>j=WooCFiF4z|K27g1LVQ~Y?JyQ+9aSr$p8CTX&E5T z|3{l73|J%nvrY29Z6zrUd}aF2ezMU2)B`4okoYfrGBC-1jP>?JIXGiHPc~1v;cDv! t><#k!QPp>IBSz=RV{70Lhyk4Z`|B^t)BD8jFc}F1jE0~8s;)ZC{{u9schUd= literal 0 HcmV?d00001