/*** generating random variables***/
set seed 19560713
set obs 1000

gen my= invnorm(uniform())
gen rx=uniform()
gen lx=uniform()

/***Summary statistics***/
sum my rx lx

/***Graph***/

graph my, bin(50)
graph rx, bin(50)
graph lx, bin(50)

/***Random Walk ***/

gen rpos=rx+rx[_n-1]
gen lpos=lx+lx[_n-1]
gen cumrpos=rpos+rpos[_n-1]
gen cumlpos=lpos+lpos[_n-1]

graph rpos lpos, connect(l)
graph cumrpos cumlpos, connect(l)