From ad5e7edd8b758b68c725961ffbe4e4338375670f Mon Sep 17 00:00:00 2001 From: nin13001 <32068635+nin13001@users.noreply.github.com> Date: Tue, 5 Dec 2017 23:44:00 -0500 Subject: [PATCH] problem 3 started --- README.md | 20 ++++++++++++++++++++ assets/README-4e4bb548.png | Bin 0 -> 32001 bytes prob3_ode.m | 5 +++++ prob3ptB.m | 11 +++++++++++ 4 files changed, 36 insertions(+) create mode 100644 assets/README-4e4bb548.png create mode 100644 prob3_ode.m create mode 100644 prob3ptB.m diff --git a/README.md b/README.md index a5c0450..c997aa6 100644 --- a/README.md +++ b/README.md @@ -47,3 +47,23 @@ y(1)=y0; end ``` ![](assets/README-23f54e97.png) + +# Problem 3 +## Part A +The analytical solution is y = cos(3t) + +## Part B +``` +dt=[.1 .001]; +figure(); +hold on +x=[0:dt(2):3]'; +y_n=zeros(length(x),2); +y_n(1,:)=[1,0]; +for i=2:length(x) + dy=prob3_ode(x,y_n(i-1,:)); + y_n(i,:)=y_n(i-1,:)+ dt(2)*dy; +end +plot(x,cos(3*x),'k-',x,y_n(:,1),'o'); +``` +![](assets/README-4e4bb548.png) diff --git a/assets/README-4e4bb548.png b/assets/README-4e4bb548.png new file mode 100644 index 0000000000000000000000000000000000000000..5dccf9e9d14aab5167e2cc81d870bd3ed81051bc GIT binary patch literal 32001 zcmeFZc{G*n+dq6EQ>dh5=B|t-Y_o_WN~VmNm3c~%OtFbX85&F*naVsxWGF*-#*NH# zsCLMlVVl|SanqVuA`;KaFFdFf*=f6uc+uE z$N^gfq0pkCf?q26N_Y{32f3(*Lv&L_U50YQt7dDY0dSKDC*B^S-d!x z1U9N8y*TeiZm*V;w{-dN|nEH`0|rz>}<1a=SK39PmNu|FB(pNMO?#* z+EVALPm1Mmh1*0Om)t!8RW&sN6LyT5k@3UQuR{)$ikqJf-)ekve0es%n7Jk;`AME3`5`)Q-e;VEg&|!}LNMrG%FW}0@#sHRa%0Hgk3fAg!GDEXar5H3 z9Qv~f-d!1Y{EVSZF$MfzS|*2djM&j5=s#=M{x%5IL3b(*xx@;80?BD$_YmYD&;LIB zzu5%Q>|c~+PyCi@q`9!W%VA-^=H~KuR_`x3WO6~URRM_}hxz=;d`Qu60?Qd5u z;SMNnwQB!(thjq~R(!D7Qjzd`S~`Fr{pXZ*?2APDhtCv$JnA*(X|(N;>bF?6W={?+ z_4s3QUsF&YYM8#TQN_Z=rRBfK?)H)K)TvXsod@7~Rx>eEZfm$RfZ1}pC;mr$5XDUA zc&)!*|J-C-e6wx5xzyx_*ZungzE|Pp`s*PW&)+Slq8n28XHuI+-RjI+RyP(#jfmA1 z719@0o<|=mWk#pF8G?~4bDOB$Y?7{JG=1+Pc#reE*JARn_m^Z@_vOYDCo0*{*Sz7i zJHNku8me@+-kr4ExBa-Y_`!g*7-xDr{DkELu@s*vTp97>?|e_5Nxrpijv6X+)HC4e9Uf+#*qxFlM4ex}V&3-s^W*#?5n;1RQQ{Ul>%N?r+~A4* zg^8$6?EdOo)xAPE4m_18zEd zZ%xMMTNnw&wnlF(j97j;}3uX)McnItU|cg~zGI5?PpcU$@Kp+8+K?k`6Z+AR0)KA6=@#oyqA zq08jcHNqnzBO4dTKJZ&cM?^HVw6y%3d;k8u%-&ir;rI8?%^Kly_kQU-Sa7y{no8Q8 zcAqP(7_d>`rwiu0=jAo{(-?Z`atXokVnjvnXbc@?Tn^WMaPq?5jCONoyoKDbq8-er z#AkKN`}ZemLWlH3^6tpegX!}nR-YbL?r+c8m&bkj^r>=dDoK0GGvBCKD^-UQL3HjH z*KRv?U$Z(0J@-sZSQ32~zg8_bV|Sa~iTy_0tL1deXTwb0Y9ud!XgAIF=M!MG#9gnA zF)zypy9Dm&Dv#+ezMAaPx8)s++5BB#mJSA&a7TSeI}_Z}6AwH-ZqMeIU`hKoYJ5G- zJo^p!{X5Yu^!^q^eQUVXc6r3P__9iab8!uGwEW7I2kY-rUl9|hI#Wfw!4w7(j^JZj zc;n5E#+%)3Et~LKteS4hGsyp8IMbb`)+?gidu!i&y4q&t`6ySLINa-!@g zv=C5Cr++yfu(Oai$i4fGn$853otRYBI$r3M*{z0XI6q!{b>Q3d@~8?_8q+Exw(AEftfHI3_*rP53znC)4;+`-1G`#?oj-ah`6LZmx=H z(j_5yGv`E|XM96a^$E>KVIx5BO^LpQrbs>C ziBG}Q*1cJpFN;Sw417{e$Q+_S}T%qq?1Cw?G zAzCq@XNfdm>egsfyAtcYP(H|~xM5mHG)3iPScD6@koiZ!mhCHODB+srwr>CQ{v39% z*nKkAh-on!tRVb^vYrWdhg{Bo>e$`xyI~&e>Ev_Mmg8J6< zmsjrJ>Vt;f*=yhP^lW@cNq6h}zmMEMQJ!e`c7-15+i7i0>Z_(x*rAo6FY1x&86z6d z!x8vjKmEU(@P8YW)hB11P!e0L*%xhHek~@tYq+)-{X8tp z+tFKsDs+5%vW}ZeVM-zP%la~DK9QY~@tUt(Mvvrv)p*JmUlGpE+860AMV3kS+gR7L#n;)600~gxBCbP{Bw6@i5s*xlk(Wz4sT%7?j8YY$s~{ zNxr|j(pfk;Bei#l-V+mJUYm`TE&-%jyUJJh(K$RQ&X zfvZtB@++v!Fuf{_=*T%sErnL=O1yq{xx{eX_GTY> zB_Yb#e@U~G`>e#XaF&T*1^zDDW~IPHg_Y|GqW{e~xFn0|qa=_8McY4wpd_j;$ zAP={zp55AezwDo|k441!V+2beMvw|v+6WN&qDb0Jd{77TDDbnMz{y~eO9T-FqpAP9 zXhQ17n(jfd zmKza0Uv#FUz2Q3x(##)o6Mw-FXS8+pc%8jL9XOXk!Czt-W)BCqZE zv|P5gA3wD?h;%8zKZwz}OrPkLB-818?@qtKGSSI>5%6V|6c?3X%jX)mwMGz#t+ri* z$#Dmqxks7cH-j|{auDWd!hA28(d*onH1wL4JOnz=2v2zhu9O9gTAls4vGmk}N+~@A zqcLhPB71((tjBx}g*bBDcuLfbo?hAgWlGNOou|>!5{xby^q8Tyu5ShRcQg-(uUjhRGvY{xfh*2yN?ASqTxW`YUAT#OWNfIY6wI^LIn~}> z-ozfFk5siJ-BErl0X>~Liq~O!$q*PHQOCu>5y4K1Cb6*BRfvuIXb=9(?daV*B__IN z`1v+flnw0YWPeLgJHBxYRR-s3QVys$zG(4ARpogQO)zPnZtUD3Fvz|ok()&7gjUP$ zbp&6Zaf;qrJ%7r#NzrS~QS=Rzb~%U6%Ff1&AzMZYN2l2rxl8b`8x8N!kk zs=T&5h_;J;`oVgv^y!b>QITi`FVeP+Rj-3YDn;OmLxWOQMfR!nuX}SPQ`;kk&B=Hb znRnBRmeGDirU~m`bI}<*(g_XKfnGe@&N;PcUJQvOLt;MVJ@%5+=h1_ovv1d(O{uOK zdh$$_!Jcp360zlrNl%xU&EmNvN`Z8^fEty6Jao4u-0Sot?7YgA`rRQ~a@e$DJ|nbB zO5(4Wp&~UqWwgF9x>5Z25a$!f?lv1NqwOU=EIG70pO%RJuf&GX%AXps1&Q6}P6X-l z$9KzG+k^0r}j{L>et2CN^rw{5jvq+@9w&9h)79uUKT-BT_VTe z@Y30jvZMyj%+iF)FZyr7-xmHNxso|u7Esn><4*T?S56+Hu|v{I%vS5y|9yNwV$L&D zZ)Ra80uI(U-QR9)g95K2pxi50N@m_2(?hOgedC>vc6v~|3Ms*^w&*}?_V=Wp9m0A`yqVLvA{v7VjPU-~Vbo$J0n!nw#fE_R%dS@#34A&*_DZaYCm0*dqz(JmU=7>ZdY0ojFN_?&p zs+)%@N8{KdN^v$4xzgV8VKr4AkAjR~RS@!*Kg#g5bk_jFLDbzU%0dbD8~?REX2-7RyR=o~v^bQX$ zv6_VHnjj7Vmx<4fL5U5|EuQ;1+~Jb(&s|22CuAK;0?O9}t78mdfP^Y=pU3XH=2m3C zeruIZ|Ag(MsD*)nnYcq>?h$;~-6v-=$9@+;X9?l8``^I~?0;8v(^ixB^b@0&vn`$4 ze)Q;oxc}a%LydufKf(s+f z>y=Fh5H-P!jPGEfdWWR1U;q1tiU5TgY7}8VwfwW3H4vS`uv%wOR~>E~`7w8aDK@6A zDWkY+%jw$9ufdI78Frqmi0LY&(R`Rqypr4i6?s z>``LMF3TVH+E{A<5chk=Ah{<04mu_x@vc3;680~biK!KwruLX0^wwu;+92ufhQ;o# zbr&D2VRlf5P9W6jiYWDa4Sib7_=&{gl{hsmEVoBJEk*8fc9lJ;Yzf^@tAjVkY1 zO_^h^ey@LP_|J<9k&iJHB)45f(nKbu2xX#?o3_1E_&rO~+a?c0BC=gQ_>T*}u;feE z*fQ>yt4}qXDSR?x>HJ2MN9z8AZs(w^I3J;G(D;swu&LUz(DreRP;8K|<5O=TF09Ot z-qgzKU*$@Az21kAs92w{w|CvW)?m*SqHj%FTI=T0_gRdlEbUG#|F-pNX2a39k8Ik# zJ+6zPWgLk(wGrpU_2w9Z1TB(WGkFPq`(Ok)ng=-8uDlR}Xufb*A<*qBto__h6-qXy?mHGB5MOfekWv6~R$h>1!h zuAJ@gm&LWK+xqVX4m6A9P!UKWrwQOM3L}8_lov&pxjfBo0p|EQq@?Oz}4Ld)~&cfyyNgH?!x7QA@l&FfZ zW?d<#@1hUEJad5jw5(SlavEuQJNx~WW9MDh@ELr8zEn7E0&P0|X?#)K@Y{vh7<;m= z#HPbYM17XKZFGwiB~!vVc-nS0)4a}A2lmJF<*?dk8W%gEcmLvv?pD>%wFaCCp67#l z6Z-Q9UE^|kQO--O0)eQco;Scqyll^R%d-ufIf7TZi0+Vaq-AtR?h9~$rm^(z;Hc4v zQQc<%tVO}~TzZY?I#dT=jCpIpXyv8|Qpb)U z+mAt#!-QB(F5WMXP!@y26sv>Mur(XmH-|hM2VPu1hW&lbi8-XnVYK#gqtxc=SIzIl^>elr%PN(=4-JVF2XX_->IbeV(L;c&xc-v_{tA{#SRPWs;xSR)z3b(&P zVT~EPriMKRrG|YDCb}KI|?M?y~ zuBM{$ad)L{d2A^_8bw+XiyU{>7YvHyMJ&1E&F>G#Ojve08)#VC2^NlHkenCkv8e1r z=BAsaSAF4-3XQGZLi%H_C)llcEbc+gy_p_savUvSH-8}XmB}{#Vp)Z?F6CvYBuag` zurr@n=m@ZwL9w5oAK?`D)|HN0#VU|wo~J^rVP$Xcz<~v^lt~T7+szd3&ldA3UH}0h zydS-NerB+4$$u$-I}xBI>44wC+%S zis7Jue9DUdpt+a10Tx%cl~k8IQ*Gh>{ZVT)c406;cw}U^^XojSQ{(!r^&&IE6gJDF z+h<+Is>|1sVxf-kmkTgUY=%OcU`;U;U5Y-g;9CbpYSwzhjH2vjp`2zg;r#vI+8G}brRjCe8CiA4zB2oYW%$V=MG~xKkEf7sVQ!EzYHkGjleXaOs82} z`{T(Orxky}6em_>_Hrye$C-B(5zx{cTFR^>cCp$gGfuyfCTceAZms59a1T@4Nm&GJ zP38^8s@tv*YWy~uP`OKzc8lUq@jjUX)!sbhF$k^j%kCR2flX=ZR8g^Grc@4^St1?( zDw0+ewF4g5L?`4XNV{3|&jnK-9kafeQ=!PaQLLtKI&`b2YBqR<;}wFe-gLp>I-& zU_(@KgUyn=j|NW1;L-;zu{E>Ane5aHg&V0*G&LxmfMSbrv7?JiL7sF+OUsr1dpT&G z6QsxbOR_cHn?sy($1Hb727jEN*!og1=>BGJXCrluJF9aF2{bVV+!-)2-IPrAP`gPb z`W$``KnE+Uv{_>JKW(|m_eUHJiYxt6WS8E@e`^V54C-S?@E zB*G_Z?t;9YtuJ1+OlSdfuDC8>U@knYax%_TYkSG4*n-pRtGk+JO20blOYoK>kmuqx zW9gF<=Lkv)w6=D?K!)D=l2l|c-@DS4-)#plGVbl|85DyFWQXD?mcz)>pD1Se0j;hQ z2TTcQ1+KuNrdsG&^-}H54+1yIyZzO9uL0+gN_WClzy4p=;^I=&NV_3Kb3&{MC#i2K z)l}B%pF~fTFqPyMH9kIG2~~80_Wn=pDE1bP84B)ACHdd=6b6w`Ed3VB6N&LSMeP1n z(#{O)6 z4neQxpqd87*`qldzlveD6N+<^9ZmmTXf^1~bkIOoreZi$q#lX^G+NEAbN;Yy_R0Jg z3<4=ScA+DqqxlwfuuHwZq?RkhnLeC-iY(g6|4UK!mJlZ27BbQ^ld?#k(Nwv<0Vbvc zLiB9;>|$&K7$h2)(JNF6j1ulIB#w&|h}3T0>YqDB@jS!q&L>I`4$#7fhX3Bn0Y5GW z2BaN)rIV(BI34$P*JD1#079g6mgm?sanMrn`-Tu6Wk?XWgSmIDE}pq5=W+;1IAnkk zJr5q&K!v_UI}hBYa!2nH<#e_U3DvWa3iifc!rdpa^Jc-xHgeGc`e1gHRHm`U#d05b^E^a zU#mRuOW666<~$4Ryft@Hzi&%?xhxCIfGJfPt1Tsl@%caq*OQ*nUff+CDN%_$6DoBr z*v08&4=pqI#^W!GD{ha;WuDPEtgxnIK*CMzc1#qs`WDr!kK_$P@gccq2Q5Q@)k=IY zHZZ`9mYAGN`BsKFkUJP-hJZog4+4|2I%VxnxWD1qRtVpuT| z&~3Ps7swz>nzf6Ii&06pt7k$i6=Cx`t1(%rYF;&I28)OXhBT7dcR`a z_Y4%bi}ee)Dgx#z^0p_t+D9?dd zd0ZTPLG+Kq``-@dUD#+830QtpV|s}D0td>zfs*?7M<=3o{!<7qne_lC=(SrxTAx5m z*RoS{so4Fk8|6+0y9C93!q4eyC^OSR4V-pR*;-GSqJBB&MbP$*=f^I|(Fd<2>aCIE ztrkYgOCD^kKnY3-Vt~ybPfJ{eODiE3u-@O@a^L*@-Fts)+M2Y#XR!1E_#1>o=c3#h zlmCBpbZ#aN&wKY8_BJNC?O%BsY7VLt`4w5R_#fU?9!s_MK)N{eTb1@83slX_uUW)4 z{e_q>Vw=)J^nvX%Y@f=BT|R+jtlQ3q@MUXZw4ymWVC}W{TDMv>t<42uu683>81vQg zvWkl4)V<-<%WiIN{r_E&|3ik7-~RQw`2Z$Xows&pY^iB$I@NMj@_BxbtaPCXpL3ZI zB3|*Sfj1`!4a7(r<~JYSJ&NyZ0l$ix@d%XU1jGQXO;0o$eXuuY8LhbcS+QW$^HwtAw zt4pkV_*(y+eIgimF*!N5)Xll&2)hl=ece})D1%~s*q$)0Nb;E!PQei z3xeq=+P~e})q_H8xjaRHX{LO@)+~0Wb7)EMgf?YlYjBp8mupIq6snhRGkyWy85$~g zxAAnV2)Z^+SwZ~}STGC@yZQ>*>P`K;CgzZ<;63>AY)ez9`%>qw>F9!gGch$Jg#$FA zw_h8iiYD9V-0JJ1MKQSMIMyWRHo{H3=Lj`#$_eY-n(fmqHDn zALq_OWt*-GrbAFz$r*L-w%aK}tTKj~`3k;%a`Ihfw~D_FB5cR2?dl5SZFsw9Hq#|2)P%_eIc42nEr$vSdRa$cC zzrSHSU-$Hctq~Au zf%pO?$dMMy>Xg&Ymgd=l-t}&)la=$rbf9peeZCwW8f$D6kGe3Mg%2E|7d_Cv@zoOr zFx^^^l?afhYQ=0o*kbA1bu#nKwH3dGK>qj;j9rYfA*rY!z(&Uq9J+p{3pt+z6d%e2 z04C`u1nlWgj1_Belh!s?maKQ2F7YXBlD?hb0$YtO;&d+msQceNw2e1#&&uYOI022{ z7U+&JbXHH$>s-UtoC3Sy`X9T|DMrThr-vJmok0Fm{<2^C5)VZ=S^_$5^m7;2(|b7a zF^phl^5?na%j5U#iQK`4&yU|*Cmzn$u%kdR7vXdHF@^-$z>$`Z?_-uXD9T>6(2w3PnOqfn}ND*$X0*vRzIZ^#`I($B<=IWQ$MALQ2>E1*(yPMn^%o zdiHMz8oq?p)zYzn2wXXn8QJB1$Vy7;w7EHcR+_3%$eDZ!N>!Qi)CsbXI5L8cu$e>A zZjBiIM!v0&H>VheU$^KQ3c)C;{EX`<>Ncga8YvZG`0Km6W_$RoM!YU6GvSe!rRt8- z@=1CDhjsXHdhvnVZL^~iq6H>YY&adNgg02W@z?2M)A@$Sdrn1<~3D3ZX-4 z8$L=G8+vx;42>rXRV08E`mK3!VqEN2PNhv&#b#e!>7`D$B?ig@h;*OwXaZL{y-p^` zoT2dR80K$%D@|7UFM>|mXx5wUBna6$s9a+G#$wE$3BiD1b{*0o*;cc@ulF!pv<{U_Q2|5s@({vrM=wl8T+}>GDvQZwRLTs%>n~ zA$s;+nj>NX8J?jm1>F?uA4esMz|)xUCJEU>cuJM< z0{YbQiE8ptI@2asZpcd!sTqrb$d_J}7X+#=UTg?3nB^howYZ>5Fnd9OkkSG`h$bYc z*b|ipBtL^r$Q`T-etfF9c}TQEb1GW13Tz1aeC5u?XWC-pXY5O zKi&7FO+*(g0_{)5eyI}d$G}EGl&g#92rg8o`X3LA*2}u{VCA+k#CIJmE&|lhzZUf0 z1pOyewFh97R@*qD6+|UWHl`5*tRTWso?2Ag*!t~Xy-ZCr%&1Bk$rQx%{BW;O&x@%Noqk z_YjtQ7?g#n$(kZLQEW&_*v3$hJ(_{Y8$;%F2}I>(+xBp5zZyYQ8@{}Qy;FSLmf_^J{T$nY>?e8dMTa5js|$P z{Yzu$yLtD%Xv$`Ca-W&&qjGov$4a~d`_~RPk;F@F^&a}16`bp= z6|&7iXX*jFzWV}Zs$jeDK_9Q=k9H^ymLM5^tFbi{9QC@m03ul3;v|LZmuIr#D-BQu zUp5d&6o()Hrcj7J{k1UD=Y7b*5u)tKPgJFB`)tt+<@_;ZvrVuyv6GUUpmcAS9@@b^ z%<_52ODVl>{Z%6xAa|Qy>X~t-JO?C#p4dZW2Vq9iFrz}Jj;Um@p}XX2coy{B^D5~3 zp21$_7ciM6QKUV0@s{Ej0KMU8^>%T%aEUdo2+r1|12x2%duX%#d*#_0w@Q~Dj({3_ z#`k;QcJ1*b3w+tihMwXbd1FQP`44GQU}ncaQOlkbhz&ulCC?_10H^;Rnk7J-oKuVj zWf0*LAbkuE3JMyk!(vf900k#yfOtv>WxpK;^zP#|8YcBm8pj0da$JVipK#GB-7HmY zv@_r#@0dxMJE1y3LDJcn9`gD$Piw300*o*BrJ9YVZ{JF6`+x@@4vedYEX|}_tuJ`p z0i(KXpUs`q7mZqnmc#?o@jTybgkE`o)3?~zH(WOt^(zvI?Tg0l5`f@L1o-5&0uZ{F z&jaANj6Wlr?N`PJs#*$?_U`uZWx-Ce#CZP`5{H42P}N~~%-_u=gRujqY& zT4JGFqjcR?$L83?-WbX*<|FMvr1@ig)GZp#YW%9oo<%q!j3CjCE!)F@^+0}^dK{*4 z9*D7kGJ5^Ux7>5;pAH~xU=?fk7TlU4&*diVPmp#hfulz7NZm1e06Bbg^n;(LURrJd z>ABgfZf9(sxA#&+GZh}s;*!=q$Q_zK~*-0HR0;aJ9F$|B5B; zb&u_}%|3d}-j8Z+ z6M_7-lUF+hX+tE+ul6+Mwux-!V!7``@$;`hh;rg~iXxmY3Y-GcJkW$+=anp>oPPiv z{j|4g=ZZcQ<<+&L31^?Jzv#`&m)y`4vaGvI$=Q_nyC`PbW_SJop5J*Q^SelUgie;q z``cqC#BIQ2He<#IAs02W@R|HM2au55iZdtz#%6eA5G&{##vh>gqsNf9%TX0%iTfZzm@qc1^ zl>9Fxv#x)(OyepIcn@~9?;(gA5a6zUggXER#lV4HbF3jnh<^db`{IWspzhxh0^%Kn z3R*6KzAe9Ai3U_dz#-}fKH_RW1|*9U4iVVxIllh)TJ8#Ezq)VIk8Yjgfb)|uP(jjF zQcyvYhqKi5nm=NqY$$%M#wO4NP$00h41X9|Qjk2$=WkS$8%%oueJ8PhBk2z0`N?Y& z#TGSkf1qypr&BSqS?-#;!{>AYz-(!nlRF7IPm}|LBp~3^F22S7}*Y|jPGtxQZ`Suo}kp5-qyX!@jtLI*TtkSsfwgGBqbFK zrnd4TQ8LMAk}C$muuEw1Vp0~CDGAF&0l0h%?TYJ5-IgTuBFKRA!`#`08Hp@Lhl3DM zc`DwF3IpqU=U@81M-a5!0EU0}xP`)}sltl!#DN_()ja70d$c#3`z?xg)d3a8J)1yZ zd|t!5rH7D-#I!uVu$JB$J>gA=WMOAdTJb~g3NQtTv;g@^n?2f9gt z^{rspY#?8K$qlz1IQ91KVie(`NWKLY!1{vWYN$MEg}zfACZ9yT7oX$KIR@Ft|%|(T&>|z26vY&Rez5N(R$=3NKsR_>B#J^eztU=2fC`P z+!Y_--)F7W0&v~r(WxN4wGE~v8X)kREa%`{ywfJQ!-S}UR92WIj@~5j2h8wivd5S3 zGewoRR)__p$%DKH>})Z%~uH0WeW3+&36=@70va|6#vM=L9y+cm@if?@maFEzy{ikrcTIon@P)NM98IIm*qp4h@U zUXj?O#6wij8EUUt7ec_H?x+T?A2dN(g#e`>M#ewf+EUxLALz1ygVv1Rp8%vI8o+I% znk5wBbuOyP3gcCd4Ht_f?q7DeCI8+sB#;R;dDBO9Vqk&=W?#)WQoWdhnaz4boYYUM z|C)c-8CAFWv@XnYuivt8dl$?V0ZJ9u{|iw`-R{f&foicCU$Xv4u4a#jnUV*7n-!RY zvK>-nQT+L>l^5ymoCUG#NL5tidhq)bdR3!6NqQ;7(hESP=E1ss;^t?-QyTg|<$|?2 zhUSKupNb*g=q{xarl*e{@%nV#%u=WXuBH6B{)4(A`&56afc@3f<|6YC8LwZTVm^E0 z$O#MYoBBer8S%9I_d1Dmua|1De*OwhwKvbI^8d2y;i)EE8vWI-{t`5Ma$=VWMF1OKk(LhEAUYCnP znK*VLHzvrbNQ8_bAMzl~$@^9g5-5pj5xlF~z^{QnmE}YAxe=C}ovhKXj{7Q$LwLy0 z0nrbQ@oQ6qRwVLt{)~4Gzsn_{a3Kt?w=}G4G@&93r1!UrlRJzMvqbVccY*2Lg_^2^ zw_~;x2-9r~ zT~?C(uOEbr8L+QLl$^z=*uUMN5?0!SxG?L5bXIpRL)nW6u^~eC7u0IL1`H0YpII3O zR&uXYwNXaJ}J=7OEw*Qgw%AY#1%RLp8cd)VT_VL3G0gZH!&(Y{4JX}>jquQcu$ zNbiP7H4DWpWZl7!Px3z_laUPulmcA}ew(s6mm$W~Fxb*4KG4pti@Xm~O$T)?*T7-$ zXj+cHo#~~z=D{dzqibx-We{aI#=eg~y?c^E#|OVyv%mz$oN{}8(H!GQr@J_cu@w|wss57SMs znHFd)e~}I%Wh&vWqp(UrmS_fz%XW&afu$Je3^Ih6(B+}})&gH|BfJNA834PH%itEX zLpq@eT(?tllqo0E;+yL%0%e0N&+c{?xh+gp{nTooUF0~XOaH$i9yJg<19T<-z+>LI zKDaue?wF3Vc~w?)typq%1ZC$^K^VN|_>n6c{;u9*Ep7t_)&6z!lME8^QT^axw%s_UW8Sm*N{s*PGIG}0{mtm|?J{;b| z>}JMUbTP`cJKyJO{nG~(rCJ7|2#5@YiT zqMk4sFS?}fh?GS4Cz>)P>%K6W(s_Fe}pmB2o6H zfXOY1BBvN2*@9Xc#J23@%yzaDEAVou9Ig_qT|E|Ml4uYjO11%7*ub7FZt zD<~yrHNYUm@3S&+r2>Hs>pEMyeW2j47Meuxk6WPEJPkjUs98q>85MwqwDHsCC zW_edO}>q6Q8n`&A4ah~7&{^>H|5Gte-{e$1iI4*D>g*EX8wGmu_f3>0h;I`8F4 z%VPu({le_#Uw||VHt$YMK{&$ba~$a^C{=CAe?caoI&d?=iIoi~o!?|I{N%9Ya9bkT zU%^InI{MvBUj5gjoVXxDTlM=spwQjdrK>L~6^a2p0Fp4GWd%U9UP>=;^kJ(6J0IG; zFyaM}PU*QBXLc<50?#{>rX1Cr_aH~=H*i}|flwUhhgk6K5^(#`anFDK4#e$mx(dxL#QV3S-&OzfnzT}bmN`|XsnMzb? zCD7JyiL4T)@qWO5K}6!%?1@1LSsTBOl|N%VdOTK5TXH~ZkD4b>YD^#hSS z`5_~q*Co{R+(lb{j6I2}(;>XpW8pXAJDcnK&&C$)or? zK}fNf0q?teM?G2*rPJT+SUx`de(j&%<94x1e6k0`8nDvwsSmX{R(b$d%iil{#w*b84r5|kG+52&# zt}9W}c@4t-xo$OXgJ>>=UB%*@oGz&Mj^9yD?^=C8~ z%rWJO)CiA6B`cMP)u!%wvya}_6?cX^IXY3?2;b6jM_&5={0&! z8+9cAD?9o+I^VDy?CkmfD2Ilu?U6}17RAw4Owk?SVRMzcQS>8E?$l0(Iu&vi!Gs!Jy(${Yf5Yhxq%8oI4LR`HraU5dLOW1JpclRt^C)WFwPK9<_xp*xMK4qT z;`eY|M%LD^>lE((Y1ko8f&Jjn*rP=C-}0#NR3|w zP^f;=+(F7TeN7%d7kOaD!;_d$@Z#pLqlgN1uk>d`?}mSFWbM`zJBiRX9Y(7HnLOtL z9*=*vHuLRM@P0EOsp{9R>1k<2jr-0OG6zNW9plqXkv9gp5`K6R30mfDa-3<<@9-g9 z8QN&1B+sA(uNU~|2r}ITpFO$TucZ27c5bfRbyP|-^@2HY+vu)xBO;jcd$-&TUOoR! zO~Vfs6LIZ-Al@F*!5wc&x~xQ_BiKu2-S zbdT(P0#L1G2!^M2C@wkne!4r!+;2vD^p6s1WB4s^CEG=-_wZ5tLt1&&Oly1OK8Sc1 z{8zM)^!d3Rkw{$dRlF)BBy_jz!@Zdd%l(yjw)6L0o(^WcdE)?^ortl0Qd(Ndk=_NL z#<&NcrO5&HT3CE}+ScB_!h}INlRo27E8T%S+atNC6o>rV6P~MDR0Tenwq1~kcYA#y zKUXhXNI}6Iu3Pu5b}C}~m8S0QZ4-O9iF2Gx#{n`8y%&J{`SYhC3F&U=d*w3~ckwjF zn0NK_^`PP~8YE1xa9IH{j}!g`_}pA2z#4Y zm1C-Zyneg%IGsRqLi;gB-2O{nJ5GGG%vD|eiDC@UFZ)WYn@#kS%-Y~PJ>MNFq!IKt zfTj#?+{gjF$W@ExVv};Yp|f<|%xoUZ0U8u?;G9d)x~ObLd$M$t1E|U8S`nTb`VeOS zTUce%)Ud=32QMS~AU9Gqnrp!Uh)&aFZl*VG{6bsz{-}9qU|^8BOlY19=TGU%xcqbt zAa|LZg=^QYamHdUq()y22)ciZmoHF=6$xZzx2G2?Yzs( z&kII;*>(5`j^y2+SFk&P$_exRu1N%EmRC~JJyA>A-%;*KrDfui85tfPCUZ41HqK2> z=D?ok<5Shk(xAPXiWCE#-a-|0vi`IlnLBp@;Q0*{E(KLEbU6Z_&t z&+L7KT&k zFJnDjU1?xw<@PxgeAn`dIpw|Ia(RCvYxY2ut8Ra;ZuUpxQ(+72=I@&wIB~U9|8*Rx z>dl+e!Z+T2OWY$nBH!DOX*BRuU8gK->T4824SL&k;ovJjXZqtW-P~H4#J9KWU%ve4 z^bByd6UD4URU=r^fBv+Ql#oElutI0fc(_P=&gw^UD_(fxZ~W1_L|t8dH~-Srp8Vq4 zy!m+-&}>y$`C#MI6Ch$kCuN>QM-M^}(f?X4`q?^d;Bah0+>-nkb&a=9U;(Ya^Zc2D zq*2@3qf94hV=)VT`9>mpQ_^+az-|8NeF5SvpwbqB;*Xzvh65g9+*#wJE42-cjgG*W zQ`gq^k_T5KyD(m3y7X~1>hF8YTMvBUqc(j|n0U59i}YiDGoxNb3SFtP024X*){#=% zrx*g@)gzUr8mY-ew8`}X(u&y15C@D2y3^+Nj%BO@cXHwn?6+iQ-T zvUgRMmzUMAUhM|4oLt?C*vQVvxI{@wnFkRgC`@cZLQYJ~vCOP2r}{v$o1<<%K2=#- zuksZ9)&*o`WnFg2PPT$dvQ)kPaV-LLze4viLV&2w=5!30ra@}$D4LxB1bHr4kz+y_ z#5&nnvfF>19l5D7XmXlr@!w?3Kg>)_gw)*)e7KkEuiotFrbu@C5A2A|{o6PQ85!~Q z0_ZN?{`3_ScADOj98uQj6n<%{?`mTic-s0}sdzj+!gB-+=#`5C6bKdTD^ojaLh{Wf zOmM0s+AzrKHmk}(`1+h^i%Yr2|I^-ghjaPAeP5)~kX1%j%BUonWk(uTB`Yh0vUhxw zl}d_ID5D`;Hc3cPva+{CX_zHslzqRhPv7Tve}B(?KYu-cJm2X% ze4Vd1UAx$xfK6bubF|DetkkL49*^WU!FDFi@DLM(N;aG|j?^gu=6Qw!&r-X)dTfJs zWA|jn)-D#rB2d+=&A)=MUKidjw%$V72$5U=l>+Hj+fxn})}2o%>v~Bl!2Ff%V$QO7 zVNO{jTgnY_N69=9y_Ng-No0<1T1q;^Z!?q1M>fo%z@?I5+kBkiDhIOVccT`X`-e^)Q?~xw-G+WyRyu*lHuJu z27Zdc{FBHx#<9cxGMfFQLyL{_w7u zr71xiSzeB*?W3a#-z&oapaxNPMPtLKz6QFiilOW3kr2M#n-ujGAy3CM{Wk$ldBn$}JTdAB*A=Sp)=s-AS~n!vr`} z!NJRwUW`ezS{w3x-4Ak=>^}Rmfi8UMwWaFQn)Bfxx#SKnw{Z)$EgT5oBWFlH+|!Xy zU~yV3JbXSmrGjJ87|)U4l}~V>6JM zzq&mF*`nv4-Q>@qA{_gnG3e>^TwUu|;>ke_m4b`?4L8p2dhgORw3A`iMLNRyN%Fnb z#%{03c;QDI%hjS24btQeG$EbgRsSY~*V+(xUh}^%+Uwi;j-2rgzI$5h_fr>x6SQhwnw$ zWB+CBGT2`nv7YPq%BPnq1HroNV|yb zk8iL@QA-VwGA#y6cu#R4rdFM7XF&A*DgF3Fsm)Q`1u9!s5J9W(fWDYo**||gm32G# za`hRnt`lB{CC!b31o7###n8U;s1kBTZS5g_X~9J+Hf_s1##GTF5dvGWQ`Y8_RsB1& z(!>iCyGXsm*f8~gGJs6s1G*zLwW+~hqKQVndUy<_V{EJAO%o}q@9)hXHgXwQXu zU1-sXWsNJY5|DHM?Pn6M$DPj*_jIb;T&mFcyHvUyjB(bd?S}SE)~W!+LZZiuuZAUFSPp9=cqM;#g&?6$t(t<7g&4SqNn;2}l zGCkdL$dSeZ$01FjQ$DudrR$BTiN2AHO#hk%F-CGi^BI;?lQFEl@~r$D-gZhLl6v)?KfD&D@(hQ zhGg}-!Xvk~r#ob0=Ivf}oWK|o0jKTAHFWz5+2rbug8^;%_S;@*jJpP>WoG2_sq)B4 z-zUc>|6+aI7rz&?>R;v$jzfUT66~|$pUhgw9&p}zCo?HKDcM||25Te#LFU?i`eSV) zHJ;aRjIq8VjSBsfry3pz1&~gniaS+~LBMyYF8DS%c}@foUxO`ANVfkDRq3lj(+eM# z?Yv8Afrcdhx_`qek6$2uUEBS=Ycm!(au+}_7H4iJDWK}*>m>auBp55zzY1n~>J_s` z8gESce%doxjtg^gLh7dKyZq>g{i);x2uIaO8$YHGw#)p^kqN0IPQJKsV2|55sTJ&Z zOIHzFqNFO)*08F*-|cto?Q;D-?Zo6rRcC2 zFQNIWcD$Qdz&CM=#1r=FBX!$W5c>Q$IclE23K`v>Y!f1TaXmOWx?;#LxG-a>?q4T0_Th-% zhBM@JOuveNS^aiV1KN~aBzrQ>`?4Dys8Y8bnf@wUp*M3wk?fsuJLI-QVChcG>+4uD z&2TI-r+IspI2GQlvyAEjLGcB9damVh6=5MU0|dW`!n6UaJcm2&YMOBs1EIAEVo^21 zX?i2qLo?`GH<$DxoHYw-ZIhDC&dfG^AtlU^N7>=sO?+I^e&1_$ z8rIwR^~{@6o!Tx^`mOJb=5?PbIw+}TGNGl?DBddPSNHS5jF6vHL*bqMcsOgW^u*Xz zL0T(&axuy>OpDLq@S4#$1*^h>##>+^~n8wQDI>#b}%RLY?93=6&#+6~;a z8Bh;Fwx41;Kd4!@_6|CsHIjB|$7gk56o;=SI7$>0ABV)C zm7eXT2gNz~n~)Kch#9=N?w4d`ufk4pv?d;bI@%iB(Xhmx9^={x@oywU>2fy^%^8)OTfQ*BX`SS@}bIL$=T zIK^NxC%IoXR7hS}=}XJ!_gqa9t3~fA?)ID>Ip%0k$qC^$-RE?uX2fux%a74I!~95T!mY?*2d%MZ`zmyDDZ< z-_NU=&r5zuMv}gE>8P@Cb=EVMQ1#PcrnM)1YF4~l_^Dq1P~ND&zxh}H3jYh}m!sjG z8arC-7J?k(NqPf|@3~yDU-@a8t*zOxbK5Zxt=Ty9h@B<>7$v*Ql1|wXAh)wMk z>VRkNr_&T^RuIGou?sq3Wcc4mC-cQY$T!X$Uk3+;1O=-m zj1l2q`oy!O1V-j|1eJuR*ZZ?dnr;tS9+cA;;o*7=L3`S(?!AKWuf{TL{r)DA#jU=r zqeG2}dA67{P9x&NP~6oiY-rOFk1cpkU6z2Y&|Kd&MjCP6w*2+8nK!LUKvaugE>)eC zG0y2?rn!DxKut4hkJ98oy*9E^4xTHCg;#Yjh;`zvN9&jGMw+_g8(YK}(QVM#+94#Y zcJycg#?0CJ)QZh^@Zs$c6gVl;{zS}zBg!}x?nOBrh!ml z-)Y65XJuuTxpbdcGoCz^Ps3&EDmYn3<39^VH6kA6E3a)kG7K`8Q!5q&EOQe$T|_}`^tc8aVV;{BEl}&1MRE5nmF-`4{v8!`@LWL{ zfBWf-?Y7~CU|Whhb9|a69b(?VhFgs(O34Lt7xSP?=Ug6E7c9xsS+gSefcBo6OuFO) zY1fkqGTwx{X`_BS8N?nHPLtwpVt9kUZ2=S#o@TRLiK?0|G?HJCwoCi$t*a80xF`RR z>>3vu@)_0?bEQ*U5^yNY?{ptL4Aqto%RDU`^&gZr%)2*xGJRecd&gOcTZz$wFbw0D z2h{GqyQayj>2?xBLUb`%JSYFfb&7N>$;*r?LwC(_g5XRf&j>&x$;*Fb-HtD!Ai6H- zT}k=HmYuNJj%{qDQ`k$NrX9)g(SOr_)K$Mu_K#K0C zc9DEYnD6~QFSLF8aRUQ`dj+Ei1qRY_uzco*vVzP2j&S*%V61mMz_nFF3VR*8+Xlu5 zLGDYcC|qSkpvaYGB#Wj%0K@)q45{;wF=p~F1Pro(#=i`5arcmK>Eeu0ds`b@FZRs1 zXt~QHU-zL9y$$9^S4+#>`{ViZ=WlwPIg=bn&%-@Fvu4ejxF1ZEdkfd@v{JT*J@7)E zaVCve*uS#SXv^B)jkkCYWK3}s#sGrb{eyZdTP#mLYb-D2XEi=|qEGgoR8(6c*4%Ds zIe`fHej)v=Hbl=erV}1LDqG$sC3W}KO2+hz3=PzO*iEXP9+y&3v=&h~K>4A^=j^;Z zF^8@kCFI)+o(Fs*bM?Au)m4=Z@?TilQ4&`rYJ~AMx`BuzkBg zU=OT3uDzc*F39%(vA~Sge0+RE{gI{}$O+t-Ns5bmCF;FHQnGz$i88eN`%)v;Loowc>MpL%xqCZ|maJ>`*G_FNX@6*GTLQ-5#oF<^_IbzBV>JSNkv z5Zkz1n3$+i70RP;V)Ed$g$1(rBE2qEThmHPU%vfuOM@tMmX0Ix=iZ*B-rciQZe^0TxKuLf6-R-C1!{QUVWh~25m zwa=eFvzQ0-_5L##yCC|Xx!8E5wq{MAo0qRjLUyzzQI+>wC+10c=G3?Vrk9-Be@yfK zscjM0m<}F1c+(?tx5*R^Z8!l3$zh@!(fIlYvie92X(5Yy^!~n2%kTazyV@!{R?_<8 z7p159nkP<#jqdJ~p|y%v0jG{6(j*V0-44MaLDIbisU5>!(tR;3^?h-sE-J3O+fem{ zgF_}c6-WUVq$@5Ak{$-}=-g#7t<%UZZkBPkZhAjG+-ow{f$6 zka8f}s`f6)5JdH&Ks)8=+;};CdCxXkpurUmZx4@+&91CGw9F~-XC|h8chF7$=ruBV z?lJw0gM$NkI;-5nl)tGcDzhI6TeY68&pKA_e9MVrl;4i8LfUn+g!L_)2EyEAhYZ(z z?8rZIk(8VMI^Zry*C~n?rbaatJ}pg+UT)mKv=>Y~=$&BSKQo0Cp`aCidH(Q4A;A$OLNeB_T}%DJ z7X&OfXGl_$yqA<){L{%9WERkH&VTac$%^i*yu1@gnNU_*vPLvh%ip)%rL{87KpQb) zZVN>nZ{D!8v)=>|=s9#qmfpwzm8sJ;GWzEC1H1z|$ZpNc{t-IEiu+=_JO0kUJv4ge z_yGCO-p@}-mYx;#a*v<@=3|ZWEU)1^5z*0);jhb35^b1xtir5V}q+Nu%)F%5u%x%m38#j(&5&U(kqkUyFAy$ zozT*X3kaar)zKOHJ@G|NLjz-U{e-r*q*)p9zL(;Ql|4y~xaR6JH`NT<8*KO}Pc**B zpZllV|MxD%X$McDIS%vq++(x_>+`6jIq9s zV+Ym$aBD=5YmfBmxcqz_V4B*R6Aqu78q#_oe?bpRwC$17C3ocS?Vxr>U|L5VKkszh zE-bvL?c?)eWMssB%uixd!^VbVc|F1^VE|iL!!GX>5pkM6ha)grQ8!anab%-QYwmK* zG2fN}kyAJ-C=$3VuFu@B3c%H)0#299r9ORSRKYJUs!H9IAKWq~<2DfD3sMX_ENmDH z`)O*$L4N?ls6x#^z>sJm{58Q5fB7tpNA0z37GC@Nh6Nn_Ki>kOUgiMCk?(tqT z`txsyT(Rj}e)4=9udlZcyN|T(SRsF)uxPZ42d7=xkKL$t!>>R=;fAolI(x1Q4N4zK zBUpc5S{qziJX!M<#$kAGB2a_Y&>m@7SyD|0E??pZl^xvTM5=n^Ub8u>k7%}!LQdyi z(fdAXO*zIf#6?{(&Yy2^>eR!#BBnq7as= z&Iz3SBKQ<=5~gn2fdY%=(Ptmil9L4z$F5bwwV-rLM2_3?TbJ90pUER5wkX_FJYi*( zv|QBHb+Q-ri^Yg*v}YRaRn^?IX{>blC7HnY>jh$BDuHM5NGU2OH?gvwEg7?x-OKHX zF6!v(r(|U8s>*qanwJr$r2(fYRuH3%^b|eWZCu#&IcDK&O*pCSG3@O{ZfA{ptE_0P zngr1DwnpC>@hUfz-sOt{9ng4#OI8oc$1vae_Rq;le#r;4r!I_2qCn^4Nk2vzN4c4Q zcZ=3P4(0}@Q4-5;d9DYfk6Vc&uoiQudCP{H8bU@t ztZgHMNy?vyKfaN*zK$MRSz;vOu%cs6G7G@5pv8A}(Sax0@nu6Pr zo%3r`Tz{oSYzbO1 zJRTOG9$Q`-D+UW+c9EumK{U7TZ*!AkchwWxx`_2$*vMQSB2eH$(&I60vp1` zyz`W1<3o9_>>~H!Fkh6!Ud{BQ<=M+U4-L2c1*?RPU5qbrjNVg7ShyW5Q2DO?Z|!_{ z9CmlD%_ZMYE}8^+yN(QhN2&FII6;6XpNSY3K0-s5H=JK);C7V%asN^Yq9yW6fWDz; zae6D*RT1InxV(mWVHAFfOZi)VS}TTp;FYSrWy#b?~} zz)#Hae8d3`d`sL~K&@a*4?<)*mlo%_DDj~m#ou4KP>-a+ZT6g4mf7PvVc zz*h_at=RC_dq|2*_wsz{GDn^-Z%NNz(3|I5Nmvlhp#Ypjomuec`o9Vd?K9le-tE_; zKk!?$StLLYa3brE9`7BfoIn=jtN+|{%%h~UUEt6q*`kTA!NEaU-?$J?>AYTuI@k>9 zhbO9d^fAVWE(N}mgBmZf!T-`SGTOjQH_E2W5L8?&bNF6bVD+cUFAdVw_q=}n7-eH( zO%)Q#-<@rFRsOOVZFfI^7@&p4KYmx#@(T*eb_1C)_u01JK1|QI=&=Qqm4V!NcS+RKe>uJ>H=aBX>RzD%FVQV9+koA%Ro5?=1ZrgML3i@uvkZMR;ayJyv7H zI@aC$)c97@b91)M=)YY0$-JVcKC|lOKk%H+8i%}-fOx^-X_}~=e+XW|MDTZBN#xm+ z4L7&oIicR;jsf6hgKu?i{RMVuSIgtNDM1>aqFHifBzCg2rT?C_#H2ODS}QrIYOtgW)Q05a{9aHp4kummd-x_2rzB05+GYs zzxi4NR&e>O@_e2~9Lpd$I5>OrzoS4!?qfxacI6%b zbAZc20I6C3g#ume|I8Wl45m%-#f3pdPq|-UncFgyrmb?!fF{pdAEc9ZM1vot%_7v^f&WRj!+*ht=kE`N$RXGKnZ48Sey#EM#9-Tz$bKBuwQuKv)f+fd&wcAy`67UK>(v@aKCQkR zTr=5j5_btdO~UH!M!0#iDGa2H9P|Cj?C>tH;?Q0fGL#&FUSI(p?eOaV2EXe44~17D z*0M7&Z?OpvfIUBHvX*##&F#4Rg~jR4={oVoeJ`UvVy`_H={L*D&&PN2JNxS$VK*Ua zW5b=msb)Wh)`-;MYlN-CQG&D*+<{>%F0Ur&Z1K9 z*G~q1DSr(ru<)C`P8;S9C1y0Th7kG8xlx`3qLbkYgi>w}f+>T?t}=xEuGbZ307OZR zj4cG0!_LI-+wg2i&w^Uwj)Bvf7-KXc|9ES^peke?+#(^gdv^!(F}XNoQCU&Ia^e}H(~UX+J$!>(-Vnzy{d$J3Px_Zev9hqt@l?He zV?8z0%sa;Ox=OR^<$u3}^$hY>;3-WcXhPf=_=X`#!Y`wKjkK$y);ilni}Cunb2+62 z2%UI=&t!2FyoNhcj}O4>a?LMNb>vx|WXKmP)qo2w`)Cw`f!z(S--S0pj3bc4BlCgd zuC2MHt^vjWEh?69`owv;yStmkw9ct`SOCc+=4utuXPl%L+YT*)phDG0h4VixgNrzf zo;ZC}I6RNz0O<@qHi3YJ*6^PITB--Mpw1XMX*;$wKSr;<%%I>XJM(|uYKJLMjc>qP zW`a!@l*uiW2f`v!_Fa#FtTy_->2B4ZAGVyqy2}_15${ z2MoMa^nU@|dfwj^78bUTN+s2qXC844GT=6WLp#OweH;!Q3QT~+I06!y-o6y8=y$47 zVLl_FkzV2VBP=Gv-VOeHp&7O){O()PKQJ)&9yX-~QAm)rdKr+7dPCZT{Ti_fy`nw# z5cY1IO_=^3%pET)&Ou>(%wpVj=yGN(?RwG!#==d&p6?rAoMABnYwtmI`c9{tl4k#A zC-CRb4WHWH`BIef(CfNOI^)iss+W*k=wEj^?7)lD8$F_5w>Ukh3SzbvdqN6E1Oh7Ov@7R495Gv*wOKhAVh0(!Cv5uC7#bH@ z&0Bn2vK|`ZXX=`Q7`H{-zyCd!ec>n5G1VBEb7_^KJWdI?=g)}9NY8c$aP=g|E={E*7YFUZ);TmdxD}!9r{1%!hgBzFWegt# zbq!7q!zO)>Wc2s=1mB(1K@Of&EQJ)8xB%1zN~G<4KKAmPUsiE(afM1lq&jfFa4IA% zs~MP0R0uZ301s8T?x=$!BRq&3nNFUAJ?mU2QnG&w^FjDK&(D2J)=k$Nu3$B7LtRkE zh+TMC*i@VLF&o5~je`iHE|V1bCwXDO5voXT%8PL@uAvV&fqbTsVF! z8_zz!K)dZ3*E}|hn>KDt1<;rb3ZFh54{?ZyiAh750^1lpsU~R31kfaM!wi|xz|xrm zPemIlGT`iw;A01hg&aKTJNN%9(N_K^M)d!1Vr2ZPm>+0-?PvbPrpmNW^N|n@@@i#! zNkjQyM)j{GqG*^q4wdctK?%Mez literal 0 HcmV?d00001 diff --git a/prob3_ode.m b/prob3_ode.m new file mode 100644 index 0000000..97cc9a1 --- /dev/null +++ b/prob3_ode.m @@ -0,0 +1,5 @@ +function dy=prob3_ode(t,y) +dy=zeros(size(y)); +dy(1)=y(2); +dy(2)=-9*y(1); +end \ No newline at end of file diff --git a/prob3ptB.m b/prob3ptB.m new file mode 100644 index 0000000..ee1e635 --- /dev/null +++ b/prob3ptB.m @@ -0,0 +1,11 @@ +dt=[.1 .001]; +figure(); +hold on +x=[0:dt(2):3]'; +y_n=zeros(length(x),2); +y_n(1,:)=[1,0]; +for i=2:length(x) + dy=prob3_ode(x,y_n(i-1,:)); + y_n(i,:)=y_n(i-1,:)+ dt(2)*dy; +end +plot(x,cos(3*x),'k-',x,y_n(:,1),'o'); \ No newline at end of file