rm(list=ls()) #clears memory
setwd("/rstuff") #sets path

a1<-c(rnorm(500, mean=15, sd=5)) 
#draws 500 random observations from a normal distribution with mean=15 and s=5
a2<-c(rnorm(500, mean=20, sd=5))
b1<-c(rnorm(500, mean=15, sd=.5))
b2<-c(rnorm(500, mean=20, sd=.5))

aball <- cbind(a1,a2,b1,b2) #binds the 4 columns into a single object
write.csv(aball, file = "repl1.csv")#outputs aball to a csv file