This is bivariate practice from the astroML textbook chapter 3.0 figures