{VERSION 5 0 "IBM INTEL LINUX" "5.0" }
{USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 
1 0 0 0 0 1 }{PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 
0 0 0 0 0 0 0 0 0 0 }0 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 0 }1 0 0 
0 8 4 0 0 0 0 0 0 -1 0 }}
{SECT 0 {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{SECT 0 {PARA 3 
"" 0 "" {TEXT -1 14 "Lecture Note 9" }}{PARA 0 "" 0 "" {TEXT -1 87 "Pr
oducing random numbers with Maple and computing the average of the pro
duced numebers." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> \+
" 0 "" {MPLTEXT 1 0 8 "restart;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 
1 0 25 "with(linalg, randmatrix):" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 0 
"" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "trials := 100;" }}{PARA 0 "> \+
" 0 "" {MPLTEXT 1 0 12 "result := 0;" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 25 "for k from 1 to trials do" }}{PARA 0 "> " 0 "" 
{MPLTEXT 1 0 0 "" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 15 "intresult := 0:
" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 43 "A := randmatrix(5,5, entries= r
and(-3..3));" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "for i from 1 to 5 d
o" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "for j from 1 to 5 do" }}{PARA 
0 "> " 0 "" {MPLTEXT 1 0 31 "intresult := intresult + A[i,j]" }}{PARA 
0 "> " 0 "" {MPLTEXT 1 0 3 "od;" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 3 "o
d;" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}{PARA 0 "> " 0 "" 
{MPLTEXT 1 0 0 "" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 30 "intresult := in
tresult / (25):" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 9 "print(k);" }}
{PARA 0 "> " 0 "" {MPLTEXT 1 0 29 "result := result + intresult;" }}
{PARA 0 "> " 0 "" {MPLTEXT 1 0 3 "od;" }}{PARA 0 "> " 0 "" {MPLTEXT 1 
0 26 "result := result / trials;" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 0 "" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }{TEXT -1 103 "The a
bove program, if apllied several times, creates the same sequence of m
atrices and thus is not very" }}{PARA 0 "" 0 "" {TEXT -1 116 "helpful \+
in performing randomized experiments. The reason is hidden behind Mapl
e's way to produce the  random numbers" }}{PARA 0 "" 0 "" {TEXT -1 
174 "in the matrix. A better way is to include Maple's command Seed = \+
randomize(). Read its description in the manual and test the following
 enlarged version of the above program." }}{PARA 0 "" 0 "" {TEXT -1 0 
"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "restart;" }}{PARA 0 "> \+
" 0 "" {MPLTEXT 1 0 25 "with(linalg, randmatrix):" }}}{EXCHG {PARA 0 "
> " 0 "" {MPLTEXT 1 0 20 "Seed := randomize();" }}{PARA 0 "> " 0 "" 
{MPLTEXT 1 0 14 "trials := 100;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 
1 0 12 "result := 0;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "for
 k from 1 to trials do" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 15 "intresult
 := 0:" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}{PARA 0 "> " 0 "" 
{MPLTEXT 1 0 43 "A := randmatrix(5,5, entries= rand(-3..3));" }}{PARA 
0 "> " 0 "" {MPLTEXT 1 0 20 "for i from 1 to 5 do" }}{PARA 0 "> " 0 "
" {MPLTEXT 1 0 20 "for j from 1 to 5 do" }}{PARA 0 "> " 0 "" {MPLTEXT 
1 0 31 "intresult := intresult + A[i,j]" }}{PARA 0 "> " 0 "" {MPLTEXT 
1 0 3 "od;" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 3 "od;" }}{PARA 0 "> " 0 
"" {MPLTEXT 1 0 28 "intresult := intresult / 25:" }}{PARA 0 "> " 0 "" 
{MPLTEXT 1 0 9 "print(k);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 29 "result
 := result + intresult;" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 3 "od;" }}
{PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "result := result / trials;" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}{MARK "11" 0 }
{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }