{VERSION 4 0 "IBM INTEL NT" "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 Output" 2 20 "" 0 1 0 0 255 1 0 0 0 0 0 0 0 0 0 1 } {CSTYLE "" -1 256 "" 1 14 0 0 1 1 0 0 0 0 0 0 0 0 0 0 }{CSTYLE "Maple \+ Input" -1 257 "Courier" 0 0 255 0 0 1 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 258 "" 0 1 0 0 255 1 0 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 259 "" 1 36 0 0 0 0 0 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 260 "" 1 18 235 121 7 1 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 261 "" 1 14 0 0 1 1 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 262 "" 1 12 0 0 255 1 0 0 0 0 0 0 0 0 0 0 } {CSTYLE "" -1 263 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 264 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 265 "" 1 12 0 0 255 1 0 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 266 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 267 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 } {CSTYLE "" -1 268 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 269 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 270 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 271 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 272 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 } {CSTYLE "" -1 273 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 274 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 275 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 276 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 277 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 } {CSTYLE "" -1 278 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 279 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 280 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 281 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 282 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 } {CSTYLE "" -1 283 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 284 "" 1 12 0 0 255 1 0 0 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 "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 "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 "" 11 12 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 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 "Author" -1 19 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }3 1 0 0 8 8 1 0 1 0 2 2 0 1 }{PSTYLE " Heading 4" -1 20 1 {CSTYLE "" -1 -1 "Times" 1 10 0 0 0 1 1 1 2 2 2 2 1 1 1 1 }1 1 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Normal" -1 256 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }3 1 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Title" -1 257 1 {CSTYLE "" -1 -1 "Times" 1 24 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 "Title" -1 258 1 {CSTYLE "" -1 -1 "Times" 1 18 0 0 0 1 2 1 2 2 2 2 1 1 1 1 }3 1 0 0 12 12 1 0 1 0 2 2 19 1 }{PSTYLE "Normal" -1 259 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 "Heading 2" -1 260 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 2 1 0 1 0 2 2 0 1 }{PSTYLE "Heading 1" -1 261 1 {CSTYLE "" -1 -1 "Times" 1 24 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 "Normal" -1 262 1 {CSTYLE "" -1 -1 "Times" 1 24 0 0 0 1 2 1 1 2 2 2 1 1 1 1 }3 1 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Normal" -1 263 1 {CSTYLE "" -1 -1 "Times" 1 24 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }3 1 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Normal" -1 264 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 1 2 2 2 1 1 1 1 }1 1 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Normal" -1 265 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 1 2 2 2 2 1 1 1 1 }1 1 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Bullet Item" -1 266 1 {CSTYLE "" -1 -1 " Times" 1 12 0 0 0 1 2 1 2 2 2 2 1 1 1 1 }1 1 0 0 3 3 1 0 1 0 2 2 15 2 }{PSTYLE "Heading 3" -1 267 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 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Maple Output " -1 268 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 } 1 3 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Heading 1" -1 269 1 {CSTYLE "" -1 -1 "Times" 1 18 0 0 0 1 2 1 2 2 2 2 1 1 1 1 }3 1 0 0 8 4 1 0 1 0 2 2 0 1 }{PSTYLE "Heading 2" -1 270 1 {CSTYLE "" -1 -1 "Times" 1 18 0 0 0 1 2 1 2 2 2 2 1 1 1 1 }3 1 0 0 8 2 1 0 1 0 2 2 0 1 }{PSTYLE "Heading 3" -1 271 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 1 1 2 2 2 2 1 1 1 1 }3 1 0 0 0 0 1 0 1 0 2 2 0 1 }} {SECT 0 {EXCHG {PARA 269 "" 0 "" {TEXT 259 15 "FINITE ELEMENTS" }} {PARA 270 "" 0 "" {TEXT -1 29 "SYMBOLIC PROGRAMMING IN MAPLE" }}{PARA 19 "" 0 "" {TEXT 260 16 "by Artur Portela" }{TEXT 261 0 "" }}{PARA 271 "" 0 "" {TEXT 256 26 "Potential-Flow Example 4: " }{TEXT 284 19 "S ymbolic Variables\n" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 47 "rest art:interface(verboseproc=3):printlevel:=3:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "libname := \"C:/mylib/fem\",libname:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "with(Plotter);" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#7:%+draw_labelG%,draw_symbolG%*draw_textG%+draw_titleG% )move_penG%*new_arrowG%)new_axesG%*new_colorG%.new_pen_styleG%1new_plo t_shadingG%/new_plot_styleG%.new_thicknessG%)pen_downG%'pen_upG%-plot_ elementG%2show_2D_animationG%-show_2D_plotG%2show_3D_animationG%-show_ 3D_plotG%3start_2D_animationG%.start_2D_plotG%3start_3D_animationG%.st art_3D_plotG%/upgrade_framesG" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "with(Cgt_fem);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#7&%(cgt_femG%/ compute_fluxesG%1element_geometryG%/read_save_dataG" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "with(G_cgt_fem);" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#7+%/animate_3D_defG%/animate_3D_rotG%3animate_velocitie sG%3plot_contour_linesG%,plot_fluxesG%5plot_line_potentialsG%*plot_mes hG%2plot_problem_dataG%0plot_velocitiesG" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{SECT 0 {PARA 20 "" 0 "" {TEXT 262 17 " Data Prep aration" }}{EXCHG {PARA 0 "" 0 "" {TEXT -1 31 "The best way to input d ata for " }{TEXT 264 7 "Cgt_fem" }{TEXT -1 25 " is to use the procedur e " }{TEXT 263 14 "read_save_data" }{TEXT -1 50 ", which reads a file \+ with the following structure:" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }} {PARA 0 "" 0 "" {TEXT -1 1 "*" }{TEXT 266 7 "control" }{TEXT -1 169 "* [title,point sources,element sources,boundary velocities] \ntitle Potential Flow in a Rectangular Domain\npoint sources\nelement sources\nboundary velocities\n\n*" }{TEXT 267 9 "materials" }{TEXT -1 104 "* [material,x-permeability,y-permeability,angle of t he local x-direction]\n1 2 2 45\n2 8 8 45\n\n*" }{TEXT 268 5 "nodes" } {TEXT -1 180 "* [node,x,y]\n1 0 0\n2 0 1 \n3 0 \+ 2 \n4 1 0 \n5 1 1 \n6 1 2\n7 2 0 \n8 2 1\n9 2 2\n10 \+ 3 0 \n11 3 1 \n12 3 2\n13 4 0 \n14 4 1\n15 4 2 \n\n *" }{TEXT 269 8 "elements" }{TEXT -1 399 "* [element,node1, node2,node3,material]\n1 1 4 5 1 \n2 1 5 2 1 \+ \n3 2 5 3 1 \n4 3 5 6 1 \n5 4 7 8 1 \n 6 5 4 8 1\n7 5 8 6 1 \n8 6 8 9 1\n9 7 10 \+ 11 1 \n10 7 11 8 1 \n11 8 11 9 1\n12 11 12 \+ 9 1 \n13 10 13 14 1\n14 10 14 11 1 \n15 11 14 12 1\n16 12 14 15 1 \n\n*" }{TEXT 270 11 "constraints" } {TEXT -1 59 "* [node,potential]\n1 10\n2 10\n3 10\n13 0\n14 0 \n15 0\n\n*" }{TEXT 271 19 "boundary velocities" }{TEXT -1 52 "* [node 1,node2,normal velocity]\n1 2 0\n2 3 0\n5 4 0\n\n*" }{TEXT 272 13 "poi nt sources" }{TEXT -1 25 "* [node,Q]\n1 -50\n\n*" }{TEXT 273 15 "element sources" }{TEXT -1 25 "* [element,q]\n2 30\n\n*" }{TEXT 274 3 "end" }{TEXT -1 1 "*" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 85 "The data blocks, with the respective keyword on t he top, can be defined in any order." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 83 "Alternatively, data can be given manual ly through the definition of the variables: " }{TEXT 275 5 "tcase" } {TEXT -1 2 ", " }{TEXT 276 7 "control" }{TEXT -1 2 ", " }{TEXT 277 4 " nods" }{TEXT -1 2 ", " }{TEXT 278 5 "elems" }{TEXT -1 2 ", " }{TEXT 279 9 "mat_props" }{TEXT -1 2 ", " }{TEXT 280 9 "bdr_conds" }{TEXT -1 2 ", " }{TEXT 281 7 "b_velts" }{TEXT -1 2 ", " }{TEXT 282 6 "p_srcs" } {TEXT -1 5 " and " }{TEXT 283 6 "e_srcs" }{TEXT -1 46 ". See bellow th e structure of these variables." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}{SECT 0 {PARA 20 "" 0 "" {TEXT 265 19 " Symbolic Variable s" }{TEXT -1 0 "" }}{EXCHG {PARA 0 "" 0 "" {TEXT -1 63 "For the sake o f simplicity, start reading the data of example 1" }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 17 "read_save_data();" }}}{EXCHG {PARA 0 "read d ata from a file (y/n) ? " 0 "" {MPLTEXT 1 0 2 "y;" }}}{EXCHG {PARA 0 " file name: " 0 "" {MPLTEXT 1 0 16 "\"dat_test1.txt\";" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#Q**control*6\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6# Q,*materials*6\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#Q(*nodes*6\"" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#Q+*elements*6\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#Q.*constraints*6\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#Q 6*boundary~velocities*6\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#Q0*point~ sources*6\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#Q2*element~sources*6\" " }}}{EXCHG {PARA 0 "save data into a file (y/n) ? " 0 "" {MPLTEXT 1 0 2 "n;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 67 "tcase;control;no ds;elems;mat_props;bdr_conds;b_velts;p_srcs;e_srcs;" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#QH~Potential~Flow~in~a~Rectangular~Domain6\"" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#7%Q.point~sources6\"Q0element~sourcesF %Q4boundary~velocitiesF%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#717$\"\"!F %7$F%\"\"\"7$F%\"\"#7$F'F%7$F'F'7$F'F)7$F)F%7$F)F'7$F)F)7$\"\"$F%7$F1F '7$F1F)7$\"\"%F%7$F5F'7$F5F)" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#727&\" \"\"\"\"%\"\"&F%7&F%F'\"\"#F%7&F)F'\"\"$F%7&F+F'\"\"'F%7&F&\"\"(\"\")F %7&F'F&F0F%7&F'F0F-F%7&F-F0\"\"*F%7&F/\"#5\"#6F%7&F/F7F0F%7&F0F7F4F%7& F7\"#7F4F%7&F6\"#8\"#9F%7&F6F>F7F%7&F7F>F;F%7&F;F>\"#:F%" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#7$7%\"\"#F%\"#X7%\"\")F(F&" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#7(7$\"\"\"\"#57$\"\"#F&7$\"\"$F&7$\"#8\"\"!7$\"#9F-7$\" #:F-" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#7%7%\"\"\"\"\"#\"\"!7%F&\"\"$F '7%\"\"&\"\"%F'" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#7#7$\"\"\"!#]" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#7#7$\"\"#\"#I" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 44 "Redefine this data with symbolic variables :" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 214 "mat_props:=[[kx,ky,alpha]]: elems:=[[1,4,5,1],[1,5,2,1],[2,5,3,1],[3,5,6,1],[4,7,8,1],[5,4,8,1],[5 ,8,6,1],[6,8,9,1],[7,10,11,1],[7,11,8,1],[8,11,9,1],[11,12,9,1],[10,13 ,14,1],[10,14,11,1],[11,14,12,1],[12,14,15,1]]:" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 67 "tcase;control;nods;elems;mat_props;bdr_conds;b _velts;p_srcs;e_srcs;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#QH~Potential~ Flow~in~a~Rectangular~Domain6\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#7%Q .point~sources6\"Q0element~sourcesF%Q4boundary~velocitiesF%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#717$\"\"!F%7$F%\"\"\"7$F%\"\"#7$F'F%7$F'F'7$ F'F)7$F)F%7$F)F'7$F)F)7$\"\"$F%7$F1F'7$F1F)7$\"\"%F%7$F5F'7$F5F)" }} {PARA 12 "" 1 "" {XPPMATH 20 "6#727&\"\"\"\"\"%\"\"&F%7&F%F'\"\"#F%7&F )F'\"\"$F%7&F+F'\"\"'F%7&F&\"\"(\"\")F%7&F'F&F0F%7&F'F0F-F%7&F-F0\"\"* F%7&F/\"#5\"#6F%7&F/F7F0F%7&F0F7F4F%7&F7\"#7F4F%7&F6\"#8\"#9F%7&F6F>F7 F%7&F7F>F;F%7&F;F>\"#:F%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#7#7%%#kxG% #kyG%&alphaG" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#7(7$\"\"\"\"#57$\"\"#F &7$\"\"$F&7$\"#8\"\"!7$\"#9F-7$\"#:F-" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#7%7%\"\"\"\"\"#\"\"!7%F&\"\"$F'7%\"\"&\"\"%F'" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#7#7$\"\"\"!#]" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#7#7$\" \"#\"#I" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "cgt_fem(yes);" } }{PARA 11 "" 1 "" {XPPMATH 20 "6#%>Initializing~global~matrices:G" }} {PARA 11 "" 1 "" {XPPMATH 20 "6$-%'matrixG6#7171\"\"!F(F(F(F(F(F(F(F(F (F(F(F(F(F(F'F'F'F'F'F'F'F'F'F'F'F'F'F'-%'vectorG6#F'" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#*(%3Assembling~elementG\"\"\"7&F%\"\"%\"\"&F%F%%\" :GF%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7171,$%#kxG#!\"\" \"\"#\"\"!F-,$F)#\"\"\"F,F-F-F-F-F-F-F-F-F-F-F-71F-F-F-F-F-F-F-F-F-F-F -F-F-F-F-F171F.F-F-,&F)F**&#F0F,F0%#kyGF0F+,$F6F/F-F-F-F-F-F-F-F-F-F-7 1F-F-F-F7,$F6F*F-F-F-F-F-F-F-F-F-F-F1F1F1F1F1F1F1F1F1F1" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#*(%3Assembling~elementG\"\"\"7&F%\"\"&\"\"#F%F%% \":GF%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7171,&%#kxG#!\" \"\"\"#*&#\"\"\"F,F/%#kyGF/F+,$F0#F/F,\"\"!,$F)F2F3F3F3F3F3F3F3F3F3F3F 371F1F(F3F3F4F3F3F3F3F3F3F3F3F3F371F3F3F3F3F3F3F3F3F3F3F3F3F3F3F371F4F 3F3F(F1F3F3F3F3F3F3F3F3F3F371F3F4F3F1F(F3F3F3F3F3F3F3F3F3F3F6F6F6F6F6F 6F6F6F6F6" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#*(%3Assembling~elementG\" \"\"7&\"\"#\"\"&\"\"$F%F%%\":GF%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-% 'matrixG6#7171,&%#kxG#!\"\"\"\"#*&#\"\"\"F,F/%#kyGF/F+,$F0#F/F,\"\"!,$ F)F2F3F3F3F3F3F3F3F3F3F3F371F1,&F)F+F0F+F1F3F)F3F3F3F3F3F3F3F3F3F371F3 F1,$F0F*F3F3F3F3F3F3F3F3F3F3F3F371F4F3F3F(F1F3F3F3F3F3F3F3F3F3F371F3F) F3F1,&F)F+*&#F/F,F/F0F/F+F3F3F3F3F3F3F3F3F3F371F3F3F3F3F3F3F3F3F3F3F3F 3F3F3F3F>F>F>F>F>F>F>F>F>" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#*(%3Assem bling~elementG\"\"\"7&\"\"$\"\"&\"\"'F%F%%\":GF%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7171,&%#kxG#!\"\"\"\"#*&#\"\"\"F,F/%#kyGF/F +,$F0#F/F,\"\"!,$F)F2F3F3F3F3F3F3F3F3F3F3F371F1,&F)F+F0F+F1F3F)F3F3F3F 3F3F3F3F3F3F371F3F1F(F3F3F4F3F3F3F3F3F3F3F3F371F4F3F3F(F1F3F3F3F3F3F3F 3F3F3F371F3F)F3F1F6F1F3F3F3F3F3F3F3F3F371F3F3F4F3F1F(F3F3F3F3F3F3F3F3F 371F3F3F3F3F3F3F3F3F3F3F3F3F3F3F3F;F;F;F;F;F;F;F;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#*(%3Assembling~elementG\"\"\"7&\"\"%\"\"(\"\")F%F%%\":G F%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7171,&%#kxG#!\"\"\" \"#*&#\"\"\"F,F/%#kyGF/F+,$F0#F/F,\"\"!,$F)F2F3F3F3F3F3F3F3F3F3F3F371F 1,&F)F+F0F+F1F3F)F3F3F3F3F3F3F3F3F3F371F3F1F(F3F3F4F3F3F3F3F3F3F3F3F37 1F4F3F3,&F)F+*&#F/F,F/F0F/F+F1F3F4F3F3F3F3F3F3F3F371F3F)F3F1F6F1F3F3F3 F3F3F3F3F3F371F3F3F4F3F1F(F3F3F3F3F3F3F3F3F371F3F3F3F4F3F3F(F1F3F3F3F3 F3F3F371F3F3F3F3F3F3F1,$F0F*F3F3F3F3F3F3F371F3F3F3F3F3F3F3F3F3F3F3F3F3 F3F3FAFAFAFAFAFA" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#*(%3Assembling~ele mentG\"\"\"7&\"\"&\"\"%\"\")F%F%%\":GF%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7171,&%#kxG#!\"\"\"\"#*&#\"\"\"F,F/%#kyGF/F+,$F0#F/F ,\"\"!,$F)F2F3F3F3F3F3F3F3F3F3F3F371F1,&F)F+F0F+F1F3F)F3F3F3F3F3F3F3F3 F3F371F3F1F(F3F3F4F3F3F3F3F3F3F3F3F371F4F3F3F6F0F3F4F3F3F3F3F3F3F3F371 F3F)F3F0,&F)#!\"$F,*&#\"\"$F,F/F0F/F+F1F3F4F3F3F3F3F3F3F371F3F3F4F3F1F (F3F3F3F3F3F3F3F3F371F3F3F3F4F3F3F(F1F3F3F3F3F3F3F371F3F3F3F3F4F3F1F(F 3F3F3F3F3F3F371F3F3F3F3F3F3F3F3F3F3F3F3F3F3F3FCFCFCFCFCFC" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#*(%3Assembling~elementG\"\"\"7&\"\"&\"\")\"\"'F %F%%\":GF%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7171,&%#kxG# !\"\"\"\"#*&#\"\"\"F,F/%#kyGF/F+,$F0#F/F,\"\"!,$F)F2F3F3F3F3F3F3F3F3F3 F3F371F1,&F)F+F0F+F1F3F)F3F3F3F3F3F3F3F3F3F371F3F1F(F3F3F4F3F3F3F3F3F3 F3F3F371F4F3F3F6F0F3F4F3F3F3F3F3F3F3F371F3F)F3F0,&F)!\"#*&F,F/F0F/F+F0 F3F)F3F3F3F3F3F3F371F3F3F4F3F0,&F)F*F0F+F3F3F3F3F3F3F3F3F371F3F3F3F4F3 F3F(F1F3F3F3F3F3F3F371F3F3F3F3F)F3F1,&F)F+*&#F/F,F/F0F/F+F3F3F3F3F3F3F 371F3F3F3F3F3F3F3F3F3F3F3F3F3F3F3FDFDFDFDFDFD" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#*(%3Assembling~elementG\"\"\"7&\"\"'\"\")\"\"*F%F%%\":G F%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7171,&%#kxG#!\"\"\" \"#*&#\"\"\"F,F/%#kyGF/F+,$F0#F/F,\"\"!,$F)F2F3F3F3F3F3F3F3F3F3F3F371F 1,&F)F+F0F+F1F3F)F3F3F3F3F3F3F3F3F3F371F3F1F(F3F3F4F3F3F3F3F3F3F3F3F37 1F4F3F3F6F0F3F4F3F3F3F3F3F3F3F371F3F)F3F0,&F)!\"#*&F,F/F0F/F+F0F3F)F3F 3F3F3F3F3F371F3F3F4F3F0F6F3F3F4F3F3F3F3F3F371F3F3F3F4F3F3F(F1F3F3F3F3F 3F3F371F3F3F3F3F)F3F1F6F1F3F3F3F3F3F371F3F3F3F3F3F4F3F1F(F3F3F3F3F3F37 1F3F3F3F3F3F3F3F3F3F3F3F3F3F3F3FAFAFAFAFA" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#*(%3Assembling~elementG\"\"\"7&\"\"(\"#5\"#6F%F%%\":GF% " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7171,&%#kxG#!\"\"\"\"# *&#\"\"\"F,F/%#kyGF/F+,$F0#F/F,\"\"!,$F)F2F3F3F3F3F3F3F3F3F3F3F371F1,& F)F+F0F+F1F3F)F3F3F3F3F3F3F3F3F3F371F3F1F(F3F3F4F3F3F3F3F3F3F3F3F371F4 F3F3F6F0F3F4F3F3F3F3F3F3F3F371F3F)F3F0,&F)!\"#*&F,F/F0F/F+F0F3F)F3F3F3 F3F3F3F371F3F3F4F3F0F6F3F3F4F3F3F3F3F3F371F3F3F3F4F3F3,&F)F+*&#F/F,F/F 0F/F+F1F3F4F3F3F3F3F371F3F3F3F3F)F3F1F6F1F3F3F3F3F3F371F3F3F3F3F3F4F3F 1F(F3F3F3F3F3F371F3F3F3F3F3F3F4F3F3F(F1F3F3F3F371F3F3F3F3F3F3F3F3F3F1, $F0F*F3F3F3F371F3F3F3F3F3F3F3F3F3F3F3F3F3F3F3FGFGFG" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#*(%3Assembling~elementG\"\"\"7&\"\"(\"#6\"\")F%F%%\": GF%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7171,&%#kxG#!\"\"\" \"#*&#\"\"\"F,F/%#kyGF/F+,$F0#F/F,\"\"!,$F)F2F3F3F3F3F3F3F3F3F3F3F371F 1,&F)F+F0F+F1F3F)F3F3F3F3F3F3F3F3F3F371F3F1F(F3F3F4F3F3F3F3F3F3F3F3F37 1F4F3F3F6F0F3F4F3F3F3F3F3F3F3F371F3F)F3F0,&F)!\"#*&F,F/F0F/F+F0F3F)F3F 3F3F3F3F3F371F3F3F4F3F0F6F3F3F4F3F3F3F3F3F371F3F3F3F4F3F3F6F0F3F4F3F3F 3F3F371F3F3F3F3F)F3F0,&F)#!\"$F,*&#\"\"$F,F/F0F/F+F1F3F4F3F3F3F371F3F3 F3F3F3F4F3F1F(F3F3F3F3F3F371F3F3F3F3F3F3F4F3F3F(F1F3F3F3F371F3F3F3F3F3 F3F3F4F3F1F(F3F3F3F371F3F3F3F3F3F3F3F3F3F3F3F3F3F3F3FIFIFI" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#*(%3Assembling~elementG\"\"\"7&\"\")\"#6\"\" *F%F%%\":GF%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7171,&%#kx G#!\"\"\"\"#*&#\"\"\"F,F/%#kyGF/F+,$F0#F/F,\"\"!,$F)F2F3F3F3F3F3F3F3F3 F3F3F371F1,&F)F+F0F+F1F3F)F3F3F3F3F3F3F3F3F3F371F3F1F(F3F3F4F3F3F3F3F3 F3F3F3F371F4F3F3F6F0F3F4F3F3F3F3F3F3F3F371F3F)F3F0,&F)!\"#*&F,F/F0F/F+ F0F3F)F3F3F3F3F3F3F371F3F3F4F3F0F6F3F3F4F3F3F3F3F3F371F3F3F3F4F3F3F6F0 F3F4F3F3F3F3F371F3F3F3F3F)F3F0F:F0F3F)F3F3F3F371F3F3F3F3F3F4F3F0,&F)F* F0F+F3F3F3F3F3F371F3F3F3F3F3F3F4F3F3F(F1F3F3F3F371F3F3F3F3F3F3F3F)F3F1 ,&F)F+*&#F/F,F/F0F/F+F3F3F3F371F3F3F3F3F3F3F3F3F3F3F3F3F3F3F3FGFGFG" } }{PARA 11 "" 1 "" {XPPMATH 20 "6#*(%3Assembling~elementG\"\"\"7&\"#6\" #7\"\"*F%F%%\":GF%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7171 ,&%#kxG#!\"\"\"\"#*&#\"\"\"F,F/%#kyGF/F+,$F0#F/F,\"\"!,$F)F2F3F3F3F3F3 F3F3F3F3F3F371F1,&F)F+F0F+F1F3F)F3F3F3F3F3F3F3F3F3F371F3F1F(F3F3F4F3F3 F3F3F3F3F3F3F371F4F3F3F6F0F3F4F3F3F3F3F3F3F3F371F3F)F3F0,&F)!\"#*&F,F/ F0F/F+F0F3F)F3F3F3F3F3F3F371F3F3F4F3F0F6F3F3F4F3F3F3F3F3F371F3F3F3F4F3 F3F6F0F3F4F3F3F3F3F371F3F3F3F3F)F3F0F:F0F3F)F3F3F3F371F3F3F3F3F3F4F3F0 F6F3F3F4F3F3F371F3F3F3F3F3F3F4F3F3F(F1F3F3F3F371F3F3F3F3F3F3F3F)F3F1F6 F1F3F3F371F3F3F3F3F3F3F3F3F4F3F1F(F3F3F371F3F3F3F3F3F3F3F3F3F3F3F3F3F3 F3FDFD" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#*(%3Assembling~elementG\"\" \"7&\"#5\"#8\"#9F%F%%\":GF%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matr ixG6#7171,&%#kxG#!\"\"\"\"#*&#\"\"\"F,F/%#kyGF/F+,$F0#F/F,\"\"!,$F)F2F 3F3F3F3F3F3F3F3F3F3F371F1,&F)F+F0F+F1F3F)F3F3F3F3F3F3F3F3F3F371F3F1F(F 3F3F4F3F3F3F3F3F3F3F3F371F4F3F3F6F0F3F4F3F3F3F3F3F3F3F371F3F)F3F0,&F)! \"#*&F,F/F0F/F+F0F3F)F3F3F3F3F3F3F371F3F3F4F3F0F6F3F3F4F3F3F3F3F3F371F 3F3F3F4F3F3F6F0F3F4F3F3F3F3F371F3F3F3F3F)F3F0F:F0F3F)F3F3F3F371F3F3F3F 3F3F4F3F0F6F3F3F4F3F3F371F3F3F3F3F3F3F4F3F3,&F)F+*&#F/F,F/F0F/F+F1F3F4 F3F371F3F3F3F3F3F3F3F)F3F1F6F1F3F3F371F3F3F3F3F3F3F3F3F4F3F1F(F3F3F371 F3F3F3F3F3F3F3F3F3F4F3F3F(F1F371F3F3F3F3F3F3F3F3F3F3F3F3F1,$F0F*F371F3 F3F3F3F3F3F3F3F3F3F3F3F3F3F3" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#*(%3As sembling~elementG\"\"\"7&\"#5\"#9\"#6F%F%%\":GF%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7171,&%#kxG#!\"\"\"\"#*&#\"\"\"F,F/%#kyGF/F +,$F0#F/F,\"\"!,$F)F2F3F3F3F3F3F3F3F3F3F3F371F1,&F)F+F0F+F1F3F)F3F3F3F 3F3F3F3F3F3F371F3F1F(F3F3F4F3F3F3F3F3F3F3F3F371F4F3F3F6F0F3F4F3F3F3F3F 3F3F3F371F3F)F3F0,&F)!\"#*&F,F/F0F/F+F0F3F)F3F3F3F3F3F3F371F3F3F4F3F0F 6F3F3F4F3F3F3F3F3F371F3F3F3F4F3F3F6F0F3F4F3F3F3F3F371F3F3F3F3F)F3F0F:F 0F3F)F3F3F3F371F3F3F3F3F3F4F3F0F6F3F3F4F3F3F371F3F3F3F3F3F3F4F3F3F6F0F 3F4F3F371F3F3F3F3F3F3F3F)F3F0,&F)#!\"$F,*&#\"\"$F,F/F0F/F+F1F3F4F371F3 F3F3F3F3F3F3F3F4F3F1F(F3F3F371F3F3F3F3F3F3F3F3F3F4F3F3F(F1F371F3F3F3F3 F3F3F3F3F3F3F4F3F1F(F371F3F3F3F3F3F3F3F3F3F3F3F3F3F3F3" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#*(%3Assembling~elementG\"\"\"7&\"#6\"#9\"#7F%F%%\" :GF%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'matrixG6#7171,&%#kxG#!\"\" \"\"#*&#\"\"\"F,F/%#kyGF/F+,$F0#F/F,\"\"!,$F)F2F3F3F3F3F3F3F3F3F3F3F37 1F1,&F)F+F0F+F1F3F)F3F3F3F3F3F3F3F3F3F371F3F1F(F3F3F4F3F3F3F3F3F3F3F3F 371F4F3F3F6F0F3F4F3F3F3F3F3F3F3F371F3F)F3F0,&F)!\"#*&F,F/F0F/F+F0F3F)F 3F3F3F3F3F3F371F3F3F4F3F0F6F3F3F4F3F3F3F3F3F371F3F3F3F4F3F3F6F0F3F4F3F 3F3F3F371F3F3F3F3F)F3F0F:F0F3F)F3F3F3F371F3F3F3F3F3F4F3F0F6F3F3F4F3F3F 371F3F3F3F3F3F3F4F3F3F6F0F3F4F3F371F3F3F3F3F3F3F3F)F3F0F:F0F3F)F371F3F 3F3F3F3F3F3F3F4F3F0,&F)F*F0F+F3F3F371F3F3F3F3F3F3F3F3F3F4F3F3F(F1F371F 3F3F3F3F3F3F3F3F3F3F)F3F1,&F)F+*&#F/F,F/F0F/F+F371F3F3F3F3F3F3F3F3F3F3 F3F3F3F3F3" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#*(%3Assembling~elementG \"\"\"7&\"#7\"#9\"#:F%F%%\":GF%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%' matrixG6#7171,&%#kxG#!\"\"\"\"#*&#\"\"\"F,F/%#kyGF/F+,$F0#F/F,\"\"!,$F )F2F3F3F3F3F3F3F3F3F3F3F371F1,&F)F+F0F+F1F3F)F3F3F3F3F3F3F3F3F3F371F3F 1F(F3F3F4F3F3F3F3F3F3F3F3F371F4F3F3F6F0F3F4F3F3F3F3F3F3F3F371F3F)F3F0, &F)!\"#*&F,F/F0F/F+F0F3F)F3F3F3F3F3F3F371F3F3F4F3F0F6F3F3F4F3F3F3F3F3F 371F3F3F3F4F3F3F6F0F3F4F3F3F3F3F371F3F3F3F3F)F3F0F:F0F3F)F3F3F3F371F3F 3F3F3F3F4F3F0F6F3F3F4F3F3F371F3F3F3F3F3F3F4F3F3F6F0F3F4F3F371F3F3F3F3F 3F3F3F)F3F0F:F0F3F)F371F3F3F3F3F3F3F3F3F4F3F0F6F3F3F471F3F3F3F3F3F3F3F 3F3F4F3F3F(F1F371F3F3F3F3F3F3F3F3F3F3F)F3F1F6F171F3F3F3F3F3F3F3F3F3F3F 3F4F3F1F(" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#%:Assembling~point~source s:G" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'vectorG6#71!#]\"\"!F(F(F(F(F (F(F(F(F(F(F(F(F(" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#%F,F?F,F*F,F,*&FAF,F?F,F :*&F1F,F0F,F:*(\"#BF,F5F,F2F,F:*(\"#WF,F*F,F7F,F:F,*&FCF,F*F,F:FIF%,$* &,,FGF,*(\"\")F,F5F,F2F,F,*(FgnF,F*F,F7F,F,*(F6F,F*F,F2F,F:*&F6F,F7F,F :F,*&F*F,,(*$F5F,F,*(FgnF,F*F,F2F,F,*&FgnF,F7F,F,F,F:FJ,$*&,.FGF6F]oF: *(FFF,F5F,F2F,F,*(FFF,F*F,F7F,F,*(F+F,F*F,F2F,F:*&F+F,F7F,F:F,*&F*F,F \\oF,F:#FJF6FY,$*&,0F(F6*(\"#?F,F0F,F2F,F,*(FJF,F5F,F2F,F:*(\"#[F,F5F, F7F,F,*(\"#7F,F*F,F7F,F:*(\"#KF,F?F,F*F,F,*&FgnF,F?F,F:F,*&F*F,FCF,F:F I,$*&,2F(F6FGF:*(F]pF,F0F,F2F,F,*(FJF,F5F,F2F,F:*(F`pF,F5F,F7F,F,*(Fbp F,F*F,F7F,F:*(FdpF,F?F,F*F,F,*&FgnF,F?F,F:F,*&F*F,FCF,F:FIFio\"\"!FaqF aq" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#*&%9Boundary~fluxes~(bflux):G\" \"\"7$%%nodeG7#%%fluxGF%" }}{PARA 12 "" 1 "" {XPPMATH 20 "6-7$\"\"\"7# ,$*&,0*&)%#kxG\"\"$F$%#kyGF$\"#?*(\"#KF$)F-F,F$F+F$F$*&\"\"#F$)F+\"\"% F$F$*(\"#[F$)F+F3F$)F-F3F$F$*(\"\"(F$F8F$F-F$F$*(\"#GF$F+F$F9F$F$*&\"# CF$F1F$F$F$,**&F8F$F-F$\"#5*&\"#;F$F1F$F$*$F*F$F$*(F?F$F+F$F9F$F$!\"\" #!\"&\"\")7$F37#,$*&,2F)F.*(F0F$F1F$F+F$F$*&F3F$F4F$F$*(F7F$F8F$F9F$F$ *&F?F$F1F$F$*&F,F$F*F$F$*(\"#BF$F8F$F-F$F$*(\"#WF$F+F$F9F$F$F$F@FG#FIF 57$F,F%7$\"#87#,$*&,0*$F4F$F3*(F.F$F*F$F-F$F$*(\"\"&F$F8F$F-F$FG*(F7F$ F8F$F9F$F$*(\"#7F$F+F$F9F$FG*(F0F$F1F$F+F$F$*&FJF$F1F$FGF$F@FG#F^oFJ7$ \"#97#,$*&,2F[oF3FEFG*(F.F$F*F$F-F$F$*(F^oF$F8F$F-F$FG*(F7F$F8F$F9F$F$ *(FaoF$F+F$F9F$FG*(F0F$F1F$F+F$F$*&FJF$F1F$FGF$F@FG#F^oF57$\"#:Fgn7$F$ 7#!#X7$F37#F^o7$F,7#\"\"!7$F5Fjp7$F^oFhp" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#**%4Flux~equilibrium:~~G\"\"\"%&SigmaGF%%'fluxesGF%%&=~-75GF%" } }{PARA 11 "" 1 "" {XPPMATH 20 "6#*&%CElement~flow~gradients~(e_grads): ~G\"\"\"7$&%%gradG6#%\"xG&F(6#%\"yGF%" }}{PARA 12 "" 1 "" {XPPMATH 20 "627$,$*&,0*&)%#kxG\"\"$\"\"\"%#kyGF+\"#?*(\"#KF+)F,F*F+F)F+F+*&\"\"#F +)F)\"\"%F+F+*(\"#[F+)F)F2F+)F,F2F+F+*(\"\"(F+F7F+F,F+F+*(\"#GF+F)F+F8 F+F+*&\"#CF+F0F+F+F+*&F)F+,**&F7F+F,F+\"#5*&\"#;F+F0F+F+*$F(F+F+*(F>F+ F)F+F8F+F+F+!\"\"#!\"&F4,$*&,(*&F)F+F,F+FD*&FDF+F8F+F+*&F*F+F7F+F+F+F@ FGFH7$,$*&,2F'F-*(F/F+F0F+F)F+F+*&F2F+F3F+F+*(F6F+F7F+F8F+F+*&F>F+F0F+ F+*&F*F+F(F+F+*(\"#BF+F7F+F,F+F+*(\"#WF+F)F+F8F+F+F+*&F)F+F@F+FGFH\"\" !FP7$F$,$FK#\"\"&F47$,$*&,0*&F0F+F)F+F/*&\"\")F+F0F+FG*(F6F+F7F+F8F+F+ *(F4F+F)F+F8F+FG*(F-F+F(F+F,F+F+FAF+*&F2F+F3F+F+F+*(,(*$F7F+F+*(FcoF+F )F+F,F+F+*&FcoF+F8F+F+F+,&F,F2F)F+F+F)F+FGFH,$*&F)F+FioFG#FIF27$,$*&,2 FaoF/*&FcoF+F0F+FG*(F6F+F7F+F8F+F+*(F-F+F)F+F8F+FG*(F-F+F(F+F,F+F+*(\" #6F+F7F+F,F+FG*&F2F+F3F+F+FEFGF+*(FioF+F]pF+F)F+FGFHFJ7$FbpFjn7$F^o,$F _p#F\\oF27$,$*&,0FaoF/*&FcoF+F0F+FG*(\"#7F+F)F+F8F+FG*(F6F+F7F+F8F+F+* (F-F+F(F+F,F+F+*(F*F+F7F+F,F+FG*&F2F+F3F+F+F+*(FioF+F]pF+F)F+FGFH,$*&* $F7F+F+F@FGFH7$,$*&,2FaoF/*&FcoF+F0F+FG*(FgqF+F)F+F8F+FG*(F6F+F7F+F8F+ F+*(F-F+F(F+F,F+F+*(F:F+F7F+F,F+FG*&F2F+F3F+F+FEFGF+*(FioF+F]pF+F)F+FG FHF^p7$FarF_q7$Fbq,$F^rF[o7$,$*&,0*$F3F+F2*(F-F+F(F+F,F+F+*(F\\oF+F7F+ F,F+FG*(F6F+F7F+F8F+F+*(FgqF+F)F+F8F+FG*(F/F+F0F+F)F+F+*&FcoF+F0F+FGF+ *&F)F+F@F+FGFHFhn7$,$*&,2FbsF2FEFG*(F-F+F(F+F,F+F+*(F\\oF+F7F+F,F+FG*( F6F+F7F+F8F+F+*(FgqF+F)F+F8F+FG*(F/F+F0F+F)F+F+*&FcoF+F0F+FGF+*&F)F+F@ F+FGFHF]r7$F[tF]sF^s" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#*&%DElement~fl ow~velocities~(e_velts):~G\"\"\"7$&%\"vG6#%\"xG&F(6#%\"yGF%" }}{PARA 12 "" 1 "" {XPPMATH 20 "627$,$*&,H*()-%$cosG6#,$*&%&alphaG\"\"\"%#PiGF /#F/\"$!=\"\"#F/)%#kxG\"\"%F/%#kyGF/!#=**\"#;F/F(F/)F5F3F/)F7\"\"$F/F/ *(F3F/F(F/)F5\"\"&F/!\"\"**\"#GF/F(F/)F5F=F/)F7F3F/FA**\"\"(F/F(F/FDF/ F7F/FA**\"#@F/F(F/F;F/FEF/FA**F6F/F(F/F5F/F**FMF/F;F/F)F/F**FMF/FF+F3F+F(F+F+*(F]wF+F3F+F6F+F+*(\"#_F+F6F+F-F+F+**FjwF +FQF+F-F+F6F+FD*(FNF+FIF+F-F+F+**FNF+FIF+F-F+FQF+FD*&F/F+FIF+F+*(F/F+F IF+FQF+FDF+*&F`oF+F-F+FDFiw,$*&*(F9F+FBF+,6*$FhnF+F_v*(FdwF+F)F+FLF+F+ *&F*F+FLF+F+*(FjyF+F,F+F)F+F+*(FgwF+F(F+F,F+F+*(F_zF+F3F+F(F+F+*(FHF+F 3F+F6F+FD*(FjwF+F6F+F-F+FD*(FNF+FIF+F-F+FD*&F/F+FIF+FDF+F+*&F-F+F`oF+F DF`y7$,$*&,NF'FgvFdqFD*(F.F+F3F+F(F+FD*(FJF+F6F+F-F+F+*(FgoF+FIF+FQF+F D**FdwF+FQF+FLF+F)F+F+*(F_vF+FQF+FhnF+F+*,FbwF+F9F+F,F+FBF+F)F+F+*,FJF +FBF+F3F+F9F+F(F+FD*,FSF+FBF+F6F+F9F+F-F+FDF`vFD**FgoF+F9F+FLF+FBF+F+* *F*F+FQF+F,F+F)F+FD**FdvF+FQF+F3F+F(F+F+**FJF+FQF+F-F+F6F+F+**FgwF+FQF +F,F+F(F+F+**FHF+FQF+F3F+F6F+FD**FNF+FIF+F-F+FQF+FD*&FgoF+FIF+F+*(F_vF +F)F+FLF+F+*(F]wF+F3F+F6F+F+*(FNF+FIF+F-F+F+F+*&F-F+F`oF+FDFiw,$*&,DF^ xF_v*&FgoF+FLF+F+*(FgoF+FLF+FQF+FDFbxFD*,FdwF+FBF+FLF+F9F+F)F+F+*,F*F+ F9F+F,F+FBF+F)F+FD**FbwF+FQF+F,F+F)F+FD*,FgwF+FBF+F,F+F9F+F(F+F+*(FPF+ F,F+F)F+F+*,FHF+FBF+F3F+F9F+F6F+FD*,FdvF+FBF+F3F+F9F+F(F+F+**FJF+FQF+F 3F+F(F+F+*(FSF+F3F+F(F+F+*,FNF+F9F+FIF+FBF+F-F+FD*,FJF+FBF+F6F+F9F+F-F +F+**FSF+FQF+F-F+F6F+F+**FgoF+F9F+FIF+FBF+FDF+*&F-F+F`oF+FDF`y7$,$*&,L F`vFD*(FgvF+F(F+F,F+F+FdqFD*(\"\"*F+F3F+F(F+FD*(FbwF+F6F+F-F+FD*(FenF+ FIF+FQF+F+*&FenF+FIF+FD*(F_vF+F)F+FLF+F+*(F]wF+F3F+F6F+F+*,F.F+F9F+F,F +FBF+F)F+FD*(FNF+FIF+F-F+F+**FJF+FQF+F-F+F6F+F+*,FgoF+FBF+F3F+F9F+F(F+ F+**FbzF+FQF+F3F+F(F+FD**FgoF+FQF+F,F+F)F+FD**FinF+F9F+FLF+FBF+FD*(F_v F+FQF+FhnF+F+**FNF+FIF+F-F+FQF+FD**FHF+FQF+F3F+F6F+FD**FgwF+FQF+F,F+F( F+F+**FdwF+FQF+FLF+F)F+F+F+*&F-F+F`oF+FDFiw,$*&,@F^xF_v*,FdwF+FBF+FLF+ F9F+F)F+F+*&FinF+FLF+FD*(FinF+FLF+FQF+F+FbxFD*,FgoF+F9F+F,F+FBF+F)F+FD *(FgoF+F,F+F)F+FD**F.F+FQF+F,F+F)F+F+*,FgwF+FBF+F,F+F9F+F(F+F+**FgoF+F QF+F3F+F(F+FD*,FHF+FBF+F3F+F9F+F6F+FD*,FbzF+FBF+F3F+F9F+F(F+FD*,FNF+F9 F+FIF+FBF+F-F+FD*,FJF+FBF+F6F+F9F+F-F+F+**FenF+F9F+FIF+FBF+F+F+*&F-F+F `oF+FDF`y7$,$*&,FFeyF_v*(F*F+FLF+FQF+FD*(F_vF+F)F+FLF+F+**FdwF+FQF+FLF +F)F+F+*(F*F+F,F+F)F+FD**\"#:F+FQF+F,F+F)F+FD*(FgvF+F(F+F,F+F+**FgwF+F QF+F,F+F(F+F+**FbzF+FQF+F3F+F(F+FD*(FjyF+F3F+F(F+FD*(F]wF+F3F+F6F+F+** FHF+FQF+F3F+F6F+FD*(FNF+FIF+F-F+F+**FNF+FIF+F-F+FQF+FD**FjwF+FQF+F-F+F 6F+F+*(FPF+F6F+F-F+FD*&FenF+FIF+FD*(FenF+FIF+FQF+F+F+*&F`oF+F-F+FDFiw, $*&*(FBF+F9F+,6F`[lF_v*(FdwF+F)F+FLF+F+*&F*F+FLF+FD*(FialF+F,F+F)F+FD* (FgwF+F(F+F,F+F+*(FHF+F3F+F6F+FD*(FbzF+F3F+F(F+FD*(FNF+FIF+F-F+FD*(Fjw F+F6F+F-F+F+*&FenF+FIF+F+F+F+*&F`oF+F-F+FDF`y7$,$*&,LF'FgvFdqFD*(F\\_l F+F3F+F(F+FD*(FbwF+F6F+F-F+FD*(FenF+FIF+FQF+F+**FdwF+FQF+FLF+F)F+F+**F gwF+FQF+F,F+F(F+F+**FgoF+FQF+F,F+F)F+FDF`vFD**FbzF+FQF+F3F+F(F+FD**FJF +FQF+F-F+F6F+F+**FinF+F9F+FLF+FBF+F+**FHF+FQF+F3F+F6F+FD*(F_vF+FQF+Fhn F+F+**FNF+FIF+F-F+FQF+FD*&FenF+FIF+FD*(F_vF+F)F+FLF+F+*(F]wF+F3F+F6F+F +*(FNF+FIF+F-F+F+*,FgoF+FBF+F3F+F9F+F(F+FD*,F.F+F9F+F,F+FBF+F)F+F+F+*& F-F+F`oF+FDFiw,$*&,@F^xF_v*&FinF+FLF+F+*(FinF+FLF+FQF+FDFbxFD*,FdwF+FB F+FLF+F9F+F)F+F+*(FgoF+F,F+F)F+F+**F.F+FQF+F,F+F)F+FD*,FgoF+F9F+F,F+FB F+F)F+FD*,FgwF+FBF+F,F+F9F+F(F+F+*,FbzF+FBF+F3F+F9F+F(F+FD**FgoF+FQF+F 3F+F(F+F+*,FHF+FBF+F3F+F9F+F6F+FD*,FNF+F9F+FIF+FBF+F-F+FD*,FJF+FBF+F6F +F9F+F-F+F+**FenF+F9F+FIF+FBF+F+F+*&F-F+F`oF+FDF`y7$,$*&,JF'FgvFdqFD*( F8F+F3F+F(F+FD*(FPF+F6F+F-F+FDF`vFD*(F_vF+FQF+FhnF+F+*,F*F+F9F+F,F+FBF +F)F+F+*&FenF+FIF+FD*(F_vF+F)F+FLF+F+*(F]wF+F3F+F6F+F+*(FNF+FIF+F-F+F+ **FNF+FIF+F-F+FQF+FD*(FenF+FIF+FQF+F+**F*F+F9F+FLF+FBF+FD**FjwF+FQF+F- F+F6F+F+**FHF+FQF+F3F+F6F+FD**F\\rF+FQF+F3F+F(F+FD**FgwF+FQF+F,F+F(F+F +**FjwF+FQF+F,F+F)F+FD**FdwF+FQF+FLF+F)F+F+F+*&F-F+F`oF+FDFiw,$*&,F/F)F/F:F/FEF/F %F/F:F/FF/F(F/F:F/FKF/FH**\"#pF/F(F/FCF/FGF/FH**F\\uF/F)F/F@F/FEF/F/**F_i lF/F(F/FgoF/F=F/F/**F\\alF/F(F/FTF/FenF/FH**F][lF/F(F/FOF/FYF/FH**F[oF /F(F/F:F/FYF/F/**F9F/F(F/FOF/FenF/F/**\"%!3#F/F(F/F:F/F\"\"'F.)F4\"\"$F.F.*(F8F.F>F.) F4F8F.F.*(\"%39F.F?F.)F>\"\"&F.F.*(\"$Y&F.)F4\"\"%F.)F>FEF.F.*(\"%%=\" F.FPF.)F4FCF.F.*(\"%orF.FJF.FNF.F.*(\"$g\"F.)F>F2F.F3F.F.*(\"%'4%F.)F> FOF.)F4FKF.F.*(\"$I#F.FfnF.FXF.F.*(\"%35F.FDF.FenF.F.*&F5F.F3F.F.*&\"$ %QF.F=F.F.*(\"%_6F.FBF.F4F.F.*(F]oF.F=F.F'F.!\"\"*(\"#fF.FSF.F>F.F.** \"$_\"F.F'F.FGF.F>F.F.**\"%sIF.F'F.FfnF.FenF.F.**\"%C5F.F'F.F3F.FXF.F. **FF.Fao**\"$o(F.F'F.F4F.FBF.Fao**\" $G\"F.F'F.F?F.FJF.Fao**\"$c'F.F'F.FDF.FenF.F.**\"$=&F.F'F.FNF.FPF.F.** \"$;\"F.F'F.FfnF.FXF.F.**FioF.F'F.FNF.FJF.Fao**FZF.F'F.FDF.FBF.Fao**\" %7HF.F'F.FSF.FPF.F.F.*&F?F.),**&F?F.F>F.\"#5*&\"#;F.FPF.F.*$FDF.F.*(\" #CF.F4F.FXF.F.F2F.Fao#\"#DF8" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#*&%?No dal~flow~fluxes~(n_fluxes):~G\"\"\"7$&%\"qG6#%\"xG&F(6#%\"yGF%" }} {PARA 12 "" 1 "" {XPPMATH 20 "617$,$*&,H*()-%$cosG6#,$*&%&alphaG\"\"\" %#PiGF/#F/\"$!=\"\"#F/)%#kxG\"\"%F/%#kyGF/!#=**\"#;F/F(F/)F5F3F/)F7\" \"$F/F/*(F3F/F(F/)F5\"\"&F/!\"\"**\"#GF/F(F/)F5F=F/)F7F3F/FA**\"\"(F/F (F/FDF/F7F/FA**\"#@F/F(F/F;F/FEF/FA**F6F/F(F/F5F/F " 0 "" {MPLTEXT 1 0 16 "n_potls;e_grads;" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 66 "subs(kx=2,ky=2,alpha=45,e_velts);subs(kx= 2,ky=2,alpha=45,n_velts);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "subs(k x=2,ky=2,alpha=45,n_fluxes);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 79 "sub s(kx=2,ky=2,alpha=45,n_kined);subs(kx=2,ky=2,alpha=45,total_kinetic_en ergy);" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#71\"#5F$F$,$*&,0*$)%#kxG\"\" %\"\"\"\"\"'*(\"#gF,)F*\"\"$F,%#kyGF,F,*(\"$W\"F,)F*\"\"#F,)F2F6F,F,*( \"\"(F,F5F,F2F,!\"\"*(\"#GF,F*F,F7F,F:*(\"#'*F,)F2F1F,F*F,F,*&\"#CF,F? F,F:F,*&F*F,,**&F5F,F2F,F$*&\"#;F,F?F,F,*$F0F,F,*(FAF,F*F,F7F,F,F,F:# \"\"&F+,$*&,2*&F5F,F7F,F4*(F/F,F0F,F2F,F,*&F-F,F)F,F,*(F>F,F?F,F*F,F,* &FAF,F?F,F:*&F1F,F0F,F:*(\"#BF,F5F,F2F,F:*(\"#WF,F*F,F7F,F:F,*&FCF,F*F ,F:FIF%,$*&,,FGF,*(\"\")F,F5F,F2F,F,*(FgnF,F*F,F7F,F,*(F6F,F*F,F2F,F:* &F6F,F7F,F:F,*&F*F,,(*$F5F,F,*(FgnF,F*F,F2F,F,*&FgnF,F7F,F,F,F:FJ,$*&, .FGF6F]oF:*(FFF,F5F,F2F,F,*(FFF,F*F,F7F,F,*(F+F,F*F,F2F,F:*&F+F,F7F,F: F,*&F*F,F\\oF,F:#FJF6FY,$*&,0F(F6*(\"#?F,F0F,F2F,F,*(FJF,F5F,F2F,F:*( \"#[F,F5F,F7F,F,*(\"#7F,F*F,F7F,F:*(\"#KF,F?F,F*F,F,*&FgnF,F?F,F:F,*&F *F,FCF,F:FI,$*&,2F(F6FGF:*(F]pF,F0F,F2F,F,*(FJF,F5F,F2F,F:*(F`pF,F5F,F 7F,F,*(FbpF,F*F,F7F,F:*(FdpF,F?F,F*F,F,*&FgnF,F?F,F:F,*&F*F,FCF,F:FIFi o\"\"!FaqFaq" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#727$,$*&,0*&)%#kxG\"\" $\"\"\"%#kyGF,\"#?*(\"#KF,)F-F+F,F*F,F,*&\"\"#F,)F*\"\"%F,F,*(\"#[F,)F *F3F,)F-F3F,F,*(\"\"(F,F8F,F-F,F,*(\"#GF,F*F,F9F,F,*&\"#CF,F1F,F,F,*&F *F,,**&F8F,F-F,\"#5*&\"#;F,F1F,F,*$F)F,F,*(F?F,F*F,F9F,F,F,!\"\"#!\"&F 5,$*&,(*&F*F,F-F,FE*&FEF,F9F,F,*&F+F,F8F,F,F,FAFHFI7$,$*&,2F(F.*(F0F,F 1F,F*F,F,*&F3F,F4F,F,*(F7F,F8F,F9F,F,*&F?F,F1F,F,*&F+F,F)F,F,*(\"#BF,F 8F,F-F,F,*(\"#WF,F*F,F9F,F,F,*&F*F,FAF,FHFI\"\"!FQ7$F%,$FL#\"\"&F57$,$ *&,0*&F1F,F*F,F0*&\"\")F,F1F,FH*(F7F,F8F,F9F,F,*(F5F,F*F,F9F,FH*(F.F,F )F,F-F,F,FBF,*&F3F,F4F,F,F,*(,(*$F8F,F,*(FdoF,F*F,F-F,F,*&FdoF,F9F,F,F ,,&F-F3F*F,F,F*F,FHFI,$*&F*F,FjoFH#FJF37$,$*&,2FboF0*&FdoF,F1F,FH*(F7F ,F8F,F9F,F,*(F.F,F*F,F9F,FH*(F.F,F)F,F-F,F,*(\"#6F,F8F,F-F,FH*&F3F,F4F ,F,FFFHF,*(FjoF,F^pF,F*F,FHFIFK7$FcpF[o7$F_o,$F`p#F]oF37$,$*&,0FboF0*& FdoF,F1F,FH*(\"#7F,F*F,F9F,FH*(F7F,F8F,F9F,F,*(F.F,F)F,F-F,F,*(F+F,F8F ,F-F,FH*&F3F,F4F,F,F,*(FjoF,F^pF,F*F,FHFI,$*&*$F8F,F,FAFHFI7$,$*&,2Fbo F0*&FdoF,F1F,FH*(FhqF,F*F,F9F,FH*(F7F,F8F,F9F,F,*(F.F,F)F,F-F,F,*(F;F, F8F,F-F,FH*&F3F,F4F,F,FFFHF,*(FjoF,F^pF,F*F,FHFIF_p7$FbrF`q7$Fcq,$F_rF \\o7$,$*&,0*$F4F,F3*(F.F,F)F,F-F,F,*(F]oF,F8F,F-F,FH*(F7F,F8F,F9F,F,*( FhqF,F*F,F9F,FH*(F0F,F1F,F*F,F,*&FdoF,F1F,FHF,*&F*F,FAF,FHFIFin7$,$*&, 2FcsF3FFFH*(F.F,F)F,F-F,F,*(F]oF,F8F,F-F,FH*(F7F,F8F,F9F,F,*(FhqF,F*F, F9F,FH*(F0F,F1F,F*F,F,*&FdoF,F1F,FHF,*&F*F,FAF,FHFIF^r7$F\\tF^sF_s" }} {PARA 12 "" 1 "" {XPPMATH 20 "6#727$#\"%:8\"$/##\"$v\"F'7$#\"$X(\"$-\" \"\"!F*7$F%#!$v\"F'7$#\"$l*F'#\"\"&\"#M7$#\"$0#\"#^F(7$F9F07$F3#!\"&F7 7$#\"$0*F'#F6F'7$#\"$?#F;F57$FEF>7$FA#F?F'7$#\"$&*)F'F.7$#\"$X%F-FC7$F NFIFJ" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#717$#\"#b\"\")#\"$v\"\"$3%7$# \"$X(\"$-\"\"\"!7$F%#!$v\"F*7$#\"$v(\"$`\"#\"#&*F67$#\"%DO\"$7'F/7$F4# !#&*F67$#\"%v8\"$1$#\"#lF<7$#\"%lEF " 0 "" {MPLTEXT 1 0 0 "" }}}}{PARA 4 "" 0 "" {TEXT 258 1 " " }}} {MARK "8 13 0 0" 0 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }