<SYSTEM_TASK:>
Solve the following problem using Python, implementing the functions described below, one line at a time
<END_TASK>
<USER_TASK:>
Description:
 def anscombe():
 """ Creates 2x2 grid plot of the 4 anscombe datasets for illustration. """

    _, ((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)