From 0908032904c495c467c92ab655c614469b778707 Mon Sep 17 00:00:00 2001 From: nin13001 <32068635+nin13001@users.noreply.github.com> Date: Tue, 5 Dec 2017 23:35:27 -0500 Subject: [PATCH] part c --- Prob2ptC.asv | 22 ++++++++++++++++++++++ Prob2ptC.m | 22 ++++++++++++++++++++++ README.md | 28 ++++++++++++++++++++++++++++ assets/README-23f54e97.png | Bin 0 -> 53314 bytes prob2plot.m | 17 +++++++++++++++++ 5 files changed, 89 insertions(+) create mode 100644 Prob2ptC.asv create mode 100644 Prob2ptC.m create mode 100644 assets/README-23f54e97.png create mode 100644 prob2plot.m diff --git a/Prob2ptC.asv b/Prob2ptC.asv new file mode 100644 index 0000000..83b6c4a --- /dev/null +++ b/Prob2ptC.asv @@ -0,0 +1,22 @@ +function y=Prob2ptC(ode,dt,y0,tlimits) +maxit=100; +t=[tlimits(1):dt:tlimits(2)]; +y=zeros(length(t),1); +y(1)=y0; + for i=2:length(t) + dy=ode(t(i-1),y(i-1)); + y(i)=y(i-1)+dt*dy; + dyc=(dy+ode(t(i),y(i)))/2; + y(i)=y(i-1)+dt*dyc; + nits=1; + + while(1) + nits=nits+1; + yold=y(i); + dyc=(dy+ode(t(i),y(i)))/2; + y(i)=y(i-1)+dt*dyc; + ea=abs(yold-y(i))/y(i)*100; + if ea<.0001|nits>maxit;break;end; + end; + end +end \ No newline at end of file diff --git a/Prob2ptC.m b/Prob2ptC.m new file mode 100644 index 0000000..83b6c4a --- /dev/null +++ b/Prob2ptC.m @@ -0,0 +1,22 @@ +function y=Prob2ptC(ode,dt,y0,tlimits) +maxit=100; +t=[tlimits(1):dt:tlimits(2)]; +y=zeros(length(t),1); +y(1)=y0; + for i=2:length(t) + dy=ode(t(i-1),y(i-1)); + y(i)=y(i-1)+dt*dy; + dyc=(dy+ode(t(i),y(i)))/2; + y(i)=y(i-1)+dt*dyc; + nits=1; + + while(1) + nits=nits+1; + yold=y(i); + dyc=(dy+ode(t(i),y(i)))/2; + y(i)=y(i-1)+dt*dyc; + ea=abs(yold-y(i))/y(i)*100; + if ea<.0001|nits>maxit;break;end; + end; + end +end \ No newline at end of file diff --git a/README.md b/README.md index 3ce5658..a5c0450 100644 --- a/README.md +++ b/README.md @@ -19,3 +19,31 @@ end plot(x,exp(-x),'k-',x,y_n,'o'); ``` ![](assets/README-66741291.png) + +## Part C + +``` +function y=Prob2ptC(ode,dt,y0,tlimits) +maxit=100; +t=[tlimits(1):dt:tlimits(2)]; +y=zeros(length(t),1); +y(1)=y0; + for i=2:length(t) + dy=ode(t(i-1),y(i-1)); + y(i)=y(i-1)+dt*dy; + dyc=(dy+ode(t(i),y(i)))/2; + y(i)=y(i-1)+dt*dyc; + nits=1; + + while(1) + nits=nits+1; + yold=y(i); + dyc=(dy+ode(t(i),y(i)))/2; + y(i)=y(i-1)+dt*dyc; + ea=abs(yold-y(i))/y(i)*100; + if ea<.0001|nits>maxit;break;end; + end; + end +end +``` +![](assets/README-23f54e97.png) diff --git a/assets/README-23f54e97.png b/assets/README-23f54e97.png new file mode 100644 index 0000000000000000000000000000000000000000..a16cda432a292feaa143df042b19514abc9c7f3c GIT binary patch literal 53314 zcmeGEbyQSs7(NPPU?2t}pdg5Z2+ANzw&C@J9LJ}t5q?Y?#a$%mIq#$^QRk%<}SS68y)&gz!NSPcBC5-=D;=*RB73Xnoo14B<1+np4m_ z}v#_k-YBK$)1o~-!aFMJ=XhzY+{=^jV@`{*3w|9dMUqW=fS>Wbd(#=LEb zCj8fkBFu(oerWLG(Y6g9^S;>J++2fD!rY40(OaYKN?)$@bfqdr7%XQj4R)vHD`iSc zGe(OJ7iG6>ETsoe)^RuP5vDt6p=ExvAtX;@1sA<2&~|)tvHzRPa*1TA9`7~{U1n)D?{d+kMCIpx5Id zRKRiZ>mOH`2mJufv#EWzqU=-_yZ-SPbson%frE1M@{Li?+1S;FaCI$Gv3J; znN?3xITe*ha~Sbcnr*QQ`}--$r?0%I`Qa;g!EHIs>sZ0v3;GxCK85bhv_wUGX0|on z**?4uI$YwEPL zG$yo@lhb~zSMln>n}@5yZwdeZf#2U%z7JCGuM8ejCDy^P~(oN5z#QAgWXUqS@A6 zFI)Al-eMXd{E3};H1;|LnwQ?m8A4fGdnaK+Z+7MwazA=vvSRwQg?(3Kwa?zz*o%rF ze8uyU#MtYv$X6JpLzskQ+i->}l!OL6qPO;9Pu|4U`-XnBq#h_#Rh5L zkd2O9yx8{o|Lowggcv)laPIVC&mNYg^c+_lt74Mz(m+8Eeix6U5mtp=WU8#5Rv0Ku zX_vvuO5dnD$zM6jQ8YL_w{;SE=0g}f1(<|izkYR<3s#OMlL%;D2>la{de39gk)xJp zZ#%>XJ2f5puAHgf_*v<8BV+1&h8jxr+K#Cxp`T*AjcCUs#$jo&x8Z7&AwG0RmR2d3 z-p^Chv$JNkz7%IijKY&%|A^IGoF+sJDSi9mD*0NT0N13Ze+z8+l;+sYnmc;S?H5T; zP>`O8NM%4;M-+a|>&(?zyt95l+<#ZLz?FQ2E5$8oI*YbvV@hKu;}>jPC-xMUsJxmzVQ{7;Yq%Jei5nEgs+9u z!6PE0P*gmvguW$^AEn~rueoFHP;&gJ*NZg5y*!h?K;R9!qUTU7g=w7)9$d2YHifT8 z$}1~>S-eP;P+V44W)SMh1P#l(L;AqIQim=u`3Gl`3EO4T0TLs)6qA{j_W0eqca5lH zV!j$Bw#;$T-VAVc4tc^WAv$Fh6`lv&a9PV&!I!AYCr47K;>=$#!rwVbA6OO%E@-nO zL*f~7gM%Zmidc;HHvjR6nRqp)@{zClt|B`3ln8?-)`>KLrT;YtZ5qvAUD5{< z5)$(DUz1YB=-Zx95K>AIAKWdfX&R$@S-$>z24T?KFjGvWL{)c4oyxVHgrZN>zb(}q zvV3_@1#p>Xi23?V^FWDRvV&WYsfAzLGF2u zQET*z)F#%(UCs8joP<-XKkkt392}8m3kwTO%AcbT z_h&_M7->{{=}JY!qJrn=2hc9vFqv>82Ee#xy{eO%`2o%=wBe2R zt8rr7Dhyj!=}RlOzk8!a(R=GLe4bhu^(>)oHtS!*J*{Fl}?_g>l&|l@31z7mYgwf ztd4n4{?~1_KfVOQ`1ECO`IQq`ZL~Z33&!`&{nXg8h-^Axra;NG9DZcq-<93o)<*wo zZSBSP>TsaL@94h9?|)*HZca2u@M?H?w||qwEftLxx~+H|&U&bCm2Y%NJI_YF%)mW(e@<;O+I0<$+Zm9w z%6^tp?VOpCa&yK7rBV2#Xvo&DMR-dSw{F~<6@^#JfIq7&HLvy`ZAZ3BtcZ4YuR70^ z%1z@BzqlUgO4HNtxGKnhQqmmiGKBc>uvTI2XR#(qgGJ3-SyfdsQrsS%sUFCIonxS& zq%@uW^#vo|AO>y?h-G_k9_2d&Hi zgJ2Q>=&9dV;!;*t4(HVOb`Z80;9Ba-so51W9r&*M1|Tjy+5szBXfgO=*lxsWmV4V> z&~dS?*k&^SMRT~bPWF!4$1C?9L8P@9EX;gqZs%4g#&sMeg&zNFM|**kZMp9!S0V<} zrjMY|poc%62c7*(eo1yL_tmPN8@^k$$6?9R_cZy9XI@QvE0x(VN8MCKciRKaTPy8m zpAK8soK>-I&(ac1uXuhZMZZ2kO@FpIB6+0T#SEedzJ(iFd$3tlS&)^Nm)D&vdrj-r z_`6??jqN%5_4&Fr-agM~VPW-tW>wr6vTNZeO>9VCf)*Cgi`PBqh`h!kG0^K zIhWBr?VmX4jk2X>mILj|Oo<;Qd&lx{X?Lau@0XB}g+Jpp>V@WEA?zASHSf>kQj?<< z(p54uR`0#|jo@qC?>v*vBp4;vazzZ11{cZ^*B__t~|tI?dSFU z@p}gYc}B9D4${7SUnhQkAXhK8*8S~b>Nxjy<(hzf3Q~7K?eNe(U1DZSm+ADi!>4cs!k1t!Gjankx z2MWw4)6_612*d~5^`eDYW~lPCR`0Vk8>{u#q)8lnx`fmRO(0HVME3tY!0o~y0%S<< z^y;i|bAfP~Vl`UPW5-=X!O9=Q-gW`PVD$60{$?xb$Ktdc5SpYkI4ZW|TO#=x(U2>? z)cMn>t`vty(Ir2P^6fL|4S*To-ab>tBooes?(gsK!B@ELEqqKZIr%FA9jv!Zv(_1N z*rlGu`0jON-Je%NJ(Ij6uNX~vC7YK+e!MxE4rAo+EGua0NR70QF(PhCnwlA*I__##<@gZit!0bykc?K!l=EI^6r|~& zo3VZZP|W(APY@X>_2>e#-i5qYk-RK3N87IC8QX8IB709@!RF#cN=i!POEs#UHy2HW zM{@Kf%k^7C|GJ)~72uk2tq;7`J9@auODpO=q~hP64@IONTKR!j6NXtj!cqbRLU=xody7CErY%onXLhOa;rDR(T%b98I3_E^W78 z9_sHY4zJi>leXylet&OkrT64DSNcw_0^OPu5ZRj5*`|*d4QBG&%}8N-GBuKAMGx{@ zW~#?~6WiJcxVO@5)SH=AVI5$%C;VtV0%vy3y#`?J)C(>Ac5SQJcdA{++z(_sPYKVe z1MuO_xd050=L^c#Pjx9xEwe8-=GuLxIV1zF%Ktq)m~@V*^`lA$pbF)ge0}$Ff>f8} zHMa!V6XFokALlL6D#*#j?(e%=j8=r|Z=m z2u?9@EflkWfXvPra_L=Qq=XR2cI#qi1BVfUO5cV?&Y zk9KrBbz;8d-c3pC0aLEz(_n`7G3+u*CH>u`h zY3a@C&Y@S1{a=$8CLnH4joMEdBmmwuuVU-ztelvpT|64<~Re$QLIVPE^o+k z7+XtmAww;%Zo&N(JR1L}<~;BGTG)3p?Kqdz5+_xIcL$Rm&{v;OH~B4!a5EWokU+B@+l7`^+gzw7<{ zrAS&TxCkjYmbtIYDR0CX$Bx^js@#4@=R9|PaBI}#?}77LBMUPc0k4OUS4?1D0(#Xt zA^bB4JVwTWUAv5LA|Yy*uOtTtsOWZRQ{q_f@H9kC1j1y)Qc;UBa&&YIj(&Q{5depJ zplp;t;$(f`pt`z}ynNglje=pFm=qntywIG+4yxU z`&lazBDG2#9X5n2VYNyVKg8Btzqd(-`ew3>0Eubr8{@2XZl$xo-K- z^2J6C&DM3VHdO{YsO(~=Q?jxGo*houmgreJYYbATYR28lF6+Q*=jgo?`MYmdZ<^5x zyW_dACBUDo;tAicTW^nj`kq~}%Qh}h7eH$kh0S-_vMgV35p;&6nb8~ZtyBVos@%7Q z0M(4mQhV?KT5nD}?FK99RzCVmE9EJuu}nYYq(TsZ zhzJbs`~gL!%!|=3%O8pl5%?;nmD43OKL??cGc<2;v{I+G_LIltGbjbEX@R%u>C4g2 z#JxL1r7mgm5>Q~#h;w=-f3+#y%GK--kfLOL;V$eo&*F`a9R!zj$= z5mzHym$XJ@Zi3`-ZAF`ESMY~wO;pGwn4pADvoRhc&O>R+w}##Kmi)tK0|b}lko1fL z93Ca&CNREBeb;(S47VC~&uAlJbU*1w4#ae7F5@7Fw{~9@rJtLkZtceW-MMz_=T`u{ znayv^98Vrhb9k6+Xy2^duH~=?oP{Pi=s$6my6gE|uJw4ePtcUtNxJWwuIR*8${-+S zQtR$8W%JIgw%M4~F&@GgwIN(vqG3~rS}%0e?B`7S>ReZ9uw~_eBVb)wOg>~#grcJh zjG|oF_NV54Or>e<<+wb=ibilB41^+5g+vLy47s#GsfM*r1%kwCHUwu`@!R_rRO< z6B81qAe)pd7f+r;Rw!hJfg-Dznbw%>zXzj#3zkb4($X1t&`*)x!O;6^X0?SX#`z64 z+cH3yd$)S!yn+`T*KsrY~Qaw;Zf>l|Iyp5;-+hR zv}#8Uyb=iYgB>N8I-;qC5_F9JFO&K>3mf88uog z%d=f_abjnR%z9D&mpR36sb+n2dAWTeEJ;c3!mbhHab_pd2S9n=?W=^`z?bolC22X_ zfAcX9zJly!0oy5`wDk1uLJJKo$9_XBAkK*CgDt?109|pAJbtybs6ZX@E#S=vMAt&A z(X7bGt6DH#yFEn}U*xosWn+4NeYsU&)Jg|+PdWqxAR&6SUPdkPV-EUxX8%5VJDQ%6 z@npDkea|iJIhsr%R99+xS^=q_Y=+fuH?p9USYukLR*UFi#fzoQjg*%@3b-1UzxB6C z=TMKmm0IYP^|xwDG}1KMvv+N6Z7(isKHfAs^F85IYBNGIrsm$AsVhnnzb}2>jkeI0 z{W2-r+{TSEk9%KP6S;R&}_T=_6=20Y0C?2lF6^h_0%}2S}qyDjSaPR z2$kh&c(}jg{P$ps8BL}ADh`1^!@e^I6+sCyNP{RV=d!h& zxxQ~0s+9xr97t({w!He#TV? zTcb!So8MlY++Azp;5O^Ycy#U9S$ZCT$rS&wRjHm;#3t`h|HK70AhHL5beLFu6&`mH z9!MmX1X&;XNH6@xf9%x_U{vy_w)XxfE!n07n>~cI=6m-vZ1M1a?4{TZATS1UJEyBZ zNFX=zOcH~y8D91V(y5E{KiVt&C^R@cirkMtv%rFdAvE#-uQ7Ll?OJO8M|K%o0`_rn zEY&1%jSx+sW!NTQUd8PGIG1yVo`osBfQ3)<@7z+!ur|*ENHPyyeraCeZT^3wv2J+C z5;X(ycZEPGxfAHHuSu+d_4+tI6Szr&_x(2~_<4H2a_Vz=9JOBt?o*_QB|U-4&H4g7EnpZGioswAd~57H zMM@Qk+LN%ghHrWl>0-n^m&HRv#$iniRC-hf&(44&{q6w%rfWIDyvrMc?#5DQK!7!vdyuJGAru-uwf87YYiup^= zvdsZA0sFZiCX=od`P1cpAuij41_hBQW29F$k}ob%I@A>E)`ii^O3K^MVu7Xa5p;yj z(-ULZ_gyy}h=gU}Iej=|+&04iG0P7Fj{Nm&VrACr;UB9DcPRh zLH8Zy^Yriw%2`~TX81tv2i@|U4i!z^V&v4X``8Y2K?Y6-;Nd{&%WQ2oyzthZx0dkdg~gw#S<-=Zp9wSKGU-Tgu74WICl5SBSM2>aOlU~loiL%?HWNQj%!)$M zk3=shF>n)p9EH&==0_qSRD47f3>1wX2Q1h&g)4WbLUrhz%tt31cTYpt|9JsID{J4e zD;S>_>Exnu^I#)J-7{**r`uW@Ec)@C*&Nh(>}<5l24Vk2yYJ4Yc=<6MET~rcp3)%V z^Mmfaa@@4Na)I>M#;b&h39E{}K}df2ME8j&Ni>)mS#-W_x!<1vtp~WI6FBF!sm8+P zQhbj_X8kqy?KsE9J~D@W;A&`8>4`MvtdA&pg!)}il~SN) z4s2#trj%F`r3+hbNl%=nOl+c12$rhyI!5ELf>cb_-2Pio*h^>*F6CzlFAu*Q?Uta@ zQd;SQvNZ)LPybsNJlT(Af(Hsam7jAzr2QN51HB7 ztV@=$xQ&fYS(eg|Z-I;Fp4eS3T`6?4a&#;t5E=tOfxMXn2DPtvf~XJt2d>R4Tj{Zq z=XVV6nWwV5hLgpCNF+rzUamo%C)B`9m3)aAy%cKzn*i|cdaJ0_A8xFMXHq>Ks;ZZ? zQz?=9O5t1JI0Lt35&GDeWc;k%Zp4o+KzdjM>`7Lw5z9ozl4Y62|)Na`qK(`-Ltvu=y5Uy-+1yJorW= z1%ZbC9nVjPwXV45L5d$?6Z}fyuxC{*8L~}R_|jwdmc2zP|3Tu+i3$i#9|uy8(mQ~i z0%JGKt~=gr5A9QU*xUxoS7lqTK`xQ)v<7t>V0qjsYO~K2^13@2e$I)nZYRb{kLDK^Gktkq zO{E_(b?)ca!7ptkk}&z-O=er7)N}MC#O*YpbQk!$zn+PMv3*n~0=(;&(9l~P92^rn zz-tfN)YGP=q)5#j)FAaYQl`n_mPqa`_gK6G!^!2f#Z_;(ol&a ztP&*^l@P#1_;1$^d{%y3>WchJL*||h1p12|?YF|WS$#YAch(uJLk2BNGeuzi&Drf(NkJ$t_WWf(orVhijoeW`XsJR(%X;kTpBNhLiT|*JHB>>_B?LGl6TX z|KzNmaO7D03DmNuf%fek)2SQYtUJOy0-AdEZ7?aZhK(S#fEpOQHUg1op*LGcVT49( z|Cz$24EB1~-p}^0{ubOOml%1fh(LmzQ_O zP4r-1mT}Wgt-`fLpSwFF8Dy2hg;Y0z_R2#%NBx4yClK0x_yDa$sd3+rsbtumV677a z(f`APX3L9c1l48k!5|{~j(YE4p+zSE%d43(Bo3{8Ho)LL(u@&y=~-1*zf^N-sbGT3 z;}}gP(i)-;Ixp#~g_gtAHf}(*^;c}|jsA7-Dxt=`K7;yBca7inrvM5S4hR*J`n)gN z1*Mw8Sj+mB9rWAVle0~##%JqaS+Y<6r}x(F5fM4CocCtuO+18vEwRI{|m}gfImXV zQ9Mz6HsY`+dmKP3oWxdcZi5ahpr8ty5D(e}8w-nl^-Bkb0x+I2WcQ|f>`z_;kty&= zXV>W8!Y2g_yMoC$j0*@b%PFMJJ{Q=u%)8wYqk67Qw3!-5STnS>I@*gXFkM^dF zX=N!2-m6BQ+8~wV-s+1Nz`Uudz0t|Vm}PeIW{!c8k))>8U?KZ>xGBOS0789;bb(A# zqIWheJ*xxQM%wclOePQ-#UbZ5sBEos|->ycBya>@dq8Og;sv-J9a}Jv^f8P+K!6?n$@INplF3aOs##( z{8NIh@-N1~*TnezimnT!Ck(@rNO{@;*%Gx${immY3t+u{+dRqiIhB|6nIed*2Qa>+ zfYEh@djXaVJIzM1-M*b4A_AdmzL|FvY(F0y?Cl{*9PYGwEP$A1yJ$QIZW~I9w4Z>Z zN`BF}j!bSo98L$81lH0&j(|%&V>IWP!sqZl_U@zsX`hLuSE<5>gQJIz)!ejv)?(*H zwx2)jv+b1OO-)F62m+b@Pmz0|c_2Ojq|%Z`mHc!*gvf+Mr*LcbuwOWQ-Bt2a$1jmS zT6v2=O`-&#N+L|XFDdLtjP%frLaFxFEc~8Pz#%@{L#iS3=P*QV^#xQ4rfHO>~5Nbd#k=DXtzou;)r+wW&=9`3HB|nSY=*ruw+N8zB zMIS>6jY9JexI?f#`QTu0@WwfxyUG^}?iICOqzV9_Z?Ja=aIA0!UNR(z4ha0d1IC8(LZ3$8Igr!Vzr;n9srM&DWJz z-%6&$M3rk+BOC(!bkIPA%PQX-@WFu^)V5g;v8sU1$k4Dm?3OmrUb;W=03IizCjlPy z=(B8oZnXE-ZCG-Yq$L-*LGX12jliSj1eu)Hn|UpUiy%e;-+K__a8PsY@XkcTOkMF@ zyk8eorZ{K76G%O}v9Q!D{ppXd3I*2Fpi`b3O%NV?^*8W~FA z-&`5BYR=KCCE4UXK6VPSi~luIk%`?OC$D`RP|JlN)EsuN9=i-YBd|9;W809RGqp;2 zCU8)=;R9%SEj54Q{WvO80FBk|-@E3#?xH!WQ0MpKfw1&-e$imqlkz% zxFvfh$WI&TItOeMB?u8X!8bDi#oRZrEQ^)*Itcz`H=esT@LsFOP-baYHp<6)0>*BJL0X~`y|K3jt1A5+CVKw6qmZY}7 z8VgDX!OEkq<$7@avC!nPBdnJ{|7v2cF~$xM(Kok`cbq~b2H*~Gk(crzZo%R-gVIbb z>boJ}IRrya7C;ckA&4`NX$}Eb1*TtswC5&VS!Aj8un+U)OCk)R4v;)Clib&TE_1f$ zXB;K+yg^1N@l`Sq(gvq=5pP(22>c^LU+eHBq~vY5D~&3-nmTjR+tlk+9v@%r>;|mh0B9dspwb_J zPTK=wnL_eNZtmUTSCjtK;EK=OI{6%I6m6v7i2bSV3`p*RvheUZGE&A%xo<;HUh$rC zd>`N$P!MyK7KK16*hXNHxr-|~M=M|qE8Mnp7}39J~@cHY1r@>Lf^-R{fkIgZ>(jF7DxL}61%%qKv(A* zHc{fB#fJr8drDW;aZ`p}3#zs%$f-yHed{-@FAgr-sRn~D84`Y!CSV?0=*f(LGRpxb z0gRnbnnKt+VFK}8in6K7=_)I@5S>mX#c~`ac5nPEU>p>>VtCZ#eHF3sBe3nx$d40jy9rWnsuEe@ zj6tXKP?41mSa(<#wm}$K08u$73r}0> z14F~bSV}XOi&%J#KR*-?enzv59bWAXzVZKqBpiQbvDzaa1F~p3r2HN{sAUmYCVEE1 zD@fm2JZuRP#KR@M)7QwCTE?!KL$YTiy?ntXLZh8B0{t#@5HRNh#~rW08lWdF-{lM1c|9Snw&^vJw-6|vOO z?B2oWV|r3-n5oTI;-Z%%%8aEXBt|1|!{(L?ZlVz8MhBd70EU)hc@Iz`86FRe%OlNEI7(7D+XD3)Xx&VyBaD#|>adkFZD2K8PCWMf z{b1}UlNTr$5Q$jc63HkI-W|sFlJDvu*=c(w?rpwb4UzD;g>q^TcUD+K>dHwq1Xw^9) zqORwcF{3N{WmgVEg`Uy|YiemNH3$)@pEwHNsiv6=kv82p8fUF>M;c)azcx*q^ZM!< z?6*>L2f8l9iJ5yyb#LB(*9}nVS#K^V!gJ|X28%dmA_08n@r)A(frH4?z-D)@ z_D==k|5mtpGd06}Xylgn{XT^Uy1KeQO+@pIS_>UUtKXgBHt)NGU+9)boDX1^04ti) zlMMvt*|qfQcntTr9Fr1(vAA;d$e^_6#K6Vg`3rEAA>Vv_?C0{uZ+~nPem$-Zp+@Rp zpdZwF5snf(wPg5cs2EJ480o=7Rgw)J$DS*yc`h@MfDorFb>un_;D9V12_+trVn80C zlIhFOL_dEOEo{1f@5Q!WIN|sL!#m#c%nZ_{`U%JP(dj$WcyI0WjNjX5ZE4op9-Ifh z?r9I6y=VxLFna|zfm5IscFqSuG-H;|eOnh0=K($)9D8KFV5NbWq#R}Q&r>1yGm=Pw z=6&w5#tV#>PCrg7XY_;PA3OnpGf1$JyH<1R9lPbr+*>A}{!bF_=Kr50T>1LU^z=M( z;s`7-4Cvj(yjFy_MR-6ckh-7^Yw_Qha`mdMBV|nt7f2Lb}&9IPlmmLF3vD3=LCtD%}HG!h(ap!8R2XVDX#+UqKKl zCsuAhTJ8_S12SE|TQ^INeP||>Y{t$&4>jI-F4h(7OA3dcf{sdbq-Ae`MD*F$jc3NO zdt@Ieqlwyvkez8=%B-`o@$*jQxo3-SZ|kC{XbGn@f{s2@I1U=5ly>&|I!|pR@&(Ug z&>$Ko90dpcCmi!4Pa&TvfV>-Y`x%FYyuj4qa&sj)w(&ec}zti8q|}gc>sv0(ib0bQwA9^VRy_QJq*Pm^}8cy*PiT{XCO! zHdKJa6oHVeLIbQh^Ww61nYPt$8)jza+w2I1=K10w+{GFJ9;>b39+~xXPy2U0KU`1! zz0}K+{(>>#qs<=HGnb(e2URDRvZm$`X*bdku@|u5IzA64N0_Et)eR zj6kxz!ZRh($aqU$C94C3oO=#Cw=mRO4wmQ?hYRbxlxxx9K8aR6J-|<%;uPWjOddngW z%Nn^vYO=z&VM$Z#WQg&DNbC}{{RABF%4xlw+c+D{*H74s3oxJtn28B&@&#swx<*9X zA__`{XFQp(7rQv}1Y}nL?_T>-$hG1MtA zC@Od*&|zFLLOMEw0(ZmzO$`HVrZ9#_V-ug01bH}_zjUp;mLc{_aL&N6-<9|%mk=jc z`*AIN81aPyQbew&ufxPV^ZHfzY~9L_AL%vGT@sk#(;w$FMyWAVs>Y#1;hGq>h$JTD zZ&F6kvm19lro{POZ)BJ7MxoHtaa83w0c_NFs=OEow{`p1j8MN`7N=t*?@-!k@S8lI zCXVxXG@Uj~P2hpTAlx=wf|v<0aTEQ87~>H!wI%gMm?<=hJVyU5Kk}>~6J%({9S5mj ze8mV}*v+lQS-8p&{*r|JMVA22B5P|~T23RMnt*+%fe97*Z7=h9Unya+{u^=YP`_g? zS(l{OjvpmCzL+KaKi8ZAbxqK(CUO!j5~iIB_r>)cgYcYF2dN?d!wV3=`R`I#-$3fg zeIH7n;M*$YKATz9m+MAbEq>*w<`#gw`meh1nAIvM+tirIpT6 zUE6-3|7IFv>VNjc8@3ng1S8IO!a{CtHDW(S;*h!R`#}0F>AekQ#VGo^AfefQdA_LSc+AW5^^fPjZIBCsT9CD_ zqBEzEWm__>l31ZeJ%zcoQ8s8{7DyRbUnAd+Ty>WkZZmG`jhbdFs-}A2&8BQ3%z*}m zZY6v0Kq0~;d~L_|U)&SJ!Y-tL@aSQm!rqS0thqPD8Zj@is&)ZnxMCCOt&d#bEtU3j z6}k3rzZ<}@Ln8l)tv$ob?X5vPGK$2+@&W(0As6hMSow+m*2LoC;{QTM7le%1c1=`h z*+`$HMCj+w)=EgI^ZW>NMStCK(~eU{5kgtmzW_4@uQ&(_dKE(XZ+|wx{v^5xQ8qJC z0z%+FPV+7pz=$OKgxLlrSO2>iQl3?jt=?vbZZ*ph0xlDj+C=qzXSEJRUAE|$3>}hf z_I14nk3aMNc<~uSWYTB3`-)x(eJS!I-gl*+@La+&d`aoFh?kcxyMALz@1-PpPeG5l zatag96f4==gyq!siMNS!=}|a6P3`|sLTi`)sg2$1Q$96@L!#l9Ny!{5gP~!|VLR>_ zs9m?w^=Oe7ydHn1B1btW{o2_OKC~(p%oed}kO$#ETybGxVOM*+E8+0@4TJoMQRl@HS;TErtieUH+MQhE0mgIvVE7AErsDs5CXvrGL1!;NIbsnU{nmWIa|_)#oq zAlz0XLI-0fMcY6*+_e`zdG2Z(91y`bFb1ZSmI}@~xww=L3=TR#g={un{Z1WJ!Mea} zqseB~ZIq{uca4moBBXAnU91U-?ue5%OMjR{ctgiAF`1*GEiF&WVKee?flW2*%PrM~ z;F$dcgw+IMzWceuh|ua!@@S30!lEKvL9ccvICL`6!8h=Y&ftYQGLO&KuDo{xJ^*8vs^CN%PP31wBSL@4!TWn^TO5krZ~#*&qX z2D~=n*ANZ398hI{#F4S5=G1$O#bOcflhKlUe}!AAW%2llT0PE%k=WL+55<-DtQ@h0 zz>h&wJOxyPtET&nVw$mtyd(0wIs`d)231cB_T{8MuW2KOS#7ID0C$^RE%#wlC`WIS zOav1e=`c8rP(c1;K{qs$32jxj>g~QaI69MY zErmKz@Q08zlfpG{9>R}9#*IPhj#3e()YMco6hAN6L*zrix=nv;=oy%;e|*;i26`9T zmIcMRQTb?pzghvl>`(O>#7H0W8xJ%S-cM(rfyEph0SASv%T1r}!@D0;yL*lqoPK}i zCBwpx$B(b>y}C%M^8F_Az&RK0)0)nd_8@fXv4_YNzQ5+qx9$c68+Ps2Ex&+kqVBFe zU*TZsT)X78%mHV;oM%0FS@$(^|10;B<^(~!!-}kjwHjn0Vg*Ah=WI)*+ z0CT2nC`?OByW1Ia__;*;`*W`szrT*pYgLYBBh=_(U(PO^QqD%`PD2?n*4G&4Ykawn zi8-?Vy-cPnC;T`HXt8Ej(N}+-L$^ZUzAv=Es?r}lVrFdIHdJDFueWg;(H6(^NZ<(A&n5*mLo`aC@$ZC0kE$2#(68o{fO z7jU^d4>{6?bB1y+dQ^BRztguxIUqW-K+>N&5zMKoLM(FJC@&uk&D?H_!^|v3WGrG| zm|iTt4QFU@&Y)HmJbTm#LQLe22m}mZ;~LjX6_cJcV%ur%?RV1cf*fLI8Yv)R;<_#19vYo zih2Gm`O^mHL-1hFK^Hq>{_Y!%JS;73Fq-zjLb|o<7OU_GLyoz)^lM~! zLdy4Xrby`NnPyyG0bvvpO$(VPz;@NN@YxQ0l}$Z1WWPah$@nG<=EK|b2P2cB`+sVD zd)7VvZh;0l-_wbZ#T9Uvask;2aB#7^y*(}~3k4@4}?ro1;R^461H3%UtkoL;%>FOy#M{*x!8F) zkQxpA^3?7jG!BnQ?1-?guCC6eE-eeoQ6gqbMEzY{^D|bu&w2s*-Jg=0gO65pxow5= z=K?tQjQg>)qGs7DPXh~m#C2T;(IA{+z|5opvf+0iaK>|Itze)*MEm6~8Z%OA<$Us{ zrE^i`h`R2N){(ehz>PWpxo)%^18<~Q+dJ4@G!jj(t9xY8pDVHYN(eVvqn?h?0+ld` zkIw~>lmHp+CsGOa9U#zh)AlJ(c*->tyvoW^;2X|rFq@1Mi;Xu#ChJ(^iu;xzxUjH2 zO=2HR3u_WPsBauB!D68(A1#f|4Ar_>yD@T|-2c(N(Mo)#)jVkYzfQ{|{Jk)ktFWLt z+}(mDnqkb92D)llJU#(&v`D_Ul}kbDqgR3QZ^&GdC6W9`MuXh7 ze+&9rDthcL^0FVD$_h$RtwP18{LDCPSA0=c|VA zD2FqE%;l9J!DrxPT_?yXGLyHiR3K3itQgUw*Esbas?6m<@bm3C_Uh9em5O5JJ8aC6 z6Z)B-s;@Uvg)$J=6&4h*j&~3f$`4bo_!4H&J-oqR%Y+s{EXbv$cG&Y2xgLKHBTLvg z;bb{C7Z-qv(nWsfqO=pc3gV}MFZp=sOoB+=s!umA0-z&sQ3A`OMKN56r)>h9{AatZ)Mcyrwi*f5IBlh)+&OtaBT zv>BK5mhMJicfiM!&NjDi=gQQ}-vo=i}1I}Dc@G*1( zIPHR$PvncBjf|4P$T@EL(jf3<_bG3pBqC z3~t)i?y(#jq*@ccM%7e|-x4rty3*m3e2@yW?2`QmyGe}wY=nIi?pkGzxozd&Kk%h! z<=%%F38Cn5X)6nP#Yxs zo%LB;0aH!T;3b1a}Gx?hqytkDO+Q5pag>mRB7D>wI44 zFJ348Vjr74cr@3ri3<*6{KgMjRcdW& zY9lsV#2$ehc}9+Rj1?pAohTf%?1Te+d-lS}?7@o}%ns%xMcH;k+Y!wjY>y!22rGT* z=@Rp}S9g;xury%O7^%ZO)v2}puEDG9>hF$1J(~(HEft|z+LuFdekklsyK?>dbAbL* zC!G9DOieqGH*g@fV&we`h$jMx{HB(c$ZPc+ez7LM+SQqGL8=?55ua^_SFHM5&C!f* zduJ3l!&?yYEe3fU932rWGV+}(Mv-xKb#=M3Cl~HJC;Mr{Vo>oo=2!$8{9lZ{2{@GP z`#(IC$|RM>`k-h=hODKcELmbKHO9WnT9&eA%hG~EVzNZ`X6$=NBvi7Llq@BZ?1U%@ zE&ubHp6~DXJ>KJekN5u^9X-dPanF6<*L9uOd7hu;3UetWeOv;nkHaI7XjhAQzq~C2 z*(^{Tg9U485c?f2M5lox->;=R=VD|Ck%~@t`ex`)``!sf5EM_W`+>~1uyboxNTPWc;oX3k z`;gc&yUj%Bts@MPUzRo`*D<@5tsz{1=ob7~#JtbJJ^ChzrG)S%Mr4z7>B+C0+#vl0Hpu1A zFFC1W)+6*VE%l}Kq6$X4Ta65LFbJ5YyffD=`x;GK1CLrlQ zq7I?(;u{LH0&f6j)eMp${(=jH5by4tL~bXEz%}n_{lqdGo<~i1qydMZ+hE z--OYOw1bpWg9G&pkA%K(Pm=vVI?&LZKv-|&1sIV zE_g*aHOuH-PnAeACszS1I|BjS>QxqWoW8yL2|XHD*Vj!N+Tr3E z%t;x{m%L067Gb;Plh)3#<`&O2r+wP_uTC?1?QPKZwDpk5of*;6rM#t&Gv@X;!{xmi z2(k9|zBhYB%Dw|o`f{Z^|1@ECY=Pxhsu}&ToO+oGT{!CF5$g86X!-M8E)=Z&2^3!+ zfdi0~Ba}4%`>730p0D%mWEz=wT71X#Y5f)~yI)=1?#VCjBhGjWm-lL5Htn;=ZF$Nm z`PE1Nar{M*>L-Lw^CDYmQP=0&8C_$xZRoX1!H4#7veTJyQpG4WmLdh5#GAQ}udK5& z7BtQX?%5hRO;TTj2Wve!yY7Qt8}ZVR;0Benv%l9D%KT4aHe|-fvuci;ke_Ln$h4nL z2D>SIpIG(L*2^hbOA!)hr(Y3Lbia>xqlaR%{ZDD=^IC&SOG@Gb9%aoZ`LYKSe4uHO z_ig}x0w7;L`V_##=Eh3E=_WQL%2Dw9H3qG*MBk(~f)C(ykbW7N6X9}U%?_{li7qZK z9d%03>9}FcoMnwqh*S1Gl7H|7#gO`3hGx3|M87b66cey5wP_sUF?`f#ip=<^`rt|> z&Iw-~_~VAwRR$WP63eu^*2m9p7QB=0eE_aj@j5J-;@44_p6L7>ufduE+{ZS3!w<1r ztD%U97xns$FrucoH~Zu)&pjCDQbf(cXW;uWqhqs1dDqxGG?c!gv}}PvIV;MQqc|O= zTTlu~*<(~h;2=IS5V7ha2;O|oV=!IGctgy=l~2Y&WK{u$NOVY`*FUqq8C76ZG1240 zBz=&smVEK3)zX}8`&Eh7t4~gcs|xiGll&yn;A0B5-^BeI+owETRAAL&SCRVbP}tfn zDWB`&rD4q4GyblVA5QC-@7WPITDx1Xvae(EtC2GUAlOu;wAGSJ(lQ^u=6mlq`P9_Z zgHZRNrZ1y)Ip@d3#-^wS`DMzNmNwA>hh|Wf#K@p5ICd-jCO2cc4y@54bOsQPkLJwe zZ?uw7JCxef-JKID7?@wY`VSvp4=gH35+HI`Py2MFi^nh(k0cIazU6&e z6qXto7)ZF*Lh~?PTJ=(H&Qd>$f|xCN=Bt3@>WkYlYT(x6BN4)f5rZbM<$H}tJWFvR z(YTAir4t+zZ*F1G%m$wiJDx29`cM>Wjg_}yICXY&c5v6J6IirUUYHX(;%m+&9-GfV z&rBXZLfMFMoEefF(26R0T29Y1-S6AeMi)hXNMUOx%L0;MlGgU%Fc8R2jTRDz&oisY z(Ptmfj&f@d%KELe+fv&h~VpACZ^{ddlr1@Ng3_-Bo zWSYSzm@$A@iIhU6_RxpSMioNu@3WZvZOM0)Lr98ucNH1_=q@!hxtsIVT%6Zn8jD`e zw*#>_CKZC^g=|yDNC~qtu7C>{nP(WPn+-dm-E8=nu)}k!q2&hdScd832Jjw$ACA}F ze1CdR{}a`}VSN(mIjr) zOT$4-+iASU4Si`4LB_p>odsr*v2Vqbai9|gf;ArWcY=#P-bRF`*Et$UnYawDaA%=B zMVTY=(CS!N5$C_)UXI1|I>(IRhVwYejXmiX?%lS}%s*BYezI#jT+_CxV3j=?c-DH* z-SU`NL;3_Lf4rZi2TY{Q%+E_6>AH1FsQ?L+0wLvHb4bUK$*x7No)(&dvD=rqu!;4I zz$9nH zq2M`2YSCULtApj-ECFdUf0K8rtL_g}5fcyVLD8-F?UmXZlDHwuDe$;d z>-)+nZfqFZj3LL2YF$z>%7tdS76Z3wQZWMT0IHSOzzvh_hyLwrX!H5B6OzYbZnW5$ z^&=QZ&ToFOI0zAH9|nAt?S~o`VH-scT`(Ryu|%%I;x_lkq`7m`*BQZ>#~ z@?wjN)Z4hHdHLzBmf zP;dXemam7Qin%-rTC!B%V8F}#f8WT+W{22zl>;j#N~#p8zS81*Uq_3YIQfdeccWwa zig6Fz52o<0$Wco^zSuvXA>dN|UO-Y=1?{WiAf(NB9dMIhp1S(mFLcF^vv>!l;}%)< z_vJ*>C~0aA^_YX`^s5$gi?s!Gohfd3ut33IFn^twVpm6Z|GCYW8GvtbrhNtJfC6l? z6Ki0y;uT)BF|1nyNBrZsS6*N9S*@P}Y$yvWyBH7TTL`rPEJfL281EFn7s$H6(JIrY zo1`)Md1(3!bm|;S*)bj5MuxG;?qUH)yl{s-p=-v#fNaA6UGM@Wo#UfM7Y-kgmFwf>XhTSaFuHZl07%Rc_7LPf>YyR>Z3R!v&a-j`%wV zKB74w`OP1!CTFkj%XkH?7sT>LEOYsV7wJb2FIdja&8dZJch1-ugBQnI+9Dw<2}@^zC)ha2NM8G|Ks@0I*WQC(1l{;y@CNin~+Om2a7pI`;t%)y;o!P zf19_-vHbGFe+!-%-n*#a0T|_(_{OjS)}tqvH&9{cuWIugY(=*(T3E`Q3_H|W5E%tR z`vkrhfZ-FvZ?fZTrgiS|aY~_aRZEg%%?`kRZ^rAGg{l3EMy~ z(#Zfg_WbM7g~mURkf*7yL1O;CPt?r@hU#FENCam^MTKSoLE%8lR2T9 zS^|FROq%ygBGGMnE?!gcH2^uiVx?Hj%mOqnEmqg<{a0rLPU7C&vjE+ z-A77K67@b@KxM8M?;#KUd~5K_Dh@)9&=Z|^OqMJ`LemY7eQCZ60C1I^3C587bzfKB zmld<*={LA9T!+lv;On6$1Y`mSt6^9W;5qH|?kmzZfCAx#Wh2~?BQ$b0r0f4ar+$Lc z^JZhq1_s_S4^ieAn5d*IK-JLGP1h06fL&|8t?mOTQd+#3Rjg#6`EqM3_*@rUe@Y;X zkJRcn`Tjv;fY}bedhtv1=9C4%1WG#6z~i;8u9Vh!5Yy8@Zu8Yz`$1OuFkAOj%oK5Y z!$#rbKEA$=k3vVb+T;&BRNOkJ1a_fDcepV#v7?3!ecHMdX|7wt3uQAXg=-=k-jmcg z`y_Hj&v^zE(PhE1Qkb3_`FTE#`YC9 zl2ccq((w>i<%uOty80t$h)k~T8F}zcm?!|8yujPl=@E#C$E&VOwnAq%#dqkQ-a!nT zlhCe2Fi&38Swm+A^!}9Z7r`*+<%~91$c`1; z>tLjB^YXUnQ64#nsa)rr{_J)G$4>g~?^67;qp^WLd6bc{sIw}Y(^Ac8d zhG1Cg*6xjrCy-xE#Iyhq}ms_kp#gIy6Ds z1C?C+_6LK<#oU|f4Jr*J(d?|tDXy-AH+eL0`Ulfyd~ny*qqg@gzCu;vZP{o>3|CX{ zWKKq+D3D2`aM1S)zz!5{s1%92Na+Wlsjr`35p?^ZeR#TA9!f7ttoPy1N?Uv;2e+YHfw!^ zKS}09vu$~kMXQ-BC+j9he$8Cde2!>R{I#(o6j8h<4#9tFDiCYW!h-9)vRz4(gc@={ z00G-eNd_Aq@t2_MZQR+}2{Hs9iR9GR&__g^{4T|nwBqWSsVOobbWShR3ZN~Hh(;(< zlq3^LjHDTF;mzPnj@wG-s5cZRtDhD2LfNsw|QZQV!-rVptm$I;5GKU^rMwGuce zuC6YJQ2{D_t=H&dw}FmMG(hhz4KX?YzbvTysker76v2_$?4{-8EQZC__BnUmQ_f`X z&$}}(QlkJ)1u9?cXxrr1@nFH^j-yW<++FTRtD_n^?0Ct9oaLrMl;p!b`+g{=-~ayo zf)%uA1|V)J7w_+fj1o5?b*L#Wj@_y|Tf43%I$obbIt59LBuZku=aU6$zTvJjVz(`M zc11uEp*$c<6eafZtsw>#z$fcLyqp1G+T}Hn2QnKAkd0)(2$ms=5n@sSW#5^o9QZd?M(cSH9l~&k^ z+A%#hLRc*Uq5`}C7#8s>t6++*qTx z5s+H8&h;KMllFW*=_+MJ`@AvLb|KxP&lA)$!6eeHFUeed$Q{+IEDm|u1C(I+(!qzf&Gln0Msw(D5$ z!;Fo0XSbDv#=C~zeKZw@YK$H-w03~93FR3;yBWiCR=XU*lSl!8Z=*h=0*HjLjX0M= zoom*;-q0HQV9lLqQZxjey>q5n%K za^>2@lpiKjV2jtU{aKWFO@Ky!<**N&@)k+UMm#phB>&jJE{RjkR@F$T9qtHQK5mrZc}jk$Q^@F^7S<09V= zR~(@jW)@l>^edkUzM9`7Al1+cCORe$b5WBo5d5eJ^b4EN%IEB2#b|`98r=~SDsc9> z#cuWZbbp}~yaSb?y&j8G2?i%&_Ek^0FDyc4Tzp0m`Rs!e%A*&{U|#s-B?l(Bi-(v} zTw6p|r7Cv?yOuzQmXR$=vc#}SF3!aTn0_^8K7@cX&W}P13BouH+L?O z>n2|qNg(3GYC5op0W5PWJq2C3!jjERADta=he@lCm8x0r`-D1OlbIW4rgPPjE@-9* z$Xqg2sknOoaealU?sfIccBsd|;{ii-8Gtv)js5&4KbsRLkV{tIGK4L|YF{Q534|!< zg^PhPLSjU4+OI`|!BO+>(>jLw8XNQpw=`bjmZ_4tt|fDKLwCXTvew>?P%c04o4g7^~0ZV?d48>fB_4crAV)oU` zDbHShP{+l`$J_H0ZjriyIc`t+ypxMItE#RXUI4Fay z^78t6B91OW^V34c8laebyWH?E=Q?E;eD%*jJa8;lm9fvcCi{0`3K8Q;1&Obi}M zb6nt^=a^{WXnJIf3>Ne+;i?+L(K=$THr5LMr%jMfU~T6mo$Du+e|;;ytI=nUe3%2a zPL%>X1a^$JI@pQ(Ip^GbzD?;{fgPt=rdQn#CKVfP?iUhrK5*qy|M`p0fVjX7gd0lA z;!7DuK{Vxz>y!X4)$4yrGdfgcr3OH}3mOct^8l~cLXeJX02S110l6Q9gigzChkOB0 zsjSP|qVdc7gB$ZyPwSk5UduBcmMg8p+wva+EK*ZrC}V}4fw9QCX{6`@{ZflUMaQB4 z@LQnHp9X?Z=U^yghv2p(H6HBzZ=WMnVkc%zFSQpC4nl#w1Yk%nb9)K=<&> z?B0pLmb#xH&E{c;~!={1xQc+DWGSoJ^)hPz%)mJ$2m0fDDXDI)_OizfAh46 zcQ&TZ{^r?!^bGnHFHj3&HWB(2_kU&_2(=rDkbCmZzKK%Gm|UuJ-Tt@vk+mItZswNG zo=^E6`JW05UhkL|hAvOo0LXtZJ^PYXo@zyjah862$W6svNExVPsJ4R95vAnJ7|)?` zJLnEk6_w7Nj>fleu!cXeK7khDI6C{Zr;OYd97Ac2?%$v~UGo|^+n_8>`h?U@QyUmJ zK7NFO^e^jO(KuRkbg0 zy9KT1zm4}c5hi_nd zuq;&x^ofg#N*QZAtE(T}M5#d*fs{?lv9VxF1%PZurfe59V#qq{Z-1&~s|NudALK{Vv>9jPL&?!{A%+ z$fAG9=eWNnCCd;tNNb@IIPQ?U`RmpajSOD2k}x0PG(9lUN*lcJ#_T2`9DB^$eczlc z>a}fR7`%5M`>evM6ZOwt&DIUK$0`bxR_QSjkdjy0RjEd5;Hr5u-j0uCntu4^7nqP)FFj~mFkz}O`g_t0rNvr~jq&zrjl_-VXbbHtef6;0 zrPKM3t_tInY_J!n6i-o15nIGlzB? z(6txL%vd-LRST>np7De~ISu(x3UmD_)k0~v_Mf$&^+La(hvouVMV(;$e6AU-Dzg~eCqn}7>l_UHX*QoPcE5t^*+MK$OoRdHlWFl-Fxe> zA<<`eS2FpjglvVjue_K-=p_b(IdzmVph=c4wY5h6o;N69%#FdgT!z3Erobh-w14b- zrGd4MEBw$f7X8rNi&Z0rM#JazpG=IS3`Iv>Cvbq?v$6;6Y*T&^Z!kU7$i2!ZD1+uL zgcd^Z{uSV?z&6a-RM>+8P^YmMAucKsyn8k7W<)5oDp{M7D!%!-a$Z+wz)-JlGFG=- zYe8T`qHRVHV_TR~W?VZ}hSFBOvDEpxedfi$`3ofcQcR^WHj>b2!*^)^`si1BMuiTn zjwVj@E|;4;3pwiW_*KQ2=kMmI?B=jbaKvN7W#YA5V!AZHpqB(D5b!}^e&Zw7zhCD7 z_PkFgXdlP*woG=i@XiG(^E(49PZj~$3VX=7sl815#Goyn5La)yH(ZS(Ee4eM#aFAM zo%)3vC}!P9tr-c77m7w?&6t=^BE(S3C4#~6D2cdCe41g~Cmrqha-4}J(=PHcL6bJP zv}VDfH=NPxG}}kAXtdP0(GB@NYZEQZFD9`Nr@zoAo!_<=tXHVaqyrTyR_?|jltPN9 z=ynjhf8_@{1||ZVv}_)Jt7!)#CZt0z+xF(OmwjU)fzhHl>SZ6)pWVJ(Z^cY49Jr>8j+jRUTILA?(1}YB+LYS=xN=>I5O7Z?~@3Ng(hE1G`cge zq@=U^01%L{0iuPj1Qx<-;KgH!&j=`5$Kkh?EyF8q^3>PZI&bJX+NR&*xS+kqqj)M> zC*H!Q%vHQVufSeQZ?D{vQTK0a=I^8E{0|SGni50+9{w2RNlU|Hcq^De z#xvyo_l~p#J5uW5vZ>KN;V$NL-a(>m@Ih%Dd!u&1Z)G`GC2>u8lEWreai(oTJg?(<}6o^ z!%?fFtuY*|;CF49&q4fBQ(aI_;L*MDwbB6IW)?eoN@+ulyGBPVM^}yjmAVaRH#^z1 zL(0wz*t?6}W|a6FP~wMCOO&A;le9QJ0Xq@>8y;&Wo=4rQ8m~2d8vI&Rlu>45D|+CB zNQPl;3{?_%k+nTD>f`>$Rz{vV{D%vWful@^`+wEOQ(95pID=zcI9&$2nTkC@zZyCC zl#j?yy+|UCJ>`^U0y+ZhTfE78wT%?n60~p&RS+A=YE|KI>r|Ok?rR;pw{(~m$cdzs zyi@ypus|Z~#ZDTm3!h`ztf)NEgknDan1;vQUTg2VT%E+|KQJW?H=yq%tFp2rD&*MT zz%7<@5^}D=I8WKd=mVev3JI)r%){LT?FF+%++CXA19T=ym;yrp7EybZV<~9VYFywKY$1tC9n)E|3_gBn5XU2KzT-*;0p?uaqm>kL?az>Bj~f6yBGu(Hy{z@mBb z84o;GxaEZ`e8PM40qrO)uQ}(Q5A%=6@FgK2Nr#otU(cmO0+^F(g(D+Y)0!=`%vb>jCPxs!}E$@dyWW(;l{q9+3I9FK?ATT3Q9 zHXB$8!!S{PqcV;e3k1Rm>-N*pV|=w$Tmf8ODGHv&Dv(dAC^9a?``=8%04O5#?ihco z*h?G}&Eu}xyK)lhT@B%#U@nU)ZSr9>@l?3>-r}tJ_Q-}SSyJN!<%m#A2~HYzWBv9w zW%M9iBJQE7c=!`)*S>B4B>h@krF7+QgZM87zBI7iWsKHFNK(g4i7p29=FEb*)gA8t zSRM4Q0B)W^wlP4XC0ouU1St{^fVe_RiPB&W}V)mXQtvT(MGWl6tx(y_7`;H_>F~mjgBM? zR8y|WQI6nON2xb756l(bkikaY+$_}~ztlIarR{rkX)szeUbN&3d;<>^SrR}HIioEq z)X1uRxK{i#v5$ZfosM3ZL1VnCe=Fd7#kmF;WRoZ?i25Ig?p58o*-1< zH!NIaC}X|%`Vp8_hKSgV$r{!;w!+f^Hek6mm%9!$N9*Vhsfw3tWH?scufKtIEq*lV zr~FngRbQt@-+YFKDYdpS))bH zkZ;f5YO&bV$S1eBWKnO??;pMREUG@Y6HB)gfF z`dZ`oj)>viazj3yK7K}i{kki<%Dd%*WsR(ieqG|*I?*!uHRQf1$MxFXw;XO1OxChJ zfng3iZr?_G-zQ+&SexBF$`c2|hXSMq0`x_~-x6B*a4L7M(l&W?2!Fe6ec0Me&QS%EV;ty&9dboN#7SV$lOG^w5Tj^nx( z3wAP<%6z)rExGcm;FCBC`z}*q1=GnEX|~CaEeGlTmM*T(1(@j@l9RZIQ#C0tkCK2B z)(J2oas}t(`?31SyquP=jvCOgnF)ZF3W?ZICkwQa`KO*$!qF}06dcs^w+Esna{B~= z1}8SW2afC9XAiNji0YOYJ(%FGEqEoCqkW8xKT!${FL@f%CY|v47a7Km0_ul^gS>W%yCg!E(ATbqH zZo>o!qIv@)z!WCWSD@0 z#P$R$p(-?KAHBoFfB2wYwZ`2g&AVbnVmb%K+y_NlIu_a>Os?(LiHfQM=R<^Nw!3xR zIb3G)mI`3^5oiVm*Ie~{SLs+2M&n1?cpx^~tv?K1O%846t)T&iqK4#mZtEeyz?%i@ zq-+{tLxZT-HBOnWS)ZNji~=m zpeVYvetv)d6oT{j*#`hT;Qk)bh=Ca`&~+4OOu&y$d-nnSSAPR;=PxS)T+7(V1L{l3 z9~IlAZaz+`a-bdp3MQM%MoA4*&_b4I&4c~j0XYg*9?0{)Z%;=$)&y-L(9P>Fl(?@0 zDeFS;Jo5B?e%l!j3An8{HDV|60Jb;XD@Be^_yTx~Zn2P2R zLIe(igp~Fg;Xq>CPBt-)>!OL;#I$|>tQ|w2xeO=|*T-GvwU1c`E&_*zyc0b}V@QMh zI*g`o@I08!*nBakn}U6cHK-5Ka$lPSn}VE|b7uGGrLpH(Fim@$c_{mH9rowbq79!O zq-mMT?yAr&a+#HlknoqObb5JkvW~ocCrOEvJ#Uf~l=(Hx6ta^h71iV%VPs~XQU-qi zrNHNUpmg84XbkjXGe8>?bv*`9LFy6D960{m^xFgIYL~miv4{=D6_7LTz;83SWNV4l zjv!V)CD`Bn#q;$HF*PxaC#oav^j0g{m1gkXG=#dw!s3p-D2K^{g0$0$QtN;eBYy!m zXN|XTrr5;n=hXe@VgwZjYUYgcawtBa+5jHHy}6@Edwq}(Ca5Nb#8R26a3@01w7oT< zAWk)k*autG=jQ00KwLb5HL@c`N(NNlKe`=&NW*P?@wLEG5QIB-tbtP=jeAxVZw^+YeDP=Q1d#~ASkuR4RFLwIr&9`Y7bIm(Pq}l6ESjREB#e{>{pZ9Kv{1Jd(6h2 zthE05!29oD;Si+|`s2$;%^dI}5*2eP4}sbc2s&c-b#%B)SLyKU*ZN<6W#|;a-X{ee zMT!=Q)8=BIRHL2OYwH0aADIMYQ%^B*d+n$QRiH%154$hgys5&@*W6zbOdgi1u1YY} z;=#OXObdDG6}UU`_{UGwccW57#SOT9AR4cBhxsPc6(0;pH#8+prQWL$pSk2>Xl3=OiWAlTW7R0(IgqdjfbZ zIJQmRyVKasvDLR?TT7Bpn`Z00(8z+8_O(AhyZ#^|g_t56?s`)T9dNZ|-;gdzm3KOU zYfnH<>;vo66&!C;f`P7X;_6zIm-!Va>s)Cq@Xnxb%Zat>Ucbl5S))bq_HO>(8Df^L zQ`P#iTxfMPTq@nF8RSS}@Q3e4A)T%fuZ3UNBkq{;4Y8i~5J{;rcXsv#H$uZmFhj9H zCI`*-V8t_Z>h2|rjgoS<;OJQg*qOVpJ))C?d&p#i(}#qjii&tL1mRwJqtXidc&YQ|fqvS_qoN7ds7vX;8F?Y`H}rb6-R;_3Gn`|A%-940A>90o~aFTkkcZf3w-d7YjF zDb;>(F;Xxj`R$f<5TzTQJss`i5fD2%X_2j?ZXD#Y4<2myn}fbSu*~0jeNm-|8hRwL z!%#)S+XctZw>8lcf$zWm0A=69&oc4>pRVsiuJiAcq|$*NzXrj%0C*nY`?>)`?quVw zDW^~V?v`}X8Y$m=_n!BzyR}mub0-xL9_t7;;?r?z>=o{CUPyiZSLM~QUGKkxd3J}O zff71X2PV%zQGPezD|5@ZJY?=gG3n+{qE_0xlNDRd66m z&P$|_cMk$?eMU@C9zcZ%?P@Sr?n*yu0YRZZ3DV^RJL<`KyOR%%0ym zW-`BohaM9DC7qqQC(~0`*I-1!gdus>TNxtp!~UfXAd|o}L;}gT-+}6=+z0IHsntn{ zyvNwoP6*Srh2zO$i5>5m>)y4{sZ?fagw|5Ki{6~6!n1<~aK0oe3Zh>1 z^dzj-0pqT}I8hbpmEV8Y*Ay6kSBoz)wFYVs$g>Zup8(4mWZJ}SVO`TK+FbGSKJu99%LOP z(DlB;J`#wI>F+YYIRj)K6YaMV16SM)kN>aGA=@?EUHH#h2yAfM7|1ZI5a3o zm!Gt@JwhfTN7munOh!VdG#fMuTvNcKoZY^yY-pO>wr-x%96HJ3cPqC5_zEBv2RVRN2K}V09uZgh7?)=zGi@ufqpjWt3my6}UP%oE-UA(c zKb%WPp)N})qr=I%;G{I`16GI!h+M~Nf`QCeJu|a3X^EB|5HjrF-$&dtIUD4q{Qf&+ zIm>X6L}9+osR4L|DVG2o1ciOMw^gx3GY|b!=R#IF4PLCOn2g3<40+ejwMy8#UxQz4 z(n(4lj-z;e#?al}eQEvIX6FX{En|O5)R~^towu|TFDL=~KBkGj5k7(J;Q;j_kdO`u z+0X%RA|F)K=!T%EaJQJ*OY)w0p`X>$827E*Yh2s?bA=?32TE!rHCj*`yy<8}3{{o8 z@u#-v$nPtuceY;GC5xw1DLrx2-0HjZck*uWUQ#64dgq3xqn<3=)G8>m5BVQv{q^+R zB#?wl+%v>91v9d0&ITO-0q+L6yEQycLR<)NHv8_zhA`q|^$?Dx{(;KGL?X%ZfN2!F z1>QvC%r=sYfi_7yh4Qq;WuHO#)_HA8Eoh6ro*vSnLflo~VyTVj=m;YR_9r5Vei%iW zGGY!Q!dYHsW`-54zCw#qRP$WMiqWs!Q$Vf5qF9~XKEFMql4QslbZh&fr}2P*}pGEGUATr+ZH)p zlBdA_hwp{kJq1q&CT*R|VZvO&=tRL<+jPa2jIF<}ij`)A;E6}UlT!YauAPg|q~#aC zlO>tyB;up#P&0A45_;4k%_^r8E>PPPnhL{GO}Cv#!^S|BE@!#)DyYttU8A%$r#?Tq>o zokrG;+DxYn|Bkv!S9z=;vV$Zx{dVI3)<&59c}yCg#?cl5(i_7)JIKf8rMaISq^mTE zaedx_X+KzE&&=e?wF7DnSg7dn)KM0>wlyV!v*xLQb_MMyVcE?l{nitfM&!%dp=2F{fP+$VZ6exak%MRLe9SQ8#c<6PhWd3`PI}@Yj)B|oq-gwAw z^HHCjQPt0}Cf$*%1VXX2$Y<#mLeO#jKvn}$9JOl_{mpwbd6@Yq-*WL?)R914R}33R zsZ;*gNV1L0cvIoLbd{^@*Hu6!-2+LX?hut5yfu#1dMsj$3ZX6+KyP1M#2H@od^d5%p=8yF#BT z*Ql)fb!5LdwfNoslLyQzIaBNHw3Z&uM4hdy+JbrIjCwJcO(~YyI?yELG}JD5@rVj# z78vzwN)ia1<{x)m2uRmj5ifl_SFk;#yyz}`1=yRSp-T{`!zZ*|m5N}aCx}V8#TP3` zvPvvob~yqEE1g0FcyRU2g1luR`Pe9K==R(k zyQTwccrh4&p81_6BG9&xxv!#SlLEA1z^NGCg(K<2PndxYJ8#MBs%hui1vrQCL}UZ7 zsxF+p;rioyo@va9H`o{Uje0^yG+Z_N=K)vUH&}9%J0V?3>Qg|w4MrsqwZ}0k?Nn{C ztgxJ&#JISNL?YGsar9Z!g>CY-0^Cq89HzV9c>L9ZjTBX|QYX}yV6?Y(&zpmsH(NQk zY?K|*IVMngJ6Fzt-m|tp_pF4nun$vMDgw_iVN?yrvSB+uki&(GDcU?rShk)pO>r@} z^EHwO_O+HcT*BwB-VCmx_Z;vA62Zz9*Nv(Q{`OI?6_o;4H!ci)F;Ix z?63djBNKH_XQ-Rv&!3CCu3Zp@dw+&Sf|xnxCzDMh3y!V|XZakr`WG$7gjM3OR$iU9 zLnb0S)UEDNKkuED(p`7KJTC%<=1FQyno{ma7F&!-aAHuk253}6n9h; z!X>DjH3OaK2;Ll?SLcO z6EZEx7>X2_{>FPKw%23dd`ka2%QSA0cq6ib`2^_BhkNeKm%5W#3Ag&(lGPvI3BzK@ z2c1xSv9*<9BV36P)=2V^+19T66f3DC%ccfsRJ`d3^+AM#QSV#;v5}VGQ}xEwM|u#5 zgFtfwa0!@hQF+Om6D@AF983Y~2q&YH9>zo{pfR7&dVOt7hm?Ig;l-`&+Yz9^xpdXn zt!cMWw`WMy0X6DyoTAU7j<)IfO#&i`W2sE`jDmBZVSkuwY9AGkQE9YioFBODEhaPe=#9AT`frZ5fzU7_&6PD7b(vqslo+8X*lh_i?|8!k?2bP$vv3Ts#< zCAKF+L@c#EuiqW+qqEIY_kiCQ99z<1Nes^zBiu=6(cYh$1V3wHorLdxC(aCsgIoq! zg=+jvHu;Uy;hiVcJ6HfQ~ezXmuXr69qs z;x_nKi-dU=mVc`9Y9|U`k5ok?zmXa6MeXlwM%{DyI*zFz&&=M6@%i^PJGGUEODL`+ zOYkI63YMJeu0Ao=yZ!AgV@iQd;I5{sBO8xKIGk5YpGdA0`Gy0R@4x#f^QACwn^(<_ zqfXP0k84LtT31(Z3S&ZKHOO`Cc(g?x$6CdC$?TG8HQS9WvB~~@(l9+c*{yz3IjC%t zBI^Hm-{?HZxgzE!@2_8!C5h=_(0MRT$6;HKiDPDkc>0N5FF;of!nDw4TA1=o=^g45 zHohkbc0y(bR(6$vT7IWKvTEWRpJoIk|?5$eZAxE#XQ0 zWT#4cM162htmn3ZshWu5>u36+4iYV1;zP^%B@&^ncl`&kx#s4_q|UHmtQjv8!~Ya# zpe6KPm-CHB#EFkx3LCGv+v>vZ729m^Q{#&EH!BP@~+>~_2vN$oqHTTf7l z*E%~CKNYnxtAjy=!btubatjF6$?o~}Ex&&2&Q52fiA+6B3?)*&%ig^xoK60HsE5i~ zm0)0dSn+-MD;5<_R0lw%+`$R^JDGf(UKw~G8CT9Z$sN>*~47K2T&G;eWYN zTJvbyIc*jT{mHHN?ABR+|IFx$&uBfZ7p{FR%HqppDByyk6GlTTQc@@!XMAMWBOciV zindTV3d7wJInPsmjjZTgzcMP#_QXYbNTB9+F7be%)DPZ}M%=aztf^|Hxd7clUq|+` z*pW`Uv$)H>dY*I8Xw->s7{bIM$8o6i@M;Ygs=6e7I3FYTEa#}is}!*HT}_Wy6)K>I&ud|@ ze@_4=hhwX=C23Raq==CKcPj1lnANK>rwcldUTQgsg<9RN`Wy@96e^(hUi2RndwWMW ziMxiDVtzwv?LSO=B7A&Y+_^pEyt1t&k98K}|8qIe|P-$Q;}CNiwjx#ZlTF5LGyCqv5VCFTk; zrH8?9`A@b$6mYv@>Gf9F?4BiuppOZ)c-Emhi+$4-Z9E6E9EGU(lYJz%QBhF@T9u|9 zza8Zz+&d0W(oQMtZ5f5CbAl>Rseae^OZQjK8|Dj6Cs@>A#I=8jCY?{F~9f2uBCHYpKXA0bI8&6S}%!Hm_5}JRUbnzZ#Fk$0M z-+%*MA?wFc@0c3?7?)Z0lYad5nX@IuOonXk*=P}HRo+Ke<%=J<{2+PpTUuu~cS-)X zHvz2rH58qt! z|5@3`LGUd+smJJ?6%=lp)2LcJX|D+JZsaxBZ;= zwTV6dY#|@m^o2o;!C>x$xN2CPG@*T-eQKw8!OXqjN^=%x-BYaQ?wWZ`!3!t6?tR%O zd{%4a`M%_#_`&@KIbcBSz_`Dgs+oQnL7^81Ux(W|O)@F!mw6Cxl%+J5ra>nH`6Di~S? z!3T-+aX76uD`)%@k6vEl&T{K_0@Spz=!0ppv|5^jZo@6!#-0hooR))f4sAOaRf`ZX z7J~BdQuVym)ro?Gw{tn)vhmt=*0hJLUhI4RL_tZVnpbwZ#&>g{)Vcg`N2HsOAP4$J zNUjO?Kc>?oa1*v+(vD71-}rLnifD8RzaiY4xK9$sEBJyof?hjq`O2f#iyey33AvBw z^3`OBf<4V=B^0sCBFy2gts4`LerDy1u(R$GTwJHC?_{#5XzY(T?Mt#_HE35*P(UJH zMeqr7IrDLH-K7~&2KbNjR)bIb!QvglsGeGPZDX4KqXUf{@q?n(>0W|P(xRd8Cw5jr zLr~MXKc-zj|GuF_aE0N%2DY|2DDf<&Q)R4QzF;c+JC0#d9Y(anL$L~OPXl}Xi4Q_z z)D{Qr9Q{uG8z3Ngg3O66(xY`aBUorzZ17VLZ^58y1C!Z3Z&j93$6YrRcJI2La5!32 z@O~0+!2hZ3yW^?u|G!U5S1!@8i?)iAkrbhn>kpE{}9|eCGT8+Rx{cM=&C3YQeBHDvxW9)w7BC=;`OKUBd0Wslhff38LRg$Lv=xdCbRPuC-Ga zvs_ePI!Aa*#frcmG4Bx2*}Z+5hUK19gJL3KxhHpL2$%)F(v_;| z6-eGC6Je;}#?BPvj8HGmYdyfLj!OIq73uDen-11VAJO=3k<~UFo^86tG}W4USD363 z!b?_{cuQGBP+eu6o{_|bj#~V||T1p{yEp&7gW-5ht7RR++L6 zGVqQWnI9jN`8G%!2d}pnu@ZagttD~ncfhaHnT+JyepDrE-nujb;I(0)>VjmbRq^1rk}#!HOv>^tBicc^7ru z$Mm+_EO6LjmMXA?%mG6&a{Q6B-r7UcZ~af!y0Q6}S@9?FI9ypb;&qmt`NM<$@T^wO zo>u;3K{*B6@+->mTA$k>eB9O)s2;HPteHCO*yMIl78Q_xq6gF-?cVrFhHd!zwajZANgt9 z^t&8AiwEZ>^o}&|^H2GOhQxI2^aE?qAIbGzS-$;CwY+A@j@`4Vcw7=w)#}0HrS%nY zrooRV`qk7+VBa`+52eVSIGw8R+r9Snxx`pT#qR>si~TbfcP)R&ruKfEe6xBheFoeK5KO0e=ITll-Jod$Hn{N)5mv$svyFBm%sS}A(J=bMb6)J6E?>ly{XL5`vLv9lRgy3s9VR}mggJzXM&XGn3Bo%^AtSsbcyp3`H{Lh@i==GW4G3w?za zQ&s)*y$3BSzDBq-V+0rf%?6p4t2?YG7zNt*wnH;S;h8Xdg8!ep_bw(*B5r(mOT&tr zKhDlGc8NFP3N13<*(h@^WrdmGh^fwxCl(EX>JJ_Mv@;bIYmj=jjr-i_Q@@^ym78`R zpO8UYEvaH#iPLUc+C%kcK*&!Vlz;cV@sKF1z=C?`-6Bu89}8vF+~Bw)GOI{eq5kIe zM3w4!JE9mBJB?rDgo>k#e6AuA>FR!9$rFIdoL@)iyT#t^evtEe*+jyF8=F+Qo zxeL}$>23Ocw07XUi)KXj$@7wz_ioVJykXf@t=4d-KK|qa$ig z-#bym`_F1U^Lb))J!wmbwRusPs)F>1aRBBz7R{z~x6%tF4@o~r_^LiecWk~9%-M9@ zdQ@{TY{~B{%18S6x9w~UJa;HTWLdwltcXHY(h{AXis23e*(5w%`F6MDpU-t^;~HFK z%~2#=`K7pboMu;1zUIQm&<+qk<;T-^_V>yGaXr1k-bjZJqHzorzE#Nr(*BP3lGO~@ zX`8d7Ud1})HKeqb|Gf4cM?G7gr1;xl_+`jiazNi?cV0%|!o<*%$_FIwJ4hV%Mk=sI zZS5NAW0tMq`)jicB>mkU>t^12T)v$nqO~u>qo~v2Y;nZC5XT8d(NL#67pKK+?oM&d zXz!`GqM3#jbT@N_*~`;NA+Ug1L4u)-IY&T-&gV1zx*J=kSKNIyd-grC5`VXF`u3%t zE?VlNcEG6aoQrCbi`IkJ&y$}mzkOdgQJ}O<+Lohi=T3UL1bb78q}Jx8uTw9td$_AR zV_p3;_z#>LHqBLW{JiMe4`T6C*6SxJvsZ4vJ-9LPNB@^i7~VX3;>Y`i*IPn&KdadO zXvu?5wdD`h?8md`&at;?82;J4G1k>vd$?lulAim9nY=LP_b$JTlvN<>3itKORhoBK z(bSIUF5@^ZZLINqR6by9ca3PLpJg$o%b(6}VG5hybn1!(wFeI*)(ebiUR|r$ zZjoh@Hxlh4xWI9O@JizAWUX3rfO-H5Mjj%|g7r3+n@wwpb)8^UXr9k0nG9a5)c(~_ z90sv8+qP6-_LnG$NR%#1sx0?LlYI4>h^yf;*Ea;;+O}a#Zv%JNyHJx$cM_a)jEuFl z$0y7Aa)VAkow_DY${?scZz{U6+58Ikju%o-ESuKf(NkMvA1aB#E1vg#;B2eD(6)JC z=+Ls+jZ4w>_cpA!!)bh4GQzv_qj=Jqsgju*>nmGlYn|w1?77t+az|~`ok6i>p#Z^4No1VJfMwvyT~C&kVIYC5j z3t(>Sd0fU>Wk(SHUQ%eq~{TUiIum75zjrx(MZHJ(HkgUB>7pQ7nZ>%+{R`wmt@ zNkP(#Gl9hBZSzMN!%zGWZ>HWmIJ-++1D;#>5-ocyA}Kn^6U2HZV_pq_;qLZ1lE;I; zxVUcK9_#;n=uJdD=gGdxx^{U_u6S5m)zSyEfgWhEiUMMMG@u4g*bQdHKHyRnl)cRw`DXE&bk3#Kincg z%2S_xxSyvGLjm~`OAWGPAOEw_ZTn7 z$RCIbXni2`ag&;Tx_dgL;xJ_{GsMJT8m4~4n>)RsJmr>2qB#oX4VINJ*k&w>`$0)B zW1Dyj!d5jjuZBaaZvT8gJ1{pw9|1a^G5!gH6;{O=X=x1Sf4t>%D>WhmlP_nFdn0`M zmL8a12J1+8vf1n=48_(yNyHzZ{U=Sw9(T%3KSRcK-M3*5=V0s&mM8T+u7ZUQwyx(i zud>CmXF98+-GWTX{)fC1YeWa5SA~izSi*;t+`sVfEsy;o``{weQzBESn_&JM)uv9$R@ae?Jm8?@D>OQixfvfT^Z|% zovi)afm43+&xXlBPX2Aa{~vaDHV39qg#Y6Re|s$7S3P116b>z8V`Fy=>@hLUf-L|{ z!Wj@wM(&_Sl2H_p(HUX0r438GUzqy6Ez|hf5m=NB(emJ?iv}gDaL^x(9Yhxaj0WQB zX!!96FCN$TY4?cP`vIyZ;6YkUR7dbT5d#63(~f@>b1Fy>cSeK#5)-|r69Wu~{Iv?GtrwZsSBNTba;*M` zkKo;%0d)e|R;T4Y@A>zqcDK!d;b`+%K&rlu!eK%P_G+DS;#AGdWKafVc;1d+5=?BT z;Bd0WCtPnXnYe14$F(~SGD3!UDhq-hWo9<}@tLF$)+ewF^(HU9zs-Q|)hS@`4D(wL z2@4Pa1fHc|v17n@#@GF%O8o1X{@yK{l@F1T=Uf_1ZuIe z?Q6?~npjOn;Kfct{0zx-gYHvtb48s8+}IkEh*BlN>SIh+$fqlm@S0iAD!YMWt~Y>O z#8+|=_dqOWw)I4Qtt6VzhTTLoQoQ*0X)^FLB_a3iMabHGnoXl9kZAN`qPzkX$7W2e z4}{+S`wEBIB#)DmYKLmd~YHFGZH6~TUW5KB~tgGIn^%q9R4ug5z@<%=)Eg>Hz`V&Iu zU$2=`bh(ouI3QmW*D?&|loZVT)9%<|I+Kv+E@Uy{V&IdowVC*AgrLlIz%d$|nvyVf zZ(#FV=KAvm8J_O_N+4n(a2B{w&ugRwJ`@?;&5jyvBDeWtWR-C*BJlgAt=qAXk_Jr_ELja%~R3&}Q zl9vU=9l8mrFou?;=y0d?&dtVQqrbs6PX^%LZc%!ORY>%dh}y>3vCofe3P;JS?gzo~ z+YtiffWBOHHpA?sE6h0&uv7-!PUV;)SotzF$WF|)!eI+KF8wJ*2j^)$HpBPdR&-3z zJ%WUh1@R*BfMP@&iqG4h+jaq;?<6|`RA$L=X&O(8lOdDXeg;%yZZ+Kq+Ww}b(Fu`R zVF|&i)UjRR%Gn{S5R_F>|i*daQ-(w;zVp9mshgjz!anbrA?awPq8(fh)ix) zb*2}(Ws@=V)^o3f+DYQ>N-?erGPf={~)1 z8cgW*jcx>(08J}P54?#M+aDZ6m6pWXd{e)IeBOF0fIk`<1EtKT2?qk7P-j4In`!k`8Ha)~P1!-T6UpLjXlY4> zFjYNs*^_a1kIEym#92%6(ix6HO4F`I@Le?9@1j!vJqCfJ-oWwD zl(vAG(r`pZQo`7wbU^{e9tj7j2G@eA^bqofs6-YQ79bw{^|d=5wFd$l^@8@!?#ji4 z+5Nr!)QkCaDVaR?>hw1;7jI6Zz^MAMo)>_P{(KV&Bd~5!Jz$08zTS}Dy<)}q5t+6L z_|gxw#bvKOi-3;9xRVzwU@BHWW=)-6EZX=pjK0?r?1LP4ek%w1eI&B9p;jFz;r1s*s%qjnnCbW>e+$Vw>a!3WpFYyf9g_ge}y% zKB*gS>bqi7oJ|2thwcqRMTPjE=Cv`}Nzd=V<>|UGEYc;Py1)&}Y&M1abg>vAgr0&E^^Dv&K}sJK5o|8!S0o z3{xuzEF5oNs=D$}_nJ6m_`4wvWl3Uq=z4d+PRJ``#63^nlx*ou=5!j>{rlac!!YYM z2&E#X2<01>&+jt!o9Qk~^W)+`U;CD3Y$rU}@k7X$4ZTo`Y8*U~6m~ zO$tlERE&-n!Z!54GuJ}qEgg2`t^zM5H^6EMW}cDa!M5l0J?lwL0~BUBthjJs_vKrf z5T)oZ96AU|Px($Ihw6K>Iy0R!ldD2Q4}&yOhw$DtDI@hUhw3wFom}7n0#G_9J#Jde60X6WVsTgi0QG!L#!adZ_0`v1g#wbK+ z=xAPqO6z^Lisa*ACj91=6VuSY?KCb>f-+9@xe9?T8&X%F7Uw2RPfw}`A|CJ>eHOZ_ zA2>QqbX=$&-(Ck7AQpChC3N8oqC;AZ9nx|Hn8T2WI!@>@`w9tJ6o8$8R8egHtXZ9( zg#?gOP%*FB2L-Co+SF8(o5*ESCfomp&PhlqR4k)Fmtjkj1>+;^|P)sm=mlcM7gETCCD zCFE7O=phqnB}v!eX4biiM~eo!Ve{T!;}Od_UCvH;$;jnzhv2|lNiZKbL3j((yoXAI zn?ypnl^mPba zQ6{R{2$boTP6sh&3`*3^1@)`gZcek_e-%u;6+H$49S{|%giXtaN`1G5F@ViMRFUy2 zbC#8*8@_eeyB9Z-MddU4ozwXX{GP5K&6H>O?MA18fum6~4YMnn{1#Ehm(|3T_7T^K ztp$%g0dd}1l2;=y&A0?xEvUxw5WBwowJyiJav5DKjCfx%ur0bXwPjz&rAwTj-M$Q|j8$z%9eAFr#Gj=f{j-D2nuBL_fXT(j zzmjc`75Tv%QVERXtlt5>dnupisuEICIF|M5u|CsJkt3fNJy#RQq00H4vAL9822 zQe(Hu?OU+Q0`XOzKQB`6`A`N-Cy%nJD){JnL(=?2bRJX(Aq*^+tpHnsVE1i(R=oY; z`l12D-ah0`^@w+jV`u2?(@iw57VG4@GQyNxKdU7K?}8IUT~rT(jO$NTtv<4qb^AmU zP2^xy_{TPIer6uO>Q@NzRa#%k#IuqiJ39yDe5(Ij;>=t%2GvEoz8?_$h*A_-HGVL+ zbpm7gbQ(;y`q9q54_W)*!VD~*0sTpoo)rX^`?i!yr_fqm=v0*8W}zoYOd26KD!Nm2 zIM=2C#@Z&VeISoYxcD=w8Q7mCK!wsrfZ%hB3H-DAj}skl{!b*!2jTfBcTB!g*u!KN z4fOQX5xB937qX09LD91IG1LhH0Yai^;xGwg*#u!qeR_DwQ1n0w(!vfmip@NjPC)9w z2`x z1RDtijx@b!m5j}uwZUp7fdag9`SN)K59H34CX7RJ09qon+L3@wfe>4oCnr1ESfdW! zg8Ps8qLSvW4|l89$S3`O_udG@E~A)*v^&_3`pc2we2{)B%g%GWH}I^q~heyeN)q;C1K|8Zn{Hkfarzx z5M2cs(C-$GdV_kCsC7hqG?;_viV6m4_iw?fgwY)e+#9le2-`gPQmY+*R7a;lkT$Et z|1U*7?a^sq{@xRW>zVB~_kKKs5xJz|wYaPnj?8g7ZjLT5$w)4gOc*By__C)okpPzl zib~@GVBwfh%{3oELOLRPE$UHx{N4+n&IJ^YY zx;I{T;6|;|RJ~F+L=77jZH&_W`2pC+6DS%yjqx70ei6zp_|It*3L^e)AQ$UMoCMQv z4Vgh`&eehA(VYK4VnBeR0d>?H#b$!?hCWe-foLL;yg96|e;YRV=B^*%-mAELTV6>2 zv1h|CoS*0(5j92>cwxk~QTYXxTRn_>BI41*?yfB67F3~pJCA|nk$nBvmHGdwLNwdN z`)?^W{f{0hYqtBwaRSxvAVtjDD9%KI{?-PuiJQ_WX#z;TT3O%u6%WG!4h>D0Uc?N| zuo88hu|wU-u`P4MVq!fi1gyxcnQtmeC#)y;&sGK*xQa6eH!0qgqjN88_kDyDt67$R zs$-#&gL6N+5gLN1NA5{v9VR3iNnA$I6f!?Eddzm~5O{6P@lPnfw$V@p{~3~j>R#

djtp6y9e{@FfZjY`ZrhzNS5zxCe77!)x?w8C!P#=H7o*<03}k zVNuGtk8o87IJkOLF=hTH7fFuwE+~e08CbbTOLWMhB#aP5aSvA3gNPR(J#t3jzK5Xt zpbqfH9QH_QLi8d%zv1M?@@}xF*X9G(4(|o;MGSlBgjheQgOT(2+ob z=v@?D8yF$Rx)t=cay=;tQWHlY40~pKDmp+S_akcKb@1zHQ28WGtrwUAcws(SpdzF++yI+%7m}Z5 zgm$wSS`=qYh#v!pkyT8FAnuSih&uX#8$b(5lms^TF z^7@8^-M%d)EG7P@z6sUis4Pwp%U-p*4_Rg;lihcFqfQe@8*;p-+!2Kw2L*BEO32r9 z8bh{*1-*tzxH~z4j#?ohd8F7>u=J*y;in7w3m8uD&Geh3rbFUyFCrl01x_gB1fk|z zJI9Du+KOwxCux@iNwRw6NoT}-%8hXwb*99+Z97z0P&rETbWEIDHS*)akoVo zxu{xw{W9VdAX~Co9l=l-Ar3*BuN5HVa!j%WLqR3Le4n)$PaC=rP`YGVXdML5k2(S8BHH5#OLF1w zi)m>t{?U`aW}p8Dn$_Fz&*RemsJXblMfuW>PI#*sjAUk%pMHu7_-49HoyZiHOrlkm zQQoy^FTK~*(ek;_QKc87esmDU^({Ses2yZe(QDI0x2wj;0kbK?tH@M_u09s!k?u8J zFd>ZG0rYi;zqR0@5%s5y%{Xi#3pA6m*`(TFb0dt=L~je$Xm_?P4BsfLajd%lxKu`^ za}Q&Tfz{CHs}xN3t6_$qSZ49K&EdDj_-k8-jx)@%x69;dwCl)mQd5VUpBosC~qiV<$FtczvG-s9Xt;nXJzTvp71^W_?FWV-nD}3!&@z7 zB{ng*%qqUP6*&6uDrCw>jd?55(+orf*3Y3RUVj0wVkNTJ*ub0OVO0W&_WzDeUCEuqUe^gOcR(9R$+xe&bD^}PZF{gAa z_GFtwd97<@Wrb!$;m-IJlgwAGt;U}|eJVUggQ}EDyzXC5pFXwRzkt%ly85i9rlw=2 zra=(94v&n4Fr==Ij?R{?TUD}--Kear6iJ*07Xq l`6=H!n+YKF|Mrhlb9R?W>xbBX&?FyHS2?&pL;2*T{{_b98`}T? literal 0 HcmV?d00001 diff --git a/prob2plot.m b/prob2plot.m new file mode 100644 index 0000000..be9521c --- /dev/null +++ b/prob2plot.m @@ -0,0 +1,17 @@ +set(0, 'defaultAxesFontSize', 16) +set(0,'defaultTextFontSize',14) +set(0,'defaultLineLineWidth',3) +dt = 0.1; +t = [0:dt:3]; +y = zeros(size(t)); +y(1) = 1; +for i = 2:length(t) + dy = prob2_ode(t,y(i-1)); + y(i) = y(i-1)+dt*dy; +end +a = figure(1); +y_heun = Prob2ptC (@(t,y) prob2_ode(t,y), 0.1, 1, [0 3]); +plot(t,exp(-t),'k-',t,y,'Xg',t,y_heun,'Or'); +legend('Analytical','Euler','Heun','Location','Northeast') +xlabel('Time') +ylabel('Position') \ No newline at end of file