* Computer Assignment # 1, Economics 435, Spring 2000
set nocolor
set maxcol=5010
sample 1 50
* H and W are matrices with one row and 5000 columns
dim H 1 5000
dim W 1 5000
set nowarn
do #=1,5000
genr x=10+uni(10)
genr g=dum(x-18)
stat g/mean=h:#
genr w:#=(h:#)*(1-h:#)/50
endo
* Note the next two commands have "primes" at the end
* This represents matrix transposition (interchange
* rows and columns)
matrix h=h'
matrix w=w'
sample 1 5000
file screen b:cp1.out
stat h w
plot h/histo range groups=10 alternate
plot w/histo range groups=10 alternate