*! version 1.0  STB-43 ip25
  program define lnsim
          version 5.0
          if "`1'" == "?" {
                  global S_1 "loc scale mean var"
                  exit                        
          }           
          drop _all
          set obs 100
          gen z = exp(loc + scale*invnorm(uniform()))
          summarize z
          post `1' loc scale _result(3) _result(4)
  end