{VERSION 4 0 "IBM INTEL LINUX22" "4.0" } {USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 1 0 0 0 0 1 }{CSTYLE "2D Math" -1 2 "Times" 0 1 0 0 0 0 0 0 2 0 0 0 0 0 0 1 }{CSTYLE "2D Comment" 2 18 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 } {CSTYLE "2D Output" 2 20 "" 0 1 0 0 255 1 0 0 0 0 0 0 0 0 0 1 } {CSTYLE "" -1 256 "" 1 18 0 0 0 0 0 1 1 0 0 0 0 0 0 1 }{CSTYLE "correc ted" -1 257 "" 0 0 255 255 0 0 0 0 0 0 0 0 0 0 0 1 }{CSTYLE "red" -1 258 "red" 0 0 255 0 0 1 0 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 259 "" 1 18 0 0 0 0 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 260 "" 1 18 0 0 0 0 0 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 261 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 }{PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }1 1 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Text Output" -1 2 1 {CSTYLE "" -1 -1 "Courier" 1 10 0 0 255 1 0 0 0 0 0 1 3 0 3 0 }1 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Heading 1" -1 3 1 {CSTYLE "" -1 -1 "Times" 1 18 0 0 0 1 2 1 2 2 2 2 1 1 1 1 }1 1 0 0 8 4 1 0 1 0 2 2 0 1 }{PSTYLE "Heading 2" -1 4 1 {CSTYLE "" -1 -1 "Times" 1 14 0 0 0 1 2 1 2 2 2 2 1 1 1 1 }1 1 0 0 8 2 1 0 1 0 2 2 0 1 }{PSTYLE "Warning" 2 7 1 {CSTYLE "" -1 -1 "" 0 1 0 0 255 1 0 0 0 0 0 0 1 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Maple Output" 0 11 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 }3 3 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Maple Plot" 0 13 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 }3 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Title" -1 18 1 {CSTYLE "" -1 -1 "Times" 1 18 0 0 0 1 2 1 1 2 2 2 1 1 1 1 }3 1 0 0 12 12 1 0 1 0 2 2 19 1 }{PSTYLE "R3 Font 0" -1 256 1 {CSTYLE "" -1 -1 "Courier" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }1 1 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "R3 Font 2" -1 257 1 {CSTYLE "" -1 -1 "Courier" 1 12 0 128 128 1 2 2 2 2 2 2 1 1 1 1 }1 1 0 0 0 0 1 0 1 0 2 2 0 1 }} {SECT 0 {EXCHG {PARA 18 "" 0 "" {TEXT 256 33 "The One Dimensional Heat Equation" }{TEXT 259 0 "" }{TEXT 260 23 " - The Leap-Frog Scheme" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 42 "restart; with(plots): with(L inearAlgebra):" }}{PARA 7 "" 1 "" {TEXT -1 50 "Warning, the name chang ecoords has been redefined\n" }}}{SECT 0 {PARA 4 "" 0 "" {TEXT -1 18 " Analytic Solution " }}{EXCHG {PARA 0 "" 0 "" {TEXT -1 74 "If u(x,t) is the temperature at position x and time t the one dimensional " } {TEXT 261 13 "heat equation" }{TEXT -1 13 " is given by:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 40 "heat := diff(u(x,t),t)=diff(u(x,t), x,x);" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%heatG/-%%diffG6$-%\"uG6$%\"xG%\"tGF--F'6$F)-%\"$G6$F,\"\"#" }}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 55 "In the lectures, we have found the particular solution " }{MPLTEXT 1 0 0 "" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 52 "u_part := (x,t) -> exp( -(k*Pi)^2 *t) * sin(k*Pi*x); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%'u_partGR6$%\"xG%\"tG6\"6$%)oper atorG%&arrowGF)*&-%$expG6#,$*()%\"kG\"\"#\"\"\")%#PiGF5F69%F6!\"\"F6-% $sinG6#*(F4F6F8F69$F6F6F)F)F)" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 19 " If k is an integer," }{MPLTEXT 1 0 0 "" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 7 "k := 1;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"kG\"\"\"" }}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 49 "we can see that u_part satifies th e heat equation" }{MPLTEXT 1 0 0 "" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 30 "subs(u(x,t)=u_part(x,t),heat);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "simplify(%);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%%diffG6$*&-%$ex pG6#,$*&)%#PiG\"\"#\"\"\"%\"tGF0!\"\"F0-%$sinG6#*&F.F0%\"xGF0F0F1-F%6$ F'-%\"$G6$F7F/" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/,$*()%#PiG\"\"#\"\" \"-%$expG6#,$*&F&F)%\"tGF)!\"\"F)-%$sinG6#*&F'F)%\"xGF)F)F0F$" }}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 38 "and the following boundary conditi ons:" }{MPLTEXT 1 0 0 "" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 39 "u_0 := u _part(0,t); u_1 := u_part(1,t);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$ u_0G\"\"!" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$u_1G\"\"!" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 23 "and initial conditions:" }{MPLTEXT 1 0 0 "" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "u_part(x,0);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 41 "u_init := unapply( eval(u_part(x,0)), x);" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#-%$sinG6#*&%#PiG\"\"\"%\"xGF(" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%'u_initGR6#%\"xG6\"6$%)operatorG%&ar rowGF(-%$sinG6#*&%#PiG\"\"\"9$F1F(F(F(" }}}}{SECT 0 {PARA 3 "" 0 "" {TEXT -1 46 "Second Order Implicit Scheme (Crank Nicholson)" }}{PARA 0 "" 0 "" {TEXT -1 217 "The leap frog scheme requires numerical soluti ons for at least two time steps in order to start.\nWe will use the im plicit second order Crank-Nicholson scheme to compute the numerical so lution for the first time step. " }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 73 "n := 20; h := 1/n;\ntau := 1/200; steps := 10;\nv := array(0.. n,0..steps);\n" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"nG\"#?" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"hG#\"\"\"\"#?" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$tauG#\"\"\"\"$+#" }}{PARA 11 "" 1 "" {XPPMATH 20 "6# >%&stepsG\"#5" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"vG-%&arrayG6%;\" \"!\"#?;F)\"#57\"" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 37 "We initializ e our numerical solution:" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 56 "for j \+ from 1 to n-1 do \n v[j,0] := u_init(j*h)\nend do:" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 32 "and set the boundary conditions:" }{MPLTEXT 1 0 0 "" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 64 "for m from 0 to steps do\n v[0,m] := 0;\n v[n,m] := 0;\nend do:" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 69 "The iteration matrix is similar to that of the simple imp licit scheme" }{MPLTEXT 1 0 0 "" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 116 "tridiag := (i,j) -> `if`(i=j,2,\n `if`(i= j+1 or i+1=j, -1,0)); \nA := Matrix(n-1,tridiag);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%(tridiagGR6$%\"iG%\"jG6\"6$%)operatorG%&arrowGF)-%#if G6%/9$9%\"\"#-F.6%5/F1,&F2\"\"\"F9F9/,&F1F9F9F9F2!\"\"\"\"!F)F)F)" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"AG-%'RTABLEG6+\"*!eG$Q\"%)anything G%'MatrixG%,rectangularG%.Fortran_orderG7\"\"\"#;\"\"\"\"#>F/" }}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 54 "However, we require an additional \+ weighting parameter " }{XPPEDIT 18 0 "alpha" "6#%&alphaG" }{TEXT -1 1 ":" }{MPLTEXT 1 0 0 "" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 50 "M := Matri x(n-1,shape=identity) + (1/2)*tau/h^2*A;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"MG-%'RTABLEG6+\"*7g.P\"%)anythingG%'MatrixG%,rectangularG%.F ortran_orderG7\"\"\"#;\"\"\"\"#>F/" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 88 "Compute the solution for the first time step (computing the rig ht hand side in b first):" }{MPLTEXT 1 0 0 "" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 163 "b := < seq( v[j,0] + (1/2)*tau/h^2*( v[j-1,0] - 2*v[ j,0] + v[j+1,0]) ,j=1..n-1) >;\ntmp := LinearSolve(M, b);\nfor j from \+ 1 to n-1 do\n v[j,1] := tmp[j]: \nend do:" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"bG-%'RTABLEG6*\"*!ozp8%)anythingG&%'VectorG6#%'colu mnG%,rectangularG%.Fortran_orderG7\"\"\"\";F1\"#>" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$tmpG-%'RTABLEG6*\"*/-RP\"%)anythingG&%'VectorG6#%'co lumnG%,rectangularG%.Fortran_orderG7\"\"\"\";F1\"#>" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 50 "And plot the results for the first two time ste ps:" }{MPLTEXT 1 0 0 "" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 109 "plot1 := listplot( [seq( v[j,0], j=0..n)]):\nplot2 := listplot( [seq( v[j,1], \+ j=0..n)]):\ndisplay(plot1,plot2);" }}{PARA 13 "" 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "6$-%'CURVESG6#777$$\"\"\"\"\"!$F*F*7$$\"\"#F*$\"+ ^YMk:!#57$$\"\"$F*$\"+W*p,4$F17$$\"\"%F*$\"+)*\\!*RXF17$$\"\"&F*$\"+CD &y(eF17$$\"\"'F*$\"+5y1rqF17$$\"\"(F*$\"+W*p,4)F17$$\"\")F*$\"+V_15*)F 17$$\"\"*F*$\"+l^c5&*F17$$\"#5F*$\"+1M)o()*F17$$\"#6F*F(7$$\"#7F*FX7$$ \"#8F*FS7$$\"#9F*FN7$$\"#:F*FI7$$\"#;F*FD7$$\"# F*F57$$\"#?F*F/7$$\"#@F*F+-F$6#77F'7$F-$\"+*Hd\"*[\"F17$F3$\"+/mkTHF17 $F8$\"+/Gq@VF17$F=$\"+QVM&f&F17$FB$\"+G*47t'F17$FG$\"+'yI8q(F17$FL$\"+ C$>=[)F17$FQ$\"+>tX`!*F17$FV$\"+\\\"p@S*F17$Fen$\"+p$o$>&*F17$FhnFar7$ F[oF^r7$F^oF[r7$FaoFhq7$FdoFeq7$FgoFbq7$FjoF_q7$F]pF\\q7$F`pFipFbp" 1 2 0 1 10 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 0 "Curve 1" "Curve 2" }}}}}{SECT 0 {PARA 3 "" 0 "" {TEXT -1 20 "The Leap-Frog Scheme" }} {PARA 0 "" 0 "" {TEXT -1 67 "Now we compute the remaining time steps u sing the leap-frog scheme." }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 169 "for m from 1 to steps-1 do\n for j from 1 to n-1 do\n v[j, m+1] := v[j,m-1] + 2*tau*(\n v[j-1,m] - 2*v[j,m] + \+ v[j+1,m] )/h^2: \n end do:\nend do:" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 128 "Finally, we plot the results. The size of the time step \+ has to bo chosen very carefully. The scheme is likely to explode . . . " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 62 "resplot := [seq( listplot( [s eq( v[j,m], j=0..n)] ), m=1..8)]:" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 33 "display(resplot,insequence=true);" }}{PARA 13 "" 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "6#-%(ANIMATEG6*7#-%'CURVESG6#777$$\"\"\"\"\"! $F.F.7$$\"\"#F.$\"+*Hd\"*[\"!#57$$\"\"$F.$\"+/mkTHF57$$\"\"%F.$\"+/Gq@ VF57$$\"\"&F.$\"+QVM&f&F57$$\"\"'F.$\"+G*47t'F57$$\"\"(F.$\"+'yI8q(F57 $$\"\")F.$\"+C$>=[)F57$$\"\"*F.$\"+>tX`!*F57$$\"#5F.$\"+\\\"p@S*F57$$ \"#6F.$\"+p$o$>&*F57$$\"#7F.Ffn7$$\"#8F.FW7$$\"#9F.FR7$$\"#:F.FM7$$\"# ;F.FH7$$\"#7$$\"#>F.F97$$\"#?F.F37$$\"#@F.F/7#-F(6#77F +7$F1$\"+uEn<9F57$F7$\"+QvV+GF57$F<$\"+DjC9TF57$FA$\"*v[nK&!\"*7$FF$\" +%*))33kF57$FK$\"*lS;L(F[r7$FP$\"+wImu!)F57$FU$\"+50')=')F57$FZ$\"+sH$ 3&*)F57$Fin$\"+KiSi!*F57$F^oFir7$FaoFfr7$FdoFcr7$FgoF`r7$FjoF]r7$F]pFi q7$F`pFfq7$FcpFcq7$FfpF`qFhp7#-F(6#77F+7$F1$\"+og_\\8F57$F7$\"+UA#em#F 57$F<$\"+BsZ;RF57$FA$\"+w`pq]F57$FF$\"+/l0+hF57$FK$\")K@zp!\")7$FP$\") %>lo(F]u7$FU$\")wb/#)F]u7$FZ$\")Bd?&)F]u7$Fin$\")AyE')F]u7$F^oFeu7$Fao Fbu7$FdoF_u7$FgoF[u7$FjoFht7$F]pFet7$F`pFbt7$FcpF_t7$FfpF\\tFhp7#-F(6# 77F+7$F1$\"+%[_ZG\"F57$F7$\");(y`#F]u7$F<$\"(*\\GP!\"(7$FA$\"(Ct#[F`w7 $FF$\"+;LF2eF57$FK$\"(PUk'F`w7$FP$\"(#f'>#Fbz7$F<$\"&oA$F\\[l7$FA$\"%aT!\"%7$FF$\"%@]Fg_l7$FK$\"$s&!\"$7 $FP$\"$K'F^`l7$FU$\"$u'F^`l7$FZ$\"$,(F^`l7$Fin$\"$4(F^`l7$F^oFf`l7$Fao Fc`l7$FdoF``l7$FgoF\\`l7$FjoFi_l7$F]pFe_l7$F`pFb_l7$FcpF__l7$FfpF\\_lF hp7#-F(6#77F+7$F1$\"&!f5F\\[l7$F7$\"$2#F^`l7$F<$\"#G!\"#7$FA$\"#TFabl7 $FF$\"$l%F^`l7$FK$\"#cFabl7$FP$\"#fFabl7$FU$\"$c'F^`l7$FZ$\"$`'F^`l7$F in$\"$#oF^`l7$F^oFbcl7$FaoF_cl7$FdoF\\cl7$FgoFibl7$FjoFfbl7$F]pFcbl7$F `pF_bl7$FcpF\\bl7$FfpFialFhp" 1 2 0 1 10 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 0 "Curve 1" }}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}}{MARK "0 0 2" 15 }{VIEWOPTS 1 1 0 3 2 1804 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }{RTABLE_HANDLES 138328580 137036012 136979680 137390204 }{RTABLE M6R0 I6RTABLE_SAVE/138328580X,%)anythingG6"6"][[[[[pdal"4"4""#!""""!F)F)F)F)F)F)F)F) F)F)F)F)F)F)F)F)F(F'F(F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F(F'F(F)F)F)F)F)F)F)F)F )F)F)F)F)F)F)F)F)F(F'F(F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F(F'F(F)F)F)F)F)F)F)F) F)F)F)F)F)F)F)F)F)F(F'F(F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F(F'F(F)F)F)F)F)F)F)F )F)F)F)F)F)F)F)F)F)F(F'F(F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F(F'F(F)F)F)F)F)F)F) F)F)F)F)F)F)F)F)F)F)F(F'F(F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F(F'F(F)F)F)F)F)F)F )F)F)F)F)F)F)F)F)F)F)F(F'F(F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F(F'F(F)F)F)F)F)F) F)F)F)F)F)F)F)F)F)F)F)F(F'F(F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F(F'F(F)F)F)F)F)F )F)F)F)F)F)F)F)F)F)F)F)F(F'F(F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F(F'F(F)F)F)F)F) F)F)F)F)F)F)F)F)F)F)F)F)F(F'F(F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F(F'F& } {RTABLE M6R0 I6RTABLE_SAVE/137036012X,%)anythingG6"6"][[[[[pdal"4"4""$!""""!F)F)F)F)F)F)F)F) F)F)F)F)F)F)F)F)F(F'F(F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F(F'F(F)F)F)F)F)F)F)F)F )F)F)F)F)F)F)F)F)F(F'F(F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F(F'F(F)F)F)F)F)F)F)F) F)F)F)F)F)F)F)F)F)F(F'F(F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F(F'F(F)F)F)F)F)F)F)F )F)F)F)F)F)F)F)F)F)F(F'F(F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F(F'F(F)F)F)F)F)F)F) F)F)F)F)F)F)F)F)F)F)F(F'F(F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F(F'F(F)F)F)F)F)F)F )F)F)F)F)F)F)F)F)F)F)F(F'F(F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F(F'F(F)F)F)F)F)F) F)F)F)F)F)F)F)F)F)F)F)F(F'F(F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F(F'F(F)F)F)F)F)F )F)F)F)F)F)F)F)F)F)F)F)F(F'F(F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F(F'F(F)F)F)F)F) F)F)F)F)F)F)F)F)F)F)F)F)F(F'F(F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F)F(F'F& } {RTABLE M6R0 I6RTABLE_SAVE/136979680X*%)anythingG6"6"\[[[[[t4"4,&-%$sinG6#,$%#PiG#""""#?!""- F)6#,$F,#F."#5F.,(F1F0F(F.-F)6#,$F,#""$F/F.,(F7F0F1F.-F)6#,$F,#F.""&F.,(F=F0F7F .*$""##F.FDFE,(FC#F0FDF=F.-F)6#,$F,#F;F5F.,(FHF0FCFE-F)6#,$F,#""(F/F.,(FMF0FHF. -F)6#,$F,#FDFAF.,(FSF0FMF.-F)6#,$F,#""*F/F.,(FXF0FSF.F.F.,&F0F.FXFDFgnFWFRFLFFF BFF2F(#"$Y#F2#"$5"F2F.F<#"$I$F2FC#"%!)>F 2FJ#F]rFQFY#!%(f""&a-$Fgn#"%q^F2F]o#FcrFQ,6F3F_pFipF.F<#"$k)F2F(F[oFJFhqFRF_qFC #"%%=&F2FY#"%&e#F2Fgn#!%"f"F2F]o#FjpFQ,6F3FgoFRFipF(FH#"$a(F2F.F<#"%iAF2FC#!%b: F2FJFdqFY#"$!**F2FgnFbsF]o#F[tFQ,6F(FaoF3FapFRFaqFhqF.F<#"$Y)FQFCFbtFJFiqFY#"#b FQFgnFhsF]o#!$(=FQ,6F(FAF3FeoFRFgp#"%o^F2F.F<#"$x$F2FCF\tFJFcqFY#"$l"F2FgnF`sF] oFdt,6#FjrF2F.F(#""%F2F3#"#7F2FR#"#KF2F<#"&O."F2FC#"%3:F2FJ#FfuFQFYFbrFgn#"$w&F 2F]o#"$k&FQFjtFctFisF_sF[rFbqFbpFboF'F& }