clear all
set obs 10
generate y = rnormal()
sort y
generate z = invnormal(_n/(_N+1))
scatter y z, name(byhand)
qnorm y, name(qnorm)
graph combine byhand qnorm

swilk y
sfrancia y

clear all
set obs 10
forvalues i=1/1000 {
	generate y`i' = rnormal()
	}
	
swilk y*