{VERSION 4 0 "Mac OS X" "4.0" } {USTYLETAB {PSTYLE "Heading 4" -1 20 1 {CSTYLE "" -1 -1 "Times" 1 10 0 0 0 1 1 1 2 2 2 2 1 0 0 1 }1 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Ord ered List 5" -1 200 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }1 1 0 0 3 3 2 144 2 0 2 2 -1 1 }{PSTYLE "Ordered List 1" -1 201 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }1 1 0 0 3 3 2 0 2 0 2 2 -1 1 }{PSTYLE "Text Output" -1 2 1 {CSTYLE "" -1 -1 "Courier" 1 10 0 0 255 1 2 2 2 2 2 1 1 0 0 1 }1 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Bullet Item" -1 15 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }1 1 0 0 3 3 2 0 2 0 2 2 -1 1 }{PSTYLE "Lef t Justified Maple Output" -1 12 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }1 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Help" -1 10 1 {CSTYLE "" -1 -1 "Courier" 1 9 0 0 255 1 2 2 2 2 2 1 1 0 0 1 }1 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Author" -1 19 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }3 1 0 0 8 8 2 0 2 0 2 2 -1 1 }{PSTYLE "Diagnostic" -1 9 1 {CSTYLE "" -1 -1 "Courier" 1 10 64 128 64 1 2 2 2 2 2 1 1 0 0 1 }1 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Norma l259" -1 202 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 1 1 2 2 2 2 1 0 0 1 }3 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Heading 2" -1 4 1 {CSTYLE " " -1 -1 "Times" 1 14 0 0 0 1 2 1 2 2 2 2 1 0 0 1 }3 1 0 0 4 4 2 0 2 0 2 2 -1 1 }{PSTYLE "Ordered List 3" -1 203 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }1 1 0 0 3 3 2 72 2 0 2 2 -1 1 } {PSTYLE "Maple Plot" -1 13 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }3 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }1 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Dash Item" -1 16 1 {CSTYLE "" -1 -1 "Tim es" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }1 1 0 0 3 3 2 0 2 0 2 2 -1 1 } {PSTYLE "Heading 3" -1 5 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 1 1 2 2 2 2 1 0 0 1 }3 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Ordered List 4" -1 204 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }1 1 0 0 3 3 2 108 2 0 2 2 -1 1 }{PSTYLE "Maple Output" -1 11 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }3 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Line Printed Output" -1 6 1 {CSTYLE "" -1 -1 "Cour ier" 1 10 0 0 255 1 2 2 2 2 2 1 1 0 0 1 }1 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "List Item" -1 14 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }1 1 0 0 3 3 2 0 2 0 2 2 -1 1 }{PSTYLE "Normal258" -1 205 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }3 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Warning" -1 7 1 {CSTYLE "" -1 -1 "C ourier" 1 10 0 0 255 1 2 2 2 2 2 1 1 0 0 1 }1 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Error" -1 8 1 {CSTYLE "" -1 -1 "Courier" 1 10 255 0 255 1 2 2 2 2 2 1 1 0 0 1 }1 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Fixed Wid th" -1 17 1 {CSTYLE "" -1 -1 "Courier" 1 10 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }3 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Heading 1" -1 3 1 {CSTYLE "" -1 -1 "Times" 1 14 0 0 0 1 2 1 2 2 2 2 1 0 0 1 }1 1 0 0 6 4 2 0 2 0 2 2 -1 1 }{PSTYLE "Title" -1 18 1 {CSTYLE "" -1 -1 "Times" 1 18 0 0 0 1 2 1 2 2 2 2 1 0 0 1 }3 1 0 0 12 12 2 0 2 0 2 2 -1 1 }{PSTYLE "Exercis e" -1 206 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }1 1 0 -12 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Normal257" -1 207 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 1 2 2 2 2 2 1 0 0 1 }3 1 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Ordered List 2" -1 208 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 0 0 1 }1 1 0 0 3 3 2 36 2 0 2 2 -1 1 } {CSTYLE "Help Variable" -1 25 "Courier" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Text" -1 200 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 } {CSTYLE "Help Bold" -1 39 "Times" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 } {CSTYLE "Page Number" -1 33 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 } {CSTYLE "2D Math Italic Small" -1 201 "Times" 1 1 0 0 0 1 1 2 2 2 2 2 0 0 0 1 }{CSTYLE "Help Nonterminal" -1 24 "Courier" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "Default" -1 38 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Maple Comment" -1 21 "Courier" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "Maple Input" -1 0 "Courier" 1 12 255 0 0 1 2 1 2 2 1 2 0 0 0 1 }{CSTYLE "2D Math Small" -1 7 "Times" 1 1 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Help Fixed" -1 23 "Courier" 1 10 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Popup" -1 31 "Times" 1 12 0 128 128 1 1 2 1 2 2 2 0 0 0 1 }{CSTYLE "Plot Title" -1 27 "Times" 1 10 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "2D Input" -1 19 "Times" 1 12 255 0 0 1 2 2 2 2 1 2 0 0 0 1 }{CSTYLE "Copyright" -1 34 "Times" 1 10 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Maple Input Placeholder" -1 202 "Courier" 1 12 200 0 200 1 2 1 2 2 1 2 0 0 0 1 }{CSTYLE "2D Math Bold Small" -1 10 "T imes" 1 1 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "2D Math" -1 2 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Help Notes" -1 37 "Times" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "Help Underlined Bold" -1 41 "Times" 1 12 0 0 0 1 1 1 2 2 2 2 0 0 0 1 }{CSTYLE "Hyperlink" -1 17 "T imes" 1 12 0 128 128 1 2 2 1 2 2 2 0 0 0 1 }{CSTYLE "2D Math Symbol 2" -1 16 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Plot Text" -1 28 "Times" 1 8 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "2D Math Italic Small201" -1 203 "Times" 1 1 0 0 0 1 1 2 2 2 2 2 0 0 0 1 }{CSTYLE "He lp Italic" -1 42 "Times" 1 12 0 0 0 1 1 2 2 2 2 2 0 0 0 1 }{CSTYLE "Ou tput Labels" -1 29 "Times" 1 8 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "H elp Heading" -1 26 "Times" 1 14 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE " Help Normal" -1 30 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE " 2D Comment" -1 18 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "H elp Maple Name" -1 35 "Times" 1 12 104 64 92 1 2 1 2 2 2 2 0 0 0 1 } {CSTYLE "2D Output" -1 20 "Times" 1 12 0 0 255 1 2 2 2 2 2 2 0 0 0 1 } {CSTYLE "Dictionary Hyperlink" -1 45 "Times" 1 12 147 0 15 1 2 2 1 2 2 2 0 0 0 1 }{CSTYLE "Help Emphasized" -1 22 "Times" 1 12 0 0 0 1 1 2 2 2 2 2 0 0 0 1 }{CSTYLE "Help Italic Bold" -1 40 "Times" 1 12 0 0 0 1 1 2 2 2 2 2 0 0 0 1 }{CSTYLE "LaTeX" -1 32 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "Help Menus" -1 36 "Times" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "Prompt" -1 1 "Courier" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "Help Underlined" -1 44 "Times" 1 12 0 0 0 1 2 2 1 2 2 2 0 0 0 1 }{CSTYLE "Help Underlined Italic" -1 43 "Times" 1 12 0 0 0 1 1 2 1 2 2 2 0 0 0 1 }{CSTYLE "2D Math Bold" -1 5 "Times" 1 12 0 0 0 1 2 1 2 2 2 2 0 0 0 1 }{CSTYLE "2D Math Italic" -1 3 "Times" 1 12 0 0 0 1 1 2 2 2 2 2 0 0 0 1 }{PSTYLE "" -1 209 1 {CSTYLE "" -1 -1 "Times " 1 14 0 0 0 1 2 1 2 2 2 2 1 0 0 1 }1 1 0 0 4 4 2 0 2 0 2 2 -1 1 }} {SECT 0 {EXCHG {PARA 18 "" 0 "" {TEXT 204 29 "Determinants by Row Redu ction" }}{PARA 0 "" 0 "" {TEXT 205 0 "" }}{PARA 0 "" 0 "" {TEXT 205 35 "Worksheet written by Russell Blyth." }}{PARA 0 "" 0 "" {TEXT 205 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 29 "restart: with(Linear Algebra):" }}}{SECT 1 {PARA 3 "" 0 "" {TEXT 206 7 "Outline" }}{EXCHG {PARA 0 "" 0 "" {TEXT 200 65 "Use Maple to compute determinants using \+ elementary row operations" }}}{EXCHG {PARA 0 "> " 0 "" {XPPEDIT 19 1 " Typesetting:-mo(\"⁢\", form = \"\", fence = \"false\", \+ separator = \"false\", lspace = \"0em\", rspace = \"0em\", stretchy = \+ \"false\", symmetric = \"false\", maxsize = \"infinity\", minsize = \" 1\", largeop = \"false\", movablelimits = \"false\", accent = \"false \", font_style_name = \"2D Input\", size = \"12\", foreground = \"[255 ,0,0]\", background = \"[255,255,255]\");" "-I#moG6#/I+modulenameG6\"I ,TypesettingGI(_syslibGF'63Q1⁢F'/%%formGQ!F'/%&fenceGQ& falseF'/%*separatorGF1/%'lspaceGQ$0emF'/%'rspaceGF6/%)stretchyGF1/%*sy mmetricGF1/%(maxsizeGQ)infinityF'/%(minsizeGQ\"1F'/%(largeopGF1/%.mova blelimitsGF1/%'accentGF1/%0font_style_nameGQ)2D~InputF'/%%sizeGQ#12F'/ %+foregroundGQ*[255,0,0]F'/%+backgroundGQ.[255,255,255]F'" }}}}{SECT 1 {PARA 3 "" 0 "" {TEXT 206 48 "Utilizing row operations to compute det erminants" }}{EXCHG {PARA 0 "" 0 "" {TEXT 205 145 "We compute the dete rminant of a matrix using elementary row operations. We need to keep t rack of the effect of each row operation as we go along." }}{PARA 0 "" 0 "" {TEXT 205 0 "" }}{PARA 0 "" 0 "" {TEXT 205 83 "We generate a ran dom 6x6 matrix, and check that it has full rank before proceeding." }} }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 47 "A := RandomMatrix(6,6,gener ator=rand(-5..5));\n" }{MPLTEXT 1 0 1 "R" }{MPLTEXT 1 0 7 "ank(A);" }} }{EXCHG {PARA 0 "" 0 "" {TEXT 205 203 "We reduce A to row echelon form , one step at a time. Start with the first column. For the row operati ons we note specifically that adding a multiple of one row to another \+ does not affect the determinant." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 54 "mf := -A[2,1]/A[1,1]: B := RowOperation(A,[2,1],mf):\n" } {MPLTEXT 1 0 22 "mf := -B[3,1]/B[1,1]: " }{MPLTEXT 1 0 32 "B := RowOpe ration(B,[3,1],mf):\n" }{MPLTEXT 1 0 22 "mf := -B[4,1]/B[1,1]: " } {MPLTEXT 1 0 32 "B := RowOperation(B,[4,1],mf):\n" }{MPLTEXT 1 0 22 "m f := -B[5,1]/B[1,1]: " }{MPLTEXT 1 0 32 "B := RowOperation(B,[5,1],mf) :\n" }{MPLTEXT 1 0 22 "mf := -B[6,1]/B[1,1]: " }{MPLTEXT 1 0 30 "B := \+ RowOperation(B,[6,1],mf);" }{MPLTEXT 1 0 2 "\n" }}}{EXCHG {PARA 0 "" 0 "" {TEXT 205 35 "Next move on the the second column:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 54 "mf := -B[3,2]/B[2,2]: B := RowOpera tion(B,[3,2],mf):\n" }{MPLTEXT 1 0 22 "mf := -B[4,2]/B[2,2]: " } {MPLTEXT 1 0 32 "B := RowOperation(B,[4,2],mf):\n" }{MPLTEXT 1 0 22 "m f := -B[5,2]/B[2,2]: " }{MPLTEXT 1 0 32 "B := RowOperation(B,[5,2],mf) :\n" }{MPLTEXT 1 0 22 "mf := -B[6,2]/B[2,2]: " }{MPLTEXT 1 0 30 "B := \+ RowOperation(B,[6,2],mf);" }}}{EXCHG {PARA 0 "" 0 "" {TEXT 205 17 "The third column:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 54 "mf := -B[ 4,3]/B[3,3]: B := RowOperation(B,[4,3],mf):\n" }{MPLTEXT 1 0 22 "mf := -B[5,3]/B[3,3]: " }{MPLTEXT 1 0 32 "B := RowOperation(B,[5,3],mf):\n" }{MPLTEXT 1 0 22 "mf := -B[6,3]/B[3,3]: " }{MPLTEXT 1 0 30 "B := RowO peration(B,[6,3],mf);" }}}{EXCHG {PARA 0 "" 0 "" {TEXT 205 14 "Fourth \+ column:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 54 "mf := -B[5,4]/B[ 4,4]: B := RowOperation(B,[5,4],mf):\n" }{MPLTEXT 1 0 22 "mf := -B[6,4 ]/B[4,4]: " }{MPLTEXT 1 0 30 "B := RowOperation(B,[6,4],mf);" }}} {EXCHG {PARA 0 "" 0 "" {TEXT 205 13 "Fifth column:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 52 "mf := -B[6,5]/B[5,5]: B := RowOperation(B,[ 6,5],mf);" }}}{EXCHG {PARA 0 "" 0 "" {TEXT 205 107 "Now the determinan t is the product of the diagonal entries. We compute the determinant d irectly as a check." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 44 "B[1, 1]*B[2,2]*B[3,3]*B[4,4]*B[5,5]*B[6,6];\n" }{MPLTEXT 1 0 15 "Determinan t(A);" }}}{SECT 1 {PARA 209 "" 0 "" {TEXT 200 8 "Exercise" }}{EXCHG {PARA 0 "" 0 "" {TEXT 200 245 "1) Generate a random 60 x 60 matrix an d time how long it takes to compute the determinant; compare this time with the time required for a 61x61 matrix. Based on your observations , what method do you think Maple is using to compute determinants?" }} }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 91 "st := time():A60 := RandomM atrix(60,60,generator=rand(-5..5)):Determinant(A60);time() - st;" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 91 "st := time():A61 := RandomMa trix(61,61,generator=rand(-5..5)):Determinant(A61);time() - st;" }}} {EXCHG {PARA 0 "" 0 "" {TEXT 200 0 "" }}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}} {MARK "0 0 0" 0 }{VIEWOPTS 1 1 0 15 10 1804 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }