{VERSION 3 0 "SUN SPARC SOLARIS" "3.0" } {USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 1 0 0 0 0 }{CSTYLE "2D Math" -1 2 "Times" 0 1 0 0 0 0 0 0 2 0 0 0 0 0 0 }{CSTYLE "Hyperlink" -1 17 "" 0 1 0 128 128 1 2 0 1 0 0 0 0 0 0 } {CSTYLE "2D Output" 2 20 "" 0 1 0 0 255 1 0 0 0 0 0 0 0 0 0 }{CSTYLE " " -1 256 "" 1 12 0 0 0 0 0 2 2 0 0 0 0 0 0 }{CSTYLE "" -1 257 "" 1 24 0 1 0 0 0 0 2 0 0 0 0 0 0 }{CSTYLE "" -1 258 "" 0 1 0 1 0 0 0 0 2 0 0 0 0 0 0 }{PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "Times" 0 14 0 0 0 0 2 2 2 0 0 0 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Text Ou tput" -1 2 1 {CSTYLE "" -1 -1 "Courier" 1 10 0 0 255 1 0 0 0 0 0 1 3 0 0 }1 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Heading 1" 0 3 1 {CSTYLE "" -1 -1 "" 1 18 0 0 0 0 0 1 0 0 0 0 0 0 0 }1 0 0 0 6 6 0 0 0 0 0 0 -1 0 }{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 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Maple Out put" 0 11 1 {CSTYLE "" -1 -1 "" 0 1 240 136 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 "" 11 12 1 {CSTYLE "" -1 -1 "" 0 1 58 61 32 0 0 0 0 0 0 0 0 0 0 }1 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 } {PSTYLE "Maple Plot" 0 13 1 {CSTYLE "" -1 -1 "" 0 1 9 137 40 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" 0 18 1 {CSTYLE "" -1 -1 "" 1 18 40 32 65 0 0 1 1 0 0 0 0 0 0 }3 0 0 -1 12 12 0 0 0 0 0 0 19 0 }{PSTYLE "Author" 0 19 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }3 0 0 -1 8 8 0 0 0 0 0 0 -1 0 }{PSTYLE "R3 Font 0" -1 256 1 {CSTYLE "" -1 -1 "Helvetica" 0 14 16 0 0 0 2 1 2 0 0 0 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "R3 Font 2" -1 257 1 {CSTYLE "" -1 -1 "Courier" 0 14 0 0 132 0 2 2 2 0 0 0 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "" 18 258 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }2 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "" 11 259 1 {CSTYLE "" -1 -1 "" 0 1 58 61 32 0 0 0 0 0 0 0 0 0 0 }1 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }} {SECT 0 {SECT 1 {PARA 258 "" 0 "" {TEXT 257 40 "Chapter 19\n\nLinear A lgebra: Applications" }{TEXT 258 1 "\n" }}{PARA 19 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT 256 65 "\251 Copyright 1997, Latin American Maple Center, Petr\363polis, Brasil" }}{PARA 0 "" 0 "" {HYPERLNK 17 " About the Worksheets" 1 "about.mws" "" }}}{SECT 1 {PARA 3 "" 0 "Kinema tics" {TEXT -1 43 "19.1 Kinematics of the Stanford Manipulator" }} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "restart; with(linalg):" }} {PARA 7 "" 1 "" {TEXT -1 32 "Warning, new definition for norm" }} {PARA 7 "" 1 "" {TEXT -1 33 "Warning, new definition for trace" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 44 "alias( seq( c[i] = cos(theta [i]), i=1..6 ), " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 40 " seq( s[i] = s in(theta[i]), i=1..6 ) ):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 52 "M := (a,alpha,d,theta) -> matrix( 4, 4, [cos(theta)," }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 48 " -sin(theta)*cos(alpha), sin(theta)*sin(alph a)," }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 50 " a*cos(theta), sin(theta), \+ cos(theta)*cos(alpha)," }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 54 " -cos(th eta)*sin(alpha), a*sin(theta), 0, sin(alpha)," }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 32 " cos(alpha), d, 0, 0, 0, 1 ] ):" }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 20 "M(a,alpha,d,theta); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7&7&-%$cosG6#%&thetaG,$*&-%$sinGF*\"\"\"-F) 6#%&alphaGF0!\"\"*&F.\"\"\"-F/F2F0*&%\"aGF0F(F07&F.*&F(F6F1F6,$*&F(F6F 7F6F4*&F9F6F.F67&\"\"!F7F1%\"dG7&F@F@F@F0" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 47 "a := vector([0$6]): " }} }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 46 "alpha := vector([-Pi/2,Pi/2 ,0,-Pi/2,Pi/2,0]): " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "d := vector([0$6]): " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 41 "d[2] := evaln(d[2]): d[3] := evaln(d[3]):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "theta := vector(6): theta[3] := 0:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 15 "for i to 6 do " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 47 " A[i] := M( a[i], alpha[i], d[i], theta[i] ) " }} {PARA 0 "> " 0 "" {MPLTEXT 1 0 5 "od; " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>&%\"AG6#\"\"\"-%'matrixG6#7&7&&%\"cGF&\"\"!,$&%\"sGF&!\"\"F/7&F 1F/F-F/7&F/F3F/F/7&F/F/F/F'" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>&%\"AG 6#\"\"#-%'matrixG6#7&7&&%\"cGF&\"\"!&%\"sGF&F/7&F0F/,$F-!\"\"F/7&F/\" \"\"F/&%\"dGF&7&F/F/F/F6" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>&%\"AG6# \"\"$-%'matrixG6#7&7&\"\"\"\"\"!F.F.7&F.F-F.F.7&F.F.F-&%\"dGF&7&F.F.F. F-" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>&%\"AG6#\"\"%-%'matrixG6#7&7&&% \"cGF&\"\"!,$&%\"sGF&!\"\"F/7&F1F/F-F/7&F/F3F/F/7&F/F/F/\"\"\"" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>&%\"AG6#\"\"&-%'matrixG6#7&7&&%\"cGF& \"\"!&%\"sGF&F/7&F0F/,$F-!\"\"F/7&F/\"\"\"F/F/7&F/F/F/F6" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#>&%\"AG6#\"\"'-%'matrixG6#7&7&&%\"cGF&,$&%\"sGF& !\"\"\"\"!F37&F0F-F3F37&F3F3\"\"\"F37&F3F3F3F6" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 46 "Tip := evalm( `&*`( seq( A[i], i=1..6 ) ) ): \+ " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 40 "collect( Tip[1,1], [c[1 ],c[2],s[1]] ); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,&*&,&*&,&*(&%\"c G6#\"\"%\"\"\"&F*6#\"\"&F-&F*6#\"\"'F-F-*&&%\"sGF+F-&F6F2F-!\"\"F-&F*6 #\"\"#F-F-*(&F6F:F-&F6F/F-F1\"\"\"F8F-&F*6#F-F-F-*&,&*(F5F?F.F?F1F?F8* &F)F?F7F?F8F-&F6FAF-F-" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 32 "T := subvector( Tip, 1..3, 4 ); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>% \"TG-%'vectorG6#7%,&*(&%\"cG6#\"\"\"F.&%\"sG6#\"\"#F.&%\"dG6#\"\"$F.F. *&&F0F-F.&F4F1F.!\"\",&*(F8\"\"\"F/F=F3F=F.*&F+F=F9F=F.*&&F,F1F.F3F=" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 36 "R := submatrix( Tip, 1..3 , 1..3 ): " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "Jv := jacobi an( T, " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 53 " [theta[1],theta[2],the ta[4],theta[5],theta[6]] ); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#Jv G-%'matrixG6#7%7',&*(&%\"sG6#\"\"\"F/&F-6#\"\"#F/&%\"dG6#\"\"$F/!\"\"* &&%\"cGF.F/&F4F1F/F7*(F9\"\"\"&F:F1F/F3F=\"\"!F?F?7',&*(F9F=F0F=F3F=F/ *&F,F=F;F=F7*(F,F=F>F=F3F=F?F?F?7'F?,$*&F0F=F3F=F7F?F?F?" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "R1 := map( diff, R, theta[1] ): \+ " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 31 "R2 := map( diff, R, the ta[2] ):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "R3 := map( diff , R, theta[4] ): " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 31 "R4 : = map( diff, R, theta[5] ):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 33 "R5 := map( diff, R, theta[6] ): " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 31 "Rtranspose := transpose( R ): " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 58 "evalm( R1 &* Rtranspose ): omega := map( sim plify, % ); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%&omegaG-%'matrixG6# 7%7%\"\"!!\"\"F*7%\"\"\"F*F*7%F*F*F*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 59 "Jomega[1] := vector([omega[3,2],-omega[3,1],omega[2,1 ]]); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>&%'JomegaG6#\"\"\"-%'vector G6#7%\"\"!F,F'" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "for i fro m 2 to 5 do" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 30 " evalm( R.i &* Rtra nspose ): " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 30 " omega := map( simpl ify, % ):" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 50 " Jomega[i] := vector( [ omega[3,2], -omega[3,1], " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 19 " \+ omega[2,1] ] ):" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 36 " print( Jomega [i]=eval(Jomega[i]) )" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 5 "od: " }} {PARA 11 "" 1 "" {XPPMATH 20 "6#/&%'JomegaG6#\"\"#-%'vectorG6#7%,$&%\" sG6#\"\"\"!\"\"&%\"cGF/\"\"!" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/&%'Jo megaG6#\"\"$-%'vectorG6#7%*&&%\"sG6#\"\"#\"\"\"&%\"cG6#F1F1*&F-\"\"\"& F.F4F1&F3F/" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/&%'JomegaG6#\"\"%-%'ve ctorG6#7%,&*(&%\"cG6#\"\"\"F1&F/6#\"\"#F1&%\"sGF&F1!\"\"*&&F6F0F1&F/F& F1F7,&*(F9\"\"\"F2F=F5F=F7*&F.F=F:F=F1*&&F6F3F1F5F=" }}{PARA 12 "" 1 " " {XPPMATH 20 "6#/&%'JomegaG6#\"\"&-%'vectorG6#7%,(**&%\"sGF&\"\"\"&% \"cG6#F0F0&F26#\"\"#F0&F26#\"\"%F0F0*(F1\"\"\"&F/F5F0&F2F&F0F0*(F.F;&F /F3F0&F/F8F0!\"\",(*(F.F;F1F;F@F;F0*(F?F;F " 0 "" {MPLTEXT 1 0 53 "Jomega := concat( seq( eval(Jomega[i]), i=1..5 ) ): " }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 33 "J := stackmatrix( Jv, Jomega ); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"JG-%'matrixG6#7(7',&*(&% \"sG6#\"\"\"F/&F-6#\"\"#F/&%\"dG6#\"\"$F/!\"\"*&&%\"cGF.F/&F4F1F/F7*(F 9\"\"\"&F:F1F/F3F=\"\"!F?F?7',&*(F9F=F0F=F3F=F/*&F,F=F;F=F7*(F,F=F>F=F 3F=F?F?F?7'F?,$*&F0F=F3F=F7F?F?F?7'F?,$F,F7*&F0F=F9F=,&*(F9F=F>F=&F-6# \"\"%F/F7*&F,F=&F:FNF/F7,(**&F-6#\"\"&F/F9F=F>F=FQF=F/*(F9F=F0F=&F:FUF /F/*(FTF=F,F=FMF=F77'F?F9*&F0F=F,F=,&*(F,F=F>F=FMF=F7*&F9F=FQF=F/,(*(F TF=F9F=FMF=F/*(F,F=F0F=FXF=F/**FTF=F,F=F>F=FQF=F/7'F/F?F>*&F0F=FMF=,&* (F0F=FQF=FTF=F7*&F>F=FXF=F/" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "print( Jv = eval(Jv) ); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/%#Jv G-%'matrixG6#7%7',&*(&%\"sG6#\"\"\"F/&F-6#\"\"#F/&%\"dG6#\"\"$F/!\"\"* &&%\"cGF.F/&F4F1F/F7*(F9\"\"\"&F:F1F/F3F=\"\"!F?F?7',&*(F9F=F0F=F3F=F/ *&F,F=F;F=F7*(F,F=F>F=F3F=F?F?F?7'F?,$*&F0F=F3F=F7F?F?F?" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 27 "subs( s[2]=0, eval(Jv) ); " }} {PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7%7',$*&&%\"cG6#\"\"\"F-&% \"dG6#\"\"#F-!\"\"*(F*\"\"\"&F+F0F-&F/6#\"\"$F-\"\"!F9F97',$*&&%\"sGF, F-F.F4F2*(F=F4F5F4F6F4F9F9F97'F9F9F9F9F9" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 29 "submatrix( %, 1..2, 1..2 ); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7$7$,$*&&%\"cG6#\"\"\"F-&%\"dG6#\"\"#F-!\" \"*(F*\"\"\"&F+F0F-&F/6#\"\"$F-7$,$*&&%\"sGF,F-F.F4F2*(F " 0 "" {MPLTEXT 1 0 9 "det(%); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#\"\"!" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 36 " addrow( Jv, 3, 1, c[2]*c[1]/s[2] ); " }}{PARA 11 "" 1 "" {XPPMATH 20 " 6#-%'matrixG6#7%7',&*(&%\"sG6#\"\"\"F-&F+6#\"\"#F-&%\"dG6#\"\"$F-!\"\" *&&%\"cGF,F-&F2F/F-F5\"\"!F:F:F:7',&*(F7\"\"\"F.F>F1F>F-*&F*F>F9F>F5*( F*F>&F8F/F-F1F>F:F:F:7'F:,$*&F.F>F1F>F5F:F:F:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 36 "addrow( %, 3, 2, s[1]*c[2]/s[2] ); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7%7',&*(&%\"sG6#\"\"\"F-&F+6#\" \"#F-&%\"dG6#\"\"$F-!\"\"*&&%\"cGF,F-&F2F/F-F5\"\"!F:F:F:7',&*(F7\"\" \"F.F>F1F>F-*&F*F>F9F>F5F:F:F:F:7'F:,$*&F.F>F1F>F5F:F:F:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "\{ %[1,1]=0, %[2,1]=0 \}; " }} {PARA 11 "" 1 "" {XPPMATH 20 "6#<$/,&*(&%\"cG6#\"\"\"F*&%\"sG6#\"\"#F* &%\"dG6#\"\"$F*F**&&F,F)F*&F0F-F*!\"\"\"\"!/,&*(F4\"\"\"F+F;F/F;F6*&F' F;F5F;F6F7" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "solve( %, \{s [1],c[1]\} ); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#<$/&%\"sG6#\"\"\"\" \"!/&%\"cGF'F)" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}} {SECT 1 {PARA 3 "" 0 "Cadmium" {TEXT -1 46 "19.2 A 3-Compartment Model of Cadmium Transfer" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "rest art; with(linalg):" }}{PARA 7 "" 1 "" {TEXT -1 32 "Warning, new defini tion for norm" }}{PARA 7 "" 1 "" {TEXT -1 33 "Warning, new definition \+ for trace" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "sequence := " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 46 " seq(seq( a.i.j = a[i,j], i=0..3 ), j=0..3 )," }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 46 " seq(seq( b.i.j = b[i,j], i=0..3 ), j=0..3 )," }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 57 " s eq(seq( t.i.j = t[i,j], i=0..3 ), j=0..3 ), c3 = c[3]:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 38 "(eval@subs)( S=sequence, 'macro(S)' ):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 40 "A := matrix( 3, 3, [ -a21, a12, 0, a21, " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 48 " -(a02+a12 +a32), a23, 0, a32, -(a03+a23) ] ); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"AG-%'matrixG6#7%7%,$&%\"aG6$\"\"#\"\"\"!\"\"&F,6$F/F.\"\"!7%F+, (&F,6$F3F.F0F1F0&F,6$\"\"$F.F0&F,6$F.F:7%F3F8,&&F,6$F3F:F0F;F0" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 31 "B := matrix( 3, 1, [0,1,0] ) ; " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"BG-%'matrixG6#7%7#\"\"!7#\" \"\"F)" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 42 "C := matrix( 2, 3 , [0,0,c3,0,a02,a03] ); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"CG-%' matrixG6#7$7%\"\"!F*&%\"cG6#\"\"$7%F*&%\"aG6$F*\"\"#&F16$F*F." }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 49 "AB := evalm( A &* B ): A2B := evalm( A &* AB ):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 30 "MC := augment( B, AB, A2B ); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#MCG -%'matrixG6#7%7%\"\"!&%\"aG6$\"\"\"\"\"#,&*&&F,6$F/F.F.F+F.!\"\"*&F+\" \"\",(&F,6$F*F/F4F+F4&F,6$\"\"$F/F4F.F.7%F.F7,(F1F.*$)F7F/F6F.*&&F,6$F /F " 0 "" {MPLTEXT 1 0 12 "rank( MC ); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#\"\"$" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "ffg ausselim( MC ); " }}{PARA 12 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7%7% \"\"\",(&%\"aG6$\"\"!\"\"#!\"\"&F+6$F(F.F/&F+6$\"\"$F.F/,2*&&F+6$F.F(F (F0F(F(*$)F*F.\"\"\"F(*&F0F;F*F(F.*&F*F;F2F(F.*$)F0F.F;F(*&F0F;F2F;F.* $)F2F.F;F(*&&F+6$F.F4F(F2F;F(7%F-F0,*F6F/FF/F@F/7%F-F-,(*(F0F;F2F ;&F+6$F-F4F(F/*(F0F;FDF;F2F;F/*(F2F;F7F;F0F;F(" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 19 "map( factor, % ); " }}{PARA 12 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7%7%\"\"\",(&%\"aG6$\"\"!\"\"#!\"\"&F+6$F(F .F/&F+6$\"\"$F.F/,2*&&F+6$F.F(F(F0F(F(*$)F*F.\"\"\"F(*&F0F;F*F(F.*&F*F ;F2F(F.*$)F0F.F;F(*&F0F;F2F;F.*$)F2F.F;F(*&&F+6$F.F4F(F2F;F(7%F-F0,$*& F0F;,*F0F(F7F(F*F(F2F(F(F/7%F-F-*(F0F;F2F;,(&F+6$F-F4F/FDF/F7F(F(" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "factor( det( MC ) ); " }} {PARA 11 "" 1 "" {XPPMATH 20 "6#,$*(&%\"aG6$\"\"\"\"\"#F(&F&6$\"\"$F)F (,(&F&6$\"\"!F,!\"\"&F&6$F)F,F1&F&6$F)F(F(F(F1" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 22 "CA := evalm( C &* A ):" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 29 "MO := stackmatrix( C, CA ); " }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%#MOG-%'matrixG6#7&7%\"\"!F*&%\"cG6#\"\"$7%F*&%\" aG6$F*\"\"#&F16$F*F.7%F**&F+\"\"\"&F16$F.F3F8*&F+\"\"\",&F4!\"\"&F16$F 3F.F>F87%*&F0F8&F16$F3F8F8,&*&F0F<,(F0F>&F16$F8F3F>F9F>F8F8*&F9F " 0 "" {MPLTEXT 1 0 20 "ffgausselim( MO ); " }}{PARA 12 "" 1 "" {XPPMATH 20 "6#-%'matrixG 6#7&7%*&&%\"aG6$\"\"!\"\"#\"\"\"&F*6$F-F.F.,**$)F)F-\"\"\"!\"\"*&&F*6$ F.F-F.F)F4F5*&F)F4&F*6$\"\"$F-F.F5*&F:F4&F*6$F,FF-F4F5*&F>F4FBF4F57%F,*&F3F4F/F4*(F>F4F)F4F/F47%F,F,*(&%\"cG6 #F " 0 "" {MPLTEXT 1 0 39 "fa ctor( det( transpose(MO) &* MO ) ); " }}{PARA 12 "" 1 "" {XPPMATH 20 "6#**)&%\"cG6#\"\"$\"\"#\"\"\")&%\"aG6$F)\"\"\"F)F*)&F-6$\"\"!F)F)F*,2 *$F0F*F/*&)&F-6$F(F)F)F*)&F-6$F3F(F)F*F/*(F0F*F;F/&F-6$F)F(F/F)*&F$F*F 7F*F/*&F0F*)F>F)F*F/*&F0F*F:F*F/**F1F/F8F/F;F*F>F*F)*(F1F*F8F*F:F*F)F/ " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "collect( op(4,%), a32 ) ; " }}{PARA 12 "" 1 "" {XPPMATH 20 "6#,.*&,&*$)&%\"aG6$\"\"!\"\"$\"\"# \"\"\"\"\"\"*$)&%\"cG6#F,F-F.F/F/)&F)6$F,F-F-F.F/*&,&*&&F)6$F+F-F/F'F. F-*(F;F.F(F/&F)6$F-F,F/F-F/F6F/F/*$)F;F-F.F/*&FAF.F'F.F/*(FAF.F(F.F>F. F-*&FAF.)F>F-F.F/" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "discri m(%,a32); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,$*&,,*$)&%\"cG6#\"\"$ \"\"#\"\"\"\"\"\"*(F'F-&%\"aG6$\"\"!F+F.&F16$F,F+F.F,*&F'F-)F4F,F-F.*& F'F-)F0F,F-F.*$F9F-F.F.)&F16$F3F,F,F-!\"%" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 17 "n := rowdim( A );" }}{PARA 11 "" 1 "" {XPPMATH 20 " 6#>%\"nG\"\"$" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "CA0B := ev alm( C&*B ): " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 58 "for i from 1 to 5 do CA.i.B := evalm( C &* A^i &* B ) od: " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 45 "J := stackmatrix( seq(CA.i.B, i=0..2*n-1) ): \+ " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 55 "J12 := map( diff, J, a1 2): J21 := map( diff, J, a21): " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 55 "J23 := map( diff, J, a23): J32 := map( diff, J, a32) : " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 55 "J02 := map( diff, J, \+ a02): J03 := map( diff, J, a03): " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 47 "JM := augment( J21, J12, J02, J23, J32, J03 ): " }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "ffgausselim( JM ): " }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "map( factor, % ); " }}{PARA 12 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7.7(*&&%\"aG6$\"\"\"\"\"#F,&F*6$ \"\"!F-F,,,*&F.\"\"\"&F*6$F-F,F,F,*$)F.F-F3F-F(F-*&F.F3&F*6$\"\"$F-F,F -*&F9F3&F*6$F0F;F,!\"\",4*&F4F3F)F3F,F6F;F(\"\"%F8FB*$)F)F-F3F,*&F)F3F 9F3F-*$)F9F-F3F,*&&F*6$F-F;F,F9F3F,F " 0 "" {MPLTEXT 1 0 35 "H := evalm( C &* (s-A)^(-1) &* B ):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 56 "evalm( denom(H[1,1]) * H ) / collect(denom(H[1,1]),s); " }} {PARA 12 "" 1 "" {XPPMATH 20 "6#*&-%'matrixG6#7$7#*(&%\"cG6#\"\"$\"\" \",&&%\"aG6$\"\"#F.F.%\"sGF.F.&F16$F-F3F.7#*&,H*(F0F.&F16$\"\"!F-F.&F1 6$F=F3F.F.*(F4F.F>\"\"\"&F16$F3F-F.F.*(F;FAF5FAF0FAF.*(F0FAF5FAF4FAF.* (F4FAFBFA&F16$F.F3F.F.*(F0FAF4FAF;FAF.*(F0FAF4FAFBFAF.*(F4FAF;FAF>FAF. *(F4FAF;FAFGFAF.*&F0FA)F4F3FAF.*&F>FAFNFAF.*&FGFAFNFAF.*&F5FAFNFAF.*&F NFAF;FAF.*&FNFAFBFAF.*$)F4F-FAF.*(F4FAF5FAF;FAF.*(F0FAF>FAF4FAF.*(F0FA F>FAFBFAF.F.,&*&*(F>FAF/FA,(F;F.FBF.F4F.F.FAF9!\"\"F.*&*(F;FAF/FAF5FAF AF9FgnF.F.FA,.FTF.*&,.FGF.F0F.F>F.F5F.F;F.FBF.F.FNFAF.*&,4*&F5FAF0FAF. *&FBFAFGFAF.*&F>FAFBFAF.*&FBFAF0FAF.*&F;FAF>FAF.*&F;FAFGFAF.*&F>FAF0FA F.*&F;FAF0FAF.*&F5FAF;FAF.F.F4FAF.F:F.FXF.FDF.Fgn" }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 30 "collect( numer(H[1,1]), s ); " }}{PARA 11 " " 1 "" {XPPMATH 20 "6#,&*(&%\"cG6#\"\"$\"\"\"&%\"aG6$F(\"\"#F)%\"sGF)F )*(F%\"\"\"F*F0&F+6$F-F)F)F)" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 30 "collect( numer(H[2,1]), s ); " }}{PARA 12 "" 1 "" {XPPMATH 20 " 6#,(*&&%\"aG6$\"\"!\"\"#\"\"\")%\"sGF)\"\"\"F**&,**&F%F-&F&6$F)F*F*F** &&F&6$F(\"\"$F*F%F-F**&F%F-&F&6$F)F6F*F**&&F&6$F6F)F*F4F-F*F*F,F*F**&F 1F-,(F3F*F7F*F:F*F*F*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 49 "cf s1 := \{ coeffs( expand(numer(H[1,1])), s ) \}; " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%cfs1G<$*(&%\"cG6#\"\"$\"\"\"&%\"aG6$F*\"\"#F+&F-6$F/ F+F+*&F'\"\"\"F,F3" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 49 "cfs2 \+ := \{ coeffs( expand(numer(H[2,1])), s ) \}; " }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%%cfs2G<%,**&&%\"aG6$\"\"!\"\"#\"\"\"&F)6$F,F-F-F-*&&F )6$F+\"\"$F-F(\"\"\"F-*&F(F4&F)6$F,F3F-F-*&&F)6$F3F,F-F1F4F-,(*(F.F4F1 F4F(F4F-*(F.F4F(F4F6F4F-*(F1F4F9F4F.F4F-F(" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 48 "cfs3 := \{ coeffs( expand(denom(H[1,1])), s ) \} \+ " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 14 " minus \{1\}; " }}{PARA 12 " " 1 "" {XPPMATH 20 "6#>%%cfs3G<%,4*&&%\"aG6$\"\"$\"\"#\"\"\"&F)6$F,F-F -F-*&&F)6$F,F+F-&F)6$F-F,F-F-*&&F)6$\"\"!F,F-F1\"\"\"F-*&F1F9F.F9F-*&& F)6$F8F+F-F6F9F-*&F " 0 "" {MPLTEXT 1 0 30 "cfs := `union`( cfs.(1. .3) ): " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 56 "eqns := map( x - > x = subs( a12=b12, a21=b21, a23=b23, " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 51 " a32=b32, a23=b23, a03=b03, a02=b02, x ), cfs ); " }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%%eqnsG<)/*&&%\"cG6#\"\"$\"\"\"&%\" aG6$F+\"\"#F,*&F(\"\"\"&%\"bGF/F,/&F.6$\"\"!F0&F4F7/,4*&F-F2&F.6$F0F,F ,F,*&&F.6$F0F+F,&F.6$F,F0F,F,*&F6F,F@F2F,*&F@F2F=F2F,*&&F.6$F8F+F,F6F2 F,*&FGF2FBF2F,*&F6F2F=F2F,*&FGF2F=F2F,*&F-F2FGF2F,,4*&F3F2&F4F>F,F,*&& F4FAF,&F4FCF,F,*&F9F,FQF2F,*&FQF2FOF2F,*&&F4FHF,F9F2F,*&FVF2FRF2F,*&F9 F2FOF2F,*&FVF2FOF2F,*&F3F2FVF2F,/,.FBF,F=F,F6F,F-F,FGF,F@F,,.FRF,FOF,F 9F,F3F,FVF,FQF,/,*FJF,FFF,FDF,FLF,,*FXF,FUF,FSF,FZF,/,(*(F=F2FGF2F6F2F ,*(F=F2F6F2F@F2F,*(FGF2F-F2F=F2F,,(*(FOF2FVF2F9F2F,*(FOF2F9F2FQF2F,*(F VF2F3F2FOF2F,/*(F(F2F-F2F=F2*(F(F2F3F2FOF2" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 35 "vars := \{b12,b21,b23,b32,b02,b03\}: " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "solve(eqns,vars); " }}{PARA 12 "" 1 "" {XPPMATH 20 "6$<(/&%\"bG6$\"\"!\"\"#&%\"aGF'/&F&6$F(\"\"$&F+F./&F &6$F/F)&F+F3/&F&6$\"\"\"F)&F+F7/&F&6$F)F/&F+F" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 72 "AA \+ := subs( a12=b12, a21=b21, a23=b23, a32=b32, \+ " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 40 " a23=b23, a03=b03, a02=b02, ev al(A) ); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#AAG-%'matrixG6#7%7%,$& %\"bG6$\"\"#\"\"\"!\"\"&F,6$F/F.\"\"!7%F+,(&F,6$F3F.F0F1F0&F,6$\"\"$F. F0&F,6$F.F:7%F3F8,&&F,6$F3F:F0F;F0" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "BB := copy(B):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 42 "CC := subs( a02=b02, a03=b03, eval(C) ); " }}{PARA 11 "" 1 " " {XPPMATH 20 "6#>%#CCG-%'matrixG6#7$7%\"\"!F*&%\"cG6#\"\"$7%F*&%\"bG6 $F*\"\"#&F16$F*F." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 37 "T := m atrix( 3, 3, (i,j)->t[i,j] ); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>% \"TG-%'matrixG6#7%7%&%\"tG6$\"\"\"F-&F+6$F-\"\"#&F+6$F-\"\"$7%&F+6$F0F -&F+6$F0F0&F+6$F0F37%&F+6$F3F-&F+6$F3F0&F+6$F3F3" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "evalm( T&*BB - B ); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7%7#&%\"tG6$\"\"\"\"\"#7#,&&F)6$F,F,F+!\"\" F+7#&F)6$\"\"$F," }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 30 "t12 := \+ 0: t32:= 0: t22 := 1:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 " evalm( CC - C&*T ); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7 $7%,$*&&%\"cG6#\"\"$\"\"\"&%\"tG6$F-F.F.!\"\"\"\"!,&F*F.*&F*\"\"\"&F06 $F-F-F.F27%,&*&&%\"aG6$F3\"\"#F.&F06$F?F.F.F2*&&F=6$F3F-F.F/F6F2,&FF.,(&FGFDF.*&F " 0 "" {MPLTEXT 1 0 19 "map( factor, % ); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7$7%,$*&&%\"cG6#\"\"$\"\"\"&%\"tG6$F-F.F.! \"\"\"\"!,$*&F*\"\"\",&F2F.&F06$F-F-F.F.F27%,&*&&%\"aG6$F3\"\"#F.&F06$ F@F.F.F2*&&F>6$F3F-F.F/F6F2,&F=F2&%\"bGF?F.,(&FHFEF.*&F=F6&F06$F@F-F.F 2*&FDF6F8F.F2" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 33 "t31 := 0: \+ t33 := 1: b02 := a02:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 40 " evalm( CC - C&*T ): map( factor, % ); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7$7%\"\"!F(F(7%,$*&&%\"aG6$F(\"\"#\"\"\"&%\"tG6$F/F0 F0!\"\",&F,F4&%\"bGF.F0,(&F76$F(\"\"$F0*&F,\"\"\"&F26$F/F;F0F4&F-F:F4 " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 9 "t21 := 0:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 32 "eqn1 := evalm(CC - C&*T)[2,3]; " } }{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%eqn1G,(&%\"bG6$\"\"!\"\"$\"\"\"*&& %\"aG6$F)\"\"#F+&%\"tG6$F0F*F+!\"\"&F.F(F4" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 24 "evalm( T&*AA - A&*T ); " }}{PARA 12 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7%7%,&*&&%\"tG6$\"\"\"F-F-&%\"bG6$\"\"#F-F- !\"\"*&&%\"aGF0F-F*\"\"\"F-,(*&F*F6&F/6$F-F1F-F-*&&F+6$F-\"\"$F-&F/6$F >F1F-F-&F5F:F2,(*&FF2&F/6$F1F>F2F-F-*&F4F6F " 0 "" {MPLTEXT 1 0 42 "e qn2 := map( op, convert(%,'listlist') ): " }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 39 "eqns := \{ eqn1, op(eqn2) \} minus \{0\}; " }} {PARA 12 "" 1 "" {XPPMATH 20 "6#>%%eqnsG<+,,&%\"bG6$\"\"#\"\"$\"\"\"*& &%\"tGF)F,,&&F(6$\"\"!F+!\"\"F'F4F,F,*&&%\"aG6$F*F,F,&F/6$F,F+F,F4*&,( &F76$F3F*F4&F76$F,F*F4&F76$F+F*F4F,F.\"\"\"F4&F7F)F4,&FAF4&F(FBF,,&*&& F/6$F,F,F,&F(F8F,F4*&F6FCFIFCF,,&FKF,FLF4,(*&FIFC&F(F@F,F,*&F9FCFFF,F, F?F4,(*&F9FCF0FCF,F5F,*&F?F,F.FCF4,,FPF4FFF4*&F.FCFFFCF,F?F,FAF,,,F1F4 F'F4*&FAF,F.FCF4&F7F2F,FDF,,(F1F,*&F=F,F.FCF4FYF4" }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 50 "solve( eqns, \{b12,b32,b23,b03,b21,t11,t13,t 23\} ); " }}{PARA 12 "" 1 "" {XPPMATH 20 "6$<*/&%\"tG6$\"\"\"\"\"$\"\" !/&F&6$\"\"#F)F*/&%\"bG6$F)F.&%\"aGF2/&F16$F.F(&F4F7/&F16$F*F)&F4F;/&F 16$F(F.&F4F?/&F1F-&F4F-/&F&6$F(F(F(<*/F,,$*&,*F,(F " 0 "" {MPLTEXT 1 0 22 "restart; with(linalg):" }}{PARA 7 " " 1 "" {TEXT -1 32 "Warning, new definition for norm" }}{PARA 7 "" 1 " " {TEXT -1 33 "Warning, new definition for trace" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "`number of C-atoms` := 10:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "Azulene_skeleton := " }}{PARA 0 "> " 0 " " {MPLTEXT 1 0 47 " [ [1,2], [1,9], [2,3], [3,10], [4,5], [4,10]," }} {PARA 0 "> " 0 "" {MPLTEXT 1 0 51 " [4,10], [5,6], [6,7], [7,8], [8 ,9], [9,10] ]: " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 32 "Hueckel _matrix := array( sparse," }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 51 " 1..` number of C-atoms`, 1..`number of C-atoms` ):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 43 "for i from 1 to nops( Azulene_skeleton ) do" }} {PARA 0 "> " 0 "" {MPLTEXT 1 0 53 " Hueckel_matrix[ op( Azulene_skele ton[i] ) ] := beta" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 4 "od: " }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 49 "Hueckel_matrix := evalm( alp ha + Hueckel_matrix " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 35 " + transp ose( Hueckel_matrix ) ): " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 46 "topological_matrix := subs( alpha=0, beta=1, " }}{PARA 0 "> " 0 " " {MPLTEXT 1 0 29 " eval( Hueckel_matrix ) ): " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 46 "factor( charpoly( topological_matrix, x ) ); \+ " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#*&,,*$)%\"xG\"\"%\"\"\"\"\"\"*$)F '\"\"$F)F**$)F'\"\"#F)!\"$F'!\"\"F*F*F*,0*$)F'\"\"'F)F**$)F'\"\"&F)F2F %!\"(F+F9F.\"#8F'!\"'!\"%F*F*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 35 "charmat( topological_matrix, x ); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7,7,%\"xG!\"\"\"\"!F*F*F*F*F*F)F*7,F)F(F)F* F*F*F*F*F*F*7,F*F)F(F*F*F*F*F*F*F)7,F*F*F*F(F)F*F*F*F*F)7,F*F*F*F)F(F) F*F*F*F*7,F*F*F*F*F)F(F)F*F*F*7,F*F*F*F*F*F)F(F)F*F*7,F*F*F*F*F*F*F)F( F)F*7,F)F*F*F*F*F*F*F)F(F)7,F*F*F)F)F*F*F*F*F)F(" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "factor( det(%) ); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#*&,,*$)%\"xG\"\"%\"\"\"\"\"\"*$)F'\"\"$F)F**$)F'\"\"#F) !\"$F'!\"\"F*F*F*,0*$)F'\"\"'F)F**$)F'\"\"&F)F2F%!\"(F+F9F.\"#8F'!\"'! \"%F*F*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 36 "eigenvalues( top ological_matrix ); " }}{PARA 12 "" 1 "" {XPPMATH 20 "6$-%'RootOfG6#,, *$)%#_ZG\"\"%\"\"\"\"\"\"*$)F)\"\"$F+F,*$)F)\"\"#F+!\"$F)!\"\"F,F,-F$6 #,0*$)F)\"\"'F+F,*$)F)\"\"&F+F4F'!\"(F-F=F0\"#8F)!\"'!\"%F," }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 37 "eigenvectors( topological_ma trix ); " }}{PARA 12 "" 1 "" {XPPMATH 20 "6$7%-%'RootOfG6#,0*$)%#_ZG \"\"'\"\"\"\"\"\"*$)F*\"\"&F,!\"\"*$)F*\"\"%F,!\"(*$)F*\"\"$F,F0*$)F* \"\"#F,\"#8F*!\"'!\"%F-F-<#-%'vectorG6#7,F-,.F$#F " 0 "" {MPLTEXT 1 0 7 "%[2]; " }}{PARA 12 "" 1 "" {XPPMATH 20 "6#7%-%'RootOf G6#,,*$)%#_ZG\"\"%\"\"\"\"\"\"*$)F*\"\"$F,F-*$)F*\"\"#F,!\"$F*!\"\"F-F -F-<#-%'vectorG6#7,,*F4F-*$)F$F3F,F3F$!\"#*$)F$F0F,F-\"\"!,*F0F-F$F3F? F5FF5,*F$F4F?F-FF?F-F " 0 "" {MPLTEXT 1 0 31 "allvalues (%,'dependent')[1]; " }}{PARA 12 "" 1 "" {XPPMATH 20 "6#7%,(#!\"\"\" \"%\"\"\"*$-%%sqrtG6#\"\"&\"\"\"#F(F'*$-F+6#,&\"#AF(F)!\"#F.F/F(<#-%'v ectorG6#7,,,#!\"&\"\"#F(*$)F$F>F.F>F)#F&F>F0FA*$)F$\"\"$F.F(\"\"!,,#F- F>F(F)#F(F>F0FHFBF&F?F5F&,,F%F(F)#!\"$F'F0FJFBF(F?F(FE,,F/F(F)#FDF'F0F MF?F&FBF&F(,,FAF(F)FAF0FAFBF(F?F(,,FHF(F)FHF0FHFBF&F?F&" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "eigenvalues := evalf( " }}{PARA 0 " > " 0 "" {MPLTEXT 1 0 54 " Eigenvals( topological_matrix, 'eigenvecto rs' ) ); " }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%,eigenvaluesG-%'vector G6#7,$!+&)RH&4#!\"*$!+%)R@p=F+$!+*4=#z:F+$!+HISwt!#5$!+)=BR+%F2$\"+w** fsZF2$\"+?Cvp))F2$\"+$HucN\"F+$\"+9Bd^;F+$\"+UoF5BF+" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 54 "subvector( eigenvectors, 1..`number of C- atoms`, 1 ); " }}{PARA 12 "" 1 "" {XPPMATH 20 "6#-%'vectorG6#7,$!++jy! f#!#5$\"%\")=!#7$\"+!G'y!f#F)$\"+<@(\\N$F)$!+iR>,;F)$!%!>#F,$\"+/S>,;F )$!+Q@(\\N$F)$\"+]*e%GaF)$!+B*e%GaF)" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 29 "evalm( ( topological_matrix " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "- eigenvalues[1] ) &* % ); " }}{PARA 12 "" 1 "" {XPPMATH 20 "6#-%'vectorG6#7,$\"#=!#5$\"#>F)$\"\"&F)$\"#7F)$\"%5=!#7$! \"%F)F,$\"#8F)$\"#HF)$\"\"#!\"*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 50 "map( fnormal, % ); # floating-point normalization " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'vectorG6#7,\"\"!F'F'F'F'F'F'F'F'F'" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 46 "occupation_numbers := [0,0,0 ,0,0,2,2,2,2,2,2]:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 50 "pi_el ectron_energy := sum( 'occupation_numbers[i] " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 37 " * ( alpha + eigenvalues[i]*beta )'," }}{PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 33 " 'i'=1..`number of C-atoms` ); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%3pi_electron_energyG,&%&alphaG\"#5%%betaG$\"+ P " 0 "" {MPLTEXT 1 0 39 "eigenvectors : = map( normalize, [ seq(" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 55 " subv ector( eigenvectors, 1..`number of C-atoms`, i )," }}{PARA 0 "> " 0 " " {MPLTEXT 1 0 35 " i=1..`number of C-atoms` ) ] ): " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 56 "electron_density := seq( sum( 'occu pation_numbers[i] * " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 57 " (eigenve ctors[i][j])^2', 'i'=1..`number of C-atoms` ), " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 30 " j=1..`number of C-atoms` ); " }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%1electron_densityG6,$\"+S$zG<\"!\"*$\"+3(*fY5F($\"+R$ zG<\"F($\"+cmX\\&)!#5$\"+b%oW')*F/$\"+-Q,+()F/$\"+$[oW')*F/$\"+KmX\\&) F/$\"+1wUF5F($\"+0wUF5F(" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 " " }}}}{SECT 1 {PARA 3 "" 0 "VectorAnalysis" {TEXT -1 20 "19.4 Vector A nalysis" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 50 "restart; \ninterf ace(showassumed=2): \nwith(linalg):" }}{PARA 7 "" 1 "" {TEXT -1 32 "Wa rning, new definition for norm" }}{PARA 7 "" 1 "" {TEXT -1 33 "Warning , new definition for trace" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 55 "alias( f=f(x,y,z), g=g(x,y,z), h=h(x,y,z), c=[x,y,z] ):" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 15 "grad( f, c ); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'vectorG6#7%-%%diffG6$%\"fG%\"xG-F(6$F*%\"yG- F(6$F*%\"zG" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "diverge( [f, g,h], c ); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,(-%%diffG6$%\"fG%\"xG \"\"\"-F%6$%\"gG%\"yGF)-F%6$%\"hG%\"zGF)" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "laplacian( f, c ); " }}{PARA 11 "" 1 "" {XPPMATH 20 " 6#,(-%%diffG6$%\"fG-%\"$G6$%\"xG\"\"#\"\"\"-F%6$F'-F)6$%\"yGF,F--F%6$F '-F)6$%\"zGF,F-" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "jacobian ( [f,g,h], c ); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7%7%-% %diffG6$%\"fG%\"xG-F)6$F+%\"yG-F)6$F+%\"zG7%-F)6$%\"gGF,-F)6$F6F/-F)6$ F6F27%-F)6$%\"hGF,-F)6$F>F/-F)6$F>F2" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "hessian( f, c ); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6# -%'matrixG6#7%7%-%%diffG6$%\"fG-%\"$G6$%\"xG\"\"#-F)6%F+F/%\"yG-F)6%F+ F/%\"zG7%F1-F)6$F+-F-6$F3F0-F)6%F+F3F67%F4F<-F)6$F+-F-6$F6F0" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "curl( [f,g,h], c ); " }} {PARA 11 "" 1 "" {XPPMATH 20 "6#-%'vectorG6#7%,&-%%diffG6$%\"hG%\"yG\" \"\"-F)6$%\"gG%\"zG!\"\",&-F)6$%\"fGF1F--F)6$F+%\"xGF2,&-F)6$F0F9F--F) 6$F6F,F2" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "diverge( %, c ) ; " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#\"\"!" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "curl( grad( f, c ), c ); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'vectorG6#7%\"\"!F'F'" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 40 "plots[coordplot3d]( prolatespheroidal );" }}{PARA 13 "" 1 "" {INLPLOT "6)-%%MESHG6$79797%\"\"!F)$\"1\"Q1_'fiF6!#:7%$\"1y5j7 %ffT(!#[x7&4o9!#;F)$\"1)y![!>\\>3\"F,7%$\"1 q!p583Z;#F6F)$\"1XUucYsD5F,7%$\"1RTM>RD'[*F67%$\"1z26 U&[CT$F6F)$\"1?`NK%G?_)F67%$\"1=!\\Kbv\"QRF6F)$\"15\"3>nzVQ(F67%$\"1'e $)omKPQ%F6F)$\"1]L8OiS'4'F67%$\"1(\\h]i\\+u%F6F)$\"1MZ(*\\vK%o%F67%$\" 11v/UG()**\\F6F)$\"1'y$Qj%*)o<$F67%$\"1h?TmI\"z:&F6F)$\"1M5ro!zZg\"F67 %$\"1uu$\\0`4@&F6F)$\"1%>FYknzVQ(F67%$\"1!y5@a[CT$F6F)$!1<`NK %G?_)F67%$\"1QTM>RD'[*F67%F:F)$!1XUucYsD5F,7%$\"1&>[x7 &4o9F6F)$!1)y![!>\\>3\"F,7%$\"1z5j7%ffT(F0F)$!1eHye$[h6\"F,7%$\"1*HO#z pd\"Q'F`oF)$!1\"Q1_'fiF6F,79F(7%$\"19)H!\\4c:rF0$\"1?Gv1sJ*3#F0F17%$\" 1q?nfpi39F6$\"1R%p-%>5OTF0F77%$\"1VpoTA-x?F6$\"1=_\"ec(o)4'F0F<7%$\"10 ')\\S`8.FF6$\"1KOo*\\@r$zF0FA7%$\"1Gj(*)R?UF$F6$\"1*y'eG&yRh*F0FF7%$\" 1=G?s5OTF0Ffq7%$\"1;)H!\\4c:rF0$\"1@Gv1s J*3#F0F[r7%$\"1es:s#yI7'F`o$\"1]1/**z*yz\"F`oF`r79F(7%$\"1Y&3_0-(QiF0$ \"1*Ri!f.P4SF0F17%Fdt$\"1LOo*\\@r$zF0F77%$\"1R-\"GOo5#=F6$\"1hR`t&H.< \"F6F<7%$\"1\\2`7[-qBF6$\"1-U$>YAJ_\"F6FA7%$\"12`l\"RM2(GF6$\"1M>Lp*3 \\%=F6FF7%$\"1G,3(4/IJ$F6$\"1*>#R]%Q\"H@F6FK7%$\"1sK+$fIyo$F6FjxFP7%$ \"1+\"yG\\$e()RF6$\"1ZL2/VmiDF6FU7%Fbt$\"11')\\S`8.FF6FZ7%$\"1d_:*e7\" RVF6$\"1D**eO$y&)y#F6Fin7%FNF?F^o7%$\"1e_:*e7\"RVF6FfzFdo7%FbtFcsFgo7% $\"1,\"yG\\$e()RF6F^zF\\p7%$\"1tK+$fIyo$F6FjxFap7%FdyFfyFdp7%$\"13`l\" RM2(GF6FayFip7%$\"1[2`7[-qBF6$\"1,U$>YAJ_\"F6F^q7%FexFgxFaq7%$\"19/v7- /N7F6$\"1VOo*\\@r$zF0Ffq7%$\"1Z&3_0-(QiF0$\"1+C1f.P4SF0F[r7%$\"1um\")y T_o`F`o$\"1N,r%*49]MF`oF`r79F(7%$\"1LS'Qi?k&[F0$\"1#RC&p\"3Yg&F0F17%$ \"1!z'eG&yRh*F0F_tF77%$\"1=#*[bBe<9F6$\"1_ml[s(fj\"F6F<7%FayFfyFA7%$\" 1$=PobyYB#F6$\"1Ig?Ll&*yDF6FF7%$\"1Jg?Ll&*yDF6$\"1\"H+\")\\ui(HF6FK7%F _y$\"1F,3(4/IJ$F6FP7%$\"1Y[SwB2/JF6$\"1)Gel[!H#e$F6FU7%FbvF]tFZ7%$\"1y lHOZrxLF6$\"1#=(\\/14)*QF6Fin7%FgpFIF^o7%$\"1zlHOZrxLF6$\"1$=(\\/14)*Q F6Fdo7%FhsF]tFgo7%Fi^lF[_lF\\p7%Fe[l$\"1H,3(4/IJ$F6Fap7%Fa^lFc^lFdp7%$ \"1%=PobyYB#F6$\"1Kg?Ll&*yDF6Fip7%$\"1L>Lp*3\\%=F6$\"1)>#R]%Q\"H@F6F^q 7%Ff]lFh]lFaq7%$\"1/oeG&yRh*F0$\"1fs!zE7&46F6Ffq7%$\"1MS'Qi?k&[F0$\"1$ RC&p\"3Yg&F0F[r7%$\"1=5N%=W!zTF`o$\"1mnw&3uG#[F`oF`r79F(7%$\"1,'>e(3q! 3$F0$\"1H(zz&RzXnF0F17%$\"1?_\"ec(o)4'F0FitF77%$\"1rM8kD_#**)F0$\"14Xi wx3p>F6F<7%FgxF^zFA7%Ff]l$\"1X[SwB2/JF6FF7%Fh]lF[_lFK7%FexF\\zFP7%Febl $\"1ET_)z+F6Faq7%$\"1F_\"ec(o)4'F0$\"1QQU+jUN8F6Ffq7%$\"1-'>e(3q!3$F0$\"1I(zz& RzXnF0F[r7%$\"1')>!Q)G+^EF`o$\"1u68,m)[!eF`oF`r79F(7%$\"1S)\\W1,a0\"F0 $\"1o4uQcZStF0F17%$\"1BGv1sJ*3#F0FauF77%$\"10'>e(3q!3$F0FdclF<7%$\"1/C 1f.P4SF0FfzFA7%$\"1QS'Qi?k&[F0F__lFF7%$\"1*RC&p\"3Yg&F0Fa_lFK7%$\"1`&3 _0-(QiF0FdzFP7%$\"1P(zz&RzXnF0F^dlFU7%$\"1C)H!\\4c:rF0F_uFZ7%$\"1x4uQc ZStF0$\"1*[#\\[HT0^F6Fin7%$\"1(3JETffT(F0FgnF^o7%$\"1z4uQcZStF0$\"1!\\ #\\[HT0^F6Fdo7%$\"1B)H!\\4c:rF0F_uFgo7%$\"1R(zz&RzXnF0F^dlF\\p7%$\"1c& 3_0-(QiF0$\"1f_:*e7\"RVF6Fap7%FeglFa_lFdp7%$\"1SS'Qi?k&[F0Fe_lFip7%$\" 1.C1f.P4SF0$\"1C**eO$y&)y#F6F^q7%$\"1/'>e(3q!3$F0FdclFaq7%$\"1EGv1sJ*3 #F0$\"1\"4mh._JX\"F6Ffq7%Fdfl$\"1p4uQcZStF0F[r7%$\"1Y`Gd4$>3*!#L$\"1T= 1FGv1sJ*3#F0FauF77%$!1)f>e (3q!3$F0FdclF<7%$!1&Ri!f.P4SF0FfzFA7%$!1FS'Qi?k&[F0Fe_lFF7%$!1(QC&p\"3 Yg&F0Fg_lFK7%$!1R&3_0-(QiF0FjzFP7%$!1A(zz&RzXnF0F^dlFU7%$!12)H!\\4c:rF 0FeuFZ7%$!1f4uQcZStF0FchlFin7%$!1r5j7%ffT(F0FboF^o7%$!1i4uQcZStF0F[ilF do7%$!11)H!\\4c:rF0F_uFgo7%$!1B(zz&RzXnF0$\"1ZbYaF!=p%F6F\\p7%$!1T&3_0 -(QiF0$\"1g_:*e7\"RVF6Fap7%Fc\\mFg_lFdp7%$!1HS'Qi?k&[F0$\"1!e'HOZrxLF6 Fip7%$!1$Ri!f.P4SF0F_jlF^q7%$!1(f>e(3q!3$F0FdclFaq7%$!1@Gv1sJ*3#F0Fgjl Ffq7%$!1Q)\\W1,a0\"F0$\"1q4uQcZStF0F[r7%$!1D`Gd4$>3*F_[mF`[mF`r79F(7%$ !1+'>e(3q!3$F0F[flF17%$!1<_\"ec(o)4'F0FitF77%$!1pM8kD_#**)F0FeblF<7%$! 1hR`t&H.<\"F6F^zFA7%$!1<#*[bBe<9F6Fi^lFF7%$!1_ml[s(fj\"F6F[_lFK7%$!1R- \"GOo5#=F6F^[lFP7%$!13Xiwx3p>F6FbdlFU7%$!1VpoTA-x?F6F[vFZ7%$!1mhaBXnU@ F6F^dlFin7%$!1q!p583Z;#F6FjoF^o7%FhamF]^mFdo7%FeamF[vFgo7%$!14Xiwx3p>F 6FbdlF\\p7%F_amFgdlFap7%F\\amF[_lFdp7%$!1=#*[bBe<9F6$\"1Z[SwB2/JF6Fip7 %$!1gR`t&H.<\"F6$\"1YL2/VmiDF6F^q7%$!1oM8kD_#**)F0FeblFaq7%$!1E_\"ec(o )4'F0FfelFfq7%$!1,'>e(3q!3$F0$\"1K(zz&RzXnF0F[r7%$!1&)>!Q)G+^EF`o$\"1v 68,m)[!eF`oF`r79F(7%$!1KS'Qi?k&[F0F`]lF17%$!1*y'eG&yRh*F0F_tF77%FebmFh ]lF<7%$!1M>Lp*3\\%=F6FfyFA7%$!1$=PobyYB#F6F^^lFF7%$!1Ig?Ll&*yDF6Fc^lFK 7%$!11`l\"RM2(GF6Ff^lFP7%$!1X[SwB2/JF6F[_lFU7%$!1Gj(*)R?UF$F6F]tFZ7%$! 1ylHOZrxLF6Fa_lFin7%$!1z26U&[CT$F6FIF^o7%$!1zlHOZrxLF6Fg_lFdo7%$!1Fj(* )R?UF$F6F]tFgo7%FcemF[_lF\\p7%$!12`l\"RM2(GF6F\\`lFap7%F]emFc^lFdp7%Fj dmFb`lFip7%$!1L>Lp*3\\%=F6Fg`lF^q7%Fi`mFh]lFaq7%$!1-oeG&yRh*F0F]alFfq7 %$!1LS'Qi?k&[F0FbalF[r7%$!1<5N%=W!zTF`oFgalF`r79F(7%$!1Y&3_0-(QiF0$\"1 (Ri!f.P4SF0F17%$!17/v7-/N7F6FgvF77%F_am$\"1gR`t&H.<\"F6F<7%$!1\\2`7[-q BF6Fj[lFA7%FffmFe`lFF7%$!1G,3(4/IJ$F6Fg`lFK7%$!1tK+$fIyo$F6$\"1Z2`7[-q BF6FP7%$!1,\"yG\\$e()RF6F\\elFU7%$!1W\"z?pgh?%F6FcsFZ7%$!1e_:*e7\"RVF6 $\"1B**eO$y&)y#F6Fin7%$!1(e$)omKPQ%F6F\\qF^o7%$!1f_:*e7\"RVF6F_jlFdo7% $!1V\"z?pgh?%F6$\"1/')\\S`8.FF6Fgo7%F`imF\\cmF\\p7%$!1uK+$fIyo$F6FjxFa p7%FhhmFg`lFdp7%$!13`l\"RM2(GF6FayFip7%$!1[2`7[-qBF6$\"1+U$>YAJ_\"F6F^ q7%F_amFahmFaq7%$!19/v7-/N7F6$\"1SOo*\\@r$zF0Ffq7%$!1[&3_0-(QiF0$\"1)R i!f.P4SF0F[r7%$!1vm\")yT_o`F`o$\"1M,r%*49]MF`oF`r79F(7%$!18)H!\\4c:rF0 FiflF17%$!1p?nfpi39F6F_wF77%Feam$\"1E_\"ec(o)4'F0F<7%$!10')\\S`8.FF6F` \\lFA7%Fbfm$\"1.oeG&yRh*F0FF7%$!1=G?s5OTF0Ffq7%$!19)H!\\4c:rF0$\"1CGv1sJ*3#F0F[r7% $!1ds:s#yI7'F`o$\"1_1/**z*yz\"F`oF`r79F(7%$!1y5j7%ffT(F0$\"1M`Gd4$>3*F _[mF17%$!1%>[x7&4o9F6FiwF77%F[bmF^flF<7%$!1RTM>RD5N%=W!zTF`oFip7%$!1QTM>RD[x7&4o9F6Fg`nFfq7%$!1z5j7%ffT(F0$\"1P`Gd4$>3*F_[ mF[r7%$!1*HO#zpd\"Q'F`o$\"1pT8@y<:y!#[F`r79F(7%F``nFg[mF17%$!1q?nfpi39 F6$!1N%p-%>5OTF0F77%$!1WpoTA-x?F6$!17_\"ec(o)4'F0F<7%$!11')\\S`8.FF6$! 1DOo*\\@r$zF0FA7%Ffem$!1zneG&yRh*F0FF7%Fc]n$!1ds!zE7&46F6FK7%Fajm$!16/ v7-/N7F6FP7%$!1S0$)4V/[XF6$!1MQU+jUN8F6FU7%$!1J0\"4'GM(z%F6$!1o?nfpi39 F6FZ7%Fe^n$!1(3mh._JX\"F6Fin7%Febn$!1#>[x7&4o9F6F^o7%$!1!yOdM\")*[\\F6 $!1)3mh._JX\"F6Fdo7%F\\^nFegnFgo7%F^gn$!1NQU+jUN8F6F\\p7%FcimF[gnFap7% Fc]nFhfnFdp7%Fb_n$!1#y'eG&yRh*F0Fip7%F]]n$!1AOo*\\@r$zF0F^q7%Feam$!16_ \"ec(o)4'F0Faq7%$!1s?nfpi39F6$!1S%p-%>5OTF0Ffq7%$!1;)H!\\4c:rF0Fg[mF[r 7%$!1es:s#yI7'F`o$!1[1/**z*yz\"F`oF`r79F(7%Figm$!1*Ri!f.P4SF0F17%F^hm$ !1LOo*\\@r$zF0F77%F_amFf`mF<7%Fdhm$!1-U$>YAJ_\"F6FA7%FffmFgdmFF7%Fhhm$ !1*>#R]%Q\"H@F6FK7%$!1sK+$fIyo$F6FdhmFP7%$!1+\"yG\\$e()RF6$!1ZL2/VmiDF 6FU7%FajmF`fnFZ7%$!1d_:*e7\"RVF6$!1D**eO$y&)y#F6Fin7%F_bnFdanF^o7%Ffim Fi[oFdo7%FajmF]]nFgo7%F`imFc[oF\\p7%F[imFdhmFap7%FhhmF[[oFdp7%F[[nFgdm Fip7%F^[n$!1,U$>YAJ_\"F6F^q7%F_amFf`mFaq7%Fd[n$!1VOo*\\@r$zF0Ffq7%$!1Z &3_0-(QiF0$!1+C1f.P4SF0F[r7%$!1um\")yT_o`F`o$!1N,r%*49]MF`oF`r79F(7%$! 1MS'Qi?k&[F0$!1#RC&p\"3Yg&F0F17%$!1$z'eG&yRh*F0FhfnF77%FebmF\\amF<7%Fg dmF[[oFA7%FjdmF]emFF7%$!1Jg?Ll&*yDF6$!1!H+\")\\ui(HF6FK7%Fffm$!1F,3(4/ IJ$F6FP7%$!1Y[SwB2/JF6$!1(Gel[!H#e$F6FU7%Fb_nFc]nFZ7%F_fm$!1\"=(\\/14) *QF6Fin7%F]dn$!1#R]%Q\"H@F6F^q7%Febm$!1^m l[s(fj\"F6Faq7%$!11oeG&yRh*F0$!1fs!zE7&46F6Ffq7%$!1NS'Qi?k&[F0Fg]oF[r7 %$!1>5N%=W!zTF`o$!1lnw&3uG#[F`oF`r79F(7%$!12'>e(3q!3$F0$!1E(zz&RzXnF0F 17%$!1I_\"ec(o)4'F0FdhnF77%$!1([LTcAD**)F0FbamF<7%$!1jR`t&H.<\"F6$!1XL 2/VmiDF6FA7%$!1?#*[bBe<9F6$!1W[SwB2/JF6FF7%$!1bml[s(fj\"F6Fj^oFK7%$!1U -\"GOo5#=F6$!1)4yG\\$e()RF6FP7%$!17Xiwx3p>F6$!1DT_)z+!Q)G+^EF`o$!1r68,m)[!eF`oF`r79F(7%$!1R)\\W1,a0\" F0$!1o4uQcZStF0F17%$!1AGv1sJ*3#F0$!1*3mh._JX\"F6F77%$!1-'>e(3q!3$F0Fha mF<7%$!1,C1f.P4SF0Fi[oFA7%FaaoFiemFF7%$!1&RC&p\"3Yg&F0Ff_oFK7%Fi[nFg[o FP7%$!1L(zz&RzXnF0$!1YbYaF!=p%F6FU7%$!1>)H!\\4c:rF0Fe^nFZ7%$!1s4uQcZSt F0$!1*[#\\[HT0^F6Fin7%$!1#3JETffT(F0FhbnF^o7%$!1t4uQcZStF0$!1!\\#\\[HT 0^F6Fdo7%$!1<)H!\\4c:rF0Fe^nFgo7%F^hoF`hoF\\p7%$!1^&3_0-(QiF0F^jmFap7% FjgoFf_oFdp7%$!1OS'Qi?k&[F0F_fmFip7%F\\]o$!1C**eO$y&)y#F6F^q7%FcgoFham Faq7%$!1CGv1sJ*3#F0$!1\"4mh._JX\"F6Ffq7%Fifo$!1p4uQcZStF0F[r7%$!1S`Gd4 $>3*F_[m$!1T=1Fe(3q!3$F0FhamF<7%$\"1\"Ri!f.P4SF0Fi[oFA7%$\"1AS'Qi?k&[ F0F_fmFF7%$\"1\"QC&p\"3Yg&F0$!1$=(\\/14)*QF6FK7%$\"1K&3_0-(QiF0FfimFP7 %$\"1:(zz&RzXnF0F`hoFU7%$\"1,)H!\\4c:rF0F^hnFZ7%$\"1_4uQcZStF0FhhoFin7 %$\"1i5j7%ffT(F0F`cnF^o7%$\"1a4uQcZStF0F`ioFdo7%$\"1*zH!\\4c:rF0Fe^nFg o7%$\"1;(zz&RzXnF0$!1ZbYaF!=p%F6F\\p7%$\"1N&3_0-(QiF0$!1g_:*e7\"RVF6Fa p7%F_\\pFa\\pFdp7%$\"1DS'Qi?k&[F0Fd_oFip7%$\"1*Qi!f.P4SF0F^joF^q7%$\"1 %f>e(3q!3$F0FhamFaq7%$\"1>Gv1sJ*3#F0FdjoFfq7%F`[p$!1q4uQcZStF0F[r7%$\" 1<`Gd4$>3*F_[mF\\[pF`r79F(7%FbjlF`foF17%$\"1D_\"ec(o)4'F0$!1OQU+jUN8F6 F77%$\"1zM8kD_#**)F0FbamF<7%$\"1iR`t&H.<\"F6Fc[oFA7%$\"1>#*[bBe<9F6Fce mFF7%$\"1`ml[s(fj\"F6F\\`oFK7%$\"1T-\"GOo5#=F6Fa[oFP7%$\"15Xiwx3p>F6$! 1ET_)z+F6F^apF\\p7%Fi`pF`imFap7%Ff`pF\\`oFdp7%Fc`pFh^oFip7%F``pFgboF^q7%F] `pFbamFaq7%$\"1L_\"ec(o)4'F0F]foFfq7%F\\gl$!1H(zz&RzXnF0F[r7%$\"1*)>!Q )G+^EF`o$!1u68,m)[!eF`oF`r79F(7%F`alFg]oF17%$\"1&z'eG&yRh*F0FhfnF77%Ff ]lF\\amF<7%$\"1N>Lp*3\\%=F6F[[oFA7%F``lF]emFF7%Fb`lFb^oFK7%Fe[lFe^oFP7 %FgbmFj^oFU7%$\"1Ij(*)R?UF$F6Fc]nFZ7%Fe_lF^_oFin7%$\"1\"y5@a[CT$F6Fa_o F^o7%Fh^mFf_oFdo7%FbvFi_oFgo7%FgbmF\\`oF\\p7%$\"14`l\"RM2(GF6FhhmFap7% Fb`lFb^oFdp7%F``lF`^oFip7%FayFf`oF^q7%Ff]lFi`oFaq7%$\"12oeG&yRh*F0F^ao Ffq7%$\"1OS'Qi?k&[F0Fg]oF[r7%F_dnFfaoF`r79F(7%$\"1U&3_0-(QiF0$!1/C1f.P 4SF0F17%$\"16/v7-/N7F6Fg\\oF77%$\"1Q-\"GOo5#=F6$!1iR`t&H.<\"F6F<7%F]im $!1.U$>YAJ_\"F6FA7%$\"10`l\"RM2(GF6$!1O>Lp*3\\%=F6FF7%$\"1E,3(4/IJ$F6$ !1-AR]%Q\"H@F6FK7%$\"1qK+$fIyo$F6$!1_2`7[-qBF6FP7%$\"1)4yG\\$e()RF6$!1 ]L2/VmiDF6FU7%$\"1T\"z?pgh?%F6$!14')\\S`8.FF6FZ7%$\"1b_:*e7\"RVF6$!1G* *eO$y&)y#F6Fin7%$\"1%e$)omKPQ%F6$!1VTM>RDLp*3\\%=F6Fip7%F]imFffpF^q7%FafpFcfpFaq7%$\"18/v7- /N7F6$!1`Oo*\\@r$zF0Ffq7%$\"1V&3_0-(QiF0F[fpF[r7%$\"1rm\")yT_o`F`o$!1S ,r%*49]MF`oF`r79F(7%FdrFa_mF17%FirFdinF77%F^s$!1@_\"ec(o)4'F0F<7%Fcs$! 1OOo*\\@r$zF0FA7%FhsFj]oFF7%F]t$!1es!zE7&46F6FK7%FbtF^hmFP7%FgtFj_pFU7 %F\\uFfenFZ7%F_uF`goFin7%FXF`anF^o7%Feu$!1!4mh._JX\"F6Fdo7%F\\uFfenFgo 7%F[vFj_pF\\p7%F^v$!18/v7-/N7F6Fap7%F]tFc[qFdp7%Fbv$!1'z'eG&yRh*F0Fip7 %Fcs$!1KOo*\\@r$zF0F^q7%F^s$!1?_\"ec(o)4'F0Faq7%F]w$!1Y%p-%>5OTF0Ffq7% FbwF^goF[r7%Fgw$!1]1/**z*yz\"F`oF`r79F(7%F.$!1nqX\">'Q;=F`oF17%F4$!1*H \"3)*fz&f$F`oF77%F:$!1sRgnd+-`F`oF<7%F?$!1s-U*)>G+pF`oFA7%FD$!1N?qo$)3 e$)F`oFF7%FI$!1KN`r\"[dk*F`oFK7%FN$!1NLwN[qt5!#JFP7%FS$!1NiA?t(4;\"Fi^ qFU7%FX$!1_9VachC7Fi^qFZ7%Fgn$!1oBTXVKj7Fi^qFin7%F\\o$!1gs%eR:jF\"Fi^q F^o7%FboFa_qFdo7%FXF^_qFgo7%FjoF[_qF\\p7%F_pFg^qFap7%FIFd^qFdp7%Fgp$!1 R?qo$)3e$)F`oFip7%F\\q$!1p-U*)>G+pF`oF^q7%F:$!1rRgnd+-`F`oFaq7%Fdq$!1/ 83)*fz&f$F`oFfq7%FiqFe]qF[r7%F^r$!1MoAkb.j:!#ZF`r-%'COLOURG6&%$RGBG$\" #5!\"\"F)F)-F$6$79797%F)F)$\"1wa'=\"y1rqF67%$\"1V*R#Q'G_@$F0F)$\"1.\"* y\"*QPyqF67%$\"1F$3ZX,rV'F0F)$\"1[cOHsI+rF67%$\"1%HQwNwAn*F0F)$\"158@ \\J\"p8(F67%$\"12XOuQu#H\"F6F)$\"1I[T'Hn#)=(F67%$\"1F2HJ:$4i\"F6F)$\"1 >f&GzvWD(F67%$\"1'3qpzoC&>F6F)$\"1!4.lXvcL(F67%$\"1&3,PyS!)G#F6F)$\"1! 4WT3M?V(F67%$\"1N\"=9%4MGEF6F)$\"1mJ8)z]Pa(F67%$\"1ogo!\\sS(HF6F)$\"1L <(yXc5n(F67%$\"1b`)e()\\fK$F6F)$\"1UF(o89U\"yF67%$\"1V5fT-q%o$F6F)$\"1 xXRm'>N(zF67%$\"1di>N\\1^SF6F)$\"1?zwYAI\\\")F67%$\"1Ua))Q5!eU%F6F)$\" 1J)yw7D>M)F67%$\"1$3BT$Ho4[F6F)$\"1iQ\")ejy^&)F67%$\"174*=!R]._F6F)$\" 1x$=Ch>$z()F67%$\"17/)>wx!3cF6F)$\"1qyYz]*\\-*F67%$\"1'>GZbSU-'F6F)$\" 1#>/8W?$*G*F67%$\"12&p'oA&GX'F6F)$\"1rXu=>%Gd*F67%$\"1v*zwh)z%*oF6F)$ \"1B8$3SXh()*F67%$\"1v!y7(G*4N(F6F)$\"1'[bbw&)*>5F,7%$\"1E0(=uxB#yF6F) $\"1%=0DwkW0\"F,7%$\"1p0%GLF*4$)F6F)$\"1-bz!yA64\"F,79Fgaq7%$\"1M`&pG* )\\3$F0$\"1t7\"3&eMe!*!#=F]bq7%$\"1+jB;ONwhF0$\"1jE[^5a8=F0Fbbq7%$\"1m :,/7[!G*F0$\"1ZVr[^*\\s#F0Fgbq7%$\"1Ff8k'y.C\"F6$\"1'3e1P!3UOF0F\\cq7% $\"1TO)\\Ts_b\"F6$\"1/VvZ=pmXF0Facq7%$\"1\\V*=4!Qt=F6$\"1F\\Bj-u+bF0Ff cq7%$\"1*f/`0f`>#F6$\"1_d\\,e:YkF0F[dq7%$\"19SblY(=_#F6$\"1;uwO@*[S(F0 F`dq7%$\"1\\M$yl,O&GF6$\"1\\uKq/$*y$)F0Fedq7%$\"13g]ObA\">$F6$\"1ohfCO Gq$*F0Fjdq7%$\"1ZAo\"*RWNNF6$\"1ywD?+5Q5F6F_eq7%$\"1yo(4Mop)QF6$\"1-N7 #)oJT6F6Fdeq7%$\"1DjU**\\_YUF6$\"1^mvUA*oC\"F6Fieq7%$\"1dYg!)p&[h%F6$ \"1'3F1FW]N\"F6F^fq7%$\"1^*R?VDF*\\F6$\"1$39)ek*fY\"F6Fcfq7%$\"1@R537 \"4Q&F6$\"1/$[o3y*z:F6Fhfq7%$\"12_i$[;-y&F6$\"139^%pCsp\"F6F]gq7%$\"1h 2c6kY\">'F6$\"1!Gs&o&yz\"=F6Fbgq7%$\"1&)3:G3^:mF6$\"1MdYV#*[U>F6Fggq7% $\"1s:34gA`qF6$\"1&H[p,952#F6F\\hq7%$\"14KG![;b](F6$\"1smz\"[=Q?#F6Fah q7%$\"1(e2r)oJtzF6$\"1cwkuqt7_U;F6$\"1b/XEVeb5F6Ffcq7%$\"1i.b#RA[#> F6$\"1,>nR#3qB\"F6F[dq7%$\"1RGX-656AF6$\"1;.7&R))4U\"F6F`dq7%$\"1%o&o* )*[>]#F6$\"1@n`#)\\!zg\"F6Fedq7%$\"1$)z5$3nzz#F6$\"1-]rfU9)z\"F6Fjdq7% $\"1sBm&4n(*4$F6$\"1?Ep]$*4#*>F6F_eq7%$\"1sBY)oszS$F6$\"17-8e5#F6Fd eq7%$\"1!Hq&f2ABPF6$\"1e)e:poFR#F6Fieq7%$\"1&))G,wih/%F6$\"1Tl?84J+EF6 F^fq7%$\"1mlcSgYxVF6$\"1%p^JgEK\"GF6Fcfq7%$\"11!e'H_\"yr%F6$\"1[;fXd&> .$F6Fhfq7%$\"1q**G\\O\"z1&F6$\"19u0Q.&pD$F6F]gq7%$\"1Cf&yw%[GaF6$\"1!p `aKv')[$F6Fbgq7%$\"1CGG&pt-!eF6$\"1DU&3c4ws$F6Fggq7%$\"14&zFs[S='F6$\" 1[!Q]zYU(RF6F\\hq7%$\"1y&)[5Hg!e'F6$\"1_8\"op'4HUF6Fahq7%$\"1CDKEdv!*p F6$\"1RNA1fo#\\%F6Ffhq79Fgaq7%$\"14(Re)p_0@F0$\"1f9![n2*HCF0F]bq7%$\"1 d06$)\\S:UF0$\"1&*4Z%oO['[F0Fbbq7%$\"1UY[%*R*RL'F0$\"1(*=KS(=)4tF0Fgbq 7%$\"1EhpB?nl%)F0$\"160+Gj!*p(*F0F\\cq7%$\"1+hD7W[h5F6$\"1j.$[J=]A\"F6 Facq7%$\"1)4*Q\\^fy7F6$\"18,Diudv9F6Ffcq7%$\"1Q7\"p1[$)\\\"F6$\"1zsQ[e =H@F6F`dq7%$\"1gr0%H.w%>F6$\"1/Z$* =SlZAF6Fedq7%$\"18vywR.y@F6$\"1bR-._e8DF6Fjdq7%$\"1!owX]lHT#F6$\"1ka?R 1r%y#F6F_eq7%$\"1IEscL)Gl#F6$\"1KgX+1fhIF6Fdeq7%$\"1E.i;LG)*GF6$\"12VJ ^szWLF6Fieq7%$\"1bRG'\\s'\\JF6$\"1gCCIe\"\\j$F6F^fq7%$\"1P%zUQqvS$F6$ \"1zZ+feaKRF6Fcfq7%$\"1DFB>*4Dn$F6$\"1yAc\"Q-$QUF6Fhfq7%$\"1eA8\"fQ]%R F6$\"1&\\V]BdA%F6$\"1!*yTb.uw[F6Fbgq7%$\"1Y)yJ)G 8:XF6$\"1p*p07T2@&F6Fggq7%$\"1><\"pewQ\"[F6$\"1#f\\Nt4bb&F6F\\hq7%$\"1 %fQ#H!oD7&F6$\"1uo&)['e<\"fF6Fahq7%$\"1hf@7^%=W&F6$\"1w!4Z/C-G'F6Ffhq7 9Fgaq7%$\"1%[elCacL\"F0$\"1!>B<%[nCHF0F]bq7%$\"1A@]%eoSn#F0$\"14eaSMRb eF0Fbbq7%$\"1^h')3)3!=SF0$\"1b)=G/7#)z)F0Fgbq7%$\"13CxL@Dq`F0$\"1=0#** =@f<\"F6F\\cq7%$\"1]IYJHfLnF0$\"1R59$=^WZ\"F6Facq7%$\"1g5z1&[36)F0$\"1 O`l(3Ggx\"F6Ffcq7%$\"1F.5>\\'[]*F0$\"1M0`.^F\"3#F6F[dq7%$\"10EN'G_=4\" F6$\"1_T,:I#3R#F6F`dq7%$\"11pQiVZN7F6$\"1]<.%\\6`q#F6Fedq7%$\"1PX$4^\\ ;Q\"F6$\"1320A/RDIF6Fjdq7%$\"1(fo()zz1`\"F6$\"1PtGo\"F6$\"1%H8.Gy\\o$F6Fdeq7%$\"1fB6'>W&Q=F6$\"1lM/I-&e-%F6Fieq7%$\" 1#p8')\\9!)*>F6$\"1a*)Q%[T]P%F6F^fq7%$\"1!Qf5k8;;#F6$\"1`?6POFLZF6Fcfq 7%$\"1$o3kpz'HBF6$\"1VNtkpG,^F6Fhfq7%$\"10*o$o*fD]#F6$\"1?IHc>%)zaF6F] gq7%$\"1'f)34*z7P0$F6$\"1:vL8$)p'o'F6F\\hq7%$\"17[M^J`\\KF6$\"1eWvt] [:rF6Fahq7%$\"1(p;5do?X$F6$\"1=G-\"yv*evF6Ffhq79Fgaq7%$\"1IviOuuvXF^iq $\"1G&zmC-D=$F0F]bq7%$\"1aavC0&4;*F^iq$\"1YKVV5erjF0Fbbq7%$\"1:=dX%3lP \"F0$\"1#R)eIl#Qd*F0Fgbq7%$\"1_'*3ajwR=F0$\"1?iq0cez7F6F\\cq7%$\"1'3xs 3EoI#F0$\"1cC\\&yKWg\"F6Facq7%$\"1QVu5GlyFF0$\"1hW%)zafK>F6Ffcq7%$\"14 5&4a@iD$F0$\"1(Q0w#=vkAF6F[dq7%$\"1k(\\+;>0u$F0$\"1(=*Q?#)e,EF6F`dq7%$ \"1Y,DbkaKUF0$\"1+$*)4s+Q%HF6Fedq7%$\"1yn\"p=?Lt%F0$\"1^Dx,l4#H$F6Fjdq 7%$\"1Zcn*=vQC&F0$\"1o1.1`>ZOF6F_eq7%$\"1D.(3^m_w&F0$\"1Z#4`$4$)4SF6Fd eq7%$\"1+sy!fr&)H'F0$\"1=xylFv!Q%F6Fieq7%$\"1q]:*[#*[%oF0$\"1DQ%oHF2w% F6F^fq7%$\"1W,_i\"e`S(F0$\"1BF+N(R0:&F6Fcfq7%$\"1/D?(zE6)zF0$\"1YQa;c* 4b&F6Fhfq7%$\"1*4bE?)Qt&)F0$\"1a'=@ZAH'fF6F]gq7%$\"1L5 o`rrQ'F6Fbgq7%$\"1p*>?]@B\")*F0$\"1a1B4&>Y#oF6Fggq7%$\"18Sih`:Y5F6$\"1 @\\5j.0u$F0F\\gsF`dq7%$!1O,DbkaKUF0$\"1,$*)4s +Q%HF6Fedq7%$!1nn\"p=?Lt%F0$\"1_Dx,l4#H$F6Fjdq7%$!1Ocn*=vQC&F0F[hsF_eq 7%$!16.(3^m_w&F0F`hsFdeq7%$!1%=(y!fr&)H'F0$\"1>xylFv!Q%F6Fieq7%$!1b]:* [#*[%oF0FjhsF^fq7%$!1F,_i\"e`S(F0F_isFcfq7%$!1&[-szE6)zF0$\"1ZQa;c*4b& F6Fhfq7%$!1!3bE?)Qt&)F0FiisF]gq7%$!155o`rrQ'F6Fbgq7%$ !1Z*>?]@B\")*F0$\"1b1B4&>Y#oF6Fggq7%$!15Sih`:Y5F6$\"1A\\5j.B<%[nCHF0F]bq7%$!1@@]%eoSn#F0F_^sFbbq7%$!1]h')3)3!=S F0Fd^sFgbq7%$!12CxL@Dq`F0Fi^sF\\cq7%$!1\\IYJHfLnF0F^_sFacq7%$!1d5z1&[3 6)F0Fc_sFfcq7%$!1C.5>\\'[]*F0Fh_sF[dq7%$!10EN'G_=4\"F6F]`sF`dq7%$!10pQ iVZN7F6Fb`sFedq7%$!1PX$4^\\;Q\"F6Fg`sFjdq7%$!1(fo()zz1`\"F6F\\asF_eq7% $!1Rgs#>tGo\"F6$\"1&H8.Gy\\o$F6Fdeq7%$!1eB6'>W&Q=F6$\"1mM/I-&e-%F6Fieq 7%$!1\"p8')\\9!)*>F6$\"1b*)Q%[T]P%F6F^fq7%$!1z$f5k8;;#F6F`bsFcfq7%$!1# o3kpz'HBF6$\"1XNtkpG,^F6Fhfq7%$!10*o$o*fD]#F6FjbsF]gq7%$!1&f)34*z7P0$F 6$\"1;vL8$)p'o'F6F\\hq7%$!16[M^J`\\KF6$\"1fWvt][:rF6Fahq7%$!1&p;5do?X$ F6FcdsFfhq79Fgaq7%$!13(Re)p_0@F0F[grF]bq7%$!1d06$)\\S:UF0F`grFbbq7%$!1 TY[%*R*RL'F0FegrFgbq7%$!1DhpB?nl%)F0FjgrF\\cq7%$!1+hD7W[h5F6F_hrFacq7% $!1)4*Q\\^fy7F6FdhrFfcq7%$!1Q7\"p1[$)\\\"F6FihrF[dq7%$!1@b&zF(>@F6FcirFedq7%$!18vywR.y@F6FhirFjdq7%$!1!owX]lHT#F6F] jrF_eq7%$!1IEscL)Gl#F6FbjrFdeq7%$!1D.i;LG)*GF6FgjrFieq7%$!1aRG'\\s'\\J F6F\\[sF^fq7%$!1P%zUQqvS$F6Fa[sFcfq7%$!1CFB>*4Dn$F6Ff[sFhfq7%$!1dA8\"f Q]%RF6F[\\sF]gq7%$!169qs&>dA%F6F`\\sFbgq7%$!1X)yJ)G8:XF6Fe\\sFggq7%$!1 =<\"pewQ\"[F6Fj\\sF\\hq7%$!1%fQ#H!oD7&F6F_]sFahq7%$!1ff@7^%=W&F6Fd]sFf hq79Fgaq7%$!1T$z2XA[q#F0$\"1!p%GRQGQt7_U;F6$\" 1a/XEVeb5F6Ffcq7%$!1i.b#RA[#>F6$\"1+>nR#3qB\"F6F[dq7%$!1RGX-656AF6$\"1 :.7&R))4U\"F6F`dq7%$!1&o&o*)*[>]#F6$\"1?n`#)\\!zg\"F6Fedq7%$!1%)z5$3nz z#F6$\"1,]rfU9)z\"F6Fjdq7%$!1sBm&4n(*4$F6$\"1>Ep]$*4#*>F6F_eq7%$!1sBY) oszS$F6$\"16-8e5#F6Fdeq7%$!1\"Hq&f2ABPF6FhcrFieq7%$!1&))G,wih/%F6$ \"1Rl?84J+EF6F^fq7%$!1nlcSgYxVF6$\"1$p^JgEK\"GF6Fcfq7%$!12!e'H_\"yr%F6 $\"1Z;fXd&>.$F6Fhfq7%$!1q**G\\O\"z1&F6$\"18u0Q.&pD$F6F]gq7%$!1Ef&yw%[G aF6$\"1)o`aKv')[$F6Fbgq7%$!1CGG&pt-!eF6$\"1CU&3c4ws$F6Fggq7%$!14&zFs[S ='F6$\"1Z!Q]zYU(RF6F\\hq7%$!1z&)[5Hg!e'F6$\"1^8\"op'4HUF6Fahq7%$!1CDKE dv!*pF6$\"1PNA1fo#\\%F6Ffhq79Fgaq7%$!1M`&pG*)\\3$F0$\"1(G63&eMe!*F^iqF ]bq7%$!1*HOih`j<'F0$\"1mE[^5a8=F0Fbbq7%$!1m:,/7[!G*F0$\"1^Vr[^*\\s#F0F gbq7%$!1Ff8k'y.C\"F6$\"1\"4e1P!3UOF0F\\cq7%$!1SO)\\Ts_b\"F6$\"15VvZ=pm XF0Facq7%$!1[V*=4!Qt=F6$\"1M\\Bj-u+bF0Ffcq7%$!1*f/`0f`>#F6$\"1hd\\,e:Y kF0F[dq7%$!19SblY(=_#F6$\"1EuwO@*[S(F0F`dq7%$!1\\M$yl,O&GF6$\"1guKq/$* y$)F0Fedq7%$!12g]ObA\">$F6$\"1!='fCOGq$*F0Fjdq7%$!1ZAo\"*RWNNF6$\"1!od --+\"Q5F6F_eq7%$!1yo(4Mop)QF6$\"1/N7#)oJT6F6Fdeq7%$!1CjU**\\_YUF6$\"1` mvUA*oC\"F6Fieq7%$!1dYg!)p&[h%F6$\"1)3F1FW]N\"F6F^fq7%$!1^*R?VDF*\\F6$ \"1%39)ek*fY\"F6Fcfq7%$!1?R537\"4Q&F6$\"11$[o3y*z:F6Fhfq7%$!11_i$[;-y& F6$\"159^%pCsp\"F6F]gq7%$!1g2c6kY\">'F6$\"1$Gs&o&yz\"=F6Fbgq7%$!1&)3:G 3^:mF6$\"1PdYV#*[U>F6Fggq7%$!1r:34gA`qF6$\"1(H[p,952#F6F\\hq7%$!13KG![ ;b](F6$\"1vmz\"[=Q?#F6Fahq7%$!1&e2r)oJtzF6$\"1fwkuqv$RF_[mF]bq7%$!1F$3ZX,rV'F0$\"1C`H#pvJ)yF_[mFbb q7%$!1%HQwNwAn*F0$\"1O]=GA^%=\"F`oFgbq7%$!12XOuQu#H\"F6$\"1Uy&ykaJe\"F `oF\\cq7%$!1F2HJ:$4i\"F6$\"12AV6'o])>F`oFacq7%$!1'3qpzoC&>F6$\"1+XkiY3 \"R#F`oFfcq7%$!1&3,PyS!)G#F6$\"1N\"y?#=/-GF`oF[dq7%$!1N\"=9%4MGEF6$\"1 gc(4K*y=KF`oF`dq7%$!1ogo!\\sS(HF6$\"1*)3:c$)=UOF`oFedq7%$!1b`)e()\\fK$ F6$\"1PvipQ6tSF`oFjdq7%$!1V5fT-q%o$F6$\"1^>ocjX7XF`oF_eq7%$!1di>N\\1^S F6$\"1wh-1P7h\\F`oFdeq7%$!1Ua))Q5!eU%F6$\"1(3Mf2V+U&F`oFieq7%$!1$3BT$H o4[F6$\"1$[(=5G;!*eF`oF^fq7%$!174*=!R]._F6$\"1o9!yRaCP'F`oFcfq7%$!17/) >wx!3cF6$\"1k(y0[9z'oF`oFhfq7%$!1'>GZbSU-'F6$\"1h:k7pcxtF`oF]gq7%$!12& p'oA&GX'F6$\"1[o>w[Y-zF`oFbgq7%$!1v*zwh)z%*oF6$\"1xt!z0$pV%)F`oFggq7%$ !1v!y7(G*4N(F6$\"11dO!*)pB+*F`oF\\hq7%$!1E0(=uxB#yF6$\"1TKIj)\\'z&*F`o Fahq7%$!1p0%GLF*4$)F6$\"1I5$4fsw,\"Fi^qFfhq79Fgaq7%$!1N`&pG*)\\3$F0$!1 k7\"3&eMe!*F^iqF]bq7%$!1,jB;ONwhF0$!1iE[^5a8=F0Fbbq7%$!1o:,/7[!G*F0$!1 XVr[^*\\s#F0Fgbq7%Fcbu$!1#3e1P!3UOF0F\\cq7%$!1TO)\\Ts_b\"F6$!1+VvZ=pmX F0Facq7%$!1\\V*=4!Qt=F6$!1@\\Bj-u+bF0Ffcq7%Fbcu$!1Xd\\,e:YkF0F[dq7%Fgc u$!14uwO@*[S(F0F`dq7%$!1]M$yl,O&GF6$!1TuKq/$*y$)F0Fedq7%$!13g]ObA\">$F 6$!1ehfCOGq$*F0Fjdq7%$!1[Ao\"*RWNNF6$!1xwD?+5Q5F6F_eq7%F[eu$!1,N7#)oJT 6F6Fdeq7%$!1DjU**\\_YUF6$!1]mvUA*oC\"F6Fieq7%$!1eYg!)p&[h%F6$!1%3F1FW] N\"F6F^fq7%$!1_*R?VDF*\\F6$!1\"39)ek*fY\"F6Fcfq7%$!1AR537\"4Q&F6$!1.$[ o3y*z:F6Fhfq7%$!12_i$[;-y&F6$!119^%pCsp\"F6F]gq7%$!1h2c6kY\">'F6$!1zAd o&yz\"=F6Fbgq7%$!1')3:G3^:mF6$!1KdYV#*[U>F6Fggq7%$!1s:34gA`qF6$!1#H[p, 952#F6F\\hq7%$!15KG![;b](F6$!1qmz\"[=Q?#F6Fahq7%$!1)e2r)oJtzF6$!1`wkuq t7_U;F6$!1b/XEVeb5F6Ffcq7%Fe \\u$!1,>nR#3qB\"F6F[dq7%Fj\\u$!1;.7&R))4U\"F6F`dq7%$!1%o&o*)*[>]#F6$!1 @n`#)\\!zg\"F6Fedq7%$!1$)z5$3nzz#F6$!1-]rfU9)z\"F6Fjdq7%Fi]u$!1?Ep]$*4 #*>F6F_eq7%F^^u$!17-8e5#F6Fdeq7%$!1!Hq&f2ABPF6$!1e)e:poFR#F6Fieq7%F f^u$!1Tl?84J+EF6F^fq7%$!1mlcSgYxVF6$!1%p^JgEK\"GF6Fcfq7%$!11!e'H_\"yr% F6$!1[;fXd&>.$F6Fhfq7%Fe_u$!19u0Q.&pD$F6F]gq7%$!1Cf&yw%[GaF6$!1!p`aKv' )[$F6Fbgq7%F_`u$!1DU&3c4ws$F6Fggq7%Fd`u$!1[!Q]zYU(RF6F\\hq7%$!1y&)[5Hg !e'F6$!1_8\"op'4HUF6Fahq7%F^au$!1RNA1fo#\\%F6Ffhq79Fgaq7%$!14(Re)p_0@F 0$!1e9![n2*HCF0F]bq7%$!1e06$)\\S:UF0$!1%*4Z%oO['[F0Fbbq7%$!1WY[%*R*RL' F0$!1&*=KS(=)4tF0Fgbq7%$!1FhpB?nl%)F0$!140+Gj!*p(*F0F\\cq7%$!1,hD7W[h5 F6$!1j.$[J=]A\"F6Facq7%$!1*4*Q\\^fy7F6$!18,Diudv9F6Ffcq7%Ffgt$!1ysQ[e= H@F6F`dq7%$!1hr0%H.w%>F6$!1/Z$*=SlZA F6Fedq7%$!19vywR.y@F6$!1bR-._e8DF6Fjdq7%Fbht$!1ja?R1r%y#F6F_eq7%$!1JEs cL)Gl#F6$!1KgX+1fhIF6Fdeq7%$!1E.i;LG)*GF6$!11VJ^szWLF6Fieq7%$!1bRG'\\s '\\JF6$!1fCCIe\"\\j$F6F^fq7%$!1Q%zUQqvS$F6$!1yZ+feaKRF6Fcfq7%$!1DFB>*4 Dn$F6$!1xAc\"Q-$QUF6Fhfq7%$!1fA8\"fQ]%RF6$!1%\\V]BdA%F6$!1*)yTb.uw[F6Fbgq7%$!1Y)yJ)G8:XF6$!1o*p07T2@&F6Fggq7%$!1><\"pew Q\"[F6$!1\"f\\Nt4bb&F6F\\hq7%$!1'fQ#H!oD7&F6$!1to&)['e<\"fF6Fahq7%$!1i f@7^%=W&F6$!1v!4Z/C-G'F6Ffhq79Fgaq7%$!1'[elCacL\"F0$!1*=B<%[nCHF0F]bq7 %$!1E@]%eoSn#F0$!11eaSMRbeF0Fbbq7%$!1eh')3)3!=SF0$!1^)=G/7#)z)F0Fgbq7% $!1=CxL@Dq`F0$!1<0#**=@f<\"F6F\\cq7%$!1jIYJHfLnF0$!1Q59$=^WZ\"F6Facq7% $!1u5z1&[36)F0$!1N`l(3Ggx\"F6Ffcq7%$!1V.5>\\'[]*F0$!1L0`.^F\"3#F6F[dq7 %$!12EN'G_=4\"F6$!1^T,:I#3R#F6F`dq7%$!13pQiVZN7F6$!1\\<.%\\6`q#F6Fedq7 %$!1SX$4^\\;Q\"F6$!1220A/RDIF6Fjdq7%$!1+'o()zz1`\"F6$!1OtGo\"F6$!1$H8.Gy\\o$F6Fdeq7%$!1iB6'>W&Q=F6$!1jM/I-&e-%F6Fieq7 %$!1&p8')\\9!)*>F6$!1`*)Q%[T]P%F6F^fq7%$!1%Qf5k8;;#F6$!1^?6POFLZF6Fcfq 7%$!1(o3kpz'HBF6$!1UNtkpG,^F6Fhfq7%$!15*o$o*fD]#F6$!1=IHc>%)zaF6F]gq7% $!1+')34*z7P0$F6$!17vL8$)p'o'F6F\\hq7%$!1<[M^J`\\KF6$!1cWvt][:rF6Fahq7%$ !1-n,r&o?X$F6$!19G-\"yv*evF6Ffhq79Fgaq7%$!1EviOuuvXF^iq$!1G&zmC-D=$F0F ]bq7%$!1ZavC0&4;*F^iq$!1YKVV5erjF0Fbbq7%$!19=dX%3lP\"F0$!1#R)eIl#Qd*F0 Fgbq7%$!1]'*3ajwR=F0$!1?iq0cez7F6F\\cq7%$!1%3xs3EoI#F0$!1cC\\&yKWg\"F6 Facq7%$!1OVu5GlyFF0$!1hW%)zafK>F6Ffcq7%$!125&4a@iD$F0$!1(Q0w#=vkAF6F[d q7%$!1h(\\+;>0u$F0$!1(=*Q?#)e,EF6F`dq7%$!1V,DbkaKUF0$!1+$*)4s+Q%HF6Fed q7%$!1un\"p=?Lt%F0$!1^Dx,l4#H$F6Fjdq7%$!1Wcn*=vQC&F0$!1o1.1`>ZOF6F_eq7 %$!1?.(3^m_w&F0$!1Z#4`$4$)4SF6Fdeq7%$!1%>(y!fr&)H'F0$!1=xylFv!Q%F6Fieq 7%$!1l]:*[#*[%oF0$!1DQ%oHF2w%F6F^fq7%$!1Q,_i\"e`S(F0$!1BF+N(R0:&F6Fcfq 7%$!1)\\-szE6)zF0$!1YQa;c*4b&F6Fhfq7%$!1%4bE?)Qt&)F0$!1a'=@ZAH'fF6F]gq 7%$!1E5o`rrQ'F6Fbgq7%$!1i*>?]@B\")*F0$!1a1B4&>Y#oF6Fggq 7%$!17Sih`:Y5F6$!1@\\5j.0u$F0Fa[xF`dq7%$\" 1K,DbkaKUF0$!1,$*)4s+Q%HF6Fedq7%$\"1in\"p=?Lt%F0$!1_Dx,l4#H$F6Fjdq7%$ \"1Icn*=vQC&F0F`\\xF_eq7%$\"11.(3^m_w&F0Fe\\xFdeq7%$\"1zry!fr&)H'F0$!1 >xylFv!Q%F6Fieq7%$\"1[]:*[#*[%oF0F_]xF^fq7%$\"1>,_i\"e`S(F0Fd]xFcfq7%$ \"1xC?(zE6)zF0$!1ZQa;c*4b&F6Fhfq7%$\"1s]l-#)Qt&)F0F^^xF]gq7%$\"1-5o`rrQ'F6Fbgq7%$\"1P*>?]@B\")*F0$!1b1B4&>Y#oF6Fggq7%$\"14Si h`:Y5F6$!1A\\5j.B<%[nCHF0F]bq7%$\"1C@]%eo Sn#F0$!13eaSMRbeF0Fbbq7%$\"1bh')3)3!=SF0$!1a)=G/7#)z)F0Fgbq7%$\"19CxL@ Dq`F0F^cwF\\cq7%$\"1dIYJHfLnF0$!1R59$=^WZ\"F6Facq7%$\"1n5z1&[36)F0Fhcw Ffcq7%$\"1N.5>\\'[]*F0$!1M0`.^F\"3#F6F[dq7%$\"11EN'G_=4\"F6FbdwF`dq7%$ \"12pQiVZN7F6$!1]<.%\\6`q#F6Fedq7%$\"1RX$4^\\;Q\"F6$!1320A/RDIF6Fjdq7% $\"1*fo()zz1`\"F6FaewF_eq7%$\"1Tgs#>tGo\"F6$!1%H8.Gy\\o$F6Fdeq7%$\"1gB 6'>W&Q=F6$!1lM/I-&e-%F6Fieq7%$\"1%p8')\\9!)*>F6F`fwF^fq7%$\"1#Qf5k8;;# F6$!1_?6POFLZF6Fcfq7%$\"1&o3kpz'HBF6$!1VNtkpG,^F6Fhfq7%$\"13*o$o*fD]#F 6$!1>IHc>%)zaF6F]gq7%$\"1)f)34*z7P0$F6$!1:vL8$)p'o'F6F\\hq7%$\"1:[M^J` \\KF6$!1dWvt][:rF6Fahq7%$\"1*p;5do?X$F6$!1;G-\"yv*evF6Ffhq79Fgaq7%Fifr Fd[wF]bq7%$\"1f06$)\\S:UF0Fi[wFbbq7%$\"1XY[%*R*RL'F0F^\\wFgbq7%$\"1Hhp B?nl%)F0Fc\\wF\\cq7%$\"1,hD7W[h5F6Fh\\wFacq7%$\"1*4*Q\\^fy7F6F]]wFfcq7 %FghrF`]wF[dq7%F\\irFe]wF`dq7%$\"1hr0%H.w%>F6Fj]wFedq7%$\"19vywR.y@F6F _^wFjdq7%$\"1\"owX]lHT#F6Fb^wF_eq7%$\"1JEscL)Gl#F6Fg^wFdeq7%$\"1F.i;LG )*GF6F\\_wFieq7%$\"1cRG'\\s'\\JF6Fa_wF^fq7%$\"1R%zUQqvS$F6Ff_wFcfq7%$ \"1EFB>*4Dn$F6F[`wFhfq7%$\"1fA8\"fQ]%RF6F``wF]gq7%$\"179qs&>dA%F6Fe`wF bgq7%$\"1Z)yJ)G8:XF6Fj`wFggq7%$\"1@<\"pewQ\"[F6F_awF\\hq7%$\"1(fQ#H!oD 7&F6FdawFahq7%$\"1if@7^%=W&F6FiawFfhq79Fgaq7%$\"1R$z2XA[q#F0$!1$p%GRQG QF6$!1->nR#3qB\"F6F[dq7%$\"1QG X-656AF6$!1<.7&R))4U\"F6F`dq7%$\"1$o&o*)*[>]#F6$!1Bn`#)\\!zg\"F6Fedq7% $\"1#)z5$3nzz#F6$!1/]rfU9)z\"F6Fjdq7%$\"1rBm&4n(*4$F6$!1AEp]$*4#*>F6F_ eq7%$\"1qBY)oszS$F6$!1:-8e5#F6Fdeq7%$\"1)Gq&f2ABPF6$!1h)e:poFR#F6Fi eq7%$\"1$))G,wih/%F6$!1Wl?84J+EF6F^fq7%$\"1klcSgYxVF6$!1)p^JgEK\"GF6Fc fq7%$\"1.!e'H_\"yr%F6$!1_;fXd&>.$F6Fhfq7%$\"1n**G\\O\"z1&F6$!1=u0Q.&pD $F6F]gq7%$\"1Af&yw%[GaF6$!1$p`aKv')[$F6Fbgq7%$\"1?GG&pt-!eF6$!1IU&3c4w s$F6Fggq7%$\"11&zFs[S='F6$!1`!Q]zYU(RF6F\\hq7%$\"1u&)[5Hg!e'F6$!1e8\"o p'4HUF6Fahq7%$\"1?DKEdv!*pF6$!1WNA1fo#\\%F6Ffhq79Fgaq7%Fjhq$!1y7\"3&eM e!*F^iqF]bq7%F`iq$!1kE[^5a8=F0Fbbq7%Feiq$!1[Vr[^*\\s#F0Fgbq7%Fjiq$!1(3 e1P!3UOF0F\\cq7%F_jq$!11VvZ=pmXF0Facq7%Fdjq$!1H\\Bj-u+bF0Ffcq7%Fijq$!1 bd\\,e:YkF0F[dq7%F^[r$!1>uwO@*[S(F0F`dq7%Fc[r$!1_uKq/$*y$)F0Fedq7%Fh[r $!1rhfCOGq$*F0Fjdq7%F]\\r$!1zwD?+5Q5F6F_eq7%Fb\\r$!1.N7#)oJT6F6Fdeq7%F g\\r$!1_mvUA*oC\"F6Fieq7%F\\]r$!1'3F1FW]N\"F6F^fq7%Fa]r$!1$39)ek*fY\"F 6Fcfq7%Ff]r$!10$[o3y*z:F6Fhfq7%F[^r$!139^%pCsp\"F6F]gq7%F`^r$!1\"Gs&o& yz\"=F6Fbgq7%Fe^r$!1NdYV#*[U>F6Fggq7%Fj^r$!1&H[p,952#F6F\\hq7%F__r$!1t mz\"[=Q?#F6Fahq7%Fd_r$!1dwkuqkyvV'F`oF`dq7%Fcdq$!1zFi^qFahq7%Fdhq$!1h?'==X`.#Fi^qFfhq-F]aq6&F_ aqF)F`aqF)-F$6$79797%F)F)$\"\"\"F)7%F)F)$\"1x:tNK.,5F,7%F)F)$\"1j%y!y] 8/5F,7%F)F)$\"1^$4p$>J45F,7%F)F)$\"1$Gb+^ul,\"F,7%F)F)$\"1:kT0y$f-\"F, 7%F)F)$\"16a(><@u.\"F,7%F)F)$\"1o:yR$[50\"F,7%F)F)$\"1y#))*ou%o1\"F,7% F)F)$\"1#>X!37&[3\"F,7%F)F)$\"1+`;hn406F,F(7%F)F)$\"1tc7w`[_6F,7%F)F)$ \"1Jk=ljsz6F,7%F)F)$\"1z^6D_S47F,7%F)F)$\"1(G/QN\"F,7%F)F)$\"1KC)p(yp'R \"F,7%F)F)$\"1sg2.xZU9F,7%F)F)$\"1c$z\"o$Q7\\\"F,7%F)F)$\"1W_\"[j!3V:F ,797%F)F)$\"1F$4)=W@)*)*F67%$\"1&GO9/UbB$F^iq$\"1qFz0K7/cF^iq$\"1bci2; W3**F67%$\"1O!p\"Q-xxkF^iq$\"1JmA7\"F0$\"1u%o#3V9R**F67%$\"1LQCiXQ L(*F^iq$\"1(GT\"f;(eo\"F0$\"14D*p'fQ!***F67%$\"1Ucxk7\"4I\"F0$\"1+,t1W C`AF0$\"1BBKZsA15F,7%$\"1c;LpBMFDGF0$\"1Ou$yB&\\:5F,7%$ \"1!QS]C/['>F0$\"1ZpLz29.MF0$\"14DG?<'o-\"F,7%$\"1l@lKj\\-BF0$\"1gXRK1 /))RF0$\"1_\"pJ=]./\"F,7%$\"1&RHKYY\\k#F0$\"1=Hhd;=\"e%F0$\"1,5-+&))f0 \"F,7%$\"1fW>+B'G*HF0$\"1fU*Rc*y$=&F0$\"1kN(p)*3Q2\"F,7%$\"1Egp,G'pM$F 0$\"1$f!eB'4rz&F0$\"1'ea&p%[Q4\"F,7%$\"1SbJ1(zzq$F0$\"1OGxW#4CU'F0F17% $\"13T=X!fm2%F0$\"1-TR\"f!)41(F0$\"1-MuOZvS6F,7%$\"1kXj$oiPX%F0$\"1u$4 HDVTr(F0$\"1&y8=X=x;\"F,7%$\"17Wd&*)p+%[F0$\"1p%>8(pC$Q)F0$\"1S=uB_4(> \"F,7%$\"1SJ$Q(*yjB&F0$\"1wVOIWnp!*F0$\"1P%)yed%*G7F,7%$\"1c:\")z)3Nk& F0$\"1L(*z7T%[x*F0$\"1N`bgM'*F67%$\"1gjN[D'o#>F0$\"1e[&H$QUPLF0$\"17w \\;rF%o*F67%$\"1[(=xsR`d#F0$\"1o%f^D>1Y%F0$\"1Xk`&*4'Rv*F67%$\"1%ROJxQ \"HKF0$\"1jGw>C.$f&F0$\"1y_ic6!Q%)*F67%$\"1#Hz;u5'*)QF0$\"1'>aiU.qt'F0 $\"1WK*=D$)R&**F67%$\"1^E#)z/7eXF0$\"1gI2ei*[*yF0$\"1(f\\p\\t%35F,7%$ \"1&=jQW\\gB&F0$\"1*z8xm.\"p!*F0$\"1jET]EjB5F,7%$\"1H\"fYb)zCfF0$\"1E! G@@0i-\"F6$\"1hl`6r!4/\"F,7%$\"19Kf#4\"zDmF0$\"1v3F\"o?w9\"F6$\"1WA]xD Lg5F,7%Fjjl$\"1'4-en29F\"F6F77%$\"1i@Ko!H.2)F0$\"1)R8!)*>#yR\"F6$\"1Wg e?;!e5\"F,7%$\"178i-'fo\"))F0$\"1lT\\z[7F:F6$\"1_T=e\"R>8\"F,7%$\"1>Ro D*4;e*F0$\"1!3f4_$ef;F6$\"1u1v:eTg6F,7%$\"1.HAPghO5F6$\"1z2TW;Z&z\"F6$ \"1Qe[R/H\">\"F,7%$\"1%eO$>K@<6F6$\"1b&y(e+2N>F6$\"1H&=4$oiC7F,7%$\"1$ G'o!4>,?\"F6$\"1.??Ssmy?F6$\"1E.vyQ\\g7F,7%$\"1t#*)Q(\\]&G\"F6$\"1A(e \"G*flA#F6$\"1U?@,d'*)H\"F,7%$\"1ZlV;tat8F6$\"1')zmQP0zBF6$\"1HpA*z@,M \"F,7%$\"1JlJb!GWY\"F6$\"1:YMnGF6$\" 1De#o-v0[\"F,797%F)F)$\"1%=XN&*>j4*F67%$\"1H6i]<]W%*F^iq$\"1+\\U)oNej \"F0$\"1[qM'f=d5*F67%$\"1j$f%H?&3*=F0$\"1R$=Yy^]F$F0$\"1w(zn%R$R8*F67% $\"1Rlhv7;TGF0$\"1Y]'=oN5#\\F0$\"1p)f5JC5=*F67%$\"1'*>Z+ng&*F67%$\"1M'Rtbh0s(F0$\"1Mv-)[SsL\"F6$\"1a sa>fQ/(*F67%$\"1q!fi*H7O()F0$\"1c\"p!*)398:F6$\"1y!p.p`\"o)*F67%$\"1&) )[SQP(p(*F0$\"1sxd_\"o@p\"F6$\"12'=>PJ_+\"F,7%$\"1NX`lSN#3\"F6$\"1\\XG LApu=F6F<7%$\"1QOXY5(**=\"F6$\"11O*>J!4h?F6$\"1)4-BdP$[5F,7%$\"1>Aeoq/ +8F6$\"1$Q:]dZB8JKtk\"F6$\"1\\;;gKE`GF6$\"1wj\"f\"z)4;\"F,7%$\"1X'fR7x&p\"F,7%$\"12YA<(ya*=F6$\"1*\\)*GWlIG$F6$\"1&)*z ]pj9B\"F,7%$\"1Ix;jsHD?F6$\"1.nBgx\"z]$F6$\"1]puw5[q7F,7%$\"1&f4\"G5If @F6$\"1I:_*3>+u$F6$\"1!\\.8(Q778F,7%$\"14k+Fpw(H#F6$\"11l_h!\\)zRF6$\" 1/iRK\"ykN\"F,7%$\"1hoo&4\")4W#F6$\"1!zW!yK!zA%F6$\"12*Ro^NOS\"F,797%F )F)$\"17=JG``7%)F67%$\"1#>&Q!H_\"H7F0$\"11o= b)F67%$\"1<5IF>n'>'F0$\"16@fh]Ht5F6$\"1$G;F.Q2j)F67%$\"15NZ`$3TY(F0$\" 1R:9!\\@GH\"F6$\"1#z\"efAMF()F67%$\"1_jbY\"ppu)F0$\"1)H2u\\>]^\"F6$\"1 9f%))F67%$\"1lpF6$\"1D>k8\"\\j7*F67%$\"1Uhf***z9F\"F6$\"1u\" 3,'zE-AF6$\"1[g]Wbm'H*F67%Fir$\"1+L1jM\")RCF6FA7%$\"1*Qir![o[:F6$\"1n_ \"4x+Co#F6$\"1[*>?%RK&p*F67%$\"1K@hnH%>p\"F6$\"1vIT$=J0$HF6$\"1O<([s*[ C**F67%$\"1H<#*zupQ=F6$\"1!Ht\\XF6$ \"1>CFa[[XMF6$\"1fD,Ej[W5F,7%$\"1q`'zY;R9#F6$\"1L'*p\\AP8PF6$\"1W)e\\j 9P2\"F,7%$\"1<#G+n6II#F6$\"1@H\"HAL*))RF6$\"1VybW<;06F,7%$\"1!\\')*zf' oY#F6$\"1*y*z7stsUF6$\"1TTORE*)Q6F,7%$\"1Ke4)*z\"ej#F6$\"1mk?#pq`c%F6$ \"1#G4L-x\\<\"F,7%$\"13.))eo@5GF6$\"15t')z$Qu'[F6$\"1aHuj%*[87F,7%$\"1 P*3E&HU!*HF6$\"1:)Q+\\k&z^F6$\"1*GCba4XD\"F,7%$\"1%*=.q'3o<$F6$\"1`0UA SR-bF6$\"11<[N?7)H\"F,797%F)F)$\"1%eUNu&\\dvF67%$\"1I&e'4C$))[\"F0$\"1 uP[xLtyDF0$\"1cg==WIlvF67%$\"1S]+\\9u!)HF0$\"1.K!Rj&zi^F0$\"1'y>&ylu)e (F67%$\"1e\\\\0,\")yWF0$\"1!>q,mEvv(F0$\"1T\"*4n1(yi(F67%$\"1vnD&4Mh)f F0$\"130RR)Go.\"F6$\"1!pOF`dFo(F67%$\"1b5Ig#Ge](F0$\"1uF(**eZ+I\"F6$\" 1xRd(f?Nv(F67%$\"1%e__+L5/*F0$\"1h/!=&G&fc\"F6$\"1I^Y\"41.%yF67%$\"1\" p!py?\\f5F6$\"1S!406%4N=F6$\"1'pcTN$HVzF67%$\"1w@+n-2<7F6$\"1)=H)Qv-3@ F6$\"1#4))e?&pi!)F67%$\"1w!*G'\\jrP\"F6$\"1s=!f9<`Q#F6$\"1p&yrQe()>)F6 7%$\"1a6k%f-,a\"F6$\"1MryUf`nEF6$\"1]%4(oSw^$)F67%$\"1!Rb5FCiq\"F6$\"1 1^dErEbHF6FF7%$\"1.bB6=(e(=F6$\"1Wha&G0\"\\KF6$\"1W[cCL!*4()F67%$\"1r? )py&R\\?F6$\"1\\0(yid'\\NF6$\"1>9W$)px:*)F67%$\"1Y\"R1yarA#F6$\"1;pQO_ adQF6$\"1Eh-T[2S\"*F67%$\"1SW4Eh^4CF6$\"12@0_VStTF6$\"164W-/E$Q*F67%$ \"1?V'ymcof#F6$\"1=]r)o()y\\%F6$\"1!p3P?Oek*F67%$\"1WA$)RNc*y#F6$\"1L6 $>yl;$[F6$\"1%pH-&[MG**F67%$\"1)))\\_'\\.))HF6$\"14^gu$Ga<&F6$\"11g!Q, PJ-\"F,7%$\"1HTlz5o#>$F6$\"14[KUe))HbF6$\"1q4.EQbb5F,7%$\"1I/SyZ#RS$F6 $\"1u@Lh1x&*eF6$\"1]_65>:!4\"F,7%$\"1?Sc!f-Ai$F6$\"1*[k@#*QQF'F6$\"1(G 7?w-q7\"F,7%$\"1vZ,\"elz%QF6$\"1$)Q4#*=([m'F6$\"1>f&H`#=m6F,797%F)F)$ \"1^GXRtg[lF67%$\"1[$4X>/#=%yb(e4mF67%$\"1Dl3;qO3pF0$\"1'[>pEkl>\"F6$\"18([wHZrl'F67%$\"1p E%\\d(=i')F0$\"1a#\\%)*[L+:F6$\"1)3sB%QY=nF67%$\"1#4B(H3RV5F6$\"1V=BHf ?2=F6$\"1d>!4!Rm$z'F67%$\"1ekKR!>FA\"F6$\"1L_\\*\\6y6#F6$\"12L,sG!H)oF 67%$\"1[?$H'Rd/9F6$\"1rmymMzKCF6$\"1PL\\l^O')pF67%$\"1F*\\&)RJ$*e\"F6$ \"1tHNKF!Gv#F6$\"1N!>KekU5(F67%$\"1/%[@9ttx\"F6$\"1v$4ie+&yIF6$\"1eE(4 wWoB(F67%Febl$\"1&[wO2g0T$F6FK7%$\"1b%*zu9([;#F6$\"1A#4iQn'\\PF6$\"1f' 3F>ura(F67%$\"1s%)4=)G^O#F6$\"1#)=RzK_'4%F6$\"1**egLZcDxF67%$\"1ungKOF qDF6$\"1`&H5_W=X%F6$\"1*)=cL*>*>zF67%$\"1d$HU%)H2y#F6$\"1M$GBPlj\"[F6$ \"1a50A9kI\")F67%$\"1#y@ZNKp*HF6$\"1e^A0\"R3>&F6$\"1Jg,\\Y;e$)F67%$\"1 WASSzK>KF6$\"1cF+f&Rgd&F6$\"1YDr%yfHg)F67%$\"1%3-a\\MHYC!=%F6$\"1GB4'3 'RSsF6$\"1DItAX`l(*F67%$\"1(f%Q4!*ySWF6$\"1%eZ**=s;p(F6$\"1n=0F0$\"1DohCur7LF0$\"1?9moU*>T&F67 %$\"1+V]B*[\"HQF0$\"1L\"F6$\"1?\"37u;g \\&F67%$\"1^YCI(4Ak*F0$\"1t=+:(z+n\"F6$\"1E$*yP&Qma&F67%$\"1]%=0$zVh6F 6$\"10c\\0&p;,#F6$\"10`i/Bs3cF67%$\"1w?Nm\\0h8F6$\"1Wdh`jTdBF6$\"1a&[g L'R#o&F67%$\"1aMV$e%[j:F6$\"1.uNNZ.3FF6$\"1#eUt(G\"yw&F67%$\"1z&*y'4X \"p\\U1$F6$\"1_CNR%[^'eF67%$\"1+C+)\\h%y>F6$\"1()*e`$ezEMF 6$\"17%RK;/Y(fF67%$\"1$zTMLm=>#F6$\"1\"*3#H&QU'z$F6FP7%$\"1$H5\")f+)4C F6$\"1&*o.sq*Q<%F6$\"1c?aej!3B'F67%$\"1\"oYpk9Fj#F6$\"14CMIb**fXF6$\"1 JJ[kA3yjF67%$\"18;ZD\"p5'GF6$\"1)[$G*39Hg]7nF67%$\"1YjFT\"**fL$F6$\"1unp1+7ydF6 $\"1()e-l\\M+pF67%$\"1,3)y5cNe$F6$\"1ALi]7!p?'F6$\"15g'oEVC5(F67%$\"11 )=4R=&QQF6$\"1/*e)y)3&[mF6$\"1(yN]c=#>tF67%$\"1RCniGT,TF6$\"1y^'4YbQ5( F6$\"1+G]>)=6b(F67%$\"1+a1(y#ysVF6$\"14*fw&>)Qd(F6$\"12__XKi)z(F67%$\" 11CIU*)=`YF6$\"1]aGl'f&f!)F6$\"1'QORIVA1)F67%$\"1X8qz2@V\\F6$\"1$*['HA #*=c)F6$\"1hT1eP_U$)F6797%F)F)$\"1k)=+8]T:%F67%$\"1h\"[LTd!o?F0$\"1&[) )G^!)>e$F0$\"1PvOYBWeTF67%$\"1s^K3%)QSTF0$\"1J%)G'3j8<(F0$\"1^h:lyKrTF 67%$\"1hVz(Rv7A'F0$\"16eo2lbx5F6$\"1s:%RJLG>%F67%$\"1@vLv%=]J)F0$\"1\\ e^YM?S9F6$\"1=+7LJ+BUF67%$\"1dv7IWfU5F6$\"1[#pY_Ee!=F6$\"1q=wn'**=E%F6 7%$\"1k)eV]TeD\"F6$\"1%G.=H\"=v@F6$\"1yfW'H.'4VF67%$\"1v6ARPor9F6$\"1? S+45.\\DF6$\"15(ftf7iO%F67%$\"1uL`mrc!p\"F6$\"1fN/E#[\"GHF6$\"1_z+_X%= V%F67%$\"1hF6$\"1.@$yP;LJ$F6$\"1c?&pyMm]%F67%$\"1)))oJ3u#R@F6 $\"1SQO.9L0PF6$\"1d#yS&yt!f%F67%Fjx$\"1g@]!R.]5%F6FU7%$\"1tY`IJn0EF6$ \"10J+]#eJ^%F6$\"1.t2vsf(y%F67%$\"1d^q(*fqYGF6$\"1$QOAUR1$\\F6$\"1'*eK L/w+\\F67%$\"1*)Qe,:i$4$F6$\"1<38?'4$e`F6$\"1:x:t30C]F67%$\"1ZM,'))HpM $F6$\"1+g%>h_qz&F6$\"1^$[2PBx:&F67%$\"1I;a0Y:2OF6$\"1uTbZ]xZiF6$\"1&[* \\cT0-`F67%$\"1_`72M$[(QF6$\"1Gj+K$39r'F6$\"1D6'e[TtX&F67%$\"1x$p8W>0: %F6$\"1`TD]0\"*)=(F6$\"1))zmbi!Ri&F67%$\"1hif0CyMWF6$\"1\\_+Z%o7o(F6$ \"1VOEnE4-eF67%$\"14'46s4#GZF6$\"16H!eY*\\*=)F6$\"1!y#yP*oA*fF67%$\"1T M3ZxSJ]F6$\"1wTa^Ql9()F6$\"1\"Rh51G[>'F67%$\"19heLI+X`F6$\"1F4y?o\"yD* F6$\"1Ptu>')=5kF6797%F)F)$\"1'H9%obK\"HilF0$\"1@f'e :F6$\"1&)pPf8-kGF67%$\"16tArOu*4\"F6$\"1D!oQ\">\"[!>F6$\"1=D!3\"3S!*GF 67%$\"1s$=nwzYK\"F6$\"12ZrfET%H#F6$\"1^1-$>`F#HF67%$\"1v^9vKN_:F6$\"1# RZQua()o#F6$\"12gIh`9hHF67%$\"1,eJxYB$y\"F6$\"1e'HH1`')3$F6$\"1a(y>lcc +$F67%$\"12#fI3,y,#F6$\"1agU&*R$\\\\$F6$\"1c3NY!zj0$F67%$\"1x!fr@PlD#F 6$\"11'*G4rV3RF6$\"1EG!3OsBDtOK$F67%$\"1I;>lm>jKF6$\"1H7?>C-_cF6$\"1T3Hj#)G2MF67%$\"1 ')zf4+RINF6$\"1(Gdr'[\"[6'F6$\"1@BP;GM#>*F6$\"17@F$)RI,UF67%$\"1\"*HU Rk)zj&F6$\"1n\\Yl*y_w*F6$\"1]**eE0OZVF6797%F)F)$\"1`Gt#Q[JU\"F67%$\"1= &o(>*o.D#F0$\"1iH\"[I`x*QF0$\"1HK\\I)=YU\"F67%$\"138pi\")Q0XF0$\"1Q2g0 7c.yF0$\"1$p7C@L!H9F67%$\"1KFohwrpnF0$\"1(3$H%\\\\D<\"F6$\"12bo]1SO9F6 7%$\"1Tc.*eO![!*F0$\"1$>&43f;n:F6$\"1>E\"*pjtY9F67%$\"1P(zDI0X8\"F6$\" 12W_D3-l>F6$\"1kGcG<1g9F67%$\"1)[)3p8bm8F6$\"1rST-k$pO#F6$\"1G'4JE/kZ \"F67%$\"1(R/]P@9g\"F6$\"1EW$e=VPx#F6$\"1)=M\\u(z&\\\"F67%$\"1PA\\s1gR =F6$\"1GWRI=G'=$F6$\"11L6]AG=:F67%$\"1#QtAX\"e\"3#F6$\"130VN[S0OF6$\"1 (HNAC/Ra\"F67%$\"1$)QsgP'yK#F6$\"1B@n1$y>.%F6$\"16dTomrs:F67%Fa^l$\"1Q X-YP)oY%F6Fin7%$\"1l(odky`$GF6$\"1Rt;t)>5\"\\F6$\"1:l:*pd,k\"F67%$\"1? -R()*fw4$F6$\"1FIt#[/`O&F6$\"1xBooc#*y;F67%$\"1r\"y\\bUjO$F6$\"1=M!4Mw 1$eF6$\"1H8>!4j6s\"F67%$\"1*=\"Rt:)>k$F6$\"1xJMDr43jF6$\"1(f4hCdpw\"F6 7%$\"1b\\zUm9DRF6$\"1n#zCT`&)z'F6$\"1kNvoFS;=F67%$\"1`Zs:HU;UF6$\"1wHB ;(eII(F6$\"1cZ)e$=gp=F67%$\"1)['G1BT;XF6$\"1DvT#ebE#yF6$\"1ls&=Qkm#>F6 7%$\"10]*RtMd#[F6$\"1F#HJu<%e$)F6$\"1eRzC$3x)>F67%$\"1v+#[SH]9&F6$\"1H UeNBX6*)F6$\"194E5)fG0#F67%$\"1+t3Kh&\\Z&F6$\"1#)f)**=-H[*F6$\"1IO%>Z` A7#F67%$\"1/N'**p'>;eF6$\"1B^6>[R25F,$\"1*pj)4F.'>#F6797%F)F)$\"1mnt&* RBBhF`o7%$\"1`Y;t*4NF#F0$\"1Ld.&yMy$RF0$\"1rR*)R2cHhF`o7%$\"1\\+\")*3= ?K$RoF0$\"1i*pj32Y=\"F6$ \"1hiQkJD!='F`o7%$\"1N^+*fz59*F0$\"1B3:I9G$e\"F6$\"1d*p\"*fBZA'F`o7%$ \"1oZhx;F6$\"1ltjZs0#G'F`o7%$\"1*o#RERg!Q\"F6$\"1>7 kXhF\"R#F6$\"1KF6$\"1^ZP!pcD`'F`o7%$\"1%[/f# o)H5#F6$\"1rX6I+[UOF6$\"1>%=uz(zUmF`o7%$\"1\"4D6r,=N#F6$\"1\"3'*H0SM2% F6$\"1ffacgwmnF`o7%$\"1P(ou_wag#F6$\"1!)R]Bx\"G^%F6F^o7%$\"1;Vq'[NX'GF 6$\"1;f0.5_h\\F6$\"1#>#oMw$p0(F`o7%$\"1k$yn#R^HJF6$\"17Xs(4.zF6$\"1g=#\\Vg'*G)F`o7%$\"1YrZd)e`([F6 $\"1n8GY#pVW)F6$\"1)>7@#zI_&)F`o7%$\"1)pMx!p$z>&F6$\"1aRt=34.!*F6$\"1# pn&4&GE$))F`o7%$\"1k$ziMc7`&F6$\"1[OW@qT!e*F6$\"1aGyr9?J\"*F`o7%$\"13! >#of+weF6$\"1FLD)3ax,\"F,$\"1sPE,Qk[%*F`o797%F)F)$!1]Gt#Q[JU\"F67%Ffj_ l$\"1jH\"[I`x*QF0$!1EK\\I)=YU\"F67%F][`lF_[`l$!1*o7C@L!H9F67%Fd[`lFf[` l$!1.bo]1SO9F67%$\"1Uc.*eO![!*F0F]\\`l$!1;E\"*pjtY9F67%Fb\\`lFd\\`l$!1 hGcG<1g9F67%$\"1*[)3p8bm8F6F[]`l$!1D'4JE/kZ\"F67%F`]`lFb]`l$!1&=M\\u(z &\\\"F67%Fg]`lFi]`l$!1.L6]AG=:F67%$\"1$QtAX\"e\"3#F6F`^`l$!1%HNAC/Ra\" F67%Fe^`l$\"1C@n1$y>.%F6$!12dTomrs:F67%Fb`lF\\_`lFdo7%$\"1m(odky`$GF6F a_`l$!16l:*pd,k\"F67%Ff_`l$\"1GIt#[/`O&F6$!1tBooc#*y;F67%$\"1s\"y\\bUj O$F6F_``l$!1E8>!4j6s\"F67%Fd``l$\"1yJMDr43jF6$!1$f4hCdpw\"F67%$\"1c\\z Um9DRF6$\"1o#zCT`&)z'F6$!1gNvoFS;=F67%Fba`l$\"1xHB;(eII(F6$!1_Z)e$=gp= F67%$\"1*['G1BT;XF6$\"1EvT#ebE#yF6$!1hs&=Qkm#>F67%F`b`l$\"1G#HJu<%e$)F 6$!1aRzC$3x)>F67%$\"1w+#[SH]9&F6$\"1IUeNBX6*)F6$!144E5)fG0#F67%F^c`l$ \"1$)f)**=-H[*F6$!1DO%>Z`A7#F67%$\"10N'**p'>;eF6Fgc`l$!1%pj)4F.'>#F679 7%F)F)$!1(H9%obK`F#HF67%Fdc_l$\"1\"RZQua()o#F6$!12gIh`9hHF67%$\"1+eJxYB$y\"F6F]d_l$! 1b(y>lcc+$F67%Fbd_lFdd_l$!1d3NY!zj0$F67%Fid_l$\"10'*G4rV3RF6$!1FG!3OsBDtOK$F67%Fcf_l$\"1G7?>C-_cF6$!1U3 Hj#)G2MF67%Fjf_l$\"1'Gdr'[\"[6'F6$!1ABP;GM#>*F6$!18@F$)RI,UF67%$\"1*)HURk)zj&F6$\"1m\\Yl*y_w*F6$!1^**eE0OZ VF6797%F)F)$!1j)=+8]T:%F67%F`g^lFbg^l$!1OvOYBWeTF67%$\"1t^K3%)QSTF0$\" 1K%)G'3j8<(F0$!1]h:lyKrTF67%$\"1iVz(Rv7A'F0F`h^l$!1r:%RJLG>%F67%$\"1Bv Lv%=]J)F0Fgh^l$!1<+7LJ+BUF67%F\\i^lF^i^l$!1p=wn'**=E%F67%Fci^lFei^l$!1 xfW'H.'4VF67%Fji^lF\\j^l$!14(ftf7iO%F67%Faj^lFcj^l$!1^z+_X%=V%F67%Fhj^ lFjj^l$!1b?&pyMm]%F67%F_[_lFa[_l$!1c#yS&yt!f%F67%FjxFf[_lF\\p7%$\"1uY` IJn0EF6F[\\_l$!1,t2vsf(y%F67%F`\\_lFb\\_l$!1&*eKL/w+\\F67%$\"1!*Qe,:i$ 4$F6$\"1=38?'4$e`F6$!19x:t30C]F67%F^]_lF`]_l$!1\\$[2PBx:&F67%Fe]_l$\"1 vTbZ]xZiF6$!1%[*\\cT0-`F67%$\"1``72M$[(QF6$\"1Hj+K$39r'F6$!1A6'e[TtX&F 67%$\"1y$p8W>0:%F6$\"1aTD]0\"*)=(F6$!1&)zmbi!Ri&F67%Fj^_l$\"1^_+Z%o7o( F6$!1UOEnE4-eF67%Fa__l$\"17H!eY*\\*=)F6$!1yFyP*oA*fF67%$\"1UM3ZxSJ]F6$ \"1xTa^Ql9()F6$!1!Rh51G[>'F67%$\"1:heLI+X`F6$\"1H4y?o\"yD*F6$!1Otu>')= 5kF6797%F)F)$!1sfbu\"3kS&F67%$\"1*o^oH)f7>F0$\"1FohCur7LF0$!1:9moU*>T& F67%$\"1,V]B*[\"HQF0$\"1>rD&[!GKmF0$!1.=[%4k(GaF67%$\"1-b8xBh`dF0$\"1J eDk*[b'**F0$!1QOQ1BvcaF67%$\"15Tb1b'**o(F0F[_^l$!1:\"37u;g\\&F67%$\"1b YCI(4Ak*F0$\"1u=+:(z+n\"F6$!1@$*yP&Qma&F67%Fg_^l$\"11c\\0&p;,#F6$!1,`i /Bs3cF67%$\"1x?Nm\\0h8F6$\"1Xdh`jTdBF6$!1]&[gL'R#o&F67%$\"1bMV$e%[j:F6 $\"1/uNNZ.3FF6$!1xDMxG\"yw&F67%$\"1!e*y'4X\"p\\U1$F6$!1[CN R%[^'eF67%$\"1,C+)\\h%y>F6$\"1))*e`$ezEMF6$!11%RK;/Y(fF67%$\"1%zTMLm=> #F6$\"1$*3#H&QU'z$F6Fap7%$\"1%H5\")f+)4CF6$\"1(*o.sq*Q<%F6$!1]?aej!3B' F67%$\"1#oYpk9Fj#F6$\"16CMIb**fXF6$!1FJ[kA3yjF67%$\"19;ZD\"p5'GF6$\"1* [$G*39Hg]7 nF67%$\"1ZjFT\"**fL$F6$\"1wnp1+7ydF6$!1#)e-l\\M+pF67%$\"1-3)y5cNe$F6$ \"1DLi]7!p?'F6$!1/g'oEVC5(F67%$\"13)=4R=&QQF6$\"11*e)y)3&[mF6$!1\"yN]c =#>tF67%$\"1TCniGT,TF6$\"1\"=l4YbQ5(F6$!1$z-&>)=6b(F67%$\"1-a1(y#ysVF6 $\"18*fw&>)Qd(F6$!1+__XKi)z(F67%$\"13CIU*)=`YF6$\"1`aGl'f&f!)F6$!1zj$R IVA1)F67%$\"1Z8qz2@V\\F6$\"1'*['HA#*=c)F6$!1aT1eP_U$)F6797%F)F)$!1]GXR tg[lF67%Fjc]lF\\d]l$!1V![Xft`b'F67%Fad]lFcd]l$!1&o1@M'ovlF67%Fhd]lFjd] l$!1[>%yb(e4mF67%F_e]lFae]l$!17([wHZrl'F67%Ffe]lFhe]l$!1(3sB%QY=nF67%F ]f]lF_f]l$!1c>!4!Rm$z'F67%Fdf]lFff]l$!11L,sG!H)oF67%F[g]lF]g]l$!1OL\\l ^O')pF67%Fbg]lFdg]l$!1M!>KekU5(F67%Fig]lF[h]l$!1cE(4wWoB(F67%FeblF`h]l Fdp7%Fch]lFeh]l$!1e'3F>ura(F67%Fjh]lF\\i]l$!1(*egLZcDxF67%Fai]lFci]l$! 1))=cL*>*>zF67%Fhi]lFji]l$!1`50A9kI\")F67%F_j]lFaj]l$!1Ig,\\Y;e$)F67%F fj]lFhj]l$!1XDr%yfHg)F67%F][^lF_[^l$!1QGI\"pKb'))F67%Fd[^lFf[^l$!1p`?o fUY\"*F67%F[\\^lF]\\^l$!1d5Ot+AY%*F67%Fb\\^lFd\\^l$!1BItAX`l(*F67%Fi\\ ^lF[]^l$!1n=0&ylu)e(F67% $\"1g\\\\0,\")yWF0$\"1$>q,mEvv(F0$!1R\"*4n1(yi(F67%$\"1xnD&4Mh)fF0$\"1 40RR)Go.\"F6$!1)oOF`dFo(F67%$\"1e5Ig#Ge](F0$\"1vF(**eZ+I\"F6$!1vRd(f?N v(F67%$\"1(e__+L5/*F0Fc\\]l$!1G^Y\"41.%yF67%$\"1#p!py?\\f5F6Fj\\]l$!1% pcTN$HVzF67%$\"1x@+n-2<7F6$\"1*=H)Qv-3@F6$!1!4))e?&pi!)F67%Ff]]l$\"1t= !f9<`Q#F6$!1m&yrQe()>)F67%$\"1b6k%f-,a\"F6$\"1NryUf`nEF6$!1Z%4(oSw^$)F 67%Fd^]l$\"12^dErEbHF6Fip7%Fi^]l$\"1Xha&G0\"\\KF6$!1U[cCL!*4()F67%F`_] l$\"1]0(yid'\\NF6$!1;9W$)px:*)F67%Fg_]l$\"1=pQO_adQF6$!1Bh-T[2S\"*F67% $\"1TW4Eh^4CF6$\"13@0_VStTF6$!134W-/E$Q*F67%$\"1@V'ymcof#F6$\"1>]r)o() y\\%F6$!1'o3P?Oek*F67%F\\a]l$\"1M6$>yl;$[F6$!1\"pH-&[MG**F67%$\"1*))\\ _'\\.))HF6$\"16^gu$Ga<&F6$!11g!Q,PJ-\"F,7%$\"1JTlz5o#>$F6$\"16[KUe))Hb F6$!1p4.EQbb5F,7%$\"1J/SyZ#RS$F6$\"1w@Lh1x&*eF6$!1]_65>:!4\"F,7%$\"1@S c!f-Ai$F6$\"1\"\\k@#*QQF'F6$!1(G7?w-q7\"F,7%$\"1wZ,\"elz%QF6$\"1&)Q4#* =([m'F6$!1>f&H`#=m6F,797%F)F)$!18=JG``7%)F67%$\"1\">&Q!H_\"H7F0$\"10o= b)F67%$\"1:5IF>n'>'F0Fbb\\l$!1%G;F.Q2j)F67%$\"1 1NZ`$3TY(F0Fib\\l$!1%z\"efAMF()F67%$\"1[jbY\"ppu)F0$\"1(H2u\\>]^\"F6$! 1:f%))F67%$\"1kpF6$!1E>k8\"\\j7*F67%Fcd\\l$\"1t\"3,'zE-AF6$!1 \\g]Wbm'H*F67%Fhu$\"1*HjIY8)RCF6F^q7%F]e\\l$\"1m_\"4x+Co#F6$!1\\*>?%RK &p*F67%$\"1J@hnH%>p\"F6$\"1uIT$=J0$HF6$!1Q<([s*[C**F67%$\"1G<#*zupQ=F6 $\"1*Gt\\XF6$\"1H\"HAL*))RF6$!1VybW<;06F,7%$\"1*[')*zf'oY#F6$\"1(y*z7stsUF6$!1TTOR E*)Q6F,7%$\"1Je4)*z\"ej#F6$\"1lk?#pq`c%F6$!1#G4L-x\\<\"F,7%$\"12.))eo@ 5GF6$\"13t')z$Qu'[F6$!1aHuj%*[87F,7%$\"1O*3E&HU!*HF6$\"18)Q+\\k&z^F6$! 1*GCba4XD\"F,7%$\"1$*=.q'3o<$F6$\"1]0UASR-bF6$!11<[N?7)H\"F,797%F)F)$! 1%=XN&*>j4*F67%$\"1F6i]<]W%*F^iq$\"1**[U)oNej\"F0$!1[qM'f=d5*F67%F_g[l Fag[l$!1w(zn%R$R8*F67%$\"1Qlhv7;TGF0$\"1X]'=oN5#\\F0$!1p)f5JC5=*F67%F] h[l$\"1$\\<*e()=xlF0$!1#H\"3+q3Z#*F67%$\"1FteB#p8w%F0Ffh[l$!1&p00`eAL* F67%$\"1U(G.)fBNdF0$\"1wkK6,sL**F0$!1l,(o!\\rO%*F67%$\"1mzbQW&4s'F0Fdi [l$!16cH&)>ng&*F67%$\"1K'Rtbh0s(F0$\"1Lv-)[SsL\"F6$!1asa>fQ/(*F67%$\"1 p!fi*H7O()F0Fbj[l$!1y!p.p`\"o)*F67%$\"1%))[SQP(p(*F0$\"1rxd_\"o@p\"F6$ !12'=>PJ_+\"F,7%F^[\\lF`[\\lFaq7%Fc[\\lFe[\\l$!1)4-BdP$[5F,7%Fj[\\l$\" 1#Q:]dZ\"F,7%F]^\\l$\"1) \\)*GWlIG$F6$!1&)*z]pj9B\"F,7%Fd^\\lFf^\\l$!1]puw5[q7F,7%F[_\\l$\"1H:_ *3>+u$F6$!1!\\.8(Q778F,7%$\"13k+Fpw(H#F6Fd_\\l$!1/iRK\"ykN\"F,7%Fi_\\l $\"1*yW!yK!zA%F6$!12*Ro^NOS\"F,797%F)F)$!1t\\9O(H\\f*F67%$\"1*oTux<_S' F^iq$\"1SsSiiT46F0$!1L.RgN%[g*F67%$\"1vNYvrO#G\"F0$\"10:'=5D6A#F0$!1,8 w>bgM'*F67%$\"1jjN[D'o#>F0$\"1j[&H$QUPLF0$!15w\\;rF%o*F67%$\"1^(=xsR`d #F0$\"1u%f^D>1Y%F0$!1Wk`&*4'Rv*F67%$\"1)ROJxQ\"HKF0$\"1rGw>C.$f&F0$!1x _ic6!Q%)*F67%$\"1(Hz;u5'*)QF0$\"10UDEM+PnF0$!1VK*=D$)R&**F67%$\"1eE#)z /7eXF0$\"1rI2ei*[*yF0$!1(f\\p\\t%35F,7%$\"1#>jQW\\gB&F0$\"16QrnO5p!*F0 $!1jET]EjB5F,7%$\"1P\"fYb)zCfF0$\"1F!G@@0i-\"F6$!1hl`6r!4/\"F,7%$\"1AK f#4\"zDmF0$\"1w3F\"o?w9\"F6$!1WA]xDLg5F,7%Fihl$\"1)4-en29F\"F6Ffq7%$\" 1t@Ko!H.2)F0$\"1+M,)*>#yR\"F6$!1Wge?;!e5\"F,7%$\"1C8i-'fo\"))F0$\"1nT \\z[7F:F6$!1^T=e\"R>8\"F,7%$\"1LRoD*4;e*F0$\"1#3f4_$ef;F6$!1u1v:eTg6F, 7%$\"10HAPghO5F6$\"1\"y5Wkraz\"F6$!1Qe[R/H\">\"F,7%$\"1'eO$>K@<6F6$\"1 e&y(e+2N>F6$!1H&=4$oiC7F,7%$\"1&G'o!4>,?\"F6$\"11??Ssmy?F6$!1E.vyQ\\g7 F,7%$\"1u#*)Q(\\]&G\"F6$\"1D(e\"G*flA#F6$!1U?@,d'*)H\"F,7%$\"1\\lV;tat 8F6$\"1*)zmQP0zBF6$!1HpA*z@,M\"F,7%$\"1LlJb!GWY\"F6$\"1@v$ezjk`#F6$!1r qp?s/%Q\"F,7%$\"1sE)H*\\Le:F6$\"1x%eIQ:\"*p#F6$!1TCDOF$3V\"F,7%$\"1W=W &=ial\"F6$\"1%G%>:YMnGF6$!1Ce#o-v0[\"F,797%F)F)$!1F$4)=W@)*)*F67%$\"1' GO9/UbB$F^iq$\"1rFz0K7/cF^iq$!1bci2;W3**F67%$\"1P!p\"Q-xxkF^iq$\"1KmA< F)>7\"F0$!1u%o#3V9R**F67%$\"1MQCiXQL(*F^iqFbdz$!14D*p'fQ!***F67%FgdzFi dz$!1BBKZsA15F,7%F^ez$\"1M\\e>MFDGF0$!1Ou$yB&\\:5F,7%Feez$\"1[pLz29.MF 0$!14DG?<'o-\"F,7%$\"1m@lKj\\-BF0$\"1hXRK1/))RF0$!1_\"pJ=]./\"F,7%Fcfz Fefz$!1,5-+&))f0\"F,7%$\"1gW>+B'G*HF0$\"1hU*Rc*y$=&F0$!1kN(p)*3Q2\"F,7 %$\"1Fgp,G'pM$F0$\"1%f!eB'4rz&F0$!1'ea&p%[Q4\"F,7%Fhgz$\"1QGxW#4CU'F0F [r7%$\"15T=X!fm2%F0$\"1.TR\"f!)41(F0$!1-MuOZvS6F,7%$\"1lXj$oiPX%F0$\"1 w$4HDVTr(F0$!1&y8=X=x;\"F,7%$\"19Wd&*)p+%[F0$\"1q%>8(pC$Q)F0$!1S=uB_4( >\"F,7%Fbiz$\"1xVOIWnp!*F0$!1P%)yed%*G7F,7%$\"1d:\")z)3Nk&F0$\"1M(*z7T %[x*F0$!1N`J45F,7%$\"1\"GnrQ f%>6F`o$\"1\\OeN0'*Q>F`o$!1$Gb+^ul,\"F,7%$\"1yUEwal.9F`o$\"1=?@,E?JCF` o$!1:kT0y$f-\"F,7%$\"1uqMz@v!p\"F`o$\"1Rpjx'o%GHF`o$!16a(><@u.\"F,7%$ \"1f(fNxU8)>F`o$\"1R\\,^jyJMF`o$!1o:yR$[50\"F,7%$\"1LF2cx-wAF`o$\"1c\\ &Hr&>URF`o$!1y#))*ou%o1\"F,7%$\"1zH*\\3;ad#F`o$\"1kz$)4:vgWF`o$!1#>X!3 7&[3\"F,7%$\"17'exVE,)GF`o$\"1v-UA`_))\\F`o$!1+`;hn406F,7%$\"1]\"='*[) y!>$F`o$\"1%fC-y2m_&F`oF`r7%$\"1Yk$[@W!3NF`o$\"1$z%H:36wgF`o$!1tc7w`[_ 6F,7%$\"1)=;L@\\D$QF`o$\"1m<(*f*p\"QmF`o$!1Jk=ljsz6F,7%$\"1#Qw`1u\\;%F `o$\"1!fQMpYR@(F`o$!1z^6D_S47F,7%$\"1du[(G/QN\"F,7%$\"1$3/&>EfqfF`o$\"1#oFo(p8M5Fi^q$!1KC)p(yp'R\"F,7% $\"1!>FizOcO'F`o$\"1Y&f`jgD5\"Fi^q$!1sg2.xZU9F,7%$\"12Xtz`$Qx'F`o$\"1y 2)RqiK<\"Fi^q$!1c$z\"o$Q7\\\"F,7%$\"1_,v+>.'>(F`o$\"1q&)\\'G*QY7Fi^q$! 1W_\"[j!3V:F,-F]aq6&F_aqF)F)F`aq-%*AXESSTYLEG6#%%NONEG-%&STYLEG6#%'HID DENG-%%VIEWG6%;$!#7Fbaq$\"#7FbaqF[efl;$F0Fbaq$\"# " 0 "" {MPLTEXT 1 0 55 "alias( f=f(u,v,w), g=g(u,v,w ), h=h(u,v,w), c=[u,v,w] ):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 44 "grad( f, c, coords = prolatespheroidal(a) );" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'vectorG6#7%*&-%%diffG6$%\"fG%\"uG\"\"\"*&%\"aG\"\"\" -%%sqrtG6#,&*$)-%%sinhG6#F,\"\"#F-\"\"\"*$)-%$sinG6#%\"vGF:F-F;F-!\"\" *&-F)6$F+FAF-*&F/\"\"\"-F26#F4F-FB*&-F)6$F+%\"wGF-*(F/\"\"\"F7\"\"\"F> \"\"\"FB" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 56 "interface( verb oseproc = 2 ): # make Maple communicative" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 56 "readlib( `plots/changecoords` ): # load library fun ction" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 31 "print( `plots/chan gecoords` ); " }}{PARA 12 "" 1 "" {XPPMATH 20 "6#R6$%\"pG%&coordG6&%'c _nameG%\"aG%\"bG%\"cG6#%aoCopyright~(c)~1994~by~the~University~of~Wate rloo.~All~rights~reserved.G6\"C'>8%\"\"\">8&#F2\"\"#>8'#F2\"\"$@%-%%ty peG6$9%%)functionGC$>8$-%#opG6$\"\"!F?@)/-%%nopsG6#F?F2>F1-FE6$F2F?/FJ F6C$>F1FN>F4-FE6$F6F?/FJF:C%>F1FN>F4FT>F8-FE6$F:F?-%&ERRORG6#%DInappro priate~number~of~parameters.G>FCF?@'-%'memberG6$FC--%(readlibG6#%1plot 3d/coordset2GF.-%plots/changecoords/threetransG6'FhoFCF1F4F8-Fhn6 $%DCannot~convert~to~coordinate~systemGF?F.F.F." }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 26 "setcoords := proc( coord )" }}{PARA 0 "> " 0 " " {MPLTEXT 1 0 15 " local c_name;" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 40 " global curl, diverge, grad, laplacian," }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 41 " scalefactors, `setcoords/cname`;" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 30 " if type(coord,function) then" }}{PARA 0 " > " 0 "" {MPLTEXT 1 0 30 " c_name := op(0,coord);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 28 " if nops(coord) > 3 then" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 52 " ERROR(`Inappropriate number of parameter s.`)" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 7 " fi" }}{PARA 0 "> " 0 " " {MPLTEXT 1 0 22 " else c_name := coord" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 5 " fi;" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 52 " if membe r(c_name,readlib(`plot3d/coordset2`)()) or" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 48 " member(c_name,readlib(`plot3d/coordset`)())" }} {PARA 0 "> " 0 "" {MPLTEXT 1 0 6 " then" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 32 " `setcoords/cname` := coord," }}{PARA 0 "> " 0 " " {MPLTEXT 1 0 48 " unprotect( curl, diverge, grad, laplacian, " } }{PARA 0 "> " 0 "" {MPLTEXT 1 0 31 " scalefactors ):" } }{PARA 0 "> " 0 "" {MPLTEXT 1 0 39 " curl := proc() linalg[curl]( \+ args," }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 46 " coords = `set coords/cname` ) end;" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 45 " diverg e := proc() linalg[diverge]( args," }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 49 " coords = `setcoords/cname` ) end;" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 39 " grad := proc() linalg[grad]( args," }} {PARA 0 "> " 0 "" {MPLTEXT 1 0 46 " coords = `setcoords/cn ame` ) end;" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 49 " laplacian := pr oc() linalg[laplacian]( args," }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 51 " \+ coords = `setcoords/cname` ) end;" }}{PARA 0 "> " 0 " " {MPLTEXT 1 0 57 " scalefactors := proc() readlib(`linalg/scalefc ts`)( " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 65 " args , `setcoords/cname` ) end; " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 46 " protect( curl, diverge, grad, laplacian, " }} {PARA 0 "> " 0 "" {MPLTEXT 1 0 29 " scalefactors )" }} {PARA 0 "> " 0 "" {MPLTEXT 1 0 6 " else" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 55 " ERROR(`Cannot convert to coordinate system`,coor d)" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 4 " fi" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 4 "end:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "set coords( prolatespheroidal(a) ):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "laplacian( f, c );" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#*&,,**% \"aG\"\"\"-%%coshG6#%\"uGF'-%$sinG6#%\"vGF'-%%diffG6$%\"fGF+F'F'**F&\" \"\"-%%sinhGF*F'F,F5-F16$F3-%\"$G6$F+\"\"#F'F'**F&F5F6F5-%$cosGF.F'-F1 6$F3F/F'F'**F&F5F6F5F,F5-F16$F3-F;6$F/F=F'F'*&*(F&F5,&*$)F6F=F5F'*$)F, F=F5F'F'-F16$F3-F;6$%\"wGF=F'F5*&F6\"\"\"F,\"\"\"!\"\"F'F5**)F&\"\"$F5 FJ\"\"\"F6\"\"\"F,\"\"\"FW" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "scalefactors( c );" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#7%*&%\"aG\" \"\"-%%sqrtG6#,&*$)-%%sinhG6#%\"uG\"\"#\"\"\"F&*$)-%$sinG6#%\"vGF1F2F& F2F$*(F%F2F-F&F5F&" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 45 "dV := convert( %, `*` ); # the volume element" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#dVG**)%\"aG\"\"$\"\"\",&*$)-%%sinhG6#%\"uG\"\"#F)\"\"\"*$)-%$ sinG6#%\"vGF1F)F2F2F-F2F5F2" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 51 "readlib( addcoords ): # load the library function" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 10 "addcoords(" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 51 " myprolatespheroidal, # name of coordinate system" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 43 " [xi,eta,phi], # list of variables" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 40 " [a*sqrt((xi^2-1)*(1 -eta^2))*cos(phi), " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 40 " a*sqrt((x i^2-1)*(1-eta^2))*sin(phi), " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 55 " \+ a*xi*eta ], # the cartesian coordinates x,y,z" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 55 " # in terms of the new coor dinates" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 50 " [a], # \+ list of names of constants" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 56 " [[1 .1],[sqrt(2)/2],[Pi/3], # default values for" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 51 " [1..1.2,-1..1,0..2*Pi], # coordplot or " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 49 " [-1.2..1.2,-1.2..1.2,-1.7.. 1.7]] # coordplot3d" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 2 "):" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 57 "assume( a>0, xi>=1, -1<=eta, eta<=1, 0<=phi, phi<=2*Pi );" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 45 "alias( f = f(xi,eta,phi), c = [xi,eta,phi] ):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 36 "setcoords( myprolatespheroidal(a) ):" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "grad( f, c );" }}{PARA 12 " " 1 "" {XPPMATH 20 "6#-%'vectorG6#7%*&-%%diffG6$%\"fG%$xi|irG\"\"\"*$- %%sqrtG6#,(*&)%#a|irG\"\"#F-)-%$absG6#*&*(F,\"\"\"-%$cosG6#%%phi|irGF= ,&!\"\"F=*$)%%eta|irGF6F-F=F-F-*$-F06#,$*&,&*$)F,F6F-F=FCF=F=FBF=FCF-! \"\"F6F-F=*&F4F-)-F96#*&*(-%$sinGF@F=F,F-FBF-F-*$-F06#FJF-FOF6F-F=*&F4 F-)-F96#FFF6F-F=F-FO*&-F)6$F+FFF-*$-F06#,(*&F4F-)-F96#*&*(F>F-FLF-FFF= F-*$-F06#FJF-FOF6F-F=*&F4F-)-F96#*&*(FVF-FLF-FFF-F-*$-F06#FJF-FOF6F-F= *&F4F-)-F96#F,F6F-F=F-FO*&-F)6$F+FAF-*$-F06#,&*&F4F-)-F96#*&-F06#FJF-F VF-F6F-F=*&F4F-)-F96#*&FbqF-F>F-F6F-F=F-FO" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 34 "map( simplify@combine@factor, % );" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'vectorG6#7%*&*&-%%diffG6$%\"fG%$xi|irG\"\"\"-%% sqrtG6#,&*$)F-\"\"#\"\"\"F.!\"\"F.F6F6*&%#a|irG\"\"\"-F06#,&*$)%%eta|i rGF5F6F7F3F.F6!\"\",$*&*&-F*6$F,F@F.,&F7F.F>F.F.F6*(F9\"\"\"-F06#F=F6- F06#,&F.F.F>F7F6FAF7*&-F*6$F,%%phi|irGF6*(F9\"\"\"-F06#F2F6-F06#FNF6FA " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "scalefactors( c ):" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "map( simplify@combine@factor , % );" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#7%*&*&%#a|irG\"\"\"-%%sqrtG6 #,&*$)%%eta|irG\"\"#\"\"\"!\"\"*$)%$xi|irGF/F0F'F0F0*$-F)6#,&F2F'F1F'F 0!\"\",$*&*(F&F0F(F0-F)6#,&F'F'F,F1F0F0,&F1F'F,F'F9F1*&F&F0-F)6#*&F8F' F?F'F0" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 36 "dV := simplify( c onvert( %, `*` ) );" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#dVG,$*&,&*$) %%eta|irG\"\"#\"\"\"\"\"\"*$)%$xi|irGF+F,!\"\"F-)%#a|irG\"\"$F,F1" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "restart; with(linalg):" }} {PARA 7 "" 1 "" {TEXT -1 32 "Warning, new definition for norm" }} {PARA 7 "" 1 "" {TEXT -1 33 "Warning, new definition for trace" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 57 "assume( a>0, xi>=1, -1<=eta, eta<=1, 0<=phi, phi<=2*Pi );" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 45 "alias( f = f(xi,eta,phi), c = [xi,eta,phi] ):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 43 "T := [a*sqrt((xi^2-1)*(1-eta^2))*cos(phi) , " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 56 " a*sqrt((xi^2-1)*(1-eta^ 2))*sin(phi), a*xi*eta ]; " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"TG7 %*(%#a|irG\"\"\"-%%sqrtG6#*&,&*$)%$xi|irG\"\"#\"\"\"F(!\"\"F(F(,&F(F(* $)%%eta|irGF1F2F3F(F2-%$cosG6#%%phi|irGF(*(F'F2F)F2-%$sinGF:F(*(F'F2F0 F(F7F(" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 32 "J := jacobian(T,[ xi,eta,phi]); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"JG-%'matrixG6#7 %7%*&**%#a|irG\"\"\"-%$cosG6#%%phi|irGF-%$xi|irGF-,&F-F-*$)%%eta|irG\" \"#\"\"\"!\"\"F8F8*$-%%sqrtG6#*&,&*$)F2F7F8F-F9F-F-F3F-F8!\"\",$*&**F, F8F.F8F?F8F6F-F8*$-F<6#F>F8FBF9,$*(F,F8-F<6#F>F8-%$sinGF0F-F97%*&**F,F 8FMF8F2F8F3F8F8*$-F<6#F>F8FB,$*&**F,F8FMF8F?F8F6F8F8*$-F<6#F>F8FBF9*(F ,F8FKF8F.F87%*&F,F8F6F8*&F,F8F2F8\"\"!" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 55 "g := evalm( transpose(J) &* J ): g := map(simplify,g );" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"gG-%'matrixG6#7%7%*&*&)%#a|i rG\"\"#\"\"\",&*$)%%eta|irGF.F/!\"\"*$)%$xi|irGF.F/\"\"\"F8F/,&F5F8F4F 8!\"\"\"\"!F;7%F;,$*&*&F,F/F0F/F/,&F4F8F1F8F:F4F;7%F;F;,**&F,F/F6F/F8* (F,F/F6F/F2F/F4*$F,F/F4*&F,F/F2F/F8" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "h1 := sqrt( g[1,1] );" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#h1G*&%#a|irG\"\"\"-%%sqrtG6#*&,&*$)%%eta|irG\"\"#\"\"\"!\"\"* $)%$xi|irGF0F1F'F1,&F3F'F2F'!\"\"F1" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "h2 := sqrt( g[2,2] );" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#h2G*&%#a|irG\"\"\"-%%sqrtG6#*&,&*$)%%eta|irG\"\"#\"\"\"!\"\"* $)%$xi|irGF0F1F'F1,&F'F'F-F2!\"\"F1" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 29 "h3 := sqrt( factor(g[3,3]) );" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#h3G*&%#a|irG\"\"\"-%%sqrtG6#**,&F'F'%%eta|irG!\"\"F' ,&F-F'F'F'F',&%$xi|irGF'F.F'F',&F1F'F'F'F'\"\"\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 42 "sftable := readlib( `linalg/scaletable` ):" } }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 37 "sftable[myprolatespheroida l] := subs(" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 47 " \{ a=_a, xi=_x, eta =_y, phi=_z \}, [h1,h2,h3] );" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>&%(s ftableG6#%4myprolatespheroidalG7%*&%#_aG\"\"\"-%%sqrtG6#*&,&*$)%#_xG\" \"#\"\"\"F+*$)%#_yGF4F5!\"\"F5,&F1F+F9F+!\"\"F5*&F*F5-F-6#*&F0F5,&F+F+ F6F9F;F5*&F*F5-F-6#**,&F+F+F8F9F+,&F8F+F+F+F+,&F3F+F9F+F+,&F3F+F+F+F+F 5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 41 "rtable := readlib(`plo t3d/coordSystems`):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 71 "rtab le[myprolatespheroidal] := subs( \{a=_a, eta=_x, psi=_y, z=_z\}, T ); " }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>&%'rtableG6#%4myprolatespheroidal G7%*(%#_aG\"\"\"-%%sqrtG6#*&,&*$)%$xi|irG\"\"#\"\"\"F+!\"\"F+F+,&F+F+* $)%#_xGF4F5F6F+F5-%$cosG6#%%phi|irGF+*(F*F5F,F5-%$sinGF=F+*(F*F5F3F+F: F+" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 49 "laplacian( f, c, coor ds=myprolatespheroidal(a) ):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "map( combine, % ):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 54 "collect( %, [ diff(f,xi$2), diff(f,eta$2), diff(f,xi)," }}{PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 43 " diff(f,eta), f ], distributed, simplify):" } }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "map(combine, %):" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "map(normal, %, expanded);" } }{PARA 12 "" 1 "" {XPPMATH 20 "6#,,*&,&*&-%%diffG6$%\"fG-%\"$G6$%$xi|i rG\"\"#\"\"\")F.F/\"\"\"F0F'!\"\"F2,&*&)%#a|irGF/F2)%%eta|irGF/F2F3*&F 6F2F1F2F0!\"\"F0*&,&-F(6$F*-F,6$F9F/F0*&F>F0F8F2F3F2F4F;F0*&*&F.F0-F(6 $F*F.F0F2F4F;F/*&*&F9F0-F(6$F*F9F0F2F4F;!\"#*&-F(6$F*-F,6$%%phi|irGF/F 2,**(F6F2F1F2F8F2F0F5F3F:F3*$F6F2F0F;F3" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 54 "collect( %, [ diff(f,xi$2), diff(f,eta$2), diff(f,xi) ," }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 43 " diff(f,eta), f ], distribute d, simplify);" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#,,*&*&%$xi|irG\"\"\"- %%diffG6$%\"fGF&F'\"\"\"*&)%#a|irG\"\"#F,,&*$)%%eta|irG\"\"#F,!\"\"*$) F&F5F,F'\"\"\"!\"\"F5*&*&,&F7F'F6F'F'-F)6$F+-%\"$G6$F&F5F'F,*&)F/\"\"# F,F1\"\"\"F:F'*&*&,&F6F'F2F'F'-F)6$F+-FA6$F4F5F'F,*&)F/\"\"#F,F1\"\"\" F:F6*&*&F4F'-F)6$F+F4F'F,*&)F/\"\"#F,F1\"\"\"F:!\"#*&-F)6$F+-FA6$%%phi |irGF5F,*&)F/\"\"#F,,**&F8F,F3F,F'F2F6F7F6F'F'\"\"\"F:F6" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}{SECT 1 {PARA 3 "" 0 "MoorePen rose" {TEXT -1 26 "19.5 Moore-Penrose Inverse" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "restart; with(linalg):" }}{PARA 7 "" 1 "" {TEXT -1 32 "Warning, new definition for norm" }}{PARA 7 "" 1 "" {TEXT -1 33 "Warning, new definition for trace" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 33 "MPinv := A -> map( limit, evalm( " }}{PARA 0 "> " 0 " " {MPLTEXT 1 0 46 " evalm( (linalg[transpose](A)&*A+x^2)^(-1) ) " }} {PARA 0 "> " 0 "" {MPLTEXT 1 0 37 " &* linalg[transpose](A) ), x=0 ); " }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%&MPinvGR6#%\"AG6\"6$%)operator G%&arrowGF(-%$mapG6%%&limitG-%&evalmG6#-%#&*G6$-F16#*&\"\"\"F9,&-F46$- &%'linalgG6#%*transposeG6#9$FC\"\"\"*$)%\"xG\"\"#F9FD!\"\"F=/FG\"\"!F( F(F(" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "M := randmatrix(3,2 ); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"MG-%'matrixG6#7%7$!#&)!#b7 $!#P!#N7$\"#(*\"#]" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "MPinv (M); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7$7%#\"$F%\"&d*) )#\"%zDF*#\"%vAF*7%#!&$49\"'q&*))#!&`f%F2#!&R\\\"\"'&yW%" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 39 "poly := proc() Randpoly(2,y) mod 2 \+ end:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 36 "M := randmatrix(3,2 ,entries=poly); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"MG-%'matrixG6 #7%7$,&*$)%\"yG\"\"#\"\"\"\"\"\"F-F0,&F+F0F0F07$F+F+7$F1F+" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "MPinv(M); " }}{PARA 259 "" 0 "" {TEXT -1 0 "" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7$7%*&*&)% \"yG\"\"#\"\"\",(*$)F+\"\"$F-F,*$F*F-!\"$!\"\"\"\"\"F5F-,.*$)F+\"\"'F- F,*$)F+\"\"%F-F9F2F<*$)F+\"\"&F-!\"'F/!\"#F5F5!\"\",$*&*&F*F-,(F4F5F/F 5F+F5F5F-F6FBF4,$*&,,F:F3F2F3F4F5F=F5F/F5F-F6FBF47%,$*&,,F=F,F:F3F/F5F 2F3F4F5F-F6FBF4*&*&F*F-,*F2F5F/F5F+F4F5F5F5F-F6FB*&*&F+F5,,F/FAF+F4F:F 5F2FAF4F5F5F-F6FB" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}} {MARK "0" 0 }{VIEWOPTS 1 1 0 1 1 1803 }