code: def anscombe():
      
      _, ((axa, axb), (axc, axd)) =  plt.subplots(2, 2, sharex='col', sharey='row')
      colors = get_color_cycle()
  
      for arr, ax, color in zip(ANSCOMBE, (axa, axb, axc, axd), colors):
          x = arr[0]
          y = arr[1]
  
          # Set the X and Y limits
          ax.set_xlim(0, 15)
          ax.set_ylim(0, 15)
  
          # Draw the points in the scatter plot
          ax.scatter(x, y, c=color)
  
          # Draw the linear best fit line on the plot
          draw_best_fit(x, y, ax, c=color)
  
      return (axa, axb, axc, axd)
docs: Creates 2x2 grid plot of the 4 anscombe datasets for illustration.