x1<-rnorm(1000,0,1)

eda1<-function(x){par(mfrow=c(2,2), oma=c(4,2,6,2), mar=c(5,5,4,2), hist(x,las=1, main = "Histogram of N(0,1)"), 
boxplot(x,las=1, main = "Boxplot of N(0,1)"),iqr<-summary(x)[5]-summary(x)[2],
plot(density(x, width=2*iqr), xlab="x",ylab="y", type="l", las=1, main="Normal Probability Curve N(0,1)"),
qqnorm(x,las=1, main="Normal Probability Plot for x" ),qqline(x))}

eda1(x1)