par(mfrow=c(1,2))
data0<-rnorm(100,mean=0,sd=1)
h1<-hist(data0,main="Histogram random normal", xlab="x",col="green",las=1,breaks=15)
data1<-rnorm(100,mean=0,sd=1)
h1<-hist(data1,main="Histogram random normal", xlab="x",col="blue",las=1,breaks=15)