nulldata 50
x = 100 * uniform()
# open a "progressive" loop, to be repeated N times
loop 1000 --progressive --quiet
   u = 10 * normal()
   # construct the dependent variable
   y = 10*x + u
   # run OLS regression
   ols y const x
   # grab the estimated slope coefficient
   a = $coeff[1]
   b = $coeff[2]
   r2 = $rsq
   # arrange for printing of stats
   print a
   print b
   # and save to file
   store coeff.gdt a b
endloop
# now one can open coeff.gdt and examine the distribution