myNorm2 = rnorm(10000)                 # Create 10,000 values from normal dist

png(filename = "illustrations/fig-4-1-normal dist histogram.png",
  units = "in",                        # Set measurements in inches
  res = 1200,                          # Set resolution at 1200dpi
  width = 6,                           # Width at 6 inches
  height = 4)                          # Height at 4 inches

par(mai = c(1, 1, .25, .25))           # Set margins - no Title

hist(myNorm2, main = NA)               # Display in histogram

dev.off()