scatter plot

plot(faithful$eruptions, faithful$waiting)

cor(faithful$eruptions, faithful$waiting)
## [1] 0.9008112