From a951e0a778319b31418bd02c4951363eabb44895 Mon Sep 17 00:00:00 2001 From: Ayush Rathore Date: Fri, 17 Nov 2017 01:25:01 -0500 Subject: [PATCH] Problem 2 --- README.md | 63 ++++++++++++++++++++++++++++++++++++++++++++++- data.csv | 26 +++++++++++++++++++ figure1.png | Bin 0 -> 9397 bytes figure2.png | Bin 0 -> 11826 bytes figure3.png | Bin 0 -> 12356 bytes figure4.png | Bin 0 -> 26855 bytes least_squares.m | 13 ++++++++++ problem_2_data.m | 19 ++++++++++++++ setdefaults.m | 3 +++ sphere_drag.m | 10 ++++++++ 10 files changed, 133 insertions(+), 1 deletion(-) create mode 100644 data.csv create mode 100644 figure1.png create mode 100644 figure2.png create mode 100644 figure3.png create mode 100644 figure4.png create mode 100644 least_squares.m create mode 100644 problem_2_data.m create mode 100644 setdefaults.m create mode 100644 sphere_drag.m diff --git a/README.md b/README.md index 1aea478..94e06bb 100644 --- a/README.md +++ b/README.md @@ -1 +1,62 @@ -# 05_curve_fitting \ No newline at end of file +# 05_curve_fitting +ME 3255 +Homework 5 + +##Problem 2 +Least Squares function +```matlab +function [a,fx,r2] = least_squares(Z,y); + % inputs: Z matrix and measurements y + % outputs: + % a - best fit constant + % fx - function at xi + % r2 - coefficient of determination + a = Z/y; + fx = Z*a; + e = y-fx; + st = std(y); + sr = std(e); + r2 = (st-sr)/st; +end +``` +set defaults +% Part A +xa = [1 2 3 4 5]'; +ya = [2.2 2.8 3.6 4.5 5.5]'; +Z = [xa.^0,xa.^1,1./xa]; +[a,fx,r2] = least_squares(Z,ya) +plot(xa,fx,xa,ya,'o','Linewidth',2) + +% Part B +xb = [0 2 4 6 8 10 12 14 16 18]'; +yb = [21.5 20.84 23.19 22.69 30.27 40.11 43.31 54.79 70.88 89.48]'; +Z = [xb.^0,-xb.^1,xb.^2]; +[a,fx,r2] = least_squares(Z,yb) +plot(xb,fx,xb,yb,'o','Linewidth',2) + +% Part C +xc = [0.5 1 2 3 4 5 6 7 9]'; +yc = [6 4.4 3.2 2.7 2.2 1.9 1.7 1.4 1.1]'; +Z = [exp(-1.5*xc),exp(-0.3*xc),exp(-0.05*xc)]; +[a,fx,r2] = least_squares(Z,yc) +plot(xc,fx,xc,yc,'o','Linewidth',2) + +% Part D +xd = [0.00000000e+00 1.26933037e-01 2.53866073e-01 3.80799110e-01 5.07732146e-01 6.34665183e-01 7.61598219e-01 8.88531256e-01 1.01546429e+00 1.14239733e+00 1.26933037e+00 1.39626340e+00 1.52319644e+00 1.65012947e+00 1.77706251e+00 1.90399555e+00 2.03092858e+00 2.15786162e+00 2.28479466e+00 2.41172769e+00 2.53866073e+00 2.66559377e+00 2.79252680e+00 2.91945984e+00 3.04639288e+00 3.17332591e+00 3.30025895e+00 3.42719199e+00 3.55412502e+00 3.68105806e+00 3.80799110e+00 3.93492413e+00 4.06185717e+00 4.18879020e+00 4.31572324e+00 4.44265628e+00 4.56958931e+00 4.69652235e+00 4.82345539e+00 4.95038842e+00 5.07732146e+00 5.20425450e+00 5.33118753e+00 5.45812057e+00 5.58505361e+00 5.71198664e+00 5.83891968e+00 5.96585272e+00 6.09278575e+00 6.21971879e+00 6.34665183e+00 6.47358486e+00 6.60051790e+00 6.72745093e+00 6.85438397e+00 6.98131701e+00 7.10825004e+00 7.23518308e+00 7.36211612e+00 7.48904915e+00 7.61598219e+00 7.74291523e+00 7.86984826e+00 7.99678130e+00 8.12371434e+00 8.25064737e+00 8.37758041e+00 8.50451345e+00 8.63144648e+00 8.75837952e+00 8.88531256e+00 9.01224559e+00 9.13917863e+00 9.26611167e+00 9.39304470e+00 9.51997774e+00 9.64691077e+00 9.77384381e+00 9.90077685e+00 1.00277099e+01 1.01546429e+01 1.02815760e+01 1.04085090e+01 1.05354420e+01 1.06623751e+01 1.07893081e+01 1.09162411e+01 1.10431742e+01 1.11701072e+01 1.12970402e+01 1.14239733e+01 1.15509063e+01 1.16778394e+01 1.18047724e+01 1.19317054e+01 1.20586385e+01 1.21855715e+01 1.23125045e+01 1.24394376e+01 1.25663706e+01]'; +yd=[9.15756288e-02 3.39393873e-01 6.28875306e-01 7.67713096e-01 1.05094584e+00 9.70887288e-01 9.84265740e-01 1.02589034e+00 8.53218113e-01 6.90197665e-01 5.51277193e-01 5.01564914e-01 5.25455797e-01 5.87052838e-01 5.41394658e-01 7.12365594e-01 8.14839678e-01 9.80181855e-01 9.44430709e-01 1.06728057e+00 1.15166322e+00 8.99464065e-01 7.77225453e-01 5.92618124e-01 3.08822183e-01 -1.07884730e-03 -3.46563271e-01 -5.64836023e-01 -8.11931510e-01 -1.05925186e+00 -1.13323611e+00 -1.11986890e+00 -8.88336727e-01 -9.54113139e-01 -6.81378679e-01 -6.02369117e-01 -4.78684439e-01 -5.88160325e-01 -4.93580777e-01 -5.68747320e-01 -7.51641934e-01 -8.14672884e-01 -9.53191554e-01 -9.55337518e-01 -9.85995556e-01 -9.63373597e-01 -1.01511061e+00 -7.56467517e-01 -4.17379564e-01 -1.22340361e-01 2.16273929e-01 5.16909714e-01 7.77031694e-01 1.00653798e+00 9.35718089e-01 1.00660116e+00 1.11177057e+00 9.85485116e-01 8.54344900e-01 6.26444042e-01 6.28124048e-01 4.27764254e-01 5.93991751e-01 4.79248018e-01 7.17522492e-01 7.35927848e-01 9.08802925e-01 9.38646871e-01 1.13125860e+00 1.07247935e+00 1.05198782e+00 9.41647332e-01 6.98801244e-01 4.03193543e-01 1.37009682e-01 -1.43203880e-01 -4.64369445e-01 -6.94978252e-01 -1.03483196e+00 -1.10261288e+00 -1.12892727e+00 -1.03902484e+00 -8.53573083e-01 -7.01815315e-01 -6.84745997e-01 -6.14189417e-01 -4.70090797e-01 -5.95052432e-01 -5.96497000e-01 -5.66861911e-01 -7.18239679e-01 -9.52873043e-01 -9.37512847e-01 -1.15782985e+00 -1.03858206e+00 -1.03182712e+00 -8.45121554e-01 -5.61821980e-01 -2.83427014e-01 -8.27056140e-02]'; +Z = [sin(xd),sin(3*xd)]; +[a,fx,r2] = least_squares(Z,yd) +plot(xd,fx,xd,yd,'o','Linewidth',2) +``` +Output: +Part A: 0.9801 +![Part A](./figure1.png) + +Part B: 0.9112 +![Part B](./figure2.png) + +Part C: 0.9462 +![Part C](./figure3.png) + +Part D: 0.9219 +![Part D](./figure4.png) diff --git a/data.csv b/data.csv new file mode 100644 index 0000000..7633db4 --- /dev/null +++ b/data.csv @@ -0,0 +1,26 @@ +Data +Re(e-4) C_D +2 0.52 +5.8 0.52 +16.8 0.52 +27.2 0.5 +29.9 0.49 +33.9 0.44 +36.3 0.18 +40 0.074 +46 0.067 +60 0.08 +100 0.12 +200 0.16 +400 0.19 + +mu=1.78e-5; +V=linspace(4,40,1000); +D=23.5; +rho=1.3; +Re=rho*V*D/mu*1e-4; +Cd_int = interp1(data(:,1), (data(:,2), Re, 'spline'); +Cd_int_pc = interp1(data(:,1), (data(:,2), Re, 'pchip'); + +plot(Re, Cd_int, Re, Cd_int_pc, data(:,1),data(:,2),'o','MarkerSize',10,LineWidth +legend('cubic','pchip','data') \ No newline at end of file diff --git a/figure1.png b/figure1.png new file mode 100644 index 0000000000000000000000000000000000000000..50aa6bb99a9086e4a95a28bf66ef644a4fc6a992 GIT binary patch literal 9397 zcmch7XH=72v-VAh5SkR}2qX|hr5KPRB1I4gA{_+jNS7*ADH0G6sS#9AdKVF}&?Ci0 zmo7+GKuYMn!`XSB=Y7xl)_TwS{v8%_FZaDOd-m*^Yi6#U$2wYSbkv;G5CqYwW3TB! z5DWuBWVr|!7|D!@TmcU%59}>32s-(d^pEVB&`B;ZNa?MvsY3aihMrLx9iP#08w{bn zRgAs$-R-^otUc`@&8r4Ct||#54IJ%|n%3TqM(&=@UcyLqHyhC_#-r?DJ|AhmvZtN3 zx0AaY(!j~p4ib|T(Q@|{5yyy0BgMooOG;k8C>6BIX$e6{Nd4MXgL|oqqrp06o3+QQ z@u)!MK0YQkSSu~u_R+QJJ3YD5X0Ja%wd=9(_ywGqT&ZhaZL-n6J=M||=cU7*IQr0L z1s<}v3dB9XcKr{1{kT)_b4AgpQ(lkllEV{_3TFH+7&XzuJ5s{U+wqC%F3-c!Zx?n` zz0%|BJ7fZE$9U4#Uv@f-@9&Oa#Kce#)HQsimIi_p5ilwUNk#f|;{WPT_^v-21SNab zc`d!m>I`F8clEp0C@^QPjV?(w(rJ*!tb??&r^+XB9Ti zOuuPpXy@QyIY0+NBC&lTy162AK>j%16v0^A_-^O%^;(*bp}G3LqDgIY-fm+b?>9Lla`@d~|r$J!n zNGkIfyn_mJZx%5<({jm9T3$Z=#fxg&7A6(lh)0iTsHt1q+JM&2M&0O*S|DUvdrD)* zt+jH7iMI|{zgHil!e)s?CZ!mO{*!eZxeR=N{+#kGm%7nBvoQEUKtRBuGd(Rinf>(X zkK1Db&UfyV%u^Wp`}^OcP8U0$Dz__?icxxEa4CT4<|X28**qncxdL|fzKh}VE{%E- zuTCjp`#`a_#QlCE&@3rA+1SwVS#t8(u*mrM&bBsZSJ&DS5E@VJsuOL~+?z)UhqB=< z(ix?rf4Kg3@|!93!39yzg_@ifzb08IwpPz1aA++IS=sL~ySLV9IG+$`QZ;lpt2ra7 zI-XF%7y^f;yEj@qI@${QCM63irQGV(#{Em8V>Ia?9LFVZ&xR@|G&_n6C&SD;cg5RN zQEn_HLvOE^IcWX`8Hg;}&ibqpe;d{ARyY{pAdR%wXLT5h&?YisAuEGM*H<)AiZDo7 z+*{dP$9s?rxB8TO@kh@qN4rp}9xO%)27M8`3SNO2VOQ;v4A~Gl7nb~<9I`<;zPfme zlS&mxK^%$Dw!G7g9eys@q&jD+*9=4>W!EM9PJDaH4Pn1OrSrKe1|;3zEv2%}kEqE+ zpdtK|YA#~v?bee}uFW18BJIHp#!2)~yS3MUN-72qDe{q}EZh=|g##O`os}Ij4b0V- zlkx+)+yH7ax4&x{a{(63yhtO}cR~e8hD!$~><}v0V=l#fI}`BXDj0GT9@b%F2DZwi z7n1&a*2x1d;?{J5adDRD+580vIpX-21#l}_S#F$QmtpD|6Pkkq**vBn!3%%Hcfq}KiT52q zjZX+FNa=!=CX$m_asq4xrK&vIjI@S$%?Smc4#`URKt0JOMK1aE5(P*ClWH<0+HY@jXR)-W&m$ zJy9*gv?uGA+jlj&G_;QJ%lng!xPL604d0~BheMyOyEH^iIY-D8D zJ09ryMEw~&L-yT$&mn^HXMD)iroV=U27y2bXW)r1$jY)l+}reLJJG=O_u(67xrwr` z$)X9POWv2cPpaA5+Xn;$sH;;cD8U#g$b|YfFMrr5ZnC8vwDXHy@RbFps`RI;C1ip> zjrcB=&Zjr=_?ct#Iavtt;eAf63FkYd(bHm6z>C8YgK3mE8)LGcQSJ5H+ystt<1q{I z^J8%EzxpC-njSL#XB!S`m2atp3My`zC!hJ=hCic2Qz*qFi)lAsH_3ghbyH)yucT;0-lV z@!mX_Wiok|ziV5Y@&|z^6FzU1C|4)Z!>#WN2{#3|8A%N$Qu3prex|o)bS5kQfT-Mi z4rd90Hz+y2O5^~u-efO0mOfW?c`t|pp(`8`ad{WWxe=$i6%z`!ZO{4^XpFuH-rb0f zxNL*ML&>M;d?H9XpUTeahy=cm;X0NK0nyI@qr%eyBt-f*Mtvg1H0|Sv z%j;licI|&+u~)}|6q^sNHIW~(fmNQsNmj8}LXu%9ZNk!ECYn*08ECG*10c3K%GIbk8BufP9!a4?14M%y|iTltu7bO>o zpVV|Y0@n3kOWPw**k?1HVpL4S zQC?RY&%4ho+9*ex1m~1vKsoZ7+Yb~5{x_Y@4ms-_*(QkL(}Rz8^E8sSwn;xiQ4^lk zARiy_+TFf=8`Ks4r~d^noD#r&nzItghn6j#1^eq9%n+35vb>hKlI(UPMrNG1i!#Ql zVz}dPfXd?ipbG#TNl8hi$PW&kJPHr4T)sS{fZyBS*G|83%h>o;dR$zb>a}aL`a(E{ zbN?^S@xOR>hvs2vi(HRE(dY|cCQb#=Xa_bwF` zl|^Y~S=q)CVOdKko3P|(xL9}h^u zOPgD0aPww~X=TaVw_R?PvT0h6zC3yHVsLECLQn5=JRW+urk&HmrG9PZ5;#~BZEjs# zH|ZNUZp@m0F?%<>C?N@nua9H-WTal-9bTLIc{WHxTx{3!u_`ABnPGFH}s6+6kt>hW^*KT6D{6TrXeo@5P%VcnXvn( zz3GfC!uvHDLyd-{RG0(mn*`M=R*?I{j62$Roh>=LKjWN9D{+98xZq5F0=ZDN!HTKobZO9o;?{~Ppxe3}qAl7SGf%yxibZD0Wc{jI~`zr>#}MUa4A>8-ak z3FsMuqNMFCqy2BtkZ%C^x(OD+!!_vu&X-~Es$Goy=!P0cYn%ZB-4F$yVc_WrWO<4Q zBFoFm1<&7Jha}vn+%dB`@c%T-p^1X82_^3ZvD%z6vO+HlT`m&BV(=PfwPilhLha{ zN2>>H(hsW6wcFcwje+}Nqw)yWlBDoEFFRC}rweNE9>mQC1Qq_I2<_W5Fe*HP5dWnW zfmUQK!!%5QCv_PHf-_|*9FI0d-#laB-riNoqsRtE1bulaCBTylJP*JwuIPNRr0lL0IpphdWB6dcFGEsQkiuIhE6-b{9XMp%#emt()v(EPz zbdR8242vui;AS@;WYb56VdG+nKK66w+(>R1Rr0BfXRHqv?{sSgDK80z?QgAKII1El zI9fxYxA&=#8yy2n13j;frUT~P>x|?c%pM;)jJh^$6w%;0L*QX!Ksp1+3&ZgjEk|o- z9SCfP3rE93pi>@FgU$dI7zn@ibPmT}9px`Sz<7RO0O&+q1OkT$RG4r{<(7Cm_ z3Z2&a_Y9`D(P0%x9b2;=ZabekuOxpWtEj0aGg)tT-Aiq591f=oBX;1P?a%O$5gdGMP)*!bP;Zx@1h+dvL+8Lg5074aX6i&On9 z7)_%x@hpFvh&ll>RY4nMnvBO<5 zL22pj@&OBHCnqNt7Z*1-`M@2IE`^>=dh2r>_WcF-4ttw}sVEffDg5bLWtDmkX4Mo- z>H>E=(qs#2vOphm@7_*0n+2lEt#X)dbsSog_4dr%!4~+p+m~X_g7haIhX)+S*n@r7 zMeo&#Zx0_nJd)nn-X8ez!>IspJdBKtzl-Qu_5bS8T=NWqG+?B>&dJLgyfx_XBZDJo z^&8uA{qX|}i}Zh31x3W~|InQ5_0#I>70i`ca*AnEwx8n|u3P@N-cgY5|2i|%>RS{` zadGi*g)8CsXg{3c*|TS${t^-vE*oCdD(L;~VQX(+QCXSv^yvyeI&8#10|O=VQ|z2{ zYMjczgqIy8S!Jf2oJJlC={7PlW#(iGiIFNRJB@~rBrMA)n$x0{> zria8^`H=F@_#f8|F@Fvs2NO(g z^pZ6yAwQ)i=abVuQVbo$Xi~B>hrsn$GxATpc}mR=U~=6m&!C|Tfleml;v_qxi2Sr- ze2#?#-i?hwniG(Qxq1CWdzfM5)0};gLzWXL^HO2+|xT zu*n7h#?5@Bp492GKY@(qWr5#co!OC7nY2Q<1WY)Rk;z-tT^f1@vT0FK+oEqlUY>xU z;Ji7FeWU<8V?Mu4fR|W_D^uu!hzbe{x^W|hflplAGca(3)XOuM zu=ILEp&(U)p8tY~h-+ZrM-Wx{`O~XSyz%kz^!^2-ZbtL~a6TI&Lv8X5kO%+d@AE3U zq81<7+1R!}QS!PBzE|PKoIpOB`Xs@fNdt=YU~!J}!m9}D4}%phgKK;9rD`fFwY$?P z{z(t#03o!gozv(ljDREYphqt#m|3hG@7H3xU~%QrrB^m`$*PLd=AiAjaA2MHB~Nb0 zhe*iFA1!?7+}hYcpz3xeqPNF`6)s%RUe0E?_IF1YlH{e*If-6I(myy|(-2CaM(9P~TgFk&4RCPeCEoV-D^x6?!q80X!N zY6T&Y4Z3Qzmv@4H-k(v^O=2axGt$Az1SvgcWGcywjqQ_QQUWtab#E3}9K==XBc7b* z`Z3?m7Yl{4h5&(dLHRtdoX}yQ{@nGkOign?rF3P(@FIwxWPUhPeoo2Lp^cWL%6K{| zl=you<4UfE(TTijOug7a$H{(?xS{N_Y#+xrb6I$TS4A5KWvtYQxk zs$hSrFYhb|)vp%sRu->WnAGlrQIMU3+Q_YL#{gLl@;v)9){pjLyFtHZ!l^QpFu|*k z>Fflk1V$#&>I1KV-O5AESas|0t@0}KL&x9me?*@qYzUvkonXK+00w~mZ~#})RKze* zaO`C_x1F3u&Ttc+kf#~P*N;7s$p&l5j^u;hqwvB7bBVtv__?n!LW2xgRwjw42W8?R zc5F~HPVqfF8-Y7OpcDC{Aodd>wkBkzF*K}2piF2915#A4Azx7I!W*i< z%;eK5fE4+`mGnNA2g*U>X-ldCHD-a@c|vU2!7ML+JT2}l%<2R@p%1JFia_=Aj?64UvLQiy z5U3>4+(1vs7B1}Q=qL#4ckkY{C0)okD9k=9Cnq=L1oDu#w>NuH$|)dAC{Y3r*!p;G z*^1lWqIwIo$Gv+sfBg7C!z#Q|pH){^mz}LdAaUqvX~k-JSy!3Ka4{i$gM+&rvda?E z(oJ8#B2fLMmd9fO%Pg#{d*hS;;u3sbi`fF@bXwT4Aj<<_GU&x3o+Ko!?apLSj`VhS zudaNd7m8Xdb`KH~7Z+cUjDU3E(9u76eh%bDg_QlLyKy5TUJgjm3I{JW-37B_4WlNMW_;z@)MnV$o0-ZY`v^vysNlCmUUpoSj zWTc&u!}T2DvDv)EH$a8n^X@!^VSqVRkXEIFIZczED45e3z<;hj>B!fv>nnMqUR&;O zLjlBJR7INeW+R|!M2L4Wn3=2(@A;t_mp%3K6){$0Ah(g^#mTfG5^3_^Ot%Wglm#iv zgn$#aGs4bY0-yUFqey)Xc$g?6w=Hz(oW;uwy)yBuO|PK_J_H&C1mbZ&$R7h^svAq= zN&k&vAwTzT&sS9Ir$z%Qxptow*V4S=C4N|V1nsmZlu}v-u02nE#P(N(jGJn&2;arp zEdJ8arZ77QSehJ{|6n335Cs`Cm?l5^0jRsb@I`*j?4!e&4C-s;(OQpkbR*IINEB?b zN8AueS9kARY(ndUV7~Je>cp z(`UP)ME6)GrPQ?G;pj?c1j>z7o zS#52Iu0Z!7pecvf74E2O-f=bnLxEVpaPIB-H}!Lq@Klr<^6<%j`JM)Oyu%EARy{Jx zE(Et4P$S@Cpj^!bXFm! z#LmdV;=#htcOX5d$_MS`r?|Z?DY1BK6*8tm2HAJu|91o;nzKPx)3nu9#D)e%cX#)j zH=q00fIgwDtc+XM7e3O6;Dip;>#AO($XhHy#FEor(J?hmdGaL7PXxlfxagt(km2*> zDl-d<>v&z71Psc#=z%P+tW?!$f1W1mw=`UNK|#TC z*#v@6R?W-!5x4E_`s2gH{cQ&ahaw|Dtg$V2C^D$MsR0m1&_zH1Kzf)dKkUG(>grv>_+iPIp&80iJt})+W5cn4 zbU)nP{kvGGYH4W+!1T|bKmW8P&6~Fg^gcz;*enop-(najG=n{(Pn9qsJVZq0>Jl2^%X7Y z5BBC#Q&TMuCfMlZcYft-rhK$09b+-iLsRkEsHmt|Sd1-=5(Y}l_W>Re%{!_-rbI&! zIg#{TfUjSFTu)fxxYLm;<27LJcevXb+)u!_<)2so+;rAj{pt%0zV2YH4bZgXUY&Al z@SY8U!-E-x<3;yiuA|raV2gLb$AWhMBoXJpQLO+wXsspK2E-I-XY^Cc`z@KNsy6Pf z&j?&luW)$9!Ob1C@`-Zb%`JXyq^ABq8x;TThXt(;NQ0eM)?>!P&7^|V0Ram-(5UR4 zz&Z?@6u2Nxen?6GI=YYsrKhJSe|+H8qLHvDmY<*BKQK^qLQPdQE< zR?s^z;0HDZF4$+yo9z-c61K$RVq>{wJf}ddeDI_6L7tyT(hb9evd)#Sk5AFdZFGJ8 z`nB`H)sL0hpUrKEGl=a&1#lsF9xYI$*WRnyrd4v{cv=2R8w<#5tJbY2`amQs(5fP$-6vC zsI00oUN#YZkxz0hZSBrs_xj>A&`8|#@X%@U|J#?SQ7%%$MwUT_`lKJPt%bqKrg%Ve zH=F%^-ub!PslUyh5M|MZit53K^Og#LFswgZbO)rNoY!*gBHfK3t+XN_AYB3if+9H} zF(Ccz(f7T+?|kQX{ytupFzngSu4k=#-RoZag_ef$RbpCV2n2Ff1)*>s0>R~hKyV5O zaY0K?RQNCOi^v6G=njF9e#HLANxnmR6EqTdsHiIvEMB=rCUX5vlXo6yV)jsc;GyGe z?cx2%%?hF}tE(m}cNeB>YXwt(o_yt=g)_+ktwPHmlD4YW%@{+H@UCGv~~IZ z>{@YO*gG=5LHhhM1cKyWGXd+DgCZah73>d=|AimaOmVcdvr znVE;JWw;Q9;gWiXsdbZ1SeDPvRI~Hl6q})}vx%yJiIEXri*6FQf}P!3U-F%yEXk8m zRc=oY55Engh0$g;NStm?*X&OA#gN(7@lvpKz{Dm^Vsf#!qJQ&kI?@f3P8q55-`^D@ zDO>jj#qI3u%=|GE(x;u$PZrwz*T+jmMMb$%`kfk2{MtTz_|Q0yydx=Tvi%rGM(o^e zYHh{9e(05sqsrw$3Y&OS7L4)%M{`(M*sOEABqygjAsGWhc4nrDzP{JDj-VeuESu-J zLf$;MdvtW9r>BQBAcl;t?b2<=E6Ih+luwy$t?$x}7Z%`rTr=r_Eq(v4nI%!$wtsTQ znpm2jZ_$^;vEKXfUDNN+KOSn`Dxy$1SJxW0JCGgyB3*EcrioV=Z!DhESJgSLVX=}T1BA zg`T|KsBA`}w8DstA;a+vE6vLly+}-G54yp#{Z=N+l-crCC?X|V1~*V>iqi9GLDQon zAvALqGa)l1-0hx`iDGXB!_z-UMSCOqR+Zn=dxY^$Dm-6*Apl501*2-h?NQi+ z#P3G`nbG{|K${b~#`BoceHyg$NH}uKH?hbrJm$IyhU_jfQJ4sCCAkfziW+cR&cf$4;Ko#MS=aJZfm1gEdP+N z1Uq^bI4ebU?%)V8O7I9*9_=m7SPHO9>ZlPaZS*zJUy1zkKeI5ssL^-k1QV0wo3D^9 zcMSj2ZlOh`Ew~!?COBfOIS_8=L4_AKvl={_iDUnT+i;K(aov~N*+{)e!O3J}> zVeY?JXa2yJN}IEdpO|ff0-qWl+PN4Kf(QS&I+Cx>IyPNvbB#_UE+gYf<7_CDGLl7F zy1lII-ZmeO`CD7@jl8N9S(X#RUQ50kR8-tX4NnF$L;;$fpKP};zciZn+Z*7}DJq#_ zW@Vk6n7|ylVs`F3I-(ItS0M34bI=1mN^QYm9k1z2Duk}9!?^(I3yO+Z+1P}IguJ}G z@c%1(sgmOGG~M|M)jy+6l1>GfG~154h-Xsv+MaK<9mz9k_9B9t0+jjv`?onyiRpFX z(4{4kxC)3%X`Y?@7ujBV1R9M7=RhxLpS0_?^znsh)s%hx7WT^~ECd%}{`G4PeGE--_ z`K1~gBjPE#fR(J2xrkJFU+wH#bQvCpbmUato~fcs;!dZJmr20F1cx{zTz1`ZZbqE{ zvE?6UiHr7<|59SCP%jOP@e|AXKgLM@;{9LOB7d=T^$`Kq%J1g+g6#uWg@u+ZK)tXC zfRHs&5S&EP=z`Cy+(nds4OuL*UQzc%9y~Jm_+70=B(Qi%lf|$%BLIW`pHFjlr`@Ce zPI3$k`v&jDJ^}z*pwzgcu7%~_sd9oJnEqYbs*EG1=#RT~AR-6;E~z6tV*0;*`ebM@ zslY9db#-d^-8?ZK;73Z~USt#0U;i#T3Xf9zkApG=De7W(*73$d{!1>vP4lNC{X7a_ zFC^T+hZ(_xVewvMdPo{9GIsDhJj7VuY9$F~1U=uD-26}zzODtLFa)BIE)zLrwn#2# zd{SIQLzK!XJm2@rmK_T=%&2!8`+IkizxRYZ-Q<5nI*rTeJD{K57TOz&eHX8|hGYe! zb%XS`Uo2g28n?I$CRd1ReA?#r=6ll#agXyIx_0!QXn!aabXmJ48B8q6V?_XoU*z`Y zAkZw5IJRCw%Y_5X6xRGK>AsFA?wMN_EphO~?;s?CFft9V>yb(!K2exm{QC5yK6`qt*zZFEzX*G^GM#!Z8zu=K|*38+N&A*;c@kkl@$X$ zy)_e9;H#nv!$YOVZMzVT;dS)>+hr9d`8fYR4Gds~>ZyFo2eZypavmNYYN>pR{&#qJ zt1BxjhcruLQd;5a%&U3f2B@lXq-az&xlSasMb*7t&cKm*BgoOr4a2BSpZ(a5Q(-i}*- z5q~Vmj^d-S=-Y}#o#qc%rwciz?ea-UNl}xOZHLL_z7?@xz7x1wTUntpXZ}oZqcVgS zx3s3lZR_j1b=C5ThdMg>ZBLw?hX4@eb3v*r^QwdA%$d5VgZ8n+?xm5Fw%N}gKM3Ks z@7xK6VlWt|MseYk;!udk8_OZ>v*V>Y6a1`+c6&bu z%JpB@%%TAnjTNb5$zin|q`DMDD68gt5bd9Mdn}Vc5rBt!JAe3`%D;fR@HbhtkW0uZHL>{j>aS%D82;~uv> z#|FieK&kCT03{t-rT+7^Y$ikb`Qc2!{}5QVAm!}2C#e7h-Gadjq1rr+bW_VdMVTJkkTc}_1KV)jf4&Kf@{Msu%6Zh588o?&Ic5tddNRELyz0u!O;M%nMpH1N zNL60-_vCAuM8_?jfHws~_kC0Y0qnt7_g+H4h)dwRk(JltxI{7)V9n*7Orr>8Jo3N? z#YNh^AJkYkJQQF!;dV2MG?Xg^YbF|MV6z{&(0L^LJ{vJIxX@-WRT&o1E1g-O;JH_N zZ}t$iRYru>B5cvh%6wpv9~tY1#OD&UYIt} zg%%P4Wr2m+p0xOQFy90VA{jaGDK;YC=fsx;@8a=g|M>*$>>&}|q{l5RRK?yn@=JS1 z=O^tKCut6qsiN*h$k7eV!jF%U zjkB5Q>E-UZ5KFo2MzUV1F)N{T@69heiP@)?Q5KBq=s zi7fSWq3B!U`rG$#K0T9i4sQCJ2D`Y5&;mb-e^#0c=wfum2#8Jb9Ey@rdE)sX5V&w_p0WSISEyje#*3X^yku*>+B#?)^3qWVRU5y zlXT}aOmGgW7-b}~hWtutP7eMH~85)p_*`xAD%*XFAaq6c+_RFPuzfc?nA3J8 zKL>p?umImgT#(E21K6}itdV26PPfuCGbWnrOw15Ikj!yUyX%T>Iiqz8LuJ92W1ljR zx%D{s9|(y8qX5QCZf;Npzs|I8m|?A$-Qe~aV=#}bjGPRxrntev%FHy~ zZ&;i|}_tp$fezUGeB{YVJ*boWbwL^puSg5i~!(r^y2UXFcJ@f+1)GPpVI zQd1g=6#%JMpY&!~)J@2OW7@B;tXRWhCw!EXsBWSQy z0Yw~EtgUMP7P@CQz5oUm1Owk9aRrl7;1YEpn#b#P$e2jp%rp+Iocu`C98UpUUE|(8 za=E9i^FSs_ANc7avR_WxQTPb0yWax|+W(pE>@?RTC?#d;a|A?*Ij>pA!otFv!cL$2 z`^Ei!|7!8s)311t>l_aWdseT*3%P`s(DU7}6^Y4Hj2!2GZ~o@(+b_#*(`c{kA%h1G zM4aba_V+!`?D2R+CN;V`NPdy5fgC82RUNRdAQ>4^_y3^sG))qkLn|^PWP7FLi1CB7 z1CG~@F&lse0Zj5goNIj485R>0^FQQ3ISlJFrP$awQ_9pUqd9GLt$^QV5JcdWUbgI( zo$i}M8KU1k#LQ4=;{fQkd2gFymO*-_98GIg|IoqSemoK6N7K{ONL>PR zC^ii@6xSQpfBqhz+1%V*S63In8@X_>!WMv2-V6TT8X7UsyqcQ(gxY3i+g{D??(Uvk z6f)c(iC@oWnpx*zpSEeDBkftAD9_Ybsk0lyDsepwsMCg7r|bViI`oYZtWeccPlRr= zIc#ysMa0$~;*ndCSyrLncrZw!U!E9rYWAW=a~g+-5xgWT5H=O4)#u+zfDnmzjMg<4 z0+{vjjo0l;Z2aVVAb!Ij97JE78CK^8=ALC>w+ew?Xx*807aw4m1`!Qn3ZUb+2=D}B zs0^;l5qN?OMDn+4U>%=8yfr4B)KL5goPo=*wo3dvmDBPoR%Zob@XV;~SUC6ySj;!_ zeFWGqR2D2Ex6F1E7R}=aMwibsDT@*K-~!f2qLHoyBohP^2>~#57_tD%8T>b7>%a^o zAp#%5UyLeAiOTm=Oxuy=-I)%?HUh)xdDuS(FgB(H=H#XJR>+Cfbn(*NAKfACtfby?R%R0aP}v}IIH}_;+uIv z5~}B&+D1m@-ZuhUxCEokZFf;DyCg*6iUVz1K>I1R8?O`7+>U@0Z3cdkr1d&A>GjGxZ_P|zuO!4N#nA=fAzt^p z4;x9C`F5Kv-gBiiSvS2_K=#9?zy`tKzpw4!P<{+3+16~M<}R1gw|lk#*6nqsey&w3 z&gTHej9~gS#1k|=JpAih@U?*b5#_%A(b3V+D^xF&lQ~VAJWdYK$CLOZiDTmWw+Z@U zHTg|T2EQ1-&HzUj_WRn1le`A4taUeGfecJbS69Gm^8q%gxy!{BfSEAJ27-I5@8Qd& zBmrUJjfXIvyGfxw3D#y`U|*lqV-I#)OX%ZR(FicRQ9?d`C79z}6DCG{>y?C0&F*eT z(52qqUaw{o!bjB)cK7zm%gQ`Bq2UsO&d7?se2uSHo31Lv3fwk)MYBd5gy56W=>)_} z%f#fY`=+zR@oIiw>Oc5M^{+-`2q-DtuM(oL{7no@pm)5pqt&Z`jADH%=hX5VbNfQO z^*|bxnV+^VA2msu8dbL!Ya#~)G0>LW2)nHsW`KowBEM$RdxJ2}r!%!I`t>r`2WU8*Q(Z(|;bh_BJVVqyiKwZsug|d7dSYs-d2S1{F|}d9#1oZ3oB{j|h}eJuplde9 znre#bw*@$Q&K%`Q5vm6=>wtC6fY?{2R}PqW&-dv%yTlnXQibdwc|^VCP_==SL?KZLp!yAavli)8i_er8B8V>I*@i74O%g--HZ2qtPeCgzR~1e3ybTq2&^g3%zl z7m($FY=}BRC>ppaOEmN{A#=x5(1Rt77Q$rrhn=XY=NC=*APw*oNWkMf;7eifIQ}T1 zHc(DFAZh+{m$^?1r>8i-gwc%vjmfgugfCe@jsm1CK6b@xcy=MSfJi2cZXMtL@|BPm zBw>~XnMJCN%t_ptSxzd$jm+k6$q<3YNVT!ot^E1FXK9XK;zg0*62*mkrN!t`Lmhx4 zQhJQR0fh-Xe^qn|0Q~9}peLNu?rNgeX{6OrGPwH%bB|)}#6AM5GW^_25cXc+79h7x zt2_KmMQad1RhXhC0t3&Xh>lNd)wIwtk3f)fML~?AY26!S_s0v?LrYCO3^}H z04s&@{Z1uw2_>M06`t?aJgdd3{3#P+TDX#R?IQT_-Cx2M{G%r zKx(#|D!OAPmI>Cv2t!(lKtEQpp5|i0N2vBgAeW&7{Teg_< z?TfwgYvg>F&#rX*ij$0ERZqiGM<6p_fq%Ak=bXon@7^|3 zyTFZ7{*$Sa43EP?5bEmwtGUljk!yt^At7L3z>oin2*sSTboX8wnnDhY)#uK%eA+V| zdun~gR6xd40hDh|siPnDQ>Hp*^DdoNi#4+w9=;Y76kP8;I5^ncd?HyUPqA)Vx6_rb z?-v*6!^efOG+rklN+qBV@IBYrz*Jrd52(3z_u2bz+^D1emX?!WFEnmMaWv|CH?xk1>8`-|m5$lG)j>@fq*5kWXaa-SFl-Iz9oK>eRGWlzp~o zb2#U126AVhQ>Jj4JbdtgJ><)O#kVn_hhi!z@QYsL?E#dEypW1MziiHI0;=Db%13WS z=t$FOd3f>&tL3pNuQG5`cm!k`iU=1~AcL{t67hzx<60`g;H*%^z*9nI1!;6!i$zEa zrij22S=jrMxa60;;HH8N_AUSjN)S6gwEf>e96B9fuh&u>jaj+|SgyF7uns(T6| z3F@7yPpa-WHCM`Lm#g;(^#r_Zv&lKYEY@U*YL?+8$A>~XzC{D|U1nQpZ=9;8?Qd?O ze6xe$icKt!4hlAG)bE*QBno}Cg2Y`9qA>cG-2s^wWUz!}wwNq-)3_ZujLYdw0&0>M zRfhxVR_y`D#oL6ga(-}J&UZ}EXb7_ctP2VjqLKqd?CX>08^AzYY_d>!W8X>Zkf3a< zY_qJiz96`2&VF``XQl3hRp{F*B+f9+fX%K}e7#9OOIdmkLVY=iLOw}=T#nt2IfUnL z#;&#mIDr8*LY|JqF8i9N4{26FK;c1+Ou4Yu;tjbt6&aAxKK}G>^12+M`?})w$A4%aoBP==EWx?+;~WZT-CXbz5|gLa zT+I%`uLjC0+Q0xj5W)k@&(_}+H&QXrgII}Z5AXnQXtLmkNB`=p=#m0lo8lgHl?xCe z){RpitLdAfuEjC(QLiV8v5HS2<;10Q@(?tQIN`TOC=n)E!RrR-*bt~=axChgJN;JU zh0(5mo_0z9=*K4tg(8uI^;ufaS&8e{H15;{=SoHQ_BmWP3V~)PUnoU)g8B$cM1nQH z8NtVwe|fe-lu&e?Mt7l9>TV}MV}jpL)D)9O#bMZYq*9x?zH0d!U;ZpPS0FR8s=2lP zk)6@uEMY_LB^MF8`ucjL0dHIhaREqelB#M!aEhp1J+KpQ#h?NT{Pwb;CLX=vjpM%r3 z(}Cx$L(kYR=C+W&8RR;{2oP;t`?Ia8@v#HO=MPSmN`CCiZfpNqc)IV#j#|^VSlea1 zeRdo|FaD`CScucTOx@0{`QydMtimj5+iglYb{zMeR42HCyKl*@QZ%I|GX!MKSJ(8+xSFl>9C?6y6z#S$`SjjbSxmw-dUH{qRV?`RA;%F}GWQMW zR(I+4Q?f7Atw6NmGh^f0d7fS7p(+PT?6&&3^hgD?^jecX$mteCY;2xaYD&k$Dx10C zJDtAmB2W-*`zPlrEB8MYsmAsVF&Ilf^iK`&e+GS==SGKay$Ok8K^AR>)TpgYUC9R) zQ+}RYsWs z%5M?^z0OIr_Htx}dqdEX5b}(p8DVE@YwP6H@Zj!Pv8I70^8YGEN&MjYj?BhU*iVgE z-x%_-cV|0Hm`dcIymEuu)EBa`mf+;8Q^7;Xc{ z`PrxMc^e}#se;cd}&9aKzZ84WKM9(Sm*wIas~#ohYy=vzjuRx*ZKU~P`a>l zSHz9JE+D!UOb-Hh;2L@l#H^x)V@{KY4In^$a2tDpmG@3iCMSw3o%urra}S>?FP?iJ z_Bbv~PqT(V)3YAbC7bd=ey42)d4N*!!s$<8pb?t_bsV|`*cKNjCs2%iBO@bCxKh9^ zl_S5<&`==Ft7>W*NIeEJWagDEHxYZ(lsnF|_=3t?-hc}N$hmm+^mYqU-=^GvkeBsN z?4UJdIN^+D6qLKu($aQ!2dhqwhAvKryqf(%g?eMEYM~v~j&uX1?`u-O&0~kkbbQ1f zT1uV-yP9D{rrjp@5okew(`%4UZ8D< z{d>77c`?@1F#o;T;{E0!oeagqOvPgZd`LpPWSlOzd1Kh}bf?CupO2qkl0LTA^8Ny# zTe(_@o*8`ehMCX$M)!ks4Bw~Oi9uFl?S82;|JvzcSn%C$y(YOo4DLjPLU(p{-kWRr z3sVqxb@gUd+~Zp_xeWQ#4+;9=)yxp^dlB3+<1(yK5ViO#d0m;Aovp2{^(-TWKxt|S z5~eo%H>-QHfQ(=sNFR^Hj`LrW_|=gn z=pY1iMEi-`AYt2iD&#nW&9<>+Rb$cV*)|*quN=0P|Id}5@IjefWGhJC6BnOXM}V7U z{wD`#%iO4u=D0{dkmZAmA;D(CjDR4XpPiJKmlL#mD1an=d$w^K5FR8IV+a)3MAkv! zR-<<;((-MqZ&JbV4(Yab?4O2$rdDq97_s(~xouF<2{`@EKt)0EWVS(l+wB8*x7_UP zzjtC$enXOnr{4&qSGTq%oAUtLR#sJA{1k3KDwaM|*_QjT=<+r%k5O2-v82TI{jxcA2E#>uIv=beiU4Gl(4?+QyxReF#y^peK|&H=TpttLK4t=sqTPv}pr=@!^dEVbjgs9h^mS%EH-}_r*s0xoQs*@J&$e1jU~Zi@$8L{e?W& z$McnA*YB(N`ksKQw-1=GY=QN7e!g{MQ+COT+~x zDmv}dcm;l!QZ*nTKu<3@2E)w8wmnm?YFK7W7#|TEl|3qY{_~9-k5bElG(pxx zl_Q+DatmN)$)<(*R#oz*`9*Lvhu~;pB^tf6^BhAcNbr}9z4e+!Dq^gDpPdUjOg@=w z((9o`0-yzzG;m!Gs5h~us51;!j9HL<4P62k^DMffRMgegH8h-_JlWjb94{n6XI2dV z8TVclcB9AtPPB(gWya3L8TZ53yS``XjTrKHt+-et6M6c32i&pA9$MLX>aneO`k!lW zvOUOy?FH12Pw~J9xgXZqf_j)-IMyb6kacfc30P)t4Wj*3_S@5$$`6~$7mfzKfh`d( zdNmtVJlKgeJ6#$S-yUxyhb!3IuaD#_1GR=TM8s#;(#p!}{{4hDahjrn0$VfEUsJ#0 zN9O0He0G0m=PPk=vuR|kVvrb4l`ga-`}CHJ)3$29w}D~^sDT3Ge+@i^j&y2bVh-GV z0^Tb&xwlrFJr6Epog8jUpVWhUi)#;JC$FJ?%hNP%4rWP;yZ`)rc5?XF6+pcc)cS|S zCf)Z}N79|!-c?ug(vhMLrfsmM4sy+M3&|TdZp=_}%sD$dN+x`rpGm0OZ|bbad~)%1 z*oa3TY{6%gKp|hH%+uh+h)%yJ!P4SCR1ANC5`$+28$bUZaPiWywa zR9g(j*46hWY2zWwi!fu5TUHu%WmwFJ;oBrF4>pt*_DulvWBeL z8T+2)dq(f~@AG*lFaqjCr_qorxmgn`nu3%MVxpQRnWDo?MQ;?Tdhaeaq z1QBJxVPGURENB({B7G|VzyX5Jepa3M9Hm?D)PJ3i62v$t~KLMqrA^N8t;&;t1@1bHcY)5ne$ zcD6`O3ma32Ux-`T&Y4?)k6#4IFCZ!;EGo#?W)>g?K}bkJ`i`b+B7VeMQ*)qhXHDOc z%<>*3Bf^m$VS2qcGnk&4Bbb)u<;zS?>DJFPA#4>b%=~Plg=|+N=E#KKXieD{6l4Y$ zEAdmr^=C;YXJYfv!LEI;V|3orPRXsFVVyuf>5F`vI2-XcOG9<4^;@E+*=SK;pU~{L zXS;=wOA&vTJhzX+hVj%8G{-Hy0=5U?hlWA&Cm+6qZ~r$w9%ql9lJ9UU&dHn3M;K+h zC31XCcew>2zQ8l|EV+kv0VSrk{f- zd$Lq^_5)Ko;=QYK(}xwJ2L)`6`Qu5T`bWL!>zthKYt1wo>gx0i3{p~3-bWiLf&Q*G zVv36l@!w)s9|%VecB|=T%tW_Y2Oqbo*&QMwpB7x@(nFqsP~Q62DJkAv;yZAbna(8l z@kZY>%pTwMQpVKOQAa}VU{UXncdOd(U+Royn{8$jRe;8w0uoKNwX}qtmvn?)`Xd^S z4_E45%VOZ>_o5L-I_qcgVH{rtzE&uE7AOM^bd6)N8$4XlZV>g@U7eNKsTw<&5s?Zw z`=H#iPeVn#Z$lujPt=DZVt=W8pdij4;f_yxK_i)$lp@|;<(w!5kMtH^dX!+1T2Aa^ zYiVh@$&;pSOFUHT=Af?LJ#5#2cdj|!u53_B5-qQ-t(`Ru4h|;%bC9_e(Jw}?v9l?j zv>j2ClhhIY<;qewBmfN_kun6tAMz@EkDL)iK+XnG;bnAtR%*@2u- zJ=WT};o0#S;vpzo#XsB2T`(y7?5&)~{rL|<(qm#`;4BwqICyzeMC_7jB*!kA1{;(s zB{Zkfo9c;t5sbnwb8l@bS77};uBxf5nB>WM1tj(&{ash;cIy9(tjNmA4fOX*>`g^; zY8NiKz3j7OoZ|K0Fw$6P86#t1(`BzOpYi?ZaS@ucT2)FfqFawmLKrRd=2%j+Whtt~ zPPB2nGauqP68|1-*75vd{jU=GknYJ#_ECUeE+Vx z^*%KfDI`>GNVBV~i;

7j~jmU{r@eXOlC0HSgblm6TNX|qzrEWALN zmfMhkX+krM<@@^jNJ&XU(r+ep`s$D0)`i#Yjh>P`egeTCr97({!47;+YM*{1=C;1O z_WP2Ig`Hj3w{PD%I`9LhuIKgT))nUDtQoy1b@oujW-7raMzIYvlGIRbrwe!5K-*8L zl41P({CtfwJ~uo28V3i^q_G+oy9S=`y}!ms$jX(e9Pp5Qr(8!W?-yG6;^N{{r%r)S zOwioiTv%Ax1$O1lMVMjx8|t_9;}&3tsOgAv!1WIwbZ###EnOBck#FExE&uRkQH7O1 zj076JVvM|XIXp2j5x6H7{}EeB*VXY_w+(;9t5>gRcplD8!@R5Ovy&!QH@Z&wpsyMu z9~AU2Em)? zviSd8MSFA({Px_`fEx8YhEUiSVd5yMN3^6+Z77=lrVVXyBnqMxgISE1WWH206fC2Z zhxl|zx*K!EQ7oU6fByn*;o}wDX&xQ{MMeenI7)C>r+o#lJ!1_d4Tyvq8Q>Xr9usDJ zfu_Gm7|{F202i2E$RB#d2haFuG=hXy$zT@mCuI!`rJnpg!g_KNK~k)^)0{j4OZ)T~ z8|}@)Q_$Hk^g2D52i$|=B`ouF1bQ7t7+C$s02^VT5sq5<9uWD}@<9m^blMNSejYrI zi3&P<)_02WEGtnXl2YCvIu^kKRM*G%3a-{{#(z;yjSp-;F&gQ1t_i#;t<)CJBcF>h zIeB?K=Pz8i(A3oAlrQ&-_TTz@)<=!lVyxcN&CM+>Ep6fPojZ3dEpzLylCpGDUAfy; zmeVA^fTMy^Bf{wA@87>aI5-$e&D~vt^X7N>B>ybc01k{@7D_y3H}KJzpztKRq+Eh2p;3dDiLOiCLmX{F9Enpz4=vhW}-N*@z$!j$UUW zNUXq7|1mIjA`ygM=O#!*z)}A(kYY@j11Qlgf&?#E$NwJiCM?X>3Sf<2vZkh{U_S?~ zmz@pNVyde}@qGOJ{9knVv=UW&9wG+A){_@XakJjir7a0i2=;kgD6tvX5OIcB|cE08xj?zHlGxoWi zTgsl31ZR*FhZ>_a5Hyo_DJIr?0ZFJr?#LmCd z^HlCQCy`giY+uh*RR8lsgTdHch=d~N&BxI;T6soB#>Ii6uJ-ntzEjRM;^N|-1BLpR zE?v5M_3HHWwENETmiddr>Dr&Yx@Au5!2>COmY&^V2iqsx!k?jS=qbd(LGX^n#cFCX zxw+Rvn+d`3 zcshsFmlOGDnk=!zcg)<2^^KWpdI8|Gu=;;R}w%i?>u| zNzzOC9&6vv!dN+;S8_?mKFT1nL=^{aE|qp!7RB)z)<62_x4G%aA`Q-0Q&Uq>F|29g zOIzC>0bD))L{c;2wL1QhTd!hHdxcxb&d|Lo7X#-PyJZ?&(Ypri=P%uky`M=jK01iSpft!=lsx$Fcdo1@oC8fTa>VyyKPST$I zw#IAkr|2@TwO<>qR;zHr1QqLZtIbNfZtR}YCnte||AN-%MLGpqp6s!S$;_PoaFx=R z50Z+K91W~<4R7hgix0`;8yFgbz>6J!uT!YyX5GS_%>}S{Hm%=A^mgyB)SEyHd_vIc zX#lcd;Iqw{;Zfe5R^MlCUM5dK`TunFF$b%L!`iZXC$5;^)7Xd@5VkWMuJy#v?-dUT z2HJTWzafAcF3&$PoiAWMB*56{W5uL8P%P4N7LolT0`c6xg{jSzLe(hlb-j<~f)Mf& z!e1+Ryjp*(3MOEm2AOz;i-eq+kewha?`JFq7PqGy5 z(bCb-^{3wKebM;)O4ETM6UZr;qLKkXvS$jS!3bVEZzcl4?Mb3K0bnZOJ_W&H3N8zy zo)|o;nER9uQWiNMQTGQzBjlqVtw=yfTm~){0W=VeW5t=IIYq_R>lXkFmP!s&pw?y2 zx%2G10u12~M+Ik-i@-${AU`T!R3tDrU^e~;>w82wgrJzkg&iM)Cr4&JFQ66wG%-pH$0H{ zwX`bA%I5T6><~mB884x83kp0ZgV{m0vBmS5{bvmYT9yAP@>Wu((_;TiCUM;rjK>SY zY?p}7s2doJE)AA~fCWSY38JffmPlaJOLIToTJ`789||mc_qRpne3~(Cp zj~_pZy$<(NjyAI@>q-c5XKOd6np+0ou=41bSVS6dCSlT6r>DD6QrkEK&#k_Uj0|8# zB_#qfYPq?&Fr+K?$^lpmFn=t~wA9kb86@q$60B4% z!+MX8ey0Gk5+5tkEgDhAZ$%E~D`Rre9||5@@(|-h1q=+tMrzSMq`%K|xJ%u4~uim6W#J;|0`rmN~_H zG7CSv`FnQX*#%B_%Q1UywlLTM93UqxT~=21E%BBEIDyE>$m9V4u}#g*c|lMvy7KPb z@<`Tb$?%5;37ctstTqYcdF3~EjlaLY=bvpEX=yl%B3pF(Gzdl`g*k3xp`M|oms~lTRWv>F}H)=wUG3(%1U!HvzeKh?E{j3 zqi?rn0S@sV39HsWn4mKM$g7S~K1j3FSk)G$sJm)0dD=y?GBS(+&ek0o9>v>L4)^r* zxZsEGtEs7}s{@U%cUc)t@RIL1=f&~h3fM8UX{wp(C|V$)Bb0{M3tv9)>C>m}rNQ5O zTK^64h`!?B1xNpgOW?G+zW#`wUEMpzaJn0rRaJxkJ60U!bs>S__!I-A4{+EkZ36%g zFJL^@+t^p?|JcgP1|ZT)O4dKUAjiN1nFSd7*nN-qWLa)zZa!>X0pivzhZ!!t zim*VOZ3s{{+MwFfrPE!O<@83EyRzsn${l2qu}D8OAq2nCY3TdQ1Eq8T=3Oa)2bH~p zJS9PRJvfQ_I6<=mfP$U^14v02$2Ec;?orBr^aHDs4@Z+x6L!y@L>-~FvwR`{0~z}6 zzrF*<8lH04VB8t11pZ-ZPC)~0AB?#?%-hOdOETWt8!zVG5lqfKE=6Gcvj*$=tR_Ij z!p5CO{`?_D#sl*6{Ns%%_shqLPX^Rqn1Xe8U@l z(OmY>flpXlSdP}m4!b~`GRRYrD)){e_9wqhrW~3}&v?t4NHCS`uFU88BRXRK057m6kmWh<3|M+$Mvc4Rn6$>q%kzG+>>rhxa^@*YwU^6%xo-uu)jjtn}I zMoy%Fh>YM2u;xOvMq@s><#uyHHda~Z(>h;Mwm0TlURh4o-(((NrCcF5x1EjSDP z#nyxsePlGW8^gjJb~yJoV$}8b%k}m3prD}G*x2|R<~%y@kpcqhE7dTihbFK0aZ3Tj zTjxTax0(@wTZlPWSRR8!ijtBN1CMT-vad>1o>eL;(=AvaZ!;fjS8PY&kK0g<=rS`^ zuhU+9{nPfBXlQ7tsXLR!-8Ui-3sNVPs+P>;xQSCZ` z>_ALLt4O%*jJNBAiHTNPHwLAXrEK7VH#g#f3^ z_lSh57CF1PHZ*jO;e4-^l_-QBqm*AgA8QZs$%AVF;8`F|C$XCZD<_sker*4Of{3M% zAOGyv#7J5F=r6ssv8b(BWXNAv{v@yAD+!WD%hIQ4^BG4VxWcP_-!+*(p4(V!mf~LW zL>9((9O)e`PeEGW9ry0CSkDtz`;azkO3Y2YH^5+=46N53J`xQ0V<6V>Yueikd8M}i zT?RMLw%TwpRTjI(*d!r74R#Yyy!IgEcnA4$))p*dGuLK;9B4K$l+~$=06&uYl8^ws zJ^F6LaYiV*EWs)BTBz+SA{=-T$+d@HA2<2#Ke#QS>W;w#ao(2`E7FvUc-;$;0F^)L zvEc~JICeku*w`*#I425lEqIS4ZY4k43kZ+b+vmJDJao6?w+)MqHV4Fu=#mowy<0qX zj~EtzpmaqXwVN=je*Ty)vI^km5gxvH?jgZ5(8uRYRGU9iZ^qVp=58t`AdIXB*pK;V zrUoRAU06QTm~oqI!dFl2C45f9()b`Rfn)qrtN>1IHX!_H$I2mAv35`WJmg1Gxs;kv z&ovDpT@mu%HM-tio)o#iSx?33gO)+w@dnBVobH3$qlWw29Kr>V-^ga5r>|zEuTJm-;H*0m;sE31edN()mF#`Am!pxzqg$38*A;m@+^Xbe zn~|~Z;Yixp2<-6d&47!Pl*JCS$^_PWZEbC~D`mx5$;rtHgYiqhrkvHYa_UoMZ~S?H z`;G(s+C31zLIAy6fIF<62@VDVKhI zB;M~rbEozCnp6%myD#x4bW|8MlXY!?$BqFijhwd@!N6cuIcy_5+@m6SxcfW8<C(&U=0yd!DP#?Epr64~>lE$Gv&;Mytr+a3Qa#^NmpZyrXE6a$NgWe?Sq` zKArGGe42PVH8s`S-%rM}`?zd+A(;CeMGSLO>>GGXUfP4-`M# z>@lo=+>@z<2|749XbPktXUW&Ha(1rhFl+z%_00|Q&Jxphtc1Tml-UW}seRu1MhXq# zu&zbN^;}7s}{M<-IKqtVkNZ6bnSu_!^fxBMQ zpO+6sW0388STU>fL;$g{G+GefyO-IVZBDBLoBB*klfn6gqY-m}UMS-!ysBKDQA^t;ke&M@-dZn$Y62z*K#Hwx;Co`8F@(DB38L-SMS?eRh}!Cc zTNj;jcjB6y+x@xh)vI+FIAtl*w z<-lWsW7_1p2zMl*dg}RTKF*3sm7>+D{E0>r!^OaU5hDZQaSd<@&;>WQnpxw%wmN&J z^DM8{*)%L&Z6)AUn1;84MRQU9IxaTK z-A7ohCSs33H0+2&tYZc+Eo!=g6?rnDS6Du8^O{iWhtr%{81qKVnKoT(X832i1Vl>9$o56$h%FKgI}Nx2YcbX1nKklc z_8ux(GbM(%+??m-W&U0b|KoZx2>>HDHa4J?W@0im;xgvhLq2u9O-cJ9;*1ooK|#3bSMM<6R29C*E zSP#%(($fA2qtRNo2M-9aM!-B60z3a$D6ab0Bq3nf0& z-4@LbsHP*lowxgZ_dAg~^pPF&o$Og+#mc0J+1j5b$7>F8c_2Y%E4}IwYXf3u%2ydE z@?=*4#Vht_sWd^@20QM1#@6X2cs0NNg5}>^9zHoBkXHv*vquk79Y?>FD5CNOd;SCG z+uPf_*hb=VWFODsH4e+8HQk9wb+jP? z!>!*!A(1&0e6R99=z02BfKa=z_bDI^PzA!vMJwA-bhYsX1m7GBQ4|P#Nn>=p{x?N> zJ&u3cEVKcdrXyvkI&amcsVGxo!tj=RFFbq1!!>}_ z#%pV5Xa1;4n*o*B=5z`;Y7+;!xhX%Mn;A(PVjt@C4gZYAyqIK{an;NW~oXni(}qGvD2Pqy4pG&LWp!RvDV3?Z!~;3Xtc3t+J5 zhc3H^8zSSs|D6erKwZGHpIr5oRmv1f5Fh2IMmifpCR_V-_i9!4 zyalR4rb-OaOs3B&1f9wgw+?;}!HV)F!tY$cK@Vr+6_fP_g|D68`o_Ctvjd-~{pt8_ zA!wHS_STE)i(Gt-T#McWTGSf8nE->r!4JP1sI$mEfHyK>`Pov!P7d)=zA8e3pA4L2DjQu6$8+ z65F7CrKAQ`=gh;nmr{M58eLX!_CcV7RPxa2ytq2P6K)Cv*K+Cmd&G?iwhYSP&E8#= zejk~7{C5hEcf+VUTR6KUh(2IVyNT8PyGjGkRV|2>pG7vo*_nay#E9d*wzq={^Jwzy zrGJ)~PTUv;35!h)+o!fJkmJ2OpFtq+MrtYvHFmwCh9j(jfx(j}*;pM_Re%+S@br5P zp2#m?cE!F8VxS@l#5r|!o`6FBv$Mhs{G$@_{N8Lzim1~`cG{GVDhWy zGR`pwhtcvItxvT9-%)LW6yF_x&LCn-Ald&_jOw8%)-4#aYCEQfT)RKzY%1zi*4f@S z_!TNX4jfyIAg)lqhMI~>T}x}b`yTJsSyWuRklXsy+V57*knYjZGN+vqEJWe>@J3vp;5y&mz`H{+!sf5KL)&|#3F702$j2_( zynvE^D2z4FK$63H@{1P zC|3Ee%7G0yL3Wj7kC}LLfaXY$<2e>6f7((}!x4p%la7Z1uW3h?&z~mFEtY8X&OT*` z;S^_&$@{G`8KzpJI&|gzRI~vq$#nM8x(WvoAFa8=OkRA9j>qVR$6s?#urOn+pz&NAEk@tn6isPHtTGo1{)r0hQx0 zOIQ25c|WiFKN+p>cQ~uyzEhF%ojc|8%hqInggA3M-Yqo3W1Yu#kcjWP4cT+VCzFaG z22WO&&woMKnPXNgo!^4&5zDm`FX-cYi&r?QBfe-FcQ?Utc8giWhb0@F0B8;fW1!3Y z!{zW-tQd&t1L#L66P>@8%>T8b(@yfR!s?TyV?*;_XXi7pxk#OovY~JSL?FXROZeTt zUnS`p5itlOGmO7BH31(0YetJ$atU1z4-jrG5cXC-cSN0X;uoQk)p>4-lSe|R>5&)V zz_>_J5MS8GZvt^AH#TPPJB6-H$}wwAbN@Kp`(=QThY(2H1eN82okevW=Hd_6A(9{v zu)rBdBA)e>gSMCazM972RItrNh2>AItbS;1mA6@wrkHNbF5LB$1GQshPX+(xR6*bO zI?j_8tNL!x&@`NETqG_i){&Xq@bG@xR{Yc*Alay{Ave-dw>oLwn%@~Q_Dpl+b}mTA zR~D*SaS)OZY*Bjs>Sca_{XUrw5>PT59;QH8;yBj!QMUHgx)ZT;LJ=BE|@W(q$0TdbQ>-Kf`<~RY->jLN}ES6nQ z9COkEB+z@I&T@E*|`EFbHJ|0P3*L zuYwz>St|ah9ibIXj zYOOd7CQSDhG@kjSp9oS%;-G(i@GX?tL1>>ea4KZDG&ET*+{c;^9?*eoR5)Ra zTk4-B%0c{5r70ydPi=)|;x+|$( z?He9~`u)m>-{+RKKT;%vKU2%(lQulp{<E)oWg&fc>2@20z{rU?Z}_N-FJCymx47$v#Tg4 zi3*8`Xt#pWn&iO*#SVUWr@Y7$Kh+iu>d*Rh*z)?VUs1Q;J32XyI1T87dPd|h-qyy3 zP?=s_T-@T$DJa0dpfUKcSXcv?jIDytcxP>GZSTEyLt85=4GO>{{4qSvklNo!IlfFu$q%Z?EYcw9@!T2N0h+R)r|#wD1^%Y+ zS|#>uR|R#KrC(Jyu3W*0hebq4dLPx+U%iS#f$EP6cbkbo)kB`Fo*X*`Ob~zxgQe!| zY-|hqK?sbdwl+HobuHvdlZeJi-{4R8(lML^}9y4wuRda zg*27JjMUU@{l9CJHV}B$;Qb-+Y3UZ8wy$-v8 zvvehveYsi5*t@M*vAXRcE5aHanRCvOlD4(AJ*cub#bU8CGR>e6vr@aE%(zH;?wq*G za-#FkA3s>|x~WmILXUXj+Ba&>laf|B+(Kb{U1U~XJzyOrphD*&O z!e7wI{cc2asWF##6n{k@o} zgZX$@F0@&|ZQZ1foHFn{<37+FGIDYnZrw6K>9Krfu&Yfo76+NLUEP)z;J5DF`AG2H zpo>Lk$4#-4>d8wdP0e1z<9%)XdtjgpW4i(2Tk_fI>63%zDI;YTJ;W4ZsV>|YxEe54 zHMO_#3(#X4%IRi-pzIgcXGBa~{rTrAPq4 z*WBDJ+V23sDR3?`GiH^;hev-_RFXuu0FBvKoNPd%OTmFcP02~_|Jdk~lzjZj9~9?7 zTS^$ANgy4iuA}1$n1!qqFWuY!&e=s{|ECh^09X~kV(JtdHG@iIk+FCQfc__qDgS_% zo1SpsaBRlv4{OE_mq8)O%E~ImV_ktR`DL9mUKjkOt{ye|0?zU|Lm|2zbSMK%4;V|Z z5JCa^ha59lC};>>@mz5Sim#xoEaASr1Wv-q`~AKaZ1tl~A25cdPoJi$QhY{1_nMMg zOeCw^g1h8hUC>G7aPzYI;a$Tj`{~k@LoRQ_%JidGF9HJYu3xHpSO^Nmpvd!Z1%tit u9t37S>IX;VpaH}l^ogR${vUiCGk?!aT-V$w5(4e}kb;b|bRp)E-~R!T89J%} literal 0 HcmV?d00001 diff --git a/figure4.png b/figure4.png new file mode 100644 index 0000000000000000000000000000000000000000..b3e448622cccd1283f811048ba11827c7efc36a9 GIT binary patch literal 26855 zcmd3OWl)#xzvoRiNQu%ZNOvP89g@=B-Jx{1lynIQf^>H)p`bKKH%NEKUOdnFpBFo4 zXZGE$gTnxKUiqo-FG5LC3Iml86@nlP8R^$55CqE#K`@_@VZm=cMTZ}PZzzt^S}qWT z`R(yPOcDbo5%?vNtBisK(lQzrHZPwIOuH-i6}78`rmL!hxvPhfvl*lyswOWg#z?Ma zX-2MKY`FOmzvzd{rm4iLGnw6ay#Ky&>=-|%8&dSC^ z&c@Eq#m&#fqxhUM7lOzknb)Ffo|(ICe!3ZFO%Hw;xCpYUHK+Wv9k zl`a8ao$;2~niTRdUVw?%&E2Xj^L{AVQ0$E&Y*RKCnb%xX;){b2fs3>D@*Ow5sg!cv zCRp?yR&-qJJh3(?)aY3#jA76>R`i}=6x+CGp+>2d{unrqe}7sjY7BnvhUJz9$IEMW zB1V9NQP8Hz*}%_6FhtbgCYPU7Y*D!rk>wMGiP}TU9|pq2qd@!oG1C(fhI+jm7(J#{g@88@%u{ zHWbw$PC>y|lfDn)VVF#sm4-6}&tiJNe*Jp3Gllfd2)upt>tHL);_`B4Hnt;zj?3@z zESFpPvFz;Zh)75zm{P0+y`Vv*zdbfK&Ni2P&-%}l{JgvmBm(C(Id`_!{r1Mjl(e*) z325Nr45;#*xX3}~3@jE+u&mhLHECGUTnnB6b(A~T|bC>z+>+8vK zy%x8fNgBE2)3dX`uFX>hes=*#q^Z|GSRc+UjQ>=dwH}lAFkrm`-*E4}PFXSj4E;)= zk<<5F4gn`YBF@godku@e*QbA#SRx<0Gn~&Vj#;#{t0W zGC=Ngxf~`WBm}X(dZn+YcQ-(r1BnnmZZ1Un&2)2MVrG~bjI3UVY4RbX-J$mx?UP1a ztKvfqfkk%&Gw{Yb5^0&NKsSdBGrobw3dR`K{jTdnMoQIrpGf+yVy5MQD4+sg>VH9O{WqxqX;kso`xs#a|#x%=+h=l|@D$#NJUA zl&ch%L=XH@#oJb~!EDUt+7fTzd{O$GqduK`J%I}rY7uoPZ;vgf3q_&Wr5KWtjA8Qp zM3;IRyk>1@>waHU1+I@d(wzEN^o!}>bR-r&LbFm+W!qUfYlUusI^K3~beSj9pVK&*l(|hlaUpRM3C_yP1^^q*n-`xAaCh(nlV5)pss!0!&hqV{5G(+_CsBs z60L?6j6X4NC<^}L^O}BDi~eOqPWM9Kx7A+Z>`-3!jd+V-=o4!srMs;yo(nG9t`5Dp z33x>t9;jb^#!B3U_^+85rkAt}PLpWB*n> zB|hwS)Ep8+f#aLZ6fyqi6f8B*ReOza-2#P|7B+bLOKsVdZpic3LAxL6P+hmGj6xv( zKWjZPggWv>F8Cq*zo^-A12IO4^BdZ`#V3Z`=`z#aUPxm`hUF#@f}5+n4Z}eM!vV35 zl%?|o>g0|0)3 z+wZjEHz;*eqKNtJjZGty#Jf7J-ha(1-`$i9o!`ho0uYa8?F(o{4w`ja*!xsPY9IB_;*wbm%op%Z4sWA%Ab-6THZ2259 zQz`V`d)+hau=t5KRik#tA(mhB%9h_rkr^QdKmTG^MT#+HazB@}T0XOe)n8hRu_@zr zp3PB)2bLz(H}@OA_ZHdq3_L%>%#qY3HUC(WODLKXEuC>N z&sDoL+#+V4yMuu=O=YlNXjo;!Ww3DDUs^}%mZ{@yiQT7&DMw2bZW9tO3OWaVUxM`^{XGtz#L zNsU1wtF`u|_T%#CO9Li3jokcT(tlK0nGBa74%fwsAjUKb_jd65J;aVL64qTau?%}-Q=TL|LD;@LJX;6vWn#xC`8kpJy~cH!EqeTQdDv^ru;HmcV+%o;NU!ADI?`9( z5)|yqzpRiee9u6 zgFg7Wc(@fS7zkL~c0Jf$YmXJ`o>^G@(!12Z`Ys7Zqq0d<^QQdO_ul+uCc;gTAtTiy z4GgbYuMrd6>~Y>QPu8Pv2PV}gbQaGs?}>3K2Hyp_u6#pzxTn~=VRe9qq_ltU@63h@ z%XueLW9cNtx?_AeusVPd4QJ_B77*Wb{3~K~mxtp`M~ru~Q}BQYebh7#MC69)HW^h} z{ScryIw=LZsG)5f464s(rEOZV4+wrkp}&k;)$Wbdm9O`?V}FI6MakZuZH0sOojD~Z zkYkGd;TSHep2+yHbov+&VucG-_o3u;D)hCJBg?y&rp`i@>^&-9DPdj1d@KPOD z6)_dcy~>{1@*1jh+r|yy(3}Z*U~!|@Wenz=_Q=EvTlE`h%L)JJ#>!nC`p~d2kKJjf zt>Ih*w5OMMFCxWAk_x?+*jPI8tsn5tYria#&4z@^j!rjANrft#{G?)6UkOyVaPTs< z)5ly}M>c3R`fkiA%@rn(kdZ9cKre}WB{`mVF8Wg*aAX`AGsw0T5jb>bPf%1cQ-e;h=>S)rC4CS$qLHZeLQqzwMas` z4L9D2trHNLoIiZDO>$MJ6XIt#WRza5ed_0H?D31!i4PARR%B{8^|hRFq1YaLh*q+& zZk(1W{gc;+TU@Tal&W{Xc}nWIkp8+0@0ClVlKS=(43xI7ejVj^tsC=bXiL0FEc?+{ zPfsT2R~1u+jcdB`HSF{x=RuI3etWvgq^1L`n`WM@QTg;5HgEK`8;|`AjyGIVGO$)* zFflLet6qsW(Ao`&-hK?~@TT89Wgz^*b>45Usm^(pgJ=4Uh(8QWbBxkDw5M!;-fdAY z!#{zh#}-z%4qmwVSF#U|I26YG3K~#9ZdX%vC%tE2yvXH=|0ra?{8^K#Jjp<7khKRKk<_td5wyxfrX|OeQybl8;MvBTqb9*S0XZuV0I*H_9 z9i%6^=AzFeL?^fY3YhN-iV=*0JA^8|>k1D~Y3On^EPpS(!cPB8J7`4@Qrf0>Qj69t zQz+vb*3;XrOlM-Bz;f2-{2Q1nkT-Pp*T6x*I<0iD?yt?>(EHdICeP|PTrr)-nJo=} zW1c#&#@(3tv!LlORi?8!PzT?B0kaiARpQe#VV z7UR9knGXj|xuUnS9Da4VN2e2hm#ixrHX0E%=sA-nW|UZfPT1Ay_tmRruzx4U;6i8Lj@06R1Xf?wL3s~LDxSq^YOU3ctTKhQ{eM~z1-rSbz0;6 z`KbHs&)6Klu39#NimK^d@4q^0?GHrI$&|)e-4n=kmn%_8uXp2FG}h|i+NZ{Ee?Jy-}O-28NTHS zBfEOM`6k#oUYfbN2ILYj2mi7KdzAb2dvv+coghJ%tO(AHdQt5f^O(gf!f=zUo~K$k zKib)zKMsZ?LgED1*;RGVJ|0(08SBKtn5m!Y+Boye55{2*?*sA|b$zRxP8?U~bM_B- zoz$Aw*;KJP&pd7aRu?Rl3mtpDRW8CMq>qv>7+06&T!hJZrhMn$Yb9j;=a}oew6<|= zo>G`i>ZCZVzv_q<=;|#>+JhJ&YDpE%@g>xCp5C(ii*@$fGt`#jPpug%H$-n6C52rv zJ_LNs(*!_R`BVP&RhQ2GmGGOp^z4EgSSIu!FFvR|^3Nub!cCp1{R82^FlX2-0@qJX zuiymx)P4Z{(ui_ z+M|dNn_A+9=CV`72mC?5%1?amm@cGiwRTltJ>hh&JC8oc6jZ$5MV}w1L#E2)6FS=u zsx~Pm7(fP>UThnY!74tt-1hNSDBKZCpd*u+3JNQUaWd4NR+B*F7w1h5JQZ-Ia54~( zRH2*~rxcw@kM?xjEV6j#6E`v8*GR6@SoHbUn5q}YsYT(paS8HQc`-lxky#n=jP4QZ zf2~FHkB=L&1D-;SS6PHoy^2lR45dQp@5$cCh4%&t`TfQAA+3^nnbBu1p{dR3W)ula z(%yb+oNsR2S!iH=ePi^)?WpdXwCrau>ApE-XJ?;sl;RQ^1u6=!887$mDXmHEq-3Zj z5g*?W;^Ha4M0G=Wj>4OV_vYBOZpm@9{Rsj~es3~Bpa8m#70MvYkKbH=hBb2}_Ho6| z{fIA_AtMvOTgrT4ITP}Zj>TEALri0S*NZC@B~zM5+LfjAy@*36TQx(U!2ONz1LL1& z2!3I9(ca6{nlb)(u)75dGHfA;VnjKJC6yJFq-3oA4AyqtdCpSvI#2qgQqoWEdN zMAt%*@l1OM-xAHft%{PV|wOnZ^j$k8BP1(qhe$rVMB$4czdz>{Z zsTim6m2P+lXZ1SW*!0BF1OS(5tkgKX5(?OmdB^gjg>x8j?+A0_PoCdU19P9%kZagH zY-UXCVEU0k={)oAae3-Y0>xevQ5%@ptJcA((Nl*gvf1=~!QtyPBIy}|2ur_XCSrWca&QFgOwKSi`%Q>;ccnsH2 zOXiTw{2{(E4R0xB@1g*(T`!2)sGE?**I}9mCYGg)Y4)pk^ozIBYF!_ezuCNZU|0WI zqBv>IZ8w-C^w0({7q2Ul8D<_!^wB4q{!Qv9XE0U$Y5+&M`93s^`{y_Krsu}DiwVUL zzdey{;OZXFPf3EPTad-dKN`h_pT~Kta4I3R)NdFruBOh}#3TN!SP#C7aV|;km z|6mS-##(j3aR0?gI7vmn&rMp?-&^Z|11z!hqr33_I`+~JlR_Q=p80f2EQA?Z4^zEu zmMk4^zx7`2DK_e{;CfTekYBly<(3bE&A|$~)Db#apY1xw3D@ZQK)4Uq=@{|V+58$p zf>QeaM)!53N$n5Mb@mY3@)zg0I2URSC(17qHnvJ(m6N{@My9lvK2Z}wUvS1khl8H+ zhm{~SkyV}ZM2gJ^F%|1g#lm&s&vTm8wu-u9llKh12K^*&Dw$6vfv8P8XD z4PJk(v*$Y1@R4R5m>471#r6m>zMXL(b$9!V)f%Wbt*cu;nrBI5zSYetfOm(bsBw%O zc5BPRLLEG_3K0DM#q$^MvNe|t3|O48in8F%be@AHV@kZ4RVx5*b!g{hrm@4M!1WNC z+1I`Lnj#X2ZJ)fmYTb`m;SnsxxGWVFzzUBD+rc6RK9j+g-EDo)h=iR7Vt^6X;;rJ) z`4%Ijd7QYK4=veaF#njjExLIsCXj(f{5WXW2&q*hM=I>yD+Ajaa=RbJjExZJ`~H_! z4dx;E@79$mJ5@aVP$#8ab8Wrtay;Vh_THV;XBtVr3rUNWFoy-p$ojvrvivr*#$C@# zp~Kj$9)>{6_yYrk&D@87c&;QE*jpLPj{=NGkEx_wZz7a!#5YInZ%%%e0{x_hUeHZ4lC@nWV(X>ItBZ^AB;CAC7n( zWa=<^!L2s8=WAJR|IMzL#{*X)S5djB$5-%a+E+^*v_}|#!n!V(9@lS;ypVeBIJvNiv#N$=(!9#0l{8*U0gL^IDg#& zAJ!)8KQSGCdj<<46*lA3l*SFU;UU}k>-2WATJshHujog|jWT#UI*cM@m0kv{HH5(j zGE?Z1{)7(*=aAmmzB&Blt1J^=xPOUtpNKkhR~|=DaKgjo+B94HP*ku|RAAh6#hTs4 zhPxnA_3h}a<8{`8%#GMgz02l7XRdkOv`(9iBje%1HqKD?mre(A2Qe{$P!!mBvp7q* zo4423dPg)*wG!l$B6GbDOB$12{DEh#lpG0`f{An-4CWc$W4BSqo?q?1TluvZusWhc z()J{ZvDoj(8F`6FQtI``tRc`CEzPy>=9#3JbO3pcoEVLhLrx`s62V6e@#eR|F+&iZ zDRnWc-vq^?JbO@Zr9ptR`D{i&1lr|FlFyYKfCLjrL)(Tl?#;c3*jEe;KB zjG1dYNHB%*#CIDD*Qj7~bY23F)`V!RH59=4pdc6ByNHao2i{4G6er>jsqjOfehsdl z6*c^!z5sdfd*rJreT1J{(sgG(Y1Yeu>M>3wcMQZ`zZSt}3uxQ6%R}{>Wa{11!ILj= zQkNJR+#9c4UMj86gSH@$ess)PJoqk;Vh2FfKY*vH#WwLq~AbBtcQeoH8@~Bo1tx$>Vqp+g)p4(S-Rw z-WzA62oVa3*sM8?+7Ejj{D~XL9siK)(oZR8tzC%qt(wSo}Kz^%qOy%)$$aM&=Ya7&z*c zUVq-Y(Vv@Z2((P^-k#3eO?BM9Oxl^n@TE^bue~tbZqz9}c2MW&dTu3iR1trU?y#<- z1CKfeRME;07!dJWj)J;j!mBLN&kR?^r523U)Q}S1_EuS;egCr;B9tTK*uD<7f4_t; zsB0Uaj+)68>n{Ez`|2piAZT#hWJ$%eh1ka}gL|DP>DUxK zs++1%uE;5)?~1Q}7rY2DWQbt;DL?s37DY|yZzUr10 zluH!GG`R=)0#X8coi65trI#I-(apDW0G5E%8h07|U>X;nSC`0KbI zm~u+La|Al(xow-6Wzv>9B*V2a&L%ZJtq=2-&Dyan*DM?mD=*P9Mk7Fg;V@oEe;X3u zz|~hIK1Y=5c-M60cqCmRk$k%N@Fyd222Nj%a}~jaXN%-+1MoK5Jtwd}$s0GaM8amZ zX&qFiA8_w=0DTXfus(x1(l^77e)PYI$jfZML~ve3>4*EZ2ugJ{wKRWZ?PI2a)#k|ezST|pN-(Z_%KT%m%?;rBExY+qvmTt&FP@cyfVv<)iJU>>aDSklr4E| zIlT-sL+xNX`FB>E6>==7r|YB{>6XuM7Eq_Qe4{WV{qy`DVk~$;XR&agQ7q0Ej;pVn zy`$%6C@B3g}g8d#WM$iWJy zYzTaxEm~z<^4OnP(0TpE3e16RD&LP=9_eor^n1an)WR0zL)~K@pM`yVcG9hSb`h zg;Iiht|PI8enMr9oANHFm&|IR0n+^x@l{<#$X4xl*#Vqr)HT`K)4b-kTy{UcPov#i zl($zjw2O^;W+@nCxLo5eKWEfDYH?u*Qc7n%TzaM#2JD!L%w%uh%&*96^n3EYt*i~T zeLe#bL=_pR!1x$td%#GjOuy?`yyYM*5f;GhNB{iTOgZ(HPQX!m$Jit**PgDVStn{6 z#A$B(ou}a+jMsL~;-Eyu-f=_3CO!4Zk`jB_|J|3m+4&ML!;|ll?*oR}RJQ3I_z$0> zVu%bs%U*G#tY^m&%V5HO=-`M(I$0!FE}A#uvS>Z;@$s4cpj|IGdQ9a#+3{SuhSko9 zN{gg4^Gz9FTzd#tm#kCQNJy?nEsY zl0+$Y!(SUfAv%uS`RX5v`AqY^9PU@p4J(J_>k9peG6G%+-Vy(Cz;y~EQp2%M`D7l2 zDb~v;-!OEPNu@Q#SeZzGQwg#neCM{4dln^MXQMojCpK@lQ~f^cE#r4Gj{3+gY9lK4+(zA^14Hk; zmS#LXBYE6%;)A${Vjsqd-}kNhVv3^+=h+%l>h!3t+&*PvGZ+x)Xg~<_rKF(uBG<=M z6nf&iL_ej);Ga`tM?DYuX1?iFOYVG{0i=7^@jIpR(G^22qu4|6n9r=(=Hov&1YSW> z9a?ilnYv|dfCXyJI(KF-^6y{mVVTB2p8&!StJ7Jn%~Or?(~WpVkVP%N83&*Y;6-uE zd~`N%Zj=R)#acf7tG;MP;=R7V+;>C?Qjf$h2htbQRYERTm5phC9xfOEECN|I=yX!Yp-iV?e>hj1h{xXh z{%Z3}`={b!E0BGNZNAzucz8c_f!B~4uY|Y7s&B!w1t#?u5=wt=6nJBs7?)3;Yp>6; z{z6MzJgP6-UwiG<20vS(z=ILLLD~JibKSxtXA_6YjzgD7U0UyZ8S^3JfZLYTeok;s z^Y#iwX3F33B@PwW^0VrU)YPK#^2wdEw6wIKpdbkeiHij<2W{>2LoUhus;^(I@2<}l zI?7&mP-b37RyA6c5N52r&Hwai_1nnF=C4I`Xdk7fFvc&rq;5^I;G;|dm$Y1TsA4(! z5$-h{O2#sf;(G|XE>D5YiNT;!S@pNa!TWSrgFs_3?S3V@H>X|i)m7eIOrZx6H3lhc z3MJbp?g!G*qt(9E)m07-4j6JaHm9@giGz-Z_CA|e*IW6~(a~L9qG3JDJ2z3J9b!U6 zX^Hwx;#g5v&4G-&J3xF5rcLAC!NHiW>n2cQxK)-vbtSl2m%j!f|Tm8j-e;7Lksj#Doci>C3)XV_w@AeEP8VYMMp*5T&~83u%xD@UfwJ| z+z)d6M#sd&#KKDX_)+j?w`Qz;qT1}|;nG}_3#Fo32I^C8R4Rux?&T-IC+6FaynJ~7 z7uRFnWy<-)sI_ZAtyosC;Ir2(5nmV{Wa3IETHvx$F zlM>Js$L|2aowSmDfHWcTRR0`OSt zj<(9j%S%peCW2RO;`7l@UZ;ab!z%HoJEF^TtB4Vua)&+>BYDm8$%GjWUON}9sZME8 z&_C97*A9L+L_)-bgv~qUEx{~}j>}?RUS3vKFJIlRM6r~&>^J{eOmwDAcss>_>XzU8 zl+o$z+5Ss=3#X3doA5npLn62x8Tn52c@M{%)wP6+9W_fDV_#$i7Z&e~%)Fjg8gN(0DizesD22H~%;5 z)BvJj;P#PmpMMn^a!b?QgB?&-&&;w6dnf<3X!+!1U&^lX22AH&!y3cqhu5*MS|LNY zsRA>HN!7mh*0n&PWGNf|pl2s33!VT>Nui55mp8AjOwG*;M=a;6Or%s)aG>#t37?zU zMP?=@J+Kz4t0yUb`8hDa=XJijF_h79e^%jjwee9Q9(j<*P@W9F>SJ6_S-OPiX%}t{ z>jWW6S3hey*{C%l6Zv2{HxiUHf!t?WFh*vw04c{gLny&M29#G`sSI(Nzw*1&uMfcni=puzK2^0*@77MYc^}w@of0e3xR$$lX zSROh53G2gTX4=K;NYMoQZ=0E`)3oBb`rW!Tmf8wW=vtj@`Lc;Q>JyBRofAfBMwP7| z)pk!{5m%Xf3zvKS?6ZuSni3wyvJCZu9|1{mT{`lF77+GQxi`|eH&lk@bzE4H(Vhfm zr|Fi_XtdZk+P>iVV024GiH`Fkj=Wl@Io_lHPW^5}%+Z58QY_2SmJeY^Ol2~3KRthu zkx?Q~nxuChu-lD}qOQV7OaxrA6lNIz5U5B!N%&&3 z>bN=l?&vekH{f`=*p)5@f6L6D5xC6R36L|h=c4t!+)(yZy;K9%y2dO zD*BT#@w{e?kBB3^-2LFZNTA?vFV5@uR$-Izmzw_*mg0o%Kz%?qSGP^2+P%h}Nn$;T zzOst5oC1V4WVv{0Fd9W%@VY&7WMNt5$XR8`CQu8XXycHlQ~KL2r?ja92;~ufL4<&_ zREhEozU_z5+VH~>Foz`4vkunosp~c%$p@-21N7h}weO)ErxEJZPp*uL zBlIr=amF!(;R38sKaS||K|ibro_b_@EF_1fP1=qbEC?a(=;7Qi2l*^-1M%Ga&Ojh6 zQtCZZ^=GctqCA*v*k2S-!$rqqkF?qlX&Tb4JD6lKK4Gv=(QTU?3hi zTO`Q8%i|Zned$Um(Nr(H3hQricWSxy#aJZ4!h`ha=Q<;VC^EN)L;KsKaiz4SOEa31 zs2eS!D+WI&fZyQa|6Q)=PXg$SnT}TF&&9weSh9zEHd+uuS{`He*)W@k`ns-CJ4M6? z3C-OiEV5SV!YIBNOhu%`5=ZhUJUWz(+{xJS@9g6)r`T~BKctg zK0ER_cM!l!U=x3#^pwxxiL-@iVWvITdNl)J%j19pERm+Z@uRdYEg zF8N8B6DD44W+l<-XBeYYk<{O$2~=+*kaY1Gf;%u+z-h%?NQd@OGWh^MtW{^};-i6R z)`^Gw4Ot~81nq|DLmBOFi%Zn-oLLi(a;|LSGTE4ezNuG@8E+AnW1bOq6Yn;tbrh}r z?Rf#Q3dOT*+MqztHDNs?UNk$5$Gc?RGP8^Lql_!!YAUP&gEC8jH=0HP( zT#na)z!~SK=Un9It7IoL?qlaudPCos&UhY*J~khl&7DAt#^CUONskFu|l#Q{IrLlmkfyctXB zH0GPBHB9Pl!t=)+Nvt%S!fNz{a+&cJZ9KG{<-s^^o^d5Dh(KROIvEn*wwb21^^Sm3 zX%ebI>Neue)f>q4k3tFtSr5P!A;;5}sp(f#!E z$<9A2@>|7P)F$iHn3g&c4kker8+s7hifrXlv^fe-`t?(N2E-L8bf0ldr3sTJIO;=< zJoQhtR6ww-b0w!_*@%0ghpMhnqPRr0d|JSz}SVEve)8wbg-Q0A!Ml*UDX z3RilYhPohI$;~}(lLqIgzAF##et3WDaa`wowwCgzn}zZoy`WW~VwZ14inSFv!iW9& z7)>0z@L z_Wc1|@>8f16pw(gGuYk3Tm~v{JQW%-@Tn+WxJ)np2lD67Q(rSYRcYqhY{cDNE{Bdx z?QYZP?TPsBIiup~{|Bl@*!Fm)p~GZ87e9Vho%+>u$3opoUvb-GGA9u%;Y|e3yVR?u zf{8R)DcNa-KCkx0_Vx8yn480j7#pu_Pn0bBJ@_o6d722Itq}>O>YByVxM9TfnQq3T z##z?yKk-CRn0zv?0EXS~T` zM+Vo~S`+jTB} zRM5j;)tX(SpP8J}xVX5=%8BiIYHH)_KQ)F{70&OguP?-hSmO2l-y4<##HjH zQezJPC^0c+E}O1#EfbG2?q12gnOgD`l@yak#{u!TMfRbAxOohoIaW3{Q7jgNcE-dZ znV*C>tP*QLw!U{2&IW1m8^Zb}wl^BMf$CL*j6hOzFpNn`863E2)GoW9UA&jX3JVS0 z2Q1@BxrKc$-SD4`Gt82vTVs=>OBb z-uZLSr1}kV2&B==MRJmx2&+wC!w_0HmbZ9jF|G%3r7kdHa+)a^m?(tD8Nj``-L?DI zPFdRV=g#1OeguNefXto86^3A_yiQE!Rj{)S9@8-h(8l8OCp295>`VH3UUkfGf8?Ps zJ@Is<*8S1X{=P_nk7dY#w;j`syd0PoN%_B@|Y_v6S@~F+X5_&ZqBNHDbgCss)z$ z%dEgi%89x@nNPkm3|QcXJQQ?-8*+h2dgx5wFD$kHemMGFL{HWqd-#S=thYASK4{CN{ClzPr*yfo35*VEU)F{i|_X*)D{~hi0I{uU$Ol(OF5KPS zb=4`j?^7WHha`w7NY~&(*H2mAF#LuCl^_5xgM1^l@5x>=jZmg;ffMGgR!SW=3+?z% z6pv|fV@R_S4T4=x1kvn`$=!upMU3}Q6g&Afi z;>*0{8BGf3I_`o}Okn-F^x6muTK{q0@4iR|DcYL39W#2o9z1HQzTi#1ei_FoeFwuJ zyb(vcQl#$Q6v8j~SgiZMkD817C$!bPi@bhWX2LDG*U9jcEg$hO86?L#V!j*&H^2hY zi$uWwkcR&K+T1(na#@dInkOAWBEWupI`5yLxd46^jX=Wnk%Mj|a4o&-S4hDMv0Et` z={kyI?1{S>pj%dEpk}Ou3Dfz;RpJqY=F^QEA=v=bD`*^oJ@t;XBIHn-DQ788?{kna ze!UFoQg8AB)xsybVL*96#Y2-Uh-E((xPWpV7wq#mN)l3Y4kpW7zR?v7=eoD$lkLvj zq_om*!oQtg69}sO)R-oh;1?VW#YTvdH;X^ITch{-y~U#EUJS@7vG*Io#21 zSM7sSKQlBvAqn;{kOKh9wGxXdbC-;CsEvt5Fi-tLPj5yg&klh0VH3CL5Gan(&d~2$ zpIGQ5#p!V|5i+tE!Rq%+q^4#EwL<6IYRG-v#>=c2;>P9yph&Q;D%-Jz;K5N(I)2YO zqvW6muAuxFIMIMMX|mgnO^w1!o4)|%f*KROq3$$Bunk^ESnhQLYhdRFBydg5brQs^ zF4}8;o|Oh*jB9&h4DNhS zHV{gnrRGhdBdzRRa2N4NSKI5Eo$Lov1*z|gDqt3W>uqVGpLA%%NYWMJLy%pqO*XJ) z?AWT7oJkK76uw9YBg_)4JZhxM9Yvpgl_MB#tRKVIp{`cBVe0O_xC3@Ypw0@UawWRf z&@r22K}}HKWIqwFDZz-Y?L8H#VA5H%FQ*U_+jYp-TaXU*PH;b)>FF6sli<9z5y$lq}ofI^LAZht{k|^ezlZWbMOqLXn+-x{)h4)3t`*G z-i2bQd=64fsO0&#bcyPR&~`a9K0Pm|B15aHiE8vAI#CDT`G69PS6XKDDK@us{0C6V zvrP|*)oyaGNud6TVn6SVcqJD4qByREl3D-@n?rj6848+}u2DEfDSsi^=-DaE8wM(c z0gBJn*|&<+aiiy7v3}7{SkN(juDJgO&=IE5@nu!P&Rz|`Lo12?=8}6_}rAM-BObioN`=nif(#*|BuRMoWOo8t}g-gE}QD zGhiSfSb&t-wpwQXqAkzX$o+u&47v9TqY^gG@6|(mZs}+=oLxiK6hIV3X=f zIUktNPA4#oFMj)>PsjADLu;FzYb)Fw&g)(|p#Oo?dp2#ExmV3&U178EerlEaWs~Lr zz@c9ziIMkw18Gw)QNwnhy$NI$=p1q2w);`wi8%uYPWw0lR)-Zol?r&sQopmC3Ip&I zk~+0g;seNbR9BxM7-5hIFuo~!%TA8B_*#lG@Cb$YRi9JV&SHzkU5MbiwSQnhRpvqy z*pUYWs(5xl1tRg=njXibaHxlc+(}4D-d*2%D1gWBW;FX;ADTW|yFN9OE4XQVhtgCN z*?qH2?LRy50pxFqsMUUtZH^+s%IzGAOT{pC@kr|iukkF~#0{FC0K8nQy$7s>!dC*k zyqAB!^4vX12 zIewtH7u4roM9^j*dWx#5j#c=cl@QrC!jQEecB2_QTu%kH4_-x#lj6)%3F7>6?ZLtH z0($lqAA(Uy6p>D^@PCRvX|8HGzNZlId@lo4a0DNylhq7>I|Txxr`R8((PjGHXB)2Y zEYn#7F9ldm3oRZ8ez$#o#JrA66B-74+uO&UF(a0xU%y^=V(3sQf7&ofnE5Obaf z3qXZ{3)mfkw(fi-?;x?<4p^Q{)7n+&=K|?9lx6O>%AZ`cTLySRP%Lb5i+7Sw*GH<^ z5d#OoI~Z)WBB1&T>_WGy)|qd?vkuE_&a3pj51Wv8M=mx{JD!9oDW)%i(LAo{`GRaU zQ!3?&^8wMy)AMp8!=W(D$;pY#@3!g?0T3X+hnvNmcN=M4U0ut|%K-Fe_t$`rD){$v z0W5NS+-Dk$1=eu8UP>gC)B)`HEFAZPff%Ra4v)f9y%XjGM%1`pOEJmLzg1)aoEtjcF7xo~>k&AS5QYrfvEDMN{>Nk* z&@}#sjHQIIZ0+pmw7Bn0mg%H0YK=}!eHr_!?{{~GMCwVt#gOGN+;*AhvPr_jU(gb} zR6Yq(Aly-23_f=t`2Z>TRw;9^89rBxB>J$=wby4+pokc}-0S{<66Dz3cd)Wa>EQi0 z9Okt@CAanynz_9>-H0aCQO%m%78nz5*^8&e&HLRxVzfm(Gq#YJ{6x>kqIvvhXXQ!v zujQd0ak~*Mgy+~ij(m>Y-^qo5H5U^T163M!1(H2n_JI~X1!OV>_+jstqrj6NgA)EF zr~j^Qx8j!0t9;j3b$NoDj`nseM^{1&XV1}z;cjI?vWK<@a-`?qxrklNCF#Jb!Q_Qt znVa>AGqU`uu412jq^Q@`BjfEN8*BK|wuwbx`lR72IhAO!jz0L%hl=u_vWt_jkk|ca z4XAvJ8MEfuk~tu^_aVBv=7ad4bpnd_5AVWG7nTL7bkRqz&vv3BBO#F&iF!@W>;E`h zT%4S)Q+ggK2HwlvA?N4kqagHkPGI2eH19Tr;8J&z6e30VgN1ukycCs?qh5AAP@;Tt z|JGPSV~!W%bZD+4ry?b10kv@kJx`5d=o=ZBfqI1*@x+IwbK-3DNiCZXm;DXxX)v%I z4fImZ4v9zr1rh#p6=uPojdPD{%R>EQof;1p$318PR7%SV^b&PTm|Y}Qc+&_!{_|V^ zk;L$}6!Hg#7fk2X$)=9B_VKaq%?g(AZRM9Q7T(^i&CS<|Z(hOpx3{+gXZc{iZj%>3Z`GPb^eTySCrMt%2V>FKBntm z(I0q-%^~wchISF#^5S$(-^=k5QbYYiLqhM4*TzpO+OO$JOAR`Ffj!@)H!g_u(;txQ z)!}InHUM5bST>_(0}<|6Tea4_t5&?{wn1PG2=}N4tCi_8yN8C(_|sjF!H8T1=2@$m zj|?lu!H3~*oKO@&`+9;jKZu0BP7c*~%o6Pw^#SxapqT2h_*|c-j;D!C$uq4t#s8Ch z>ZaAz>r_JO7yZM(Lw`sK%Er@KK1Bxa)4is|ml`E3ESeL}j+FAfIvvJ8&mF7dD886J zM^Xa5Ogb4Qj@1t!5<`oYX-0I>m~8d`B6>%#Q&;)scT_)PeEES z5bigNS3g*Z9um=M_h^1Jm+Jll{U?>4S)lDPi(YIS% z&@$ID5=Yub$5}FuVpjpEgEKvrfAR!SIYWRnfCdC&sbHxE$xj8ew=}+~CLA$Ku}{wy z16r?AtMCu-=n`kg7EWrGacQXfWX5T|9|jx&U6PxXEBHkj#rU`3!Hs_GoV{J-G@ZIR zmerWTjV~pv0Yl@45=iIZLjaC6bzFw2I5*RmEK)fAAV6>t6=WgarWlr%|KN&`ags3-|E;IW&}iXn z9GmCL%oy6+H^j;N zGDlu`CR|r#DsR+V9-xA&Iiccua^`pq_?JZZZ!?9oS!8#l>t!cqmeBC6FNH`9IkJS& z+Px5&#ILrrF^FOt3a@*>B@g<(*JSK+3J|%FP$r;qJLsQX6xd=Aih_obFj)>@=((A; zu?|hZ73e+@GOhHdp}#%O-?^7h&I9U3v`PeW1kn-c=QC7BZsz zLjI7_H(W0%%y&*kS$(OyJcT}3o3xkN#r|hE1T%s8C?!5pT$16e^0UW14I)ez+rNK$ zExTjUQk6wrb_gg>raacY4gLo%mb3iSZ(Bi!U=L*`j?p~x2O*RF2`zLyJ~k5Gm`zxb zPLzKtDWb51sG=Q8SCd;!R>0O)yeVS5N+uT`ewC<&{f$JO-U+Fs_gU@D;&YJ$cM$4y zoD6D6@nyYNdb7CYD70#byr8o{!*%+8f>?-oxEPAB7F+V27f`e1_I{9kwj4+4Fw(ug zQ#TbbHJZ0@NAJGn)(@RH+=l2L=H1Drwp~(<&-1dZ=b>N^eKD>n=RozK=3Kq0B=|>9 zh?ZKcz3Io{N#VoC*Lz2U;Dsn{35JmZCp)j`)d!B`LR9y_S_>5fs5Y6(2sdm3vvWY} z`n%C;$0jA>!07!wnx+#H7?V*z&E*1yC___8SxCnlaZC!k2>M@2k{=oSZy~}J5)Th0 zb`BlkpAr9!HR|HKn;wG?_*&0Tj`zegYT0PysObV?-84h;=P4F{Pm|EvU+H4Y#7`I% zn$AO%kPwyhfi8^B>_Sl>s}6?<3D484j{f%{K+FSxNNk01CASU2-{B~n0kLP82W z=efw&(>Vb|=Q3;(wfDXJJt;WY;qzB{BfzORS?oVK%E)2Qk+RpI@9G);uYIxJN#RYj zqG33WLaxt)F>;STLF(@I>p@Jo9j>eXNFwPTRC=2KwM1KiO}2T{*}_H}0vll0mo)Z1 zyK&WtXzALTPpN5(pHp4o+sEF};?uD&3dd`_?qgc><27?vZh$ zB3%8+UC_zDA6(detq1jvSj?mB$dI}3XIdmQSkRkm;-AtO4*eY~t3(=1AZeqpOZw8h zRAICd;9f*31ee_Ie1B4**Zk5GOvWwa}Y8W5@CsU z>HcK0{#ShgB~9}HCiFWX1c)YZ6Gr23XxC8a{iVI4P`DhYn2+FYEXRmIC$6chHhw%? zP?vGMVh|lGfN!>#3^N>Dm*rPcBy#`C9Be4KNue&G@)aCCNZ~w5n7jG-ne|h$1y~sM zZmI;0-%{laA$*s4!5~}R5vO;k&|u~xg9{J*!8XG9DO@enl;kkI?8=U~CCg$Hm*GbL z>s*al=lSyY8L0-Yi#u+uNdtP)e`h02TmUUZC%g&f{)RcU&?DV~(>u5@y$%3%y+wh- zdWKo@?4~h?iF?Tq#OI)PwO^V?EQ6Xw?IW7UV>we~^XTFC!izhg+Oh`lOG&7Ct%#3M zJ3;ox^c{Aw>q!LrAK>9Z(Rp!|F_!#KaCBB!ac5UV;uZr?rx1n$MwoGYnnxKF?Pd=z z8(h2*e)Zb!CcNY?X9zAS@Etg)p0o1WHyS_*=?9;TRl#l?3QQq!gm@x5PAG`1C@MvM~mUtX_ouRx#w zjV4@ydjhq??H)<%#syloMtQ=oETbvYJk-hif9Dd7e z0BeqMt@8MQyG{^KHnI=6F9e{Uvu4`=wgdCAgp!}4u+|Jaefx&^X(>n`JPUA_+B1>h z(am^Ny`r(J)qQ!vz5;uDqy$x=ni5YR}OGg{QWpSqY8$WgTPyHQXpG!u%2IfkY8IZUa=+@i%Ck`hF7Z-E?)BE z?%WSf;!FCSA6~~Y2Q`iIGXr*nXu1rlD?Z(3bRPyNy5)0bH`qt3J!%~Me=hD^DAA}0 zKCwO;#9tw)xHnH6juCvl{+T=euOfyOwJ$@ERUOBDM>&N2IM+`Ad(Nhb6AHEC%HWD; zR>r39sErMBECUL$VR;x5hFBKr@xmf}Y;gHObcN;79_VUHpZr^Jz$bQ6U_Trgbe(B1>EYT~H#CfYc4d zwdp0dBV`pDmiS^-WBRpiwjkF_*g|1@OY2*)T}qn;M2@r)q@N^7Wo`)7TIVkCn%eW8 zhrmqWuCniUzXlhR(eOWcJZl0;5!hc9czvS(6Fo8|thPrbk}9)~92VgB*5F>WewU2v zW7EUrp7HflhrQ#3k`SbSI57%@2knQ}vBA}iERI@iwP9mMvgB~CX}1BH!#?cKv*z2N zt^W__$Bc^nASuqZ#T5)b60uD*L*N5{9^aByff^3H5pm$vQ?_fFVe^kT&ik0>pn1E< zGW8~9Z~bgE*Zr;0*^4cC_Yyi9a!I$Q%F!DTzL<<7xo`=%4iqG*^;_uW$cP;N=K~b4 z@BLgdJfOvKopO%O{CS$yJ-2`Ou!ht$X3t<;?AWvG2RPQjz-q0}&|04!aDO&87Shs? zEf8P>wo{}dl4j*#$NU3fb^1PO9~FD~!t3PAk^dUA<}R)peIVlq?x#d1Ee@Kly;RUA zr@LJ$P7dXB4$O6RzkymoOy?(Y74M(zBQW*v2 z7#yZ|L=={->HrrzO*-!34U^n++C*nR-Mih>@laH#xfE~)R&&l2ietf57p_5!8# z2yQthx*k1J2s%zp7q zpSOo?lrW$b5vaiiFb(=BHSpJ=S_7X(^Y&fDJcR~9{nze~ml>!i2AD95iqvgN`3s$A zbBbTAn{J{7yzu~^tru1qKGQhm_%Xb$eys^}1voVZ{KWBOVs0AB zOb*MasnT)Qdfl!^YR8CP8jMo;2jwkHD#u~J{*n-aYlG<_1Kur1PjiJJToR7(6ukvC z2C2C;&K}aKm%pTyQS)|AxNPDix*(!!_wBQl6Z@~K>S2C^{ zC?d+)&6`ArwLnWO1*9p!r&lB#Q7KGr5EV!?Ve-gf>H`_fjJPQc6g2GOnkZ>#>Pm*o@0 zJ;q>FoV1O3F4oh7Vu5FQC)w0D`;(K_y8FW|{oRkvU(&$HP0k;+(n$i$wqkV))yrcU z>t&?4fP=&Ze_W3`(ppSbbbs=TXuJ(Epp=8oaRN_(XfaP65jM-|RRIUp;?=nAfZ$o3 z)66YZb+=Ab^c^=IubBrVYX>e`MF?Bv5Y%LH5Hc$3ic(td_lC4i^T+^!aCPka_xO1+ zHT5n>_+N%20A#g7eX0;fKw$Ae_bmQ*kY{Q!lade*ZT0t z63}6KzA#J3=(6Fn&EGH097zN1v!Yik?3Q>N)yhE&41khb%8~rH7MVqP%_2R_Lz0UH zV1f)vUWzLz^;OwSnTAa=>sCcX^+VEYL3Xb$mRkn9^BqV>OZmE=>(iWmM6NCcz{_5@ zq{IKUTn)Kcg1XQm{eZQPj~nmi65_E-BCHhlCj$3U(lQ{-@(IPYg!oKAC3G(oRwKL1 z7aZ}NlY(+@&LX6h27yZBtV$PkBS%2Mz|4`>=p>RXZ+w}NXWh-hL=HV{K-6m1`N^r7 z*0ivVc%0Re95HWM%2T(NTFV#~!v-<}0+-F7r5Tr-X63}6A#`?y5W?S-DR`PkMb)Xm zz-w=8^|kA+to1|1Wl%Mm$J^t+1vANE{(;IZhSz` z?LJ?6_rE-y7NVYn$=Ev(zPE6_fP#n&AZ+ff7NI7;)&<|ofWOQRMb}Zd%j;b$ERgdE zxj!CtTfZ&hSz6fSzi`x>)44G-LJAFn&C;4(bIuH4(1`IZ9rgP(#Jx(1(W^b|s$#bW zVebUQ5=fRF?^t|!c%8a8y41A)994bzzFrshKZ{vSb;G-FG@i)&g17xpCFkySr`8sV z0s~P!d(E=VezP3?`n>qjFBti6^uC!ucuA_{d#AI+Je{}oTlMagX-w-M=cU0xtNdCR zFU9ADUhFmU@11+k%~ACK@z}1C`xB>s>S-5!`{gDNy!YI!Ga~o- zd3ytJj|RB;s1B9aDC&u|ca2P%ea8u7&&wHnw6`uwD8KbfhG8yhx${&Sv<2c;=xRfa$Fa0ZCPxl)heMDbZD#00E7d%@vKDKF0-Ws(~ zSv0mezT#x77v9-vOxqK6b@1_+B{6Fbt35L49qIG?hIbGuni8#4vQ;r&Gu`MgTJ{6G zt;R%mL*{9~x!IrUjOia^$JXwz$GQI`i>qbr^dCOx4;Z>>k+WxbKsxz7{c0K?Kgrwe zvpco+8~dG{2hxMAv?+q0S65d@{=Jr#mO_0=>F(Bv;B*fBNJsAV&pkB5tJ8kMMo~cG zfH(E@==UPCT;fODd}Hl0s`CB)se(B*fW)TqyddaqU& zw`02%&;EAw9g5_?|8z_5?gvTHGD-P>7<}E0z;mnoof)q3;Ch)XLnR(RROEmc#u+jY zcfae~QZ>b9gBtW^H^*}DEeE7<=1=UmEa8z#WvgXBIp*pVF!;HBJpL5B4H(9!FQ#{8 zT==lRebp~<-J@O)EuIE^*WBFHYh24M4~Q|={k7`%%+8+PO}H{;Rq?*h16ho<;`4yZ zr#!Q^1E|qG+GbYsHTY}3OP42BtziVKBl9mdV~GS`_us`9sCu+rAN_t9PtLZrQF&(0 zKxPsdMWkb}?3JXA-;Nri~d6^N&W+| zI3-V{%dcYW+}sUJssL-`3@!m-0%_&$tfEtYomcHNDSLa>4`x1-Y{8VCIxiYoVe>-!$MaU`<`rbff zR&FjPy8Hu?w?B)B&(BHK{TPK`4Gp>7M(6FV6uNzDZia<1+lagf;)xNxFUD^qer4Rh zC&PBb;sP3ik)yvJHItjzE!P&jR`zpUuwaI8YUh+W@}0`0%fRo0g^xbY`x{H!x3S;&3%|=!_L&f1bL9&i3tgF2{3AAlxfbBE%jltejz*D{ zzrFU9z>^1ad^5|(aP`^xe1Q6yNYkvBVk!}W5@YQn?a^zQprEZe;-dWCJMuJD<2UD&gbEH% zhFj?^y!JaT8V40xmRk%0#Zc~+@MEWWsfOe~G{VsN?nj!yc&O&fU{)=+oM}C|t#Nh| zrD?-56%Su^Ez3eAQ5cOa?}iX=(S+GfJ&(ENXW{aRW9Q}q&1y#Uzj)9K56L?Gw`?bH z@QC0;c%3v>@b$JJXi3n9AliTfIHMV*u|nT?q@u24Z(4D zh_Bo;%I2(>Q%$J4Idw}21P1E$83E6PA{kIbzn@7DOw3-7Ek9vv)gUJL;~zAF?)Y+v zB>bpX>np#X*Dre=ii*}?E6F6V4L{WXPu>7Lrm`G2F9K6-GP zgL`LjsuKMUi5{Q(?tF#C@O?bIyAK7H_U~A1O;)zHwqg^pq_P_aspdJKCe2roAjHn{AAF}TSFbaIrE($_EL6Dz4Km2A}s>#EGm%#{pZ5jwh@+}xFo;zlr@9*V~z z5c9&BsBtM>-QDG7WsZ)HSXfwGB#LlNPELRFiyOiuVWc?HuQCZjuT6(=^S5m6ern3CsE&G#n(6_4W10{N0FFqrr_sID7$#y{82} zTN@i2E2}3v5<;Aod36y|8(+F!bLr?G)!nnNID19uv)4+bTV*}C<_s=ae$NAj9Y>+( zs;ZD4fmM9J1x~BqON7m-%<_Ze8sc4y)n1qEd!Jdu`ca(=W z;d}*ibB4g)lU0++9sIt*DFfh>5BiqrL=3rxk z4{B*>wAJB8S5{YVIDZQWK+PJ|Qd9eLb-oks5x9Llr{4Md(@;}OJbfxCAONRzcxZC@ zl|GFB4Na{bGM$(|)UBP*8l6LV%*HmbU|$jq$pPD3?Re5K`r9H&4L#Q$v9s3bUGD>a zEiS^9+rf3O%xI0-az?u3?+Y^W^7+NZZUSjA@4|^#x3{)_>nG^GfkQq%f2OTB5OvO& zv0c$mXm{P7emsdxbGVIzll}Sg>Y-=uum|K292^{gRVc_!#TEp1+DLZ?Y3bnMwauNK z{;n!f;vAkBH_2u|g$XLin3*CSxQ#q1 z%rHV+6Qkb^VkOL>!K1DNn?)9Z5X1PldZD0zFEQXPJRA{JlrjI0-u{a&Yjn)BIUUUb QVl#xSq|&n@iI?yG4^*)Vu>b%7 literal 0 HcmV?d00001 diff --git a/least_squares.m b/least_squares.m new file mode 100644 index 0000000..773c826 --- /dev/null +++ b/least_squares.m @@ -0,0 +1,13 @@ +function [a,fx,r2] = least_squares(Z,y); + % inputs: Z matrix and measurements y + % outputs: + % a - best fit constant + % fx - function at xi + % r2 - coefficient of determination + a = Z/y; + fx = Z*a; + e = y-fx; + st = std(y); + sr = std(e); + r2 = (st-sr)/st; +end diff --git a/problem_2_data.m b/problem_2_data.m new file mode 100644 index 0000000..558debc --- /dev/null +++ b/problem_2_data.m @@ -0,0 +1,19 @@ + +% part a +xa=[1 2 3 4 5]'; +ya=[2.2 2.8 3.6 4.5 5.5]'; + +% part b + +xb=[0 2 4 6 8 10 12 14 16 18]' +yb=[21.5 20.84 23.19 22.69 30.27 40.11 43.31 54.79 70.88 89.48]'; + +% part c + +xc=[0.5 1 2 3 4 5 6 7 9]'; +yc=[6 4.4 3.2 2.7 2.2 1.9 1.7 1.4 1.1]'; + +% part d + +xd=[0.00000000e+00 1.26933037e-01 2.53866073e-01 3.80799110e-01 5.07732146e-01 6.34665183e-01 7.61598219e-01 8.88531256e-01 1.01546429e+00 1.14239733e+00 1.26933037e+00 1.39626340e+00 1.52319644e+00 1.65012947e+00 1.77706251e+00 1.90399555e+00 2.03092858e+00 2.15786162e+00 2.28479466e+00 2.41172769e+00 2.53866073e+00 2.66559377e+00 2.79252680e+00 2.91945984e+00 3.04639288e+00 3.17332591e+00 3.30025895e+00 3.42719199e+00 3.55412502e+00 3.68105806e+00 3.80799110e+00 3.93492413e+00 4.06185717e+00 4.18879020e+00 4.31572324e+00 4.44265628e+00 4.56958931e+00 4.69652235e+00 4.82345539e+00 4.95038842e+00 5.07732146e+00 5.20425450e+00 5.33118753e+00 5.45812057e+00 5.58505361e+00 5.71198664e+00 5.83891968e+00 5.96585272e+00 6.09278575e+00 6.21971879e+00 6.34665183e+00 6.47358486e+00 6.60051790e+00 6.72745093e+00 6.85438397e+00 6.98131701e+00 7.10825004e+00 7.23518308e+00 7.36211612e+00 7.48904915e+00 7.61598219e+00 7.74291523e+00 7.86984826e+00 7.99678130e+00 8.12371434e+00 8.25064737e+00 8.37758041e+00 8.50451345e+00 8.63144648e+00 8.75837952e+00 8.88531256e+00 9.01224559e+00 9.13917863e+00 9.26611167e+00 9.39304470e+00 9.51997774e+00 9.64691077e+00 9.77384381e+00 9.90077685e+00 1.00277099e+01 1.01546429e+01 1.02815760e+01 1.04085090e+01 1.05354420e+01 1.06623751e+01 1.07893081e+01 1.09162411e+01 1.10431742e+01 1.11701072e+01 1.12970402e+01 1.14239733e+01 1.15509063e+01 1.16778394e+01 1.18047724e+01 1.19317054e+01 1.20586385e+01 1.21855715e+01 1.23125045e+01 1.24394376e+01 1.25663706e+01]'; +yd=[9.15756288e-02 3.39393873e-01 6.28875306e-01 7.67713096e-01 1.05094584e+00 9.70887288e-01 9.84265740e-01 1.02589034e+00 8.53218113e-01 6.90197665e-01 5.51277193e-01 5.01564914e-01 5.25455797e-01 5.87052838e-01 5.41394658e-01 7.12365594e-01 8.14839678e-01 9.80181855e-01 9.44430709e-01 1.06728057e+00 1.15166322e+00 8.99464065e-01 7.77225453e-01 5.92618124e-01 3.08822183e-01 -1.07884730e-03 -3.46563271e-01 -5.64836023e-01 -8.11931510e-01 -1.05925186e+00 -1.13323611e+00 -1.11986890e+00 -8.88336727e-01 -9.54113139e-01 -6.81378679e-01 -6.02369117e-01 -4.78684439e-01 -5.88160325e-01 -4.93580777e-01 -5.68747320e-01 -7.51641934e-01 -8.14672884e-01 -9.53191554e-01 -9.55337518e-01 -9.85995556e-01 -9.63373597e-01 -1.01511061e+00 -7.56467517e-01 -4.17379564e-01 -1.22340361e-01 2.16273929e-01 5.16909714e-01 7.77031694e-01 1.00653798e+00 9.35718089e-01 1.00660116e+00 1.11177057e+00 9.85485116e-01 8.54344900e-01 6.26444042e-01 6.28124048e-01 4.27764254e-01 5.93991751e-01 4.79248018e-01 7.17522492e-01 7.35927848e-01 9.08802925e-01 9.38646871e-01 1.13125860e+00 1.07247935e+00 1.05198782e+00 9.41647332e-01 6.98801244e-01 4.03193543e-01 1.37009682e-01 -1.43203880e-01 -4.64369445e-01 -6.94978252e-01 -1.03483196e+00 -1.10261288e+00 -1.12892727e+00 -1.03902484e+00 -8.53573083e-01 -7.01815315e-01 -6.84745997e-01 -6.14189417e-01 -4.70090797e-01 -5.95052432e-01 -5.96497000e-01 -5.66861911e-01 -7.18239679e-01 -9.52873043e-01 -9.37512847e-01 -1.15782985e+00 -1.03858206e+00 -1.03182712e+00 -8.45121554e-01 -5.61821980e-01 -2.83427014e-01 -8.27056140e-02]'; diff --git a/setdefaults.m b/setdefaults.m new file mode 100644 index 0000000..8c3c5c8 --- /dev/null +++ b/setdefaults.m @@ -0,0 +1,3 @@ +set(0, 'defaultAxesFontSize', 16) +set(0,'defaultTextFontSize',14) +set(0,'defaultLineLineWidth',3) diff --git a/sphere_drag.m b/sphere_drag.m new file mode 100644 index 0000000..e75818f --- /dev/null +++ b/sphere_drag.m @@ -0,0 +1,10 @@ +function [Cd_out] = sphere_drag(Re_in,spline_type) + % interpolation for drag coeff based upon Reynolds number + % output is + % Cd: drag coeff + % input is + % Re_in: Reynolds number + % and + % spline_type: 'linear', 'cubic', or 'spline'/'pchip' + + \ No newline at end of file