R code used to produce the graph shown below.

library(ellipse) 
data(airquality) 
complete <- airquality[complete.cases(airquality),] 
plotcorr(cor(complete), col=myblue)